CN102494622A - 一种水滴直径的测量方法 - Google Patents

一种水滴直径的测量方法 Download PDF

Info

Publication number
CN102494622A
CN102494622A CN2011104098325A CN201110409832A CN102494622A CN 102494622 A CN102494622 A CN 102494622A CN 2011104098325 A CN2011104098325 A CN 2011104098325A CN 201110409832 A CN201110409832 A CN 201110409832A CN 102494622 A CN102494622 A CN 102494622A
Authority
CN
China
Prior art keywords
light intensity
water droplet
array element
water drop
intensity
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
Application number
CN2011104098325A
Other languages
English (en)
Other versions
CN102494622B (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN 201110409832 priority Critical patent/CN102494622B/zh
Publication of CN102494622A publication Critical patent/CN102494622A/zh
Application granted granted Critical
Publication of CN102494622B publication Critical patent/CN102494622B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Length Measuring Devices By Optical Means (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开一种水滴直径的测量方法,根据水滴对光的折射和衰减,得到获得被测水滴在光强感应阵列上投影的最大长度L以及最大宽度W;获得被测水滴在光强感应阵列上投影中光强衰减最大的阵列元素;根据均匀水滴的厚度计算公式得到被测水滴厚度H。由于水滴在空气中受到大气阻力的作用,会使得水滴成椭球形,由此根据L、W、H得到椭球型被测水滴的体积,将椭球型被测水滴体积换算成球形体积,由此得到被测水滴的直径。本发明的优点在于:方法简单,水滴模型更接近于现实生活中的各式各样的水滴,测量更加精确。本方法可以应用于观测水滴大小对结冰的影响,以及观测雨水的大小对农作物生长的影响等领域,为这些领域的研究提供更精确的数据。

Description

一种水滴直径的测量方法
技术领域
本专利涉及到光学工程、数据信号处理、软件工程等领域,具体来说是一种测量速度较低,相互间隔较远,不会相互影响的低速稀疏水滴直径的测量方法。 
背景技术
测量水滴直径的方法主要有:色斑法、面粉法、激光法、摄影法、浸入法等。色斑法和摄影法是应用比较广泛的水滴直径测量方法,这两种方法都是利用到水滴在同一材料上形成的色斑大小与水滴的直径大小成正比的假定,预先率定好水滴直径与色斑直径的关系,然后通过量水滴在滤纸、晒图纸、相纸等材料上形成的色斑大小推断出水滴的直径,或者是摄影图像的成像计算出水滴的直径。这两种方法都有相同的部分,就是根据水滴的周向判断水滴的直径,没有考虑厚度方向。而实际情况是水滴在空气中运动时,由于受到空气阻力的作用,往往是成椭球形,为了提高精度就不能完全按球形考虑水滴的直径。本专利就提出一种按椭球计算水滴的体积的模型,并换算成就球形计算出水滴的当量直径的方法。 
发明内容
为了解决上述问题,本发明提出一种测量速度较低,相互间隔较远,不会相互影响的低速稀疏水滴直径的测量方法,以椭球形模型对水滴进行分析,得到更加符合实际情况的高精度水滴直径。 
本发明通过以下步骤来完成: 
步骤1:设备调试; 
开启光源与光强感应器,使光源发出的光照照向凹面镜,通过调整凹面镜角度,使光照通过凹面镜平行射出,并照射到光强感应器中的光强感应阵列上。 
步骤2:测得无水滴时光强感应阵列上各个阵列元素的光照强度; 
每间隔30秒,通过光强感应器测得光强感应阵列上各个阵列元素aij的光照强度Iij,通过分析处理装置接收三次光照强度Iij,取平均值 
Figure BDA0000118230580000021
为无水滴情况下各个阵列元素的光照强度;其中,i为光强感应阵列横坐标,j为光强感应阵列纵坐标。 
步骤3:获得被测水滴在光强感应阵列上投影的最大长度L; 
令被测水滴在光强感应阵列上投影为区域A,则区域A上的阵列元素apq上的光强发生衰减,通过光强感应器测得区域A上的阵列元素apq上衰减后的光强为 
Figure BDA0000118230580000022
根据区域A内所有阵列元素的坐标位置,选择间距最大的两个阵列元素 
Figure BDA0000118230580000023
和 
Figure BDA0000118230580000024
且得到 
Figure BDA0000118230580000025
和 
Figure BDA0000118230580000026
之间的距离L; 
步骤4:获得被测水滴在光强感应阵列上投影的最大宽度W; 
在区域A内选择两个阵列元素 
Figure BDA0000118230580000027
和 
Figure BDA0000118230580000028
使 
Figure BDA0000118230580000029
和 连线与 
Figure BDA00001182305800000211
和 
Figure BDA00001182305800000212
连线相互垂直,且 
Figure BDA00001182305800000213
和 
Figure BDA00001182305800000214
之间的距离W为 
Figure BDA00001182305800000215
和 
Figure BDA00001182305800000216
垂直方向上的最大距离。 
步骤5:获得区域A中光强衰减最大的阵列元素 
Figure BDA00001182305800000217
根据步骤3中获得的无水滴情况下光强感应阵列上每个阵列元素的光照强度,得到区域A中无水滴时阵列元素apq的对应的光照强度精准值 
Figure BDA00001182305800000218
则在区域A中,光强衰减最大的阵列元素 
Figure BDA00001182305800000219
为 
Figure BDA00001182305800000220
中的最大值所对应的阵列元素,且阵列元素 
Figure BDA00001182305800000221
在无水滴时光照强度精准值为 穿过水滴衰减后阵 列元素 
Figure BDA0000118230580000031
的光照强度 
步骤6:测得被测水滴直径; 
根据 
dI I = - c 1 dh - - - ( 1 )
可得 
I = I 0 e - c 1 h - - - ( 2 )
h = - 1 / c 1 ln I I O - - - ( 3 )
上式中,h为均匀水滴的厚度,I为初始入射光穿过厚度为h的均匀水滴衰减后的光强,I0为初始入射光光强,c1为衰减系数。 
由此可以得出被测水滴的厚度H,即: 
则: 
H = - 1 c 1 ln I a pmqm ′ I ‾ a pmqm
根据步骤4得到的W与步骤5得到的L,通过椭球形的体积公式,得到: 
Figure BDA0000118230580000037
令被测水滴的直径为D,则: 
Figure BDA0000118230580000038
令V=V椭球,可得 
D = LWH 3 .
本发明的优点在于: 
1、本发明水滴直径测量方法,操作简单,水滴模型更接近于现实生活中的各式各样的水滴,测量结果更加精确。 
2、本方法可以应用于观测水滴大小对结冰的影响,以及观测雨水的大小对 农作物生长的影响等领域,为这些领域的研究提供更精确的数据。 
附图说明
图1为本发明方法流程图; 
图2为设备调试过程中光照射入、射出示意图; 
图3为光照穿过水滴示意图; 
图4为被测水滴在光强感应阵列上投影示意图。 
图中: 
1-光源  2-凹面镜  3-光强感应器  4-光强感应阵列  5-被测水滴 
具体实施方式
下面通过附图对本发明做进一步说明。 
本发明一种稀疏水滴直径测量方法,如图1所示,通过以下步骤来完成: 
步骤1:设备调试; 
开启光源1与光强感应器3,使由光源1发出的光照照向凹面镜2,通过调整凹面镜2角度,使光照通过凹面镜2平行射出,并照射到光强感应器3中的光强感应阵列4上,如图2所示。其中,光强感应器3与凹面镜2间为被测水滴通道, 
步骤2:测得无水滴时光强感应阵列4上各个阵列元素的光照强度; 
由于光源1的供电电压可能不完全稳定,因此本发明中在无水滴的情况下,每间隔30秒,通过光强感应器3测得光强感应阵列4上各个阵列元素aij的光照强度Iij,通过分析处理装置接收三次光照强度Iij,取平均值 为无水滴情况下各个阵列元素的光照强度;其中,i为光强感应阵列4横坐标,j为光强感应阵列4纵坐标,如图3所示。 
步骤3:获得被测水滴5在光强感应阵列4上投影的最大长度L; 
被测水滴5由凹面镜2与光强感应器3间经过时,入射光穿过水滴后光强衰减,如图3所示。衰减后的入射光使光强感应阵列4上产生被测水滴5的投影,如图4所示,令投影为区域A,则区域A上的阵列元素apq上的光强发生衰减,通过光强感应器3测得区域A上的阵列元素apq上衰减后的光强为 
Figure BDA0000118230580000051
根据区域A内所有阵列元素的坐标位置,选择间距最大的两个阵列元素 
Figure BDA0000118230580000052
和 且得到 和 
Figure BDA0000118230580000055
之间的距离L。 
步骤4:获得被测水滴5在光强感应阵列4上投影的最大宽度W; 
如图3所示,在区域A内选择两个阵列元素 
Figure BDA0000118230580000056
和 
Figure BDA0000118230580000057
使 
Figure BDA0000118230580000058
和 
Figure BDA0000118230580000059
连线与 
Figure BDA00001182305800000510
和 连线相互垂直,且 
Figure BDA00001182305800000512
和 
Figure BDA00001182305800000513
之间的距离W为 
Figure BDA00001182305800000514
和 
Figure BDA00001182305800000515
垂直方向上的最大距离。 
步骤5:获得区域A中光强衰减最大的阵列元素 
Figure BDA00001182305800000516
根据步骤3中获得的无水滴情况下光强感应阵列4上每个阵列元素的光照强度的精准值,得到区域A中无水滴时阵列元素apq的对应的光照强度精准值 
Figure BDA00001182305800000517
则在区域A中,光强衰减最大的阵列元素 
Figure BDA00001182305800000518
为 
Figure BDA00001182305800000519
中的最大值所对应的阵列元素,且阵列元素 
Figure BDA00001182305800000520
在无水滴时光照强度精准值为 穿过水滴衰减后阵列元素 
Figure BDA00001182305800000522
的光照强度 
Figure BDA00001182305800000523
步骤6:测得被测水滴5直径; 
根据 
dI I = - c 1 dh - - - ( 1 )
可得 
I = I 0 e - c 1 h - - - ( 2 )
h = - 1 / c 1 ln I I O - - - ( 3 )
上式中,h为均匀水滴的厚度,I为初始入射光穿过厚度为h的均匀水滴衰减后的光强,I0为初始入射光光强,c1为衰减系数。 
由此可以得出被测水滴5的厚度H,即: 
H = - 1 c 1 ln I a pmqm ′ I ‾ a pmqm .
在本发明中在光强感应器3的光强感应阵列4上设有厚度一定的透明膜,通过透明膜可以让光线穿过,且防止水滴不慎掉在光强感应阵列4上。光照通过透明膜后同样会产生一定衰减,令c2为透明膜的衰减系数,透明膜厚度为h′,则根据公式(2)可得: 
Figure BDA0000118230580000063
k = e - c 2 h ′ , 可得: 
I透明膜=kI0透明膜
其中I0透明膜为透明膜上的入射光强,I透明膜为经过透明膜后衰减的光强。 
由于光照在透明膜上各处对光强的衰减作用相同,由此在可以得出被测水滴5的厚度H。即: 
H = - 1 c 1 ln I a pmqm ′ k I ‾ a pmqm k
则: 
H = - 1 c 1 ln I a pmqm ′ I ‾ a pmqm
根据步骤4得到的W与步骤5得到的L,通过椭球形的体积公式,得到: 
Figure BDA0000118230580000067
令被测水滴5的直径为D,则: 
Figure BDA0000118230580000071
令V=V椭球,可得 
D = LWH 3
由此通过上述方法便可以获得更加精确的被测水滴5直径,为观测水滴大小对结冰的影响,以及观测雨水的大小对农作物生长的影响等领域的研究提供更精确的数据。 

Claims (2)

1.本发明一种水滴直径测量方法,其特征在于:通过以下步骤来完成:
步骤1:设备调试;
开启光源与光强感应器,使光源发出的光照照向凹面镜,通过调整凹面镜角度,使光照通过凹面镜平行射出,并照射到光强感应器中的光强感应阵列上;被测水滴由凹面镜与光强感应器间经过;
步骤2:测得无水滴时光强感应阵列上各个阵列元素的光照强度;
每间隔30秒,通过光强感应器测得光强感应阵列上各个阵列元素aij的光照强度Iij,通过分析处理装置接收三次光照强度Iij,取平均值
Figure FDA0000118230570000011
为无水滴情况下各个阵列元素的光照强度;其中,i为光强感应阵列横坐标,j为光强感应阵列纵坐标;
步骤3:获得被测水滴在光强感应阵列上投影的最大长度L;
令被测水滴在光强感应阵列上投影为区域A,则区域A上的阵列元素apq上的光强发生衰减,通过光强感应器测得区域A上的阵列元素apq上衰减后的光强为根据区域A内所有阵列元素的坐标位置,选择间距最大的两个阵列元素
Figure FDA0000118230570000013
Figure FDA0000118230570000014
且得到
Figure FDA0000118230570000015
Figure FDA0000118230570000016
之间的距离L;
步骤4:获得被测水滴在光强感应阵列上投影的最大宽度W;
在区域A内选择两个阵列元素
Figure FDA0000118230570000017
Figure FDA0000118230570000018
使
Figure FDA0000118230570000019
Figure FDA00001182305700000110
连线与
Figure FDA00001182305700000111
Figure FDA00001182305700000112
连线相互垂直,且
Figure FDA00001182305700000113
Figure FDA00001182305700000114
之间的距离W为
Figure FDA00001182305700000115
Figure FDA00001182305700000116
垂直方向上的最大距离;
步骤5:获得区域A中光强衰减最大的阵列元素
Figure FDA00001182305700000117
根据步骤3中获得的无水滴情况下光强感应阵列上每个阵列元素的光照强度的精准值,得到区域A中无水滴时阵列元素apq的对应的光照强度精准值
Figure FDA0000118230570000021
则在区域A中,光强衰减最大的阵列元素
Figure FDA0000118230570000022
Figure FDA0000118230570000023
中的最大值所对应的阵列元素,且阵列元素在无水滴时光照强度为
Figure FDA0000118230570000025
穿过水滴衰减后阵列元素
Figure FDA0000118230570000026
的光照强度
Figure FDA0000118230570000027
步骤6:测得被测水滴直径;
根据
dI I = - c 1 dh - - - ( 1 )
可得
I = I 0 e - c 1 h - - - ( 2 )
h = - 1 / c 1 ln I I O - - - ( 3 )
上式中,h为均匀水滴的厚度,I为初始入射光穿过厚度为h的均匀水滴衰减后的光强,I0为初始入射光光强,c1为衰减系数;
由此可以得出被测水滴的厚度H,即:
则:
H = - 1 c 1 ln I a pmqm ′ I ‾ a pmqm
根据步骤4得到的W与步骤5得到的L,通过椭球形的体积公式,得到:
令被测水滴的直径为D,则:
令V=V椭球,可得
D = LWH 3 .
2.如权利要求1所述一种水滴直径测量方法,其特征在于:所述光强感应器的光强感应阵列上设置有厚度均匀的透明膜,令c2为透明膜的衰减系数,透明膜厚度为h′,则:
k = e - c 2 h ′ , 可得:
I透明膜=kI0透明膜
其中I0透明膜为透明膜上的入射光强,I透明膜为经过透明膜后衰减的光强;由此在光强感应阵列上具有透明膜时,被测水滴的厚度为:
H = - 1 c 1 ln I a pmqm ′ k I ‾ a pmqm k .
CN 201110409832 2011-12-09 2011-12-09 一种水滴直径测量方法 Expired - Fee Related CN102494622B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110409832 CN102494622B (zh) 2011-12-09 2011-12-09 一种水滴直径测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110409832 CN102494622B (zh) 2011-12-09 2011-12-09 一种水滴直径测量方法

Publications (2)

Publication Number Publication Date
CN102494622A true CN102494622A (zh) 2012-06-13
CN102494622B CN102494622B (zh) 2013-07-24

Family

ID=46186466

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110409832 Expired - Fee Related CN102494622B (zh) 2011-12-09 2011-12-09 一种水滴直径测量方法

Country Status (1)

Country Link
CN (1) CN102494622B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105728229A (zh) * 2016-03-01 2016-07-06 中国飞机强度研究所 一种基于喷嘴动力参数的喷雾水滴平均直径快速估算方法
CN105758320A (zh) * 2016-03-25 2016-07-13 中国水利水电科学研究院 一种公路车轮水珠雾化半径的监测系统和方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6078332A (ja) * 1983-10-05 1985-05-04 Mitsubishi Heavy Ind Ltd 水滴径測定方法
JPS62240840A (ja) * 1986-04-14 1987-10-21 Mitsubishi Heavy Ind Ltd 降雨粒径分布計測方法
US4771181A (en) * 1986-07-14 1988-09-13 Doryokuro Kakunenryo Kaihatsu Jigyodan Method for detecting dripping droplet with refracted and reflected light
EP0641599A1 (de) * 1993-09-08 1995-03-08 Roche Diagnostics GmbH Verfahren und Vorrichtung zur Dosierung von Flüssigkeiten
JPH08136227A (ja) * 1994-11-02 1996-05-31 Nippon Steel Corp パイプ外径測定装置
CN101819146A (zh) * 2010-06-02 2010-09-01 长沙理工大学 非接触式汽轮机蒸汽湿度和水滴直径的测量方法及其装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6078332A (ja) * 1983-10-05 1985-05-04 Mitsubishi Heavy Ind Ltd 水滴径測定方法
JPS62240840A (ja) * 1986-04-14 1987-10-21 Mitsubishi Heavy Ind Ltd 降雨粒径分布計測方法
US4771181A (en) * 1986-07-14 1988-09-13 Doryokuro Kakunenryo Kaihatsu Jigyodan Method for detecting dripping droplet with refracted and reflected light
EP0641599A1 (de) * 1993-09-08 1995-03-08 Roche Diagnostics GmbH Verfahren und Vorrichtung zur Dosierung von Flüssigkeiten
JPH08136227A (ja) * 1994-11-02 1996-05-31 Nippon Steel Corp パイプ外径測定装置
CN101819146A (zh) * 2010-06-02 2010-09-01 长沙理工大学 非接触式汽轮机蒸汽湿度和水滴直径的测量方法及其装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
黄竹青等: "基于CCD的蒸汽湿度及水滴直径测量", 《光电工程》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105728229A (zh) * 2016-03-01 2016-07-06 中国飞机强度研究所 一种基于喷嘴动力参数的喷雾水滴平均直径快速估算方法
CN105758320A (zh) * 2016-03-25 2016-07-13 中国水利水电科学研究院 一种公路车轮水珠雾化半径的监测系统和方法
CN105758320B (zh) * 2016-03-25 2018-11-20 中国水利水电科学研究院 一种公路车轮水珠雾化半径的监测系统和方法

Also Published As

Publication number Publication date
CN102494622B (zh) 2013-07-24

Similar Documents

Publication Publication Date Title
Florens et al. Defining the roughness sublayer and its turbulence statistics
Bucher et al. Experiments on light pulse communication and propagation through atmospheric clouds
Ni et al. Lagrangian acceleration measurements in convective thermal turbulence
CN202599852U (zh) 气液两相泡状流中气泡的识别及测量装置
Su et al. Planar measurements of the full three-dimensional scalar dissipation rate in gas-phase turbulent flows
CN102706529B (zh) 超声速流场密度场的校准以及测量超声速密度场的方法
CN107014741A (zh) 流式细胞仪
CN107941454B (zh) 基于丙酮示踪剂的高超声速plif风洞实验方法
Kupc et al. Laboratory characterization of a new nano-water-based CPC 3788 and performance comparison to an ultrafine butanol-based CPC 3776
CN103398925A (zh) 一种利用双光谱成像进行雾场粒径平面分布测量的装置
CN107036943B (zh) 一种基于作物孔隙率相似性的雾滴沉积测试方法及测试装置
Atkins Velocity field measurement using particle image velocimetry (PIV)
CN108169510B (zh) 基于单光场相机的微尺度流动三维速度场测量装置和方法
CN204988587U (zh) 一种高精度非接触气动摩擦阻力测量装置
CN103808958A (zh) 一种测量明渠紊流粘性底层流速的方法及装置
CN103413014A (zh) 一种基于二向性反射分布函数(brdf)原型反演地表反照率的方法
Cozzi et al. Development of a background-oriented schlieren technique with telecentric lenses for supersonic flow
CN102494622B (zh) 一种水滴直径测量方法
Gibeau et al. A modular, 3D-printed helium-filled soap bubble generator for large-scale volumetric flow measurements
CN103852813B (zh) 雨滴三维尺度检测装置及利用该装置计算雨滴体积的方法
Chang et al. Turbulent Prandtl number in neutrally buoyant turbulent round jet
CN106556579A (zh) 基于激光的团雾图像全天候自适应检测方法
CN102998096A (zh) 凸透镜焦距的测量方法
CN106225861A (zh) 一种液体流量的测定装置及航空器燃油油量测定方法
CN116046617A (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: 20130724

Termination date: 20131209