中国激光, 2017, 44 (8): 0804001, 网络出版: 2017-09-13   

激光多普勒流速仪的干涉条纹理论分析及测量 下载: 884次

Theory Analysis and Measurement for Interference Fringes of Laser Doppler Velocimeter
作者单位
1 中国计量大学计量测试工程学院, 浙江 杭州 310018
2 中国计量科学研究院, 北京 100013
3 北京理工大学, 北京 100081
摘要
激光多普勒流速仪(LDV)具有准确度高、非接触测量、动态响应快等优点, 在流速测量中得到了广泛应用。条纹的非一致性是影响激光多普勒流速仪测量精度的重要因素, 因此精确获知测量体中干涉条纹间距及梯度分布情况, 是激光多普勒流速仪准确测量流速的前提。对激光多普勒流速仪测量体中的干涉条纹间距及梯度分布进行了理论分析, 揭示了高斯光束固有传播特性与其测量体中干涉条纹分布之间的关系, 确定了影响干涉条纹分布的参数。使用光束分析仪分别获得绿光和紫光的束腰半径为114 μm和83 μm, 并确定了束腰的空间位置, 利用测得的相关参数量化了绿光和紫光测量体中任意位置的干涉条纹间距及梯度分布, 绿光和紫光测量体中归一化后的条纹梯度最大分别可达0.46%和0.60%。通过与转盘装置所测得的绿光和紫光干涉条纹间距结果进行比较, 对干涉条纹分布的理论分析及测量结果进行了验证, 最大相对误差分别为0.87%和0.78%。
Abstract
Laser Doppler velocimeter (LDV) is widely applied in the measurement of fluid velocity with the advantages of high accuracy, non-contact measurement and quick dynamic response. The fringe nonuniformity is a key factor to influence the measurement accuracy of LDV. To accurately measure fluid velocity of LDV, the interference fringe spacing and the gradient distribution in the measuring volume have to be determined exactly. The interference fringe spacing and the gradient distribution of LDV in the measuring volume are analyzed in theory. The relationship between the Gaussian beam propagation characteristics and its interference fringe distribution are revealed. The parameters influencing the interference fringe distribution are determined. The beam waist radii of the green beam and purple beam are 114 μm and 83 μm respectively, which are obtained by beam profiler. Meanwhile, the spatial location of the beam waist is determined. The interference fringe spacing and gradient distribution of green and purple beams in the measuring volume at any locations are quantified by using measured parameters. In the measuring volume of green and purple beams, normalized fringe gradient maximums can reach 0.46% and 0.60%, respectively. The interference fringe spacing results of green and purple beams measured by spinning-disk system are compared, and theoretical analysis and measuring results of the interference fringes distribution are verified. The maximum relative errors are 0.87% and 0.78%, respectively.

1 引言

激光多普勒测速技术自20世纪60年代末出现以来,随着激光技术的高速发展,由于它所具有非接触测量、准确度高、空间分辨率高和动态响应快等优点,使其成为获取准确无干扰的流体或粒子速度的首选技术[1]。目前,激光多普勒流速仪(LDV)已经被广泛应用在能源、水利、化工、医学、冶金、钢铁、航空、机械制造、汽车制造等行业中[2-4]。由于采用激光测量技术,LDV对于流场没有干扰作用,测速范围宽,而且由于多普勒频率与速度成线性关系,和该点的温度及压力都无关,所以LDV是目前世界上速度测量精度最高的仪器。

目前常见的LDV多采用双光束模式,根据光的干涉理论,两束相干光在空间相交,在相交的区域会形成测量体,测量体中会产生明暗相间的干涉条纹。当散射粒子穿越测量体时,粒子散射的光信号的多普勒频率与粒子的运动速度存在线性关系[5],其中多普勒频率和干涉条纹间距的测量精度决定了LDV的测量精度。由于多普勒频率的测量精度高,可达到10-5~10-6 Hz,条纹的非一致性是影响LDV测量精度的重要因素。在实际应用中,由于高斯光束固有的传播特性及光路布局不当等原因,导致LDV发射的两束激光并未在束腰处相交,测量体中的干涉条纹间距不再处处相等,而是存在梯度分布。影响高斯光束传播特性重要的参数之一是束腰处光斑半径,它的大小和与测量体的空间位置关系都影响着干涉条纹的梯度分布。例如两束激光在束腰后相交,则形成的测量体中的干涉条纹间距并不一致,而是沿着光传播方向由小逐渐增大。当统一采用固定值来计算散射粒子穿越测量体中任意位置的流速时,会引起LDV的多普勒频率展宽,测量信号偏置,信噪比下降等后果,从而使流速测量的准确度降低[6]

此外,基于相同的条件,在远离测量体中心位置的条纹梯度分布会呈现一定的非线性。而对一些要求测量湍流强度的场合,条纹的梯度分布及非线性会放大湍流强度的真实值。而且,实际上在同一时刻会有多个粒子穿过测量体中的干涉条纹,信号接收系统会同时接收来自不同粒子的散射信号。而多个粒子的散射信号存在相互的干扰作用,在这种作用下,条纹也会发生畸变,使速度测量的统计结果变差,从而影响LDV的测量准确度[7]。因此,研究测量体中干涉条纹间距及其梯度分布的规律,准确测量测量体内各处的干涉条纹间距,量化条纹的非一致性对测量结果的影响,对修正流速测量结果的偏置和提升LDV测量的准确度水平有重要的意义。

测量体中干涉条纹的非一致性是无法消除且难以修正的,因此研究干涉条纹的非一致性对测量结果的影响尤为重要。最初研究干涉条纹的非一致性是为了解释实际测量中多普勒频率展宽的原因。Hanson[7]在高斯光束的模型下发现多普勒频率展宽主要是由测量体中干涉条纹的梯度分布造成的,并通过线性近似得到了沿光束传播方向上的条纹梯度分布。随后,Durst等[8]验证了Hanson的计算结果,并计算了高斯光束模型下的干涉条纹分布情况。通过计算发现在垂直于光束传播方向上的干涉条纹间距也存在梯度分布,实验结果也验证了其结论的正确性。以此为基础,Durst等[9]通过实验结果分析了半导体LDV系统的精度,后期的研究重心为通过实验测量实际的干涉条纹间距及其梯度分布。Miles等[10-11]在前人的研究基础上应用两种方法对条纹间距及其分布进行了测量和评价,并探索了入射光为非高斯光束的情况。Zhang等[12-13]从湍流度测量的角度提出了条纹的非一致性对测量结果的影响。Mickan等[14-15]出于准确计量的需求,实现了对条纹间距的准确测量,并给出了条纹间距测量结果的不确定度。

前人的研究多集中在条纹的非一致性本身,且实验结果的精度都较低,对条纹间距及其梯度的具体分布情况的关注度不足[16]。而且这些研究并没有通过理论分析建立高斯光束固有传播特性与其条纹分布之间的关系,而是直接通过实验测量干涉条纹间距。在这样的研究现状下,高斯光束固有传播特性与其测量体中的干涉条纹分布之间的关系是割裂的,条纹的非一致性对测量准确度的影响机理并未得到清晰的阐述,不同情况下的条纹间距变化及其梯度分布也未被量化。因此,本文在理论分析和实验结合的基础上,揭示了高斯光束固有的传播特性与其测量体中干涉条纹分布之间的关系,通过实验准确测量高斯光束相关参数即不同光束的束腰大小及空间位置,从而获得不同光束形成的测量体中的条纹间距及其梯度分布情况,提高LDV的测量精度。

2 条纹分布的理论分析

在条纹分布的理想模型中,将激光束视为均匀的平面波,得到简化的条纹间距为:

S=λ2sinα,(1)

式中,λ为激光束的波长,α为两束激光之间夹角的半角。 而实际上,激光束是一种特殊的高斯球面波,为了更深入地研究实际应用中条纹分布的规律,需要从高斯光束的基本性质出发,建立高斯光束固有的传播特性和条纹分布之间的关系,确定影响条纹分布情况的参数。

2.1 高斯光束的基本性质

假设激光在光谐振腔里以TEM00波模传播,这种基模光束的场振幅的横向分布具有高斯函数的形式,如图1所示。因此,基模光束通常称为高斯光束。其中对高斯光束有重要意义的两个参数分别是光束横截面上的光斑半径ω(z)和波前的曲率半径R(z)。在高斯光束的束腰(z=0)处这两个参数的表达式为[17]:

ωz2=ω021+λz/πω022,(2)R(z)=z1+πω02/λz2(3)

由(2)、(3)式可以清楚地看到,ωRz的变化可由束腰半径ω0这一参数完全确定,且ωω0,Rz图2所示为高斯光束和这两个参数随z变化的一般规律,ω(z)随z按双曲线的规律而增大,这表明光束在传播过程中逐渐发散。

图 1. 高斯光束的横向场分布

Fig. 1. Transverse field distribution of Gaussian beam

下载图片 查看所有图片

图 2. 高斯光束

Fig. 2. Gaussian beam

下载图片 查看所有图片

如果高斯光束穿过薄透镜,它将转换成另一束高斯光束,两个基本参数即聚焦后束腰的位置和半径也会随之改变,如图3所示。设新的束腰位置离薄透镜的距离为z1(焦距为f),此时聚焦后的束腰半径为ω1,则新的光束参数为[18]:

z1=f+(z0-f)f2/z0-f2+πω02/λ2,(4)1/ω12=1/ω021-z0/f2+πω0/λ2/f2(5)

当束腰的位置和半径这两个参数已知时,就可以计算出经透镜聚焦后光束任意一点的半径ω和波面的曲率半径R

图 3. 高斯光束的转换

Fig. 3. Transformation of Gaussian beam

下载图片 查看所有图片

2.2 条纹分布的理论分析

在光学近轴近似(即满足sinθθ,cosθ≈1)的条件下,高斯光束可表示为[19]:

E(xi,yi,zi)=2π121ω(zi)expxi2+yi2ω2(zi)×exp-jkzi+(zi)-jk2xi2+yi2R(zi),(6)

式中,ω(zi)和R(zi)分别代表光束的光斑半径和波前的曲率半径。(xi,yi,zi)是以束腰处为原点的坐标系中的坐标。额外相位移φ(zi)为:

φ(zi)=arctanzizRi(7)

图4所示,以两束激光中心线相交处为原点建立新的直角坐标系,坐标系中的坐标为(x,y,z)。设两束激光束腰处的坐标分别为(xw1,zw1)和(xw2,zw2),相交后夹角的半角为α。光束坐标(xi,yi,zi)与任意(x,y,z)的位置之间的关系可通过坐标系转换得到:

x1=xcosα+zsinαy1=yz1=-(x-xw1)sinα+(z-zw1)cosα,(8)x2=xcosα-zsinαy2=yz2=(x-xw2)sinα+(z-zw2)cosα(9)

从(6)式可以得到每束光的相位是:

θi(x,y,z)=-kzi+φ(zi)-k2xi2+yi2R(zi)(10)

而干涉条纹在测量体中具有固定的相位差:

θ1(x,y,z)-θ2(x,y,z)=-k(z1-z2)+φ(z1)-φ(z2)-k2x12+y12R(z1)-x22+y22R(z2)=2πn(11)

将(3)、(7)、(8)、(9)式代入(11)式中化简,y轴方向上的条纹间距变化忽略不计,则可得条纹间距L为:

L=dndx-1=λ2sinα1-12tanα×x1z1z12+zR12-x2z2z22+zR22-1=λ2sinα1+x1z1z12+zR12-x2z2z22+zR222tanα-x1z1z12+zR12-x2z2z22+zR22,(12)

式中,光束瑞利长度为zRiω0i2/λ。(12)式和坐标转换公式(8)、(9)式可以用来计算测量体中任意位置的条纹间距。

图 4. 计算条纹间距分布的几何坐标系

Fig. 4. Geometric coordinate system for calculation of fringe spacing distribution

下载图片 查看所有图片

为了阐明条纹间距变化的量级及系统参数,限制特定的条件可以使(12)式的表达形式变得更加简明。对于典型的双光束模式,两束激光的直径和束腰位置均相同,在这种情况下,即:

zR1=zR2=zRzw1=zw2=zwxw1=-xw2,(13)

(12)式可化简为:

L=λ2sinα1+x1z1(z12+zR2)tanα-x1z1(14)

x=z=0处对(14)式求导可得:

1LLz=-1fDfDz=z1cosαz12+zR2=cosαR(z),(15)

归一化后的条纹梯度与波前曲率半径R的线性近似关系式(15)式同Hanson得到的结果一致。当条纹间距变化仅限纵轴(x=0)时,又因为x1=-x2,z1=z2=z,则(14)式可化为:

L=λ2sinα1+zcos2α(zcos2α-zw)zR2cos2α-zw(zcos2α-zw)(16)

(16)式表明:条纹间距的大小取决于激光波长λ和光束夹角α,而条纹在z轴上的分布梯度主要受高斯光束瑞利长度zRi=π ω0i2和束腰在空间的具体位置zw所影响,且在远离原点的位置会呈现非线性。所以需要通过实验准确测量高斯光束的两个未知参数即束腰处的光斑半径ω0和束腰在z轴的具体位置zw。根据建立的高斯光束相关参数与其测量体中干涉条纹间距之间的关系式(16)式可以计算得出测量体中的条纹分布。

3 实验及结果分析

通过第2节对条纹分布的理论分析,建立了高斯光束固有的传播特性与干涉条纹分布之间的关系,明确了条纹在z轴上的分布梯度主要受高斯光束的束腰处光斑半径ω0和束腰在z轴的具体位置zw所影响。需要通过实验准确测量不同光束的束腰大小及空间位置,从而得到量化的不同光束所形成的测量体中的条纹间距及其梯度分布情况。

3.1 实验布置及过程

本文主要的研究对象为双光束模式LDV的干涉条纹分布情况,其光路布置原理如图5所示。两束相干光由光纤平行发射至前透镜,会聚后在焦距f处形成测量体。当散射粒子穿越测量体时,散射光信号经接收透镜聚焦及视场光阑后到达光电倍增管,对倍增管接收到的光电流信号进行滤波处理和频谱分析,从而得到多普勒频率。在LDV的对面放置光束分析仪,通过调整LDV探头的水平角度及高度,使LDV探头发射出的激光束垂直入射到光束分析仪的镜头中心。使用光束分析仪可获得束腰直径,并通过测量LDV探头到光束分析仪镜头之间的距离,确定束腰的具体位置,通过(16)式计算得到测量体中任意位置条纹间距。

图 5. 双光束模式LDV的光路示意图

Fig. 5. Schematic of light path for double beam pattern LDV

下载图片 查看所有图片

实验现场布置如图6所示,测量环境保持光线强度分布均匀且恒定,以此减少外界环境对测量结果的影响。通过Beamgage软件观察光束,光束轮廓的2D显示效果如图7所示,激光横向的振幅呈高斯分布,每一小格的分辨率为10 μm×10 μm。适当调整LDV发射的激光功率以防相机过曝光,影响实验结果。控制坐标架,使LDV在平行于光束传播的方向往复移动,粗略地找到光斑直径由大变小再变大的转变区域,在此区域内确定某一点为坐标零点,之后正式开始测量。测量步长由10 mm逐渐缩小到1 mm,越靠近光腰的区域测量密度越大。在采集数据的过程中由于数据会不断波动,所以选取每50帧数据的平均值作为此点的数据,等数据平稳后记录并作图。

图 6. 实验现场布置图

Fig. 6. Experiment site layout

下载图片 查看所有图片

图 7. 2D光束轮廓显示图

Fig. 7. 2D figure for beam profile

下载图片 查看所有图片

3.2 实验结果及分析

实验分别对绿光(λ=514.5 nm)和紫光(λ=476 nm)在f=800 mm的焦距条件下的束腰直径及位置进行了测量。相关参数如表1所示。

表 1. 条纹间距理论参数

Table 1. Theory parameters of fringe spacing

Wave length λ /nmBeam diameter /mmBeam spacing /mmBeam angle sinα
514.52.239.070.024
4762.239.990.024

查看所有表

光斑直径随距离变化的曲线如图8所示。绿光的测量结果如图8(a)所示, 以LDV发射透镜到光束分析仪镜头的距离等于795 mm处作为坐标原点,以光束传播方向为正方向,测得光斑直径最小的位置在37 mm处,即LDV发射透镜到光束分析仪镜头的距离为832 mm时为此激光束的束腰处,其直径为l1=620.7 μm。已知LDV的透镜焦距及光束夹角,则绿光束腰在z轴的坐标为:

zω1=(832×cosα)-800=31.752mm,(17)

即两束绿光实际相交于束腰前。紫光的测量结果如图8(b)所示,以LDV发射透镜到光束分析仪镜头的距离等于760 mm处作为坐标原点,以光束传播方向为正方向,测得光斑直径最小的位置在32 mm处,即LDV发射透镜到光束分析仪镜头的距离为792 mm时为此激光束的束腰处,其直径为l2=455.5 μm,则紫光束腰在z轴的坐标为:

zω2=(792×cosα)-800=-8.235mm,(18)

即两束紫光实际相交于束腰后。瑞利长度zRi=πω0i2λω0i为电磁场强度降到轴向的1/e2的点的半径,而光束分析仪测量出的光斑半径定义为电磁场强度降到轴向的1/e的点的半径,所以要对束腰直径的测量结果进行处理,统一除以e=2.72后得到:

ω01=l1/(2×2.72)=114μm,(19)ω02=l2/(2×2.72)=83μm(20)

分别将上述测量得到的高斯光束的相关参数代入(16)式中,计算得到绿光和紫光在f=800 mm的焦距条件下条纹的分布情况,如图9所示。如图9(a)所示,以光束传播方向为正方向,在测量体中7 mm的范围内,绿光的条纹间距由10.714 μm变化到10.393 μm;如图9(b)所示,在测量体中5.5 mm的范围内,紫光的条纹间距由9.701 μm变化到9.862 μm。

图 8. 光斑直径变化图。(a)绿光;(b)紫光

Fig. 8. Variation diagram of the beam diameter. (a) Green beam; (b) purple beam

下载图片 查看所有图片

图 9. 条纹间距及梯度分布图。(a)绿光;(b)紫光

Fig. 9. Fringe spacing and gradient distribution. (a) Green beam; (b) purple beam

下载图片 查看所有图片

在相同的条件下,利用转盘装置分别测量绿光和紫光测量体中的条纹间距。转盘装置的测量原理如图10所示,把一个直径为5 μm的钨丝安装在已知直径为d的金属转盘的圆柱曲面上,来模拟流体中粒子的运动,产生多普勒频率,达到获取测量体中条纹信息的目的。干涉条纹间距的表达式为:

δ=π·υ·dfDwire,(21)

式中,d为转盘直径,υ为转盘的转动频率, fDwire为LDV测量钨丝运动得到的多普勒频率。

图 10. 转盘装置测量条纹间距原理示意图

Fig. 10. Principle diagram of measuring fringe spacing of spinning-disk system

下载图片 查看所有图片

将绿光和紫光测量体中干涉条纹间距的实验结果与理论计算结果进行对比,如图9所示。可以看出,理论计算值和实验测量值的变化趋势具有一致性,在测量体中间区域基本吻合,两端则呈现非线性。绿光条纹间距测量值的最大相对误差为0.87%,紫光条纹间距测量值的最大相对误差为0.78%,结果验证了本文干涉条纹分布的理论分析及测量结果的正确性。

将测量得到的束腰半径代入(14)式分别计算绿光和紫光归一化后的条纹梯度,其变化如图9所示。以激光传播方向为正方向,绿光相交于束腰前,测量体内的条纹间距由大逐渐变小,条纹梯度为负,在测量体内归一化条纹梯度最大可达0.46%;紫光相交于束腰后,测量体内的条纹间距由小逐渐变大。条纹梯度为正,在测量体内归一化条纹梯度最大可达0.60%。可见在测量体中离束腰位置越远处的条纹梯度越大,导致LDV测速的误差也越大。

4 结论

本文对激光多普勒流速仪测量体中的干涉条纹间距及梯度分布进行了理论分析,阐述了条纹非一致性对测量准确度的影响机理,揭示了高斯光束固有传播特性与其测量体中干涉条纹分布之间的关系,并量化了两种波长下的条纹间距变化及其梯度分布,对提高LDV的测量精度有重要的意义。通过理论分析确定了影响干涉条纹分布的参数,使用光束分析仪分别获得绿光和紫光的束腰半径为114 μm和83 μm,并确定了束腰的空间位置,利用实验测得的相关参数最终获得了绿光和紫光测量体中任意位置的干涉条纹间距及梯度分布,绿光和紫光测量体中归一化后的条纹梯度最大分别可达0.46%和0.60%。通过与转盘装置所测得的绿光和紫光干涉条纹间距结果进行比对,对本文干涉条纹分布的理论分析及测量结果进行了验证,最大相对误差分别为0.87%和0.78%。

参考文献

[1] Miles P C. Geometry of the fringe field formed in the intersection of two Gaussian beams[J]. Applied Physics, 1996, 35(30): 5887-5895.

[2] 孙渝生, 赵建新, 范丽娟. 激光多普勒测速技术的最新发展[J]. 传感器世界, 1998( 8): 20- 26.

[3] ToyoshimaM, OhashiT, ShinoharaS, et al. New digital displacement measuring circuit for aperiodic vibration using self-mixing type LDV[C]. Instrumentation and Measurement Technology Conference, 1994: 336- 339.

[4] Nam SH, Park SK, Kim DH, et al. Time series analysis on dysfunction of the autonomous nerve system measured non-invasively with the temperature variation using laser Doppler velocimeter[C]. International Conference of the IEEE Engineering in Medicine & Biology Society, 1997: 1701- 1702.

[5] 沈熊. 激光多普勒测速技术及应用[M]. 北京: 清华大学出版社, 2004: 17-36, 264- 272.

[6] Ruck B. Distortion of LDA fringe pattern by tracer particles[J]. Experiments in Fluids, 1991, 10(6): 349-354.

[7] Hanson S. Broadening of the measured frequency spectrum in a differential laser anemometer due to interference plane gradients[J]. Applied Physics, 1973, 6(2): 164-172.

[8] Durst F, Stevenson W H. Influence of Gaussian beam properties on laser Doppler signals[J]. Applied Optics, 1979, 18(4): 516-524.

[9] Durst F, Muller R, Naqwi A. Measurement accuracy of semiconductor LDA systems[J]. Experiments in Fluids, 1990, 10(2): 125-137.

[10] Miles P C, Witze P O. Fringe field quantification in an LDV probe volume by use of a magnified image[J]. Experiments in Fluids, 1994, 16(5): 330-335.

[11] Miles P C. Geometry of the fringe field formed in the intersection of two Gaussian beams[J]. Applied Optics, 1996, 35(30): 5887-5895.

[12] Zhang Z, Eisele K. On the broadening of the flow turbulence due to fringe distortion in LDA measurement volumes[J]. Experiments in Fluids, 1998, 25(4): 371-374.

[13] Zhang Z, Eisele K. On the overestimations of the flow turbulence due to fringe distortion in LDA measurement volumes[J]. Experiments in Fluids, 1998, 25(4): 371-374.

[14] MickanB. Uncertainty about the LDA calibration facility[R]. PTB Report, 2000.

[15] Yeh TT, Hall JM. Uncertainty of Nist airspeed calibrations[R]. Nist Report, 2000.

[16] CuiL, WangC, LiC. Establishment of new air velocity standard in NIM[C]. 16 th Flow Measurement Conference , 2013.

[17] HansonS. Visualization of alignment errors and heterodyning constraints in laser Doppler velocimeters[C]. Copenhagen: Proceedings of the LDA-Symposium, 1976: 176- 182.

[18] DurstF, Stevenson WH. Moire patterns to visually model laser-Doppler signals[C]. Copenhagen: Proceedings of the LDA-Symposium, 1976: 183- 205,.

[19] Naqwi A, Durst F. Focusing of laser diode beams: a simple mathematical model[J]. Applied Optics, 1990, 29(12): 1780-1785.

刘子君, 崔骊水, 谢代梁. 激光多普勒流速仪的干涉条纹理论分析及测量[J]. 中国激光, 2017, 44(8): 0804001. Liu Zijun, Cui Lishui, Xie Dailiang. Theory Analysis and Measurement for Interference Fringes of Laser Doppler Velocimeter[J]. Chinese Journal of Lasers, 2017, 44(8): 0804001.

本文已被 3 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!