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

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

Info

Publication number
CN102494622B
CN102494622B CN 201110409832 CN201110409832A CN102494622B CN 102494622 B CN102494622 B CN 102494622B CN 201110409832 CN201110409832 CN 201110409832 CN 201110409832 A CN201110409832 A CN 201110409832A CN 102494622 B CN102494622 B CN 102494622B
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.)
Expired - Fee Related
Application number
CN 201110409832
Other languages
English (en)
Other versions
CN102494622A (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

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

Abstract

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

Description

一种水滴直径测量方法
技术领域
本专利涉及到光学工程、数据信号处理、软件工程等领域,具体来说是一种测量速度较低,相互间隔较远,不会相互影响的低速稀疏水滴直径的测量方法。
背景技术
测量水滴直径的方法主要有:色斑法、面粉法、激光法、摄影法、浸入法等。色斑法和摄影法是应用比较广泛的水滴直径测量方法,这两种方法都是利用到水滴在同一材料上形成的色斑大小与水滴的直径大小成正比的假定,预先率定好水滴直径与色斑直径的关系,然后通过量水滴在滤纸、晒图纸、相纸等材料上形成的色斑大小推断出水滴的直径,或者是摄影图像的成像计算出水滴的直径。这两种方法都有相同的部分,就是根据水滴的周向判断水滴的直径,没有考虑厚度方向。而实际情况是水滴在空气中运动时,由于受到空气阻力的作用,往往是成椭球形,为了提高精度就不能完全按球形考虑水滴的直径。本专利就提出一种按椭球计算水滴的体积的模型,并换算成就球形计算出水滴的当量直径的方法。
发明内容
为了解决上述问题,本发明提出一种测量速度较低,相互间隔较远,不会相互影响的低速稀疏水滴直径的测量方法,以椭球形模型对水滴进行分析,得到更加符合实际情况的高精度水滴直径。
本发明通过以下步骤来完成:
步骤1:设备调试;
开启光源与光强感应器,使光源发出的光照照向凹面镜,通过调整凹面镜角度,使光照通过凹面镜平行射出,并照射到光强感应器中的光强感应阵列上。
步骤2:测得无水滴时光强感应阵列上各个阵列元素的光照强度;
每间隔30秒,通过光强感应器测得光强感应阵列上各个阵列元素aij的光照强度Iij,通过分析处理装置接收三次光照强度Iij,取平均值为无水滴情况下各个阵列元素的光照强度;其中,i为光强感应阵列横坐标,j为光强感应阵列纵坐标。
步骤3:获得被测水滴在光强感应阵列上投影的最大长度L;
令被测水滴在光强感应阵列上投影为区域A,则区域A上的阵列元素apq上的光强发生衰减,通过光强感应器测得区域A上的阵列元素apq上衰减后的光强为
Figure GDA00002976392900022
根据区域A内所有阵列元素的坐标位置,选择间距最大的两个阵列元素
Figure GDA00002976392900023
Figure GDA00002976392900024
且得到
Figure GDA00002976392900025
Figure GDA00002976392900026
之间的距离L;
步骤4:获得被测水滴在光强感应阵列上投影的最大宽度W;
在区域A内选择两个阵列元素
Figure GDA00002976392900028
使
Figure GDA00002976392900029
Figure GDA000029763929000210
连线与
Figure GDA000029763929000212
连线相互垂直,且
Figure GDA000029763929000213
Figure GDA000029763929000214
之间的距离W为
Figure GDA000029763929000215
Figure GDA000029763929000216
垂直方向上的最大距离。
步骤5:获得区域A中光强衰减最大的阵列元素
根据步骤2中获得的无水滴情况下光强感应阵列上每个阵列元素的光照强度,得到区域A中无水滴时阵列元素apq的对应的光照强度
Figure GDA000029763929000218
则在区域A中,光强衰减最大的阵列元素
Figure GDA000029763929000219
Figure GDA000029763929000220
中的最大值所对应的阵列元素,且阵列元素在无水滴时光照强度为
Figure GDA000029763929000222
穿过水滴衰减后阵列元素
Figure GDA000029763929000223
的光照强度
Figure GDA00002976392900038
步骤6:测得被测水滴直径;
根据
dI I = - c 1 dh            (1)
可得
I = I 0 e - c 1 h       (2)
h = - 1 / c 1 ln I I 0          (3)
上式中,h为均匀水滴的厚度,I为初始入射光穿过厚度为h的均匀水滴衰减后的光强,I0为初始入射光光强,c1为衰减系数。
由此可以得出被测水滴的厚度H,即:
则:
H = - 1 c 1 ln I a p m q m ′ I ‾ a p m q m
根据步骤4得到的W与步骤5得到的L,通过椭球形的体积公式,得到:
令被测水滴的直径为D,则:
Figure GDA00002976392900036
令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,取平均值
Figure GDA00002976392900041
为无水滴情况下各个阵列元素的光照强度;其中,i为光强感应阵列4横坐标,j为光强感应阵列4纵坐标,如图3所示。
步骤3:获得被测水滴5在光强感应阵列4上投影的最大长度L;
被测水滴5由凹面镜2与光强感应器3间经过时,入射光穿过水滴后光强衰减,如图3所示。衰减后的入射光使光强感应阵列4上产生被测水滴5的投影,如图4所示,令投影为区域A,则区域A上的阵列元素apq上的光强发生衰减,通过光强感应器3测得区域A上的阵列元素apq上衰减后的光强为
Figure GDA00002976392900051
根据区域A内所有阵列元素的坐标位置,选择间距最大的两个阵列元素
Figure GDA00002976392900052
且得到
Figure GDA00002976392900054
Figure GDA00002976392900055
之间的距离L。
步骤4:获得被测水滴5在光强感应阵列4上投影的最大宽度W;
如图3所示,在区域A内选择两个阵列元素
Figure GDA00002976392900056
Figure GDA00002976392900057
使
Figure GDA00002976392900058
Figure GDA00002976392900059
连线与
Figure GDA000029763929000510
Figure GDA000029763929000511
连线相互垂直,且
Figure GDA000029763929000512
Figure GDA000029763929000513
之间的距离W为
Figure GDA000029763929000514
Figure GDA000029763929000515
垂直方向上的最大距离。
步骤5:获得区域A中光强衰减最大的阵列元素
Figure GDA000029763929000516
根据步骤2中获得的无水滴情况下光强感应阵列4上每个阵列元素的光照强度的精准值,得到区域A中无水滴时阵列元素apq的对应的光照强度则在区域A中,光强衰减最大的阵列元素
Figure GDA000029763929000518
Figure GDA000029763929000519
中的最大值所对应的阵列元素,且阵列元素在无水滴时光照强度为穿过水滴衰减后阵列元素
Figure GDA000029763929000522
的光照强度
Figure GDA000029763929000523
步骤6:测得被测水滴5直径;
根据
dI I = - c 1 dh        (1)
可得
I = I 0 e - c 1 h           (2)
h = - 1 / c 1 ln I I 0          (3)
上式中,h为均匀水滴的厚度,I为初始入射光穿过厚度为h的均匀水滴衰减后的光强,I0为初始入射光光强,c1为衰减系数。
由此可以得出被测水滴5的厚度H,即:
H = - 1 c 1 ln I a p m q m ′ I ‾ a p m q m .
在本发明中在光强感应器3的光强感应阵列4上设有厚度一定的透明膜,通过透明膜可以让光线穿过,且防止水滴不慎掉在光强感应阵列4上。光照通过透明膜后同样会产生一定衰减,令c2为透明膜的衰减系数,透明膜厚度为h′,则根据公式(2)可得:
k = e - c 2 h ′ , 可得:
Figure GDA00002976392900064
其中I0透明膜为透明膜上的入射光强,I透明膜为经过透明膜后衰减的光强。
由于光照在透明膜上各处对光强的衰减作用相同,由此在可以得出被测水滴5的厚度H。即:
H = 1 c 1 ln I a p m q m ′ k I ‾ a p m q m k
则:
H = - 1 c 1 ln I a p m q m ′ I ‾ a p m q m
根据步骤4得到的W与步骤5得到的L,通过椭球形的体积公式,得到:
Figure GDA00002976392900067
令被测水滴5的直径为D,则:
令V=V椭球,可得
D = LWH 3
由此通过上述方法便可以获得更加精确的被测水滴5直径,为观测水滴大小对结冰的影响,以及观测雨水的大小对农作物生长的影响等领域的研究提供更精确的数据。

Claims (2)

1.一种水滴直径测量方法,其特征在于:通过以下步骤来完成:
步骤1:设备调试;
开启光源与光强感应器,使光源发出的光照照向凹面镜,通过调整凹面镜角度,使光照通过凹面镜平行射出,并照射到光强感应器中的光强感应阵列上;被测水滴由凹面镜与光强感应器间经过;
步骤2:测得无水滴时光强感应阵列上各个阵列元素的光照强度;
每间隔30秒,通过光强感应器测得光强感应阵列上各个阵列元素aij的光照强度Iij,通过分析处理装置接收三次光照强度Iij,取平均值
Figure FDA00002976392800011
为无水滴情况下各个阵列元素的光照强度;其中,i为光强感应阵列横坐标,j为光强感应阵列纵坐标;
步骤3:获得被测水滴在光强感应阵列上投影的最大长度L;
令被测水滴在光强感应阵列上投影为区域A,则区域A上的阵列元素apq上的光强发生衰减,通过光强感应器测得区域A上的阵列元素apq上衰减后的光强为
Figure FDA00002976392800012
根据区域A内所有阵列元素的坐标位置,选择间距最大的两个阵列元素
Figure FDA00002976392800013
且得到之间的距离L;
步骤4:获得被测水滴在光强感应阵列上投影的最大宽度W;
在区域A内选择两个阵列元素
Figure FDA00002976392800016
Figure FDA00002976392800017
使
Figure FDA00002976392800018
Figure FDA00002976392800019
连线与
Figure FDA000029763928000110
Figure FDA000029763928000111
连线相互垂直,且
Figure FDA000029763928000113
之间的距离W为
Figure FDA000029763928000114
Figure FDA000029763928000115
垂直方向上的最大距离;
步骤5:获得区域A中光强衰减最大的阵列元素
Figure FDA000029763928000116
根据步骤2中获得的无水滴情况下光强感应阵列上每个阵列元素的光照强度,得到区域A中无水滴时阵列元素apq的对应的光照强度
Figure FDA00002976392800021
则在区域A中,光强衰减最大的阵列元素
Figure FDA00002976392800022
Figure FDA00002976392800023
中的最大值所对应的阵列元素,且阵列元素
Figure FDA00002976392800024
在无水滴时光照强度为
Figure FDA000029763928000214
穿过水滴衰减后阵列元素
Figure FDA00002976392800025
的光照强度
Figure FDA00002976392800026
步骤6:测得被测水滴直径;
根据
dI I = - c 1 dh             (1)
可得
I = I 0 e - c 1 h               (2)
h = - 1 / c 1 ln I I 0              (3)
上式中,h为均匀水滴的厚度,I为初始入射光穿过厚度为h的均匀水滴衰减后的光强,I0为初始入射光光强,c1为衰减系数;
由此可以得出被测水滴的厚度H,即:
则:
H = - 1 c 1 ln I a p m q m ′ I ‾ a p m q m
根据步骤4得到的W与步骤5得到的L,通过椭球形的体积公式,得到:
Figure FDA000029763928000211
令被测水滴的直径为D,则:
Figure FDA000029763928000212
令V=V椭球,可得
D = LWH 3 .
2.如权利要求1所述一种水滴直径测量方法,其特征在于:所述光强感应器的光强感应阵列上设置有厚度均匀的透明膜,令c2为透明膜的衰减系数,透明膜厚度为h′,则:
Figure FDA00002976392800031
k = e - c 2 h ′ , 可得:
Figure FDA00002976392800032
其中I0透明膜为透明膜上的入射光强,I透明膜为经过透明膜后衰减的光强;
由此在光强感应阵列上具有透明膜时,被测水滴的厚度为:
H = 1 c 1 ln I a p m q m ′ k I ‾ a p m q m 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 CN102494622A (zh) 2012-06-13
CN102494622B true 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)

Families Citing this family (2)

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

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN101819146A (zh) * 2010-06-02 2010-09-01 长沙理工大学 非接触式汽轮机蒸汽湿度和水滴直径的测量方法及其装置

Family Cites Families (3)

* 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 降雨粒径分布計測方法
JPH08136227A (ja) * 1994-11-02 1996-05-31 Nippon Steel Corp パイプ外径測定装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN101819146A (zh) * 2010-06-02 2010-09-01 长沙理工大学 非接触式汽轮机蒸汽湿度和水滴直径的测量方法及其装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
JP昭60-78332A 1985.05.04
JP昭62-240840A 1987.10.21
JP特开平8-136227A 1996.05.31
基于CCD的蒸汽湿度及水滴直径测量;黄竹青等;《光电工程》;20100430;第37卷(第4期);39-43 *
黄竹青等.基于CCD的蒸汽湿度及水滴直径测量.《光电工程》.2010,第37卷(第4期),

Also Published As

Publication number Publication date
CN102494622A (zh) 2012-06-13

Similar Documents

Publication Publication Date Title
Ni et al. Lagrangian acceleration measurements in convective thermal turbulence
Su et al. Planar measurements of the full three-dimensional scalar dissipation rate in gas-phase turbulent flows
Bosbach et al. Large scale particle image velocimetry with helium filled soap bubbles
CN102818756B (zh) 基于激光能量陷阱法的pm2.5颗粒的测定方法及装置
CN102706529B (zh) 超声速流场密度场的校准以及测量超声速密度场的方法
CN107014741A (zh) 流式细胞仪
Atkins Velocity field measurement using particle image velocimetry (PIV)
US20050018882A1 (en) Controlled surface wave image velocimetry
CN103398925A (zh) 一种利用双光谱成像进行雾场粒径平面分布测量的装置
CN108169510A (zh) 基于单光场相机的微尺度流动三维速度场测量装置和方法
CN102494622B (zh) 一种水滴直径测量方法
CN103808958A (zh) 一种测量明渠紊流粘性底层流速的方法及装置
Gibeau et al. A modular, 3D-printed helium-filled soap bubble generator for large-scale volumetric flow measurements
Oliveira et al. Lagrangian and Eulerian statistics of pipe flows measured with 3D-PTV at moderate and high Reynolds numbers
Chang et al. Turbulent Prandtl number in neutrally buoyant turbulent round jet
Bissell et al. An approach to spray characterization by combination of measurement techniques
KR102136843B1 (ko) 무인항공기를 이용한 대기 중 메탄가스 배출량 측정방법
US20240193794A1 (en) Method and Apparatus for Tracking Motion of Objects in Three-Dimensional Space
CN104885988A (zh) 测量鱼类嗅觉灵敏度的装置及方法
CN104165829A (zh) 单分散微液滴直径测量装置及测量方法
CN106841076A (zh) 基于车载双光路差分吸收光谱的源清单获取方法与装置
CN208847854U (zh) 一种激光雨滴谱仪的标定装置
CN108489547B (zh) 一种雨滴参数测试装置
Obligado et al. Constrained dynamics of an inertial particle in a turbulent flow
Raddatz et al. Size selection of sub-and super-micron clay mineral kaolinite particles using a custom-built Maxi-DMA

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