CN104698219B - 一种基于近场散射的流动二维速度场测量装置及方法 - Google Patents

一种基于近场散射的流动二维速度场测量装置及方法 Download PDF

Info

Publication number
CN104698219B
CN104698219B CN201510118897.2A CN201510118897A CN104698219B CN 104698219 B CN104698219 B CN 104698219B CN 201510118897 A CN201510118897 A CN 201510118897A CN 104698219 B CN104698219 B CN 104698219B
Authority
CN
China
Prior art keywords
field
laser
images
measurement
light
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.)
Active
Application number
CN201510118897.2A
Other languages
English (en)
Other versions
CN104698219A (zh
Inventor
许传龙
谭浩
张彪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201510118897.2A priority Critical patent/CN104698219B/zh
Publication of CN104698219A publication Critical patent/CN104698219A/zh
Application granted granted Critical
Publication of CN104698219B publication Critical patent/CN104698219B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种基于近场散射的流动二维速度场测量方法及装置,其中测量装置包括一激光器、一空间滤波器、一准直透镜、一透镜组、一CCD相机及一计算机,与现有技术相比,本发明具有以下优点:该测量装置相当简单,对光路对中要求不高,并且在无需添加去除中心光强的复杂装置的情况下,可实现杂散光的有效去除;该测量技术所需的示踪粒子粒径可小于激光器发出的入射光波长,并且示踪粒子的浓度也可以相对较高,可应用于纳米流体流动二维速度场的测量;该测量技术的数据分析简单、快速,可实现在线测量,并且采用差分处理方法保证了该速度测量方法精度和线性度不受限于示踪粒子的粒径及浓度大小。

Description

一种基于近场散射的流动二维速度场测量装置及方法
技术领域
本发明属于速度测量技术领域,具体涉及一种基于近场散射的流动二维速度场测量方法及装置。
背景技术
在现代流体动力学实验中,如何精确、快速且简单地测量流体流动的流场分布一直是速度测量领域中一个很具挑战性的难题。通常是在流体中加入很小的示踪粒子,通过照射流体,检测散射光来实现速度分布的测量。目前,常见的测量二维速度场的技术主要有激光散斑测速(LSV)、在LSV基础上发展起来的粒子图像测速(PIV)和全息粒子图像测速(HPIV)等。
LSV和PIV(HPIV)都是通过比较已知时间间隔的两帧图像,采用空间互相关方法分析图像,可以测量出局部示踪粒子的位移,进而计算出流体流动二维速度场分布。PIV是20世纪80年代发展起来的研究流体动力学的一种非常重要的测量技术,但由于对流场中示踪粒子浓度有很大的限制,实验测量的流场空间分辨率并不理想,目前商品化的PIV仪器设备复杂、价格昂贵,这极大地限制了PIV测量仪的应用及发展。随后发展起来的HPIV是在原有的PIV技术上改良的可用来测量三维速度场分布的一种速度测量技术,但该技术的装置较为复杂,数据处理具有滞后性,并且在线测量的实现难度较大。此外,当实验投放的粒子浓度过高时,散斑噪音会严重影响干涉成像效果和粒子识别。
流体中加入的示踪粒子的粒径和浓度制约着如上现有的测量技术,无论是LSV还是PIV(HPIV),所使用的示踪粒子的粒径不能过小,浓度不能太高,否则无法保证系统的精度和线性度,而且杂散光对系统的干扰也较大。
发明内容
技术问题:本发明针对现有二维速度场测量技术如LSV、PIV等在测量上的不足,提出一种基于近场散射的流动二维速度场测量方法及装置。该测量方法可实现垂直于光轴的流动方向上二维速度场的实时测量,并且与示踪粒子的浓度和粒径无关。
技术方案:本发明基于近场散射光学原理,提出一种基于近场散射的流动二维速度场测量方法及装置,其基本思路为:所述的激光器发出相干光束,经由所述空间滤波器,滤除大部分由激光器产生的杂散光,再由所述透镜组的准直透镜,对激光进行准直,准直光束照射待测溶液发生散射,然后调节测量区域后的所述透镜焦距,使得距离测量区域Z处的透射光与散射光的叠加所成的散斑图像,经所述透镜放大,被所述CCD相机所记录,每隔一定的时间间隔由CCD相机采样一次,并且由计算机存储处理。采用自相关的方法对所得的散斑图像进行处理分析,得到流体的二维流场分布。
一种基于近场散射的流动二维速度场测量装置,其特征在于,包括:
一激光器,用于发出相干光束;
一空间滤波器,用于滤除由激光器产生的杂散光,得到高斯光束;
一准直透镜,用于将激光器发出的光束转换成直光束;
一透镜组,用于设定焦距,将距离测量区域Z处的散斑放大成像;
一CCD相机,用于设定时间间隔采集所述透镜组成像的近场散斑图像;
一计算机,对所述CCD相机采集的近场散斑图像进行计算得到流体的流场分布。
所述计算机对近场散斑图像进行计算的方法是:通过比较已知时间间隔的两幅图像,采用自相关方法,得到在该时间间隔内流体的位移,从而计算出流体的流场分布。
所述的激光器为连续发光的He-Ne激光器。
一种基于近场散射的流动二维速度场测量装置测量流动二维速度场的方法,其特征在于,包括如下步骤:
步骤一、激光器发出的相干光束照射含有示踪粒子的流体,调整透镜组焦距,使得CCD相机采集距离管道中心Z处的图像,并设定CCD采样时间间隔τ,通过CCD相机获取透射光和散射光干涉后光强分布的散斑图像:
其中,表示散斑图像,其中是对应的成像平面的位置向量,t是对应的成像时刻,e0是透射光场,|e0|2是静态透射光强,es是散射光场,Re表示取实部运算,*号是取共轭运算,其中距离Z的选取必须满足D*=2Z·NA<D,NA是透镜的数值孔径,D是相干光源的光束尺寸;
步骤二、CCD相机采集的N帧散斑图像相邻两组图像采样的时间间隔为τ,对相邻两组图像散斑强度值作差,得到N-1组差分信号图像
步骤三、将所得到的每组差分信号图像划分成S个小的计算区域,计算每个小的子区域内的平均速度,平均速度的计算方法是:
A、根据差分信号图像的自相关函数得到该时间间隔内示踪粒子的运动位移其中,差分信号的自相关函数为:
其中,表示示踪粒子的运动位移,第一项是正值,表示的是散射场的自相关函数,它正比于测量区域内示踪粒子的数量,均为散射场的互相关函数;
的定义为:
其中,<…>表示总体平均运算符;
散射场的互相函数值是其自相关函数值的一种平移,即
B、根据步骤A得到的示踪粒子的运动位移计算s示踪粒子的运动速度为
有益效果:与现有技术相比,本发明具有以下优点:该测量装置相当简单,对光路对齐的要求不高,并且在无需添加去除中心光强的复杂装置的情况下,可实现杂散光的有效去除;该测量技术所需的示踪粒子粒径可小于激光器发出的入射光波长,并且示踪粒子的浓度也可以相对较高,可应用于纳米流动二维速度场的测量;该测量技术的数据分析相当简单、快速,可实现在线测量,并且采用差分处理方法保证了该测量方法相对于示踪粒子数的线性度。
附图说明
图1为本发明所述基于近场散射的流动二维速度场测量装置的结构示意图;
图2为差分图像子区域划分示意图。
图3为子区域内自相关计算结果理论示意图。
其中,激光器1、空间滤波器2、准直透镜3、透镜组4、CCD相机5、计算机6。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明。应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明的工作原理:激光器发出相干光束,经由空间滤波器,滤除大部分由激光器产生的杂散光,再由准直透镜,对激光进行准直,准直光束照射待测溶液,因颗粒的存在而发生散射,然后调节测量区域后的透镜焦距,使得距离测量区域Z处的透射光与散射光的叠加所成的散斑图像,被CCD相机所记录,每隔一定的时间间隔由CCD相机采样一次,并且由计算机存储处理。采用自相关的方法对所得的散斑图像进行处理分析,得到流体的二维流场分布。
如图1所示,本发明所述基于近场散射的颗粒粒度测量装置主要包括激光器1、空间滤波器2、准直透镜3、透镜组4、CCD相机5、计算机6。光源可采用连续激光器线光源1(如He-Ne激光器1,波长为632.8nm,D~10mm)。透镜组2为普通光学镜头(如透镜放大倍数40x,NA~0.65,具体参数可以根据实验需求进行选择)。CCD相机5通过信号电缆与计算机6相连。在计算机6内的图像数据处理软件的控制下,可完成图像采集、存储与处理等整体协调。
结合图1,下面结合长度为2mm,壁厚为1mm的方管内待测溶液颗粒粒度测量,以He-Ne激光器(λ~632.8nm,D~10mm)为光源,放大倍数M=20的透镜(NA~0.5),12位的CCD相机(1024×768pixel,pixel sizeΔl~4.65μm,frame rate~800fps),对本发明的工作过程及图像处理流程具体说明:
步骤一、CCD相机采样平面距离方管Z=1.5mm(D*=2Z·NA=2×1.5mm×0.5=1.5mm<D=10mm,满足近场散射条件)。调整CCD相机位置,使得镜头的焦点刚好在采样平面上。激光器发出的相干光束照射管道,管道内是含有示踪粒子的纳米流体,由于光的散射作用,CCD相机采集距离测量区域Z处的由透射光和散射光干涉形成的散斑图像,并传输到计算机上进行存储和处理。设定CCD相机每隔τ=0.1s拍摄一次图像,拍摄N=600帧,总的测量时间T=N·τ=60s。
步骤二、CCD相机采集的600帧散斑图像相邻两组图像采样的时间间隔为τ,对这相邻两组图像散斑强度值作差,得到差分信号选取t1和t2这两个时刻的图像作差分处理如下:
从式(2)可以清晰地看到差分信号与静态透射光强度|e0|2无关。
将所得到的1024×768像素的差分图像按照16×16划分成S=3072个小的计算区域,记为(δf12)1,(δf12)2……(δf12)3072。对于每一个小的区域按步骤三所述计算得到该子区域内的平均速度。
步骤三、计算差分信号的自相关函数这样可以得到该时间间隔内示踪粒子的运动位移自相关函数在图像上的理论结果示意图如图3所示。
式(2)中,表示示踪粒子的运动位移,第一项是正值,表示的是散射场的自相关函数,它正比于测量区域内示踪粒子的数量,但不包含示踪粒子的运动信息。均为散射场的互相关函数。
的定义为:
式(3)中,<…>表示总体平均运算符。
当流体中所有的示踪粒子以相同的平均速度<v>垂直于光轴运动τ时间,则在t+τ时刻时,在成像平面位置的散射场仅是在和t时刻的散射场的一种平移。因此,散射场的互相函数值也是其自相关函数值的一种平移,即
在这里,示踪粒子运动位移为运动时间为τ,则运动速度为
至此,每一块区域t1时刻的速度已计算出。
步骤四、按步骤二和步骤三依次处理t2和t3时刻的图像、t3和t4时刻的图像、t5和t6时刻的图像,……,即可得到t2,t3,t4,……的速度分布。

Claims (2)

1.一种采用基于近场散射的流动二维速度场测量装置测量流动二维速度场的方法,基于近场散射的流动二维速度场测量装置包括:
一激光器,用于发出相干光束;
一空间滤波器,用于滤除由激光器产生的杂散光,得到高斯光束;
一准直透镜,用于将激光器发出的光束转换成直光束;
一透镜组,用于设定焦距,将距离测量区域Z处的散斑放大成像;
一CCD相机,用于设定时间间隔采集所述透镜组成像的近场散斑图像;
一计算机,对所述CCD相机采集的近场散斑图像进行计算得到流体的流场分布;
其特征在于,包括如下步骤:
步骤一、激光器发出的相干光束照射含有示踪粒子的流体,调整透镜组焦距,使得CCD相机采集距离管道中心Z处的图像,并设定CCD采样时间间隔τ,通过CCD相机获取透射光和散射光干涉后光强分布的散斑图像:
<mrow> <mi>f</mi> <mrow> <mo>(</mo> <mover> <mi>r</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mrow> <mo>|</mo> <msub> <mi>e</mi> <mn>0</mn> </msub> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>2</mn> <mi>Re</mi> <mo>{</mo> <msubsup> <mi>e</mi> <mn>0</mn> <mo>*</mo> </msubsup> <msub> <mi>e</mi> <mi>s</mi> </msub> <mrow> <mo>(</mo> <mover> <mi>r</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>}</mo> </mrow>
其中,表示散斑图像,其中是对应的成像平面的位置向量,t是对应的成像时刻,e0是透射光场,|e0|2是静态透射光强,es是散射光场,Re表示取实部运算,*号是取共轭运算,其中距离Z的选取必须满足D*=2Z·NA<D,NA是透镜的数值孔径,D是相干光源的光束尺寸;
步骤二、CCD相机采集的N帧散斑图像相邻两组图像采样的时间间隔为τ,对相邻两组图像散斑强度值作差,得到N-1组差分信号图像
步骤三、将所得到的每组差分信号图像划分成S个小的计算区域,计算每个小的子区域内的平均速度,平均速度的计算方法是:
A、根据差分信号图像的自相关函数得到该时间间隔内示踪粒子的运动位移其中,差分信号的自相关函数为:
<mrow> <msub> <mi>g</mi> <mi>&amp;tau;</mi> </msub> <mrow> <mo>(</mo> <mover> <mi>x</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>)</mo> </mrow> <mo>=</mo> <mn>4</mn> <msup> <mrow> <mo>|</mo> <msub> <mi>e</mi> <mn>0</mn> </msub> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mi>Re</mi> <mo>{</mo> <mn>2</mn> <mi>G</mi> <mrow> <mo>(</mo> <mover> <mi>x</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>-</mo> <mi>G</mi> <mrow> <mo>(</mo> <mover> <mi>x</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>G</mi> <mrow> <mo>(</mo> <mo>-</mo> <mover> <mi>x</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
其中,表示示踪粒子的运动位移,第一项是正值,表示的是散射场的自相关函数,它正比于测量区域内示踪粒子的数量,均为散射场的互相关函数;
的定义为:
<mrow> <mi>G</mi> <mrow> <mo>(</mo> <mo>&amp;PlusMinus;</mo> <mover> <mi>x</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>&lt;</mo> <msub> <mi>e</mi> <mi>s</mi> </msub> <mrow> <mo>(</mo> <mover> <mi>r</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <msubsup> <mi>e</mi> <mi>s</mi> <mo>*</mo> </msubsup> <mrow> <mo>(</mo> <mover> <mi>r</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;PlusMinus;</mo> <mover> <mi>x</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mi>t</mi> <mo>+</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>&gt;</mo> </mrow>
其中,<…>表示总体平均运算符;
散射场的互相关函数值是其自相关函数值的一种平移,即
<mrow> <mi>G</mi> <mrow> <mo>(</mo> <mover> <mi>x</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>G</mi> <mrow> <mo>(</mo> <mover> <mi>x</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>-</mo> <mi>&amp;Delta;</mi> <mover> <mi>r</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>,</mo> <mn>0</mn> <mo>)</mo> </mrow> </mrow>
B、根据步骤A得到的示踪粒子的运动位移计算小的子区域内示踪粒子的运动速度为
<mrow> <mo>&lt;</mo> <mover> <mi>v</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&gt;</mo> <mo>=</mo> <mi>&amp;Delta;</mi> <mover> <mi>r</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>/</mo> <mi>&amp;tau;</mi> <mo>.</mo> </mrow>
2.根据权利要求1所述的方法,其特征在于:所述的激光器为连续发光的He-Ne激光器。
CN201510118897.2A 2015-03-18 2015-03-18 一种基于近场散射的流动二维速度场测量装置及方法 Active CN104698219B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510118897.2A CN104698219B (zh) 2015-03-18 2015-03-18 一种基于近场散射的流动二维速度场测量装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510118897.2A CN104698219B (zh) 2015-03-18 2015-03-18 一种基于近场散射的流动二维速度场测量装置及方法

Publications (2)

Publication Number Publication Date
CN104698219A CN104698219A (zh) 2015-06-10
CN104698219B true CN104698219B (zh) 2018-01-02

Family

ID=53345554

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510118897.2A Active CN104698219B (zh) 2015-03-18 2015-03-18 一种基于近场散射的流动二维速度场测量装置及方法

Country Status (1)

Country Link
CN (1) CN104698219B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108918360A (zh) * 2018-07-27 2018-11-30 国家电网有限公司 泥沙颗粒沉积模拟三维试验装置和方法

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105222828B (zh) * 2015-09-30 2017-05-24 东南大学 一种贴壁射流速度场和浓度场的同步测量装置与方法
CA3006467C (en) * 2015-12-11 2022-06-21 University Of Helsinki Properties of a surface and subsurface structures with white light interferometry using photonic jets
CN105890671B (zh) * 2016-06-15 2018-04-24 浙江大学 用于高温环镜颗粒场三维测量的探针
CN106645777B (zh) * 2016-11-10 2023-10-27 华东师范大学 一种基于数字微镜的空间滤波测速装置
CN108663542B (zh) * 2018-05-15 2020-03-10 中国空气动力研究与发展中心低速空气动力研究所 一种高精度piv基准速度场装置
CN108918361B (zh) * 2018-07-27 2020-09-11 国家电网有限公司 推移质沉积模拟试验装置和方法
CN110672559B (zh) * 2019-10-25 2021-01-19 西安交通大学 一种同时测量二元系热扩散率和互扩散系数的装置及方法
CN117348266B (zh) * 2023-12-06 2024-01-30 中国航空工业集团公司沈阳空气动力研究所 一种用于piv系统的激光散斑生成装置与方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101408555A (zh) * 2007-10-08 2009-04-15 南京理工大学 激光散斑测试系统
CN102253674A (zh) * 2011-04-25 2011-11-23 武汉工业学院 一种智能稳流装置和方法
JP5932626B2 (ja) * 2012-03-16 2016-06-08 本田技研工業株式会社 3次元空間の粒子画像流速測定装置
CN103472256B (zh) * 2013-09-25 2015-09-16 东南大学 基于面阵ccd空间滤波器的流动二维速度场测量方法及装置
CN204514952U (zh) * 2015-03-18 2015-07-29 东南大学 一种基于近场散射的流动二维速度场测量装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108918360A (zh) * 2018-07-27 2018-11-30 国家电网有限公司 泥沙颗粒沉积模拟三维试验装置和方法
CN108918360B (zh) * 2018-07-27 2020-09-11 国家电网有限公司 泥沙颗粒沉积模拟三维试验装置和方法

Also Published As

Publication number Publication date
CN104698219A (zh) 2015-06-10

Similar Documents

Publication Publication Date Title
CN104698219B (zh) 一种基于近场散射的流动二维速度场测量装置及方法
JP5364203B2 (ja) 観察装置
JP6596498B2 (ja) 流体試料中のナノ粒子を光学的に検出する方法および装置
US8830476B2 (en) Methods and apparatuses for contact-free holographic imaging of aerosol particles
JP6549747B2 (ja) 粒子測定装置および粒子測定方法
US7911590B2 (en) Method and system for measuring the speed of blood flow
US20060175561A1 (en) Particle shadow velocimetry
EP2387708A1 (en) Automated real-time particle characterization and three-dimensional velocimetry with holographic video microscopy
Hyun et al. Assessment of PIV to measure mean velocity and turbulence in open-channel flow
CN103674791A (zh) 一种基于双光束照射的干涉粒子成像测量方法
Meier et al. Imaging laser Doppler velocimetry
CN103645341A (zh) 全流场3d可视化测速方法
CN104697906A (zh) 一种基于近场散射的颗粒粒度测量装置及方法
WO2018190162A1 (ja) 粒子測定装置および粒子測定方法
GB2589012A (en) Particulate observation device and particulate observation method
Wu et al. Direct particle depth displacement measurement in DHPTV using spatial correlation of focus metric curves
CN204514952U (zh) 一种基于近场散射的流动二维速度场测量装置
JP5557114B2 (ja) 空間フィルタ流速計および流速測定方法
Kuok et al. Application of particle image velocimetry (PIV) for measuring water velocity in laboratory sedimentation tank
Hain et al. Principles of a volumetric velocity measurement technique based on optical aberrations
Soria et al. 4D digital holographic PIV/PTV with 3D volume deconvolution and predictive inverse reconstruction
Yamamoto et al. 3D particle measurements by single beam two-views magnified digital in-line holography
Potenza et al. A new particle sizing technique based on near field scattering
Melzer et al. Diagnostic Methods in Dusty Plasmas
Mugikura Application of defocusing technique to bubble depth measurement

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant