CN102494622B - 一种水滴直径测量方法 - Google Patents
一种水滴直径测量方法 Download PDFInfo
- 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
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上衰减后的光强为根据区域A内所有阵列元素的坐标位置,选择间距最大的两个阵列元素和且得到和之间的距离L;
步骤4:获得被测水滴在光强感应阵列上投影的最大宽度W;
步骤5:获得区域A中光强衰减最大的阵列元素
根据步骤2中获得的无水滴情况下光强感应阵列上每个阵列元素的光照强度,得到区域A中无水滴时阵列元素apq的对应的光照强度则在区域A中,光强衰减最大的阵列元素为中的最大值所对应的阵列元素,且阵列元素在无水滴时光照强度为穿过水滴衰减后阵列元素的光照强度
步骤6:测得被测水滴直径;
根据
可得
则 (3)
上式中,h为均匀水滴的厚度,I为初始入射光穿过厚度为h的均匀水滴衰减后的光强,I0为初始入射光光强,c1为衰减系数。
由此可以得出被测水滴的厚度H,即:
则:
根据步骤4得到的W与步骤5得到的L,通过椭球形的体积公式,得到:
令被测水滴的直径为D,则:
令V球=V椭球,可得
本发明的优点在于:
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上衰减后的光强为根据区域A内所有阵列元素的坐标位置,选择间距最大的两个阵列元素和且得到和之间的距离L。
步骤4:获得被测水滴5在光强感应阵列4上投影的最大宽度W;
根据步骤2中获得的无水滴情况下光强感应阵列4上每个阵列元素的光照强度的精准值,得到区域A中无水滴时阵列元素apq的对应的光照强度则在区域A中,光强衰减最大的阵列元素为中的最大值所对应的阵列元素,且阵列元素在无水滴时光照强度为穿过水滴衰减后阵列元素的光照强度
步骤6:测得被测水滴5直径;
根据
可得
则 (3)
上式中,h为均匀水滴的厚度,I为初始入射光穿过厚度为h的均匀水滴衰减后的光强,I0为初始入射光光强,c1为衰减系数。
由此可以得出被测水滴5的厚度H,即:
在本发明中在光强感应器3的光强感应阵列4上设有厚度一定的透明膜,通过透明膜可以让光线穿过,且防止水滴不慎掉在光强感应阵列4上。光照通过透明膜后同样会产生一定衰减,令c2为透明膜的衰减系数,透明膜厚度为h′,则根据公式(2)可得:
令 可得:
其中I0透明膜为透明膜上的入射光强,I透明膜为经过透明膜后衰减的光强。
由于光照在透明膜上各处对光强的衰减作用相同,由此在可以得出被测水滴5的厚度H。即:
则:
根据步骤4得到的W与步骤5得到的L,通过椭球形的体积公式,得到:
令被测水滴5的直径为D,则:
令V球=V椭球,可得
由此通过上述方法便可以获得更加精确的被测水滴5直径,为观测水滴大小对结冰的影响,以及观测雨水的大小对农作物生长的影响等领域的研究提供更精确的数据。
Claims (2)
1.一种水滴直径测量方法,其特征在于:通过以下步骤来完成:
步骤1:设备调试;
开启光源与光强感应器,使光源发出的光照照向凹面镜,通过调整凹面镜角度,使光照通过凹面镜平行射出,并照射到光强感应器中的光强感应阵列上;被测水滴由凹面镜与光强感应器间经过;
步骤2:测得无水滴时光强感应阵列上各个阵列元素的光照强度;
每间隔30秒,通过光强感应器测得光强感应阵列上各个阵列元素aij的光照强度Iij,通过分析处理装置接收三次光照强度Iij,取平均值为无水滴情况下各个阵列元素的光照强度;其中,i为光强感应阵列横坐标,j为光强感应阵列纵坐标;
步骤3:获得被测水滴在光强感应阵列上投影的最大长度L;
令被测水滴在光强感应阵列上投影为区域A,则区域A上的阵列元素apq上的光强发生衰减,通过光强感应器测得区域A上的阵列元素apq上衰减后的光强为根据区域A内所有阵列元素的坐标位置,选择间距最大的两个阵列元素和且得到和之间的距离L;
步骤4:获得被测水滴在光强感应阵列上投影的最大宽度W;
根据步骤2中获得的无水滴情况下光强感应阵列上每个阵列元素的光照强度,得到区域A中无水滴时阵列元素apq的对应的光照强度则在区域A中,光强衰减最大的阵列元素为中的最大值所对应的阵列元素,且阵列元素在无水滴时光照强度为穿过水滴衰减后阵列元素的光照强度
步骤6:测得被测水滴直径;
根据
可得
则 (3)
上式中,h为均匀水滴的厚度,I为初始入射光穿过厚度为h的均匀水滴衰减后的光强,I0为初始入射光光强,c1为衰减系数;
由此可以得出被测水滴的厚度H,即:
则:
根据步骤4得到的W与步骤5得到的L,通过椭球形的体积公式,得到:
令被测水滴的直径为D,则:
令V球=V椭球,可得
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)
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)
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)
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 | パイプ外径測定装置 |
-
2011
- 2011-12-09 CN CN 201110409832 patent/CN102494622B/zh not_active Expired - Fee Related
Patent Citations (3)
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)
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 |