CN101576568A - 单相机全场多普勒测速装置 - Google Patents
单相机全场多普勒测速装置 Download PDFInfo
- Publication number
- CN101576568A CN101576568A CNA2009100993672A CN200910099367A CN101576568A CN 101576568 A CN101576568 A CN 101576568A CN A2009100993672 A CNA2009100993672 A CN A2009100993672A CN 200910099367 A CN200910099367 A CN 200910099367A CN 101576568 A CN101576568 A CN 101576568A
- Authority
- CN
- China
- Prior art keywords
- measuring device
- speed measuring
- lens
- light
- iodine
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 claims abstract description 17
- 229910052740 iodine Inorganic materials 0.000 claims abstract description 17
- 239000011630 iodine Substances 0.000 claims abstract description 17
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 230000001360 synchronised effect Effects 0.000 claims abstract description 7
- PNDPGZBMCMUPRI-UHFFFAOYSA-N iodine Chemical compound II PNDPGZBMCMUPRI-UHFFFAOYSA-N 0.000 claims description 11
- 239000005304 optical glass Substances 0.000 claims description 9
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 6
- 239000007787 solid Substances 0.000 claims description 5
- 238000000827 velocimetry Methods 0.000 claims 2
- 238000005259 measurement Methods 0.000 abstract description 16
- 239000002245 particle Substances 0.000 description 11
- 230000003595 spectral effect Effects 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000000034 method Methods 0.000 description 5
- 238000002834 transmittance Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000000498 cooling water Substances 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000001499 laser induced fluorescence spectroscopy Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种单相机全场多普勒测速装置。Nd:YAG激光器发出的调制激光经分光镜进行分光后,一路经圆柱透镜和第一球面透镜形成片光照亮流场测试区域,另一路经平面镜、第一碘分子滤波器、第二球面透镜、光电二极管和计算机接口接到控制计算机,计算机接口的三路控制信号分别经频率控制器、温度控制器、同步控制器与Nd:YAG激光器连接,流场测试区域的散射光经透镜、第二碘分子滤波器和CCD相机接入控制计算机。与传统DGV测速装置相比,只采用一个CCD相机,简化测速装置,降低设备成本。与PIV测速装置相比,可以满足大尺度流场全场测速的需要。针对大尺度流场,具有超过传统DGV和PIV两种测速装置的测量精度。
Description
技术领域
本发明涉及流速的测量装置,具体地说是涉及一种单相机全场多普勒测速装置。
背景技术
随着激光、计算机和传感器等技术的发展,用于流场测量的新技术和新仪器不断被创造出来,如激光多普勒风速仪(LDA)、粒子成像测速仪(PIV)、平面激光诱导荧光技术(PLIF)等应用十分广泛。
但目前还有一些流场没有很好手段来进行非接触测量或测量中存在较大误差,例如,高马赫数和超高马赫数流场测量、大尺度风洞流场实时测量、燃烧场多参数实时测量等。对于基于PIV技术的测量装置而言,由于对颗粒散射光要求很高,要求流场内颗粒较大,但对于高马赫数流场而言跟随性会随之降低,颗粒太小则散射光强变弱,粒子配对成功率下降,无法满足全场测速的要求。近些年发展起来的多普勒全场测速(DGV)技术应用在上述应用中具有一定潜力。
DGV的基本工作原理是测量运动粒子散射光的多普勒频移。当激光被移动粒子散射时,会产生多普勒频移,根据多普勒频移公式有:
DGV使用一个窄线宽激光来照亮流场中的某一平面。示踪粒子散射光经过分子滤波器后由一相机采集。同时该图像在不经过滤波器的情况下被另外一个相机-参考相机所采集。分子滤波器是一个简单的光学元件(两端开窗的玻璃圆柱),内装吸收分子,该分子具有能够和激光光谱相匹配的吸收带。这就形成了一个具有有限长度斜边的透射率曲线(如图1所示),通过滤波器的光谱强度和光的频率有关。Iv/I0v是分子滤波器的光谱透射率,Iv定义为通过滤波器后的光谱强度(在频率v处),I0v为滤波器之前的光谱强度。通过滤波器的光谱强度是散射光的光谱强度与滤波器分子吸收线的卷积。信号相机的每一个像素记录了积分光谱强度,I=∫Ivdv。参考相机采集了未通过滤波器的流场图像,即与未经滤波器的积分光谱强度I0(I0=∫I0vdv)。以积分透射率为自变量,频移(或频率函数ζ)为因变量绘制频移-透射率变化曲线(如图2所示),实际测试中,一旦每个像素点上的透射率通过两个相机信号对比得到后,多普勒频移就可以根据滤波器的频率函数计算出,示踪粒子速度(即流动速度)可以由频移计算出。
尽管DGV有其优点,但是在应用过程中还存在一些问题。为消除激光频率和强度变化对测量结果的影响,需用参考信号来进行修正,因此DGV技术中每个速度分量需用两个相机,即信号相机和参考相机,由此造成设备造价昂贵,且不可避免地存在光路对正误差。为此,研究人员提出多种改进方法,如双频激光、镜头分区等,这些方法效果都不是很理想。
发明内容
针对目前全场测速技术中存在的缺陷,本发明的目的在于提供一种单相机全场多普勒测速装置。
本发明解决其技术问题所采用的技术方案是:
Nd:YAG激光器发出的调制激光经分光镜进行分光后,一路经圆柱透镜和第一球面透镜形成片光照亮流场测试区域,另一路经平面镜、第一碘分子滤波器、第二球面透镜、光电二极管和计算机接口接到控制计算机,计算机接口的三路控制信号分别经频率控制器、温度控制器、同步控制器与Nd:YAG激光器连接,流场测试区域的散射光经透镜、第二碘分子滤波器和CCD相机接入控制计算机。
所述的第一、第二碘分子滤波器均为在石英玻璃筒体的两端加上光学玻璃窗口,石英玻璃筒体外壁分别开有与真空泵相连的连接口,与冷端盲头的连接接口,经截止阀与装有固态碘的侧臂的连接接口。
所述的第一碘分子滤波器是接受流场测试区域的散射光经透镜后通过一个光学玻璃窗口接收,再经过另一个光学玻璃窗口,由CCD相机接收;第二碘分子滤波器是接受经平面镜的反射光,输出至第二球面透镜。
本发明具有的有益的效果是:
1.与传统DGV测速装置相比,只采用一个CCD相机,简化测速装置,降低设备成本。
2.与PIV测速装置相比,可以满足大尺度流场全场测速的需要。
3.针对大尺度流场,具有超过上述两种测速装置的测量精度。
附图说明
图1是滤波器透射率曲线。
图2是滤波器多普勒频移函数曲线。
图3是频率调制DGV光强调制示意图。
图4是典型的检测信号及其频谱特性。
图5是一次谐波和二次谐波幅值A1fm(Δf)和A2fm(Δf)。
图6是比值q(Δf)=A1fm(Δf)/A2fm(Δf)。
图7是系统实施示意图。
图8是碘分子滤波器结构图。
图9是频率采样点分布原理图。
图10是一、二次谐波比值算法原理图。
图中:1、Nd:YAG激光器,2、频率控制器,3、温度控制器,4、同步控制器,5、计算机接口,6、控制计算机,7、CCD相机,8碘分子滤波器,9、平面镜,10、球面透镜,11、圆柱透镜,12、分光镜,13、光电二极管,14、流场测试区域,15、透镜,16、热电偶,17、连接口,18、冷端盲头,19、截止阀,20、侧壁,21、冷却水,22、筒体。
具体实施方式
下面结合附图和实施方式对本发明作进一步说明。
本发明解决其技术问题的基本方法是采用频率调制激光代替常规DGV中恒定频率激光(测量期间),假设其频率
fL(t)=fc+fhcos(2πfmt)
其中,fc为激光中心频率,fh为频率调制幅度,fm为调制频率。
由于分子滤波器的频率变化-光强变化的转换作用,通过分子滤波器后,相机(或探头)检测到的是光强调制信号(如图3所示)。由于分子滤波器透射率曲线是非线性的,光强信号中包含了各次谐波。如图4所示,给出的是典型光强调制信号及其能谱(仅显示一次和二次谐波)。当激光中心频率变化时,一、二次谐波的幅值随之变化,如图5所示,示意性给出了一次谐波和二次谐波幅值随激光中心频率的变化曲线。很显然,粒子运动产生的散射光频移相当于激光中心频率发生变化,而这个变化将反应在两个谐波的幅值变化上。而这两个谐波幅值A1fm(Δf)和A2fm(Δf)都正比于光强,比值q(Δf)=A1fm(Δf)/A2fm(Δf)的数值仅取决于中心频率的位置(如图6所示),因此,FM-DGV中考察多普勒频移引起的比值q(Δf)的变化,即可得到多普勒频移,不再需要常规DGV系统中的参考图像,分光过程自然也省略了。
如图7所示,Nd:YAG激光器1发出的调制激光经分光镜12进行分光后,一路经圆柱透镜11和第一球面透镜10形成片光照亮流场测试区域14,另一路经平面镜9、第一碘分子滤波器8、第二球面透镜10、光电二极管13和计算机接口5接到控制计算机6,计算机接口5的三路控制信号分别经频率控制器2、温度控制器3、同步控制器4与Nd:YAG激光器1连接,流场测试区域14的散射光经透镜15、第二碘分子滤波器8和CCD相机7接入控制计算机6。
如图7所示,Nd:YAG激光器1发出的调制激光首先进行分光,大部分用于测量,一小部分用于激光频率和能量监测;第一部分激光经过圆柱透镜11和球面透镜10形成片光,用于照亮被测流场。示踪粒子的散射光经过透镜15和碘分子滤波器8后被CCD相机7采集。相机采集的一组图像(至少三幅),传输到计算机6进行处理和计算。系统中Nd:YAG激光器1采用的是单块固体环形连续激光器(532nm),激光线宽可窄至100kHz,能量为1W。频率控制器2将激光信号调制成正弦波信号,调制频率为2~21kHz,调制幅值为±200MHz~400MHz,温度控制器3采用PID控温技术,控温精度控制到0.01度,同步控制器4用于负责激光照明和相机同步。本测速装置中,频率控制器2选用FC200,温度控制器3选用TC350,同步控制器4选用DG580,计算机接口5采用的控制卡为PCI1506。激光频率监测,激光相机同步以及激光温度监测等功能通过控制计算机6的操作软件和数据处理软件来实现。
如图8所示,所述的第一、第二碘分子滤波器8均为在石英玻璃筒体22的两端加上光学玻璃窗口,石英玻璃筒体22外壁分别开有与真空泵相连的连接口17,与冷端盲头18的连接接口,经截止阀19与装有固态碘的侧臂20的连接接口。
所述的第一碘分子滤波器8是接受流场测试区域14的散射光经透镜15后通过一个光学玻璃窗口接收,再经过另一个光学玻璃窗口,由CCD相机7接收;第二碘分子滤波器8是接受经平面镜9的反射光,输出至第二球面透镜10。
筒体部分用加热带缠绕并作绝热处理,保持一定的加热温度(340-380K)。滤波器筒体22温度控制不要求很精确,只要保证碘不粘在窗口上即可。与筒体相连的侧臂20保持在一个较低的温度,一般310-325K。由于滤波器工作时对侧臂20温度特别敏感,通常通过一个精密系统来保持温度恒定(如恒温水浴),通过冷却水21和热电偶16的反馈值来保持温度恒定。滤波器上的冷端盲头18用来调节碘蒸汽分压(摩尔浓度),以便获得理想的透射率曲线。当固态和气态碘达到平衡后,关闭截至阀19,使侧臂20与滤波器筒体22分开。
CCD相机13采用Nikon单CCD相机,采集频率为10Hz,最高可达13Hz,单幅分辨率为1600×1200,采集的信号直接通过1394或usb通讯接口存储到控制计算机6。控制计算机可采用目前主流的计算机,如DELL的OPTIPLEX330品牌机。本装置中平面镜9选用Φ30×5,球面透镜10选用f200,圆柱透镜11选用f50,分光镜12分光比针对532nm为10∶1,透镜15选用f50。
数据处理算法原理,如图9,10所示。对于周期性调制曲线,只需测量一个周期内的数据点就可以,理论上知道(0)-(4)五个点的值就可以确定这条曲线。g(1)和g(3)数值可假定为相同,g(0)和g(4)为同一值,这样最少只需采集3幅图像就可以满足计算要求,采用公式如下所示,
Ai,j(1fM)和Ai,j(2fM)分别表示一次谐波幅值和二次谐波的幅值gi,j(0),gi,j(1),gi,j(2)为对应的三次图像采集信号。
通过上述公式,可以对所拍摄流场的图像进行单个像素点分析处理,获该像素点所对应的一次和二次谐波幅值之比,并根据标定过的其商与对应的多普勒频移对应的曲线关系,获得该点的多普勒频移量,最终获得该点的速度值。同样,对每个像素点进行类似的数据处理,进而可获得整个流场的速度分布。
Claims (3)
1、一种单相机全场多普勒测速装置,其特征在于:Nd:YAG激光器(1)发出的调制激光经分光镜(12)进行分光后,一路经圆柱透镜(11)和第一球面透镜(10)形成片光照亮流场测试区域(14),另一路经平面镜(9)、第一碘分子滤波器(8)、第二球面透镜(10)、光电二极管(13)和计算机接口(5)接到控制计算机(6),计算机接口(5)的三路控制信号分别经频率控制器(2)、温度控制器(3)、同步控制器(4)与Nd:YAG激光器(1)连接,流场测试区域(14)的散射光经透镜(15)、第二碘分子滤波器(8)和CCD相机(7)接入控制计算机(6)。
2、根据权利要求1所述的一种单相机全场多普勒测速装置,其特征在于:所述的第一、第二碘分子滤波器(8)均为在石英玻璃筒体(22)的两端加上光学玻璃窗口,石英玻璃筒体(22)外壁分别开有与真空泵相连的连接口(17),与冷端盲头(18)的连接接口,经截止阀(19)与装有固态碘的侧臂(20)的连接接口。
3、根据权利要求1或2所述的一种单相机全场多普勒测速装置,其特征在于:所述的第一碘分子滤波器(8)是接受流场测试区域(14)的散射光经透镜(15)后通过一个光学玻璃窗口接收,再经过另一个光学玻璃窗口,由CCD相机(7)接收;第二碘分子滤波器(8)是接受经平面镜(9)的反射光,输出至第二球面透镜(10)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100993672A CN101576568B (zh) | 2009-06-04 | 2009-06-04 | 单相机全场多普勒测速装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100993672A CN101576568B (zh) | 2009-06-04 | 2009-06-04 | 单相机全场多普勒测速装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101576568A true CN101576568A (zh) | 2009-11-11 |
CN101576568B CN101576568B (zh) | 2010-12-01 |
Family
ID=41271552
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100993672A Expired - Fee Related CN101576568B (zh) | 2009-06-04 | 2009-06-04 | 单相机全场多普勒测速装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101576568B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103472256A (zh) * | 2013-09-25 | 2013-12-25 | 东南大学 | 基于面阵ccd空间滤波器的流动二维速度场测量方法及装置 |
CN106771344A (zh) * | 2016-12-20 | 2017-05-31 | 中国科学院力学研究所 | 高速稀薄气体流场的单束激光多维速度测量系统及方法 |
CN108761126A (zh) * | 2018-03-16 | 2018-11-06 | 天津大学 | 一种基于飞秒激光光致化学发光的测速装置及方法 |
CN111381620A (zh) * | 2019-12-12 | 2020-07-07 | 贵州电网有限责任公司 | 防冰摄像机加热策略软件实现方案 |
CN115560948A (zh) * | 2022-11-14 | 2023-01-03 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种干涉瑞利散射湍流度测量方法及装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1238704C (zh) * | 2004-03-23 | 2006-01-25 | 天津大学 | 高速复杂流场的全场光学测量方法与测量装置 |
US20060044547A1 (en) * | 2004-08-25 | 2006-03-02 | Hoffman Richard G Ii | Photonic crystal based optical doppler projectile detection system |
US7773229B2 (en) * | 2007-07-26 | 2010-08-10 | The United States Of America As Represented By The Secretary Of The Navy | Doppler asymmetric spatial heterodyne spectroscopy |
CN201413339Y (zh) * | 2009-06-04 | 2010-02-24 | 中国计量学院 | 一种单相机全场多普勒测速装置 |
-
2009
- 2009-06-04 CN CN2009100993672A patent/CN101576568B/zh not_active Expired - Fee Related
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103472256A (zh) * | 2013-09-25 | 2013-12-25 | 东南大学 | 基于面阵ccd空间滤波器的流动二维速度场测量方法及装置 |
CN103472256B (zh) * | 2013-09-25 | 2015-09-16 | 东南大学 | 基于面阵ccd空间滤波器的流动二维速度场测量方法及装置 |
CN106771344A (zh) * | 2016-12-20 | 2017-05-31 | 中国科学院力学研究所 | 高速稀薄气体流场的单束激光多维速度测量系统及方法 |
CN106771344B (zh) * | 2016-12-20 | 2019-06-28 | 中国科学院力学研究所 | 高速稀薄气体流场的单束激光多维速度测量系统及方法 |
CN108761126A (zh) * | 2018-03-16 | 2018-11-06 | 天津大学 | 一种基于飞秒激光光致化学发光的测速装置及方法 |
CN111381620A (zh) * | 2019-12-12 | 2020-07-07 | 贵州电网有限责任公司 | 防冰摄像机加热策略软件实现方案 |
CN115560948A (zh) * | 2022-11-14 | 2023-01-03 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | 一种干涉瑞利散射湍流度测量方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN101576568B (zh) | 2010-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101576568B (zh) | 单相机全场多普勒测速装置 | |
US7407817B2 (en) | Surface plasmon resonance sensors and method for detecting samples using the same | |
CN101504367B (zh) | 同时监测一氧化碳和二氧化碳浓度的装置 | |
CN102042971B (zh) | 一体化光学臭氧产量检测装置及标定方法和测量方法 | |
CN103076107A (zh) | 基于太赫兹脉冲测量的燃烧温度传感装置及方法 | |
CN201837582U (zh) | 一体化光学臭氧产量检测装置 | |
CN108760681A (zh) | 一种基于波形分解的路径平均温度测量系统与方法 | |
CN112378616B (zh) | 基于波长调制吸收光谱的高速流场多参数测量系统及方法 | |
CN102967583B (zh) | 一种用于测量液相气体折射率的测量仪及方法 | |
CN104280362A (zh) | 一种高温水汽激光光谱在线检测系统 | |
CN109991189A (zh) | 一种基于波数漂移修正的固定点波长调制气体浓度测量装置及其测量方法 | |
CN111351770A (zh) | 一种波长调制吸收光谱吸收光程与气体浓度同步反演方法 | |
CN101122555A (zh) | 基于后向光子相关光谱的高浓度超细颗粒测量装置及方法 | |
CN201413339Y (zh) | 一种单相机全场多普勒测速装置 | |
CN114384045A (zh) | 一种痕量气体浓度和路径长度实时检测系统与方法 | |
CN103344607B (zh) | Tdlas中激光波长检测控制装置及其控制方法 | |
CN204575674U (zh) | 基于旋转光栅的激光多普勒测速装置 | |
CN202974864U (zh) | 一种用于测量液相气体折射率的测量仪 | |
CN111337454B (zh) | 基于激光干涉技术的溶液浓度快速检测方法 | |
CN108776087A (zh) | 一种用于测量气泡参数的相位临界角散射测量方法及装置 | |
CN115266582A (zh) | 一种法布里-波罗谐振腔式盐度传感器 | |
CN106483161A (zh) | 一种基于反射式光纤位移传感器的液体体胀系数测量仪 | |
Sun et al. | A method for measuring and correcting errors in the thickness of semiconductor thin films based on reflection spectroscopy fitting technology | |
CN115343233B (zh) | 一种开放路径上痕量气体浓度的实时测量方法与装置 | |
CN108037143A (zh) | 一种气体折射率的测量方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20101201 Termination date: 20130604 |