CN102607436B - 一种测量固体壁面波动液膜层平均厚度的方法 - Google Patents
一种测量固体壁面波动液膜层平均厚度的方法 Download PDFInfo
- Publication number
- CN102607436B CN102607436B CN201210090079.2A CN201210090079A CN102607436B CN 102607436 B CN102607436 B CN 102607436B CN 201210090079 A CN201210090079 A CN 201210090079A CN 102607436 B CN102607436 B CN 102607436B
- Authority
- CN
- China
- Prior art keywords
- liquid film
- pixel
- wall surface
- solid wall
- thickness
- 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
- Application Of Or Painting With Fluid Materials (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种测量固体壁面上波动液膜层平均厚度的方法:a)采用放大摄影***对液滴产生、开始撞击固体壁面、壁面上液膜厚度增加直到液膜厚度达到动态稳定的整个过程进行可视化拍摄;b)定位液膜厚度在原始图片高度方向上的像素范围,对选定范围内的各像素点逐一读取,得到灰度级矩阵;计算灰度级矩阵中各列上不同行间的灰度级梯度,得到各列中灰度级骤然增加的像素点,该些像素点组成的曲线是液膜表面所对应的位置;c)采用统计方法得到液膜表面以下区域的像素点总个数a1;d)计算固体壁面Y方向的像素点个数Nh;e)液膜层平均厚度为a1与Nh的比值。本发明可广泛用于各种液滴碰撞固体壁面过程的可视化观测和壁面液膜层平均厚度的定量测量。
Description
技术领域
本发明涉及一种可用于测量固体壁面波动液膜层平均厚度的方法。
背景技术
在低表面温度喷雾冷却等液滴碰撞固体壁面过程中,冲击液滴会在固体表面形成一层薄液膜,液滴主要是与薄液膜层间发生相互作用。液膜层的厚度直接影响其与壁面间的换热效果以及液膜层内产生汽泡的行为特性。由于液滴撞击薄液膜层的过程随机性非常大,涉及到多相流、液滴碰撞、相变等多个复杂环节,因此在液膜形成过程中波动非常剧烈,要对波动液膜层的厚度进行定量测量是一项难度巨大的工作。迫切需要围绕液滴碰撞热固体壁面薄液膜层等进行***深入研究,探索和发展微尺度下流动与传热过程新的实验技术与研究方法。
发明内容
本发明的目的在于提供一种测量固体壁面上波动液膜层平均厚度的方法。
为实现上述目的,本发明提供的测量固体壁面上波动液膜层平均厚度的方法,步骤如下:
a)采用放大摄影***对液滴产生、开始撞击固体壁面、壁面上液膜厚度增加直到液膜厚度达到动态稳定的整个过程进行可视化拍摄;
b)定位液膜厚度在原始图片高度方向上的像素范围,对选定范围内的各像素点逐一读取,得到灰度级矩阵;计算灰度级矩阵中各列上不同行间的灰度级梯度,得到各列中灰度级骤然增加的像素点,该些像素点组成的曲线是液膜表面所对应的位置;
c)采用统计方法得到液膜表面以下区域的像素点总个数a1;
d)计算固体壁面Y方向的像素点个数Nh;
e)液膜层平均厚度通过下式计算得到:
式中λ为单位像素对应的实际长度,单位是mm/pixel;n为图片上行像素数。
所述的测量方法,其中,步骤b中定位液膜厚度在原始图片高度方向上的像素范围是由软件CamControl完成。
所述的测量方法,其中,步骤b中对选定范围内的各像素点逐一读取,得到灰度级矩阵,是由Matlab软件完成。
所述的测量方法,其中,步骤d中固体壁面Y方向的像素点个数Nh通过下式计算得到:
式中,Nhi为任意选取不同固体壁面高度处像素点个数。
本发明提供的方法可以用于***研究液滴喷射启动时不同壁面温度下的液膜形成过程,揭示液膜达到动态稳定后流量和壁面温度对平均液膜厚度的影响规律。
附图说明
图1为本发明的液膜层平均厚度测量示意图;其中
(a)是原始图片,(b)是灰度矩阵,(c)是液膜平均厚度的测量。
图2为利用本发明方法测得的固体壁面液膜平均厚度及其变化。
具体实施方式
本发明提供的测量固体壁面上波动液膜层平均厚度的方法,步骤如下:
a)采用高速显微放大摄影***对液滴产生、开始撞击固体壁面、壁面上液膜厚度增加直到液膜厚度达到动态稳定的整个过程进行可视化拍摄;
b)对高速摄影拍摄得到的照片进行图像处理,采用灰度级梯度法捕捉到液膜表面所对应的位置。然后确定液膜表面以下区域在图片上占据的像素点个数与固体壁面在高度方向占据的像素点个数,两者的比值即为波动液膜层的平均厚度。
本发明可广泛用于各种液滴碰撞固体壁面过程的可视化观测和壁面液膜层平均厚度的定量测量。
下面结合附图对本发明作详细说明。
首先打开冷光源和高速摄影仪,调节冷光源强度与显微放大设备与被测壁面之间的距离,使得被测壁面能在高速摄影中清晰呈现,测量拍摄视野的实际尺寸,然后将流量和加热功率调到设定值,开始拍摄以及启动液滴喷射***,最后当拍摄时间足够长以至于液膜波动达到动态稳定时,关闭液滴喷射***,停止拍摄。
本实施例中摄影仪的像素以512×512为例(也可以根据需要调整为其它数值),被测壁面的温度为22℃,流量为1.39×10-5m3/s。
通过高速摄影仪自带软件CamControl定位液膜厚度在原图Y方向上的像素范围(见图1a),然后采用Matlab软件对选定范围内的各像素点逐一读取,得到如图1b所示的灰度级矩阵,图中每个方格表示一个独立的像素单位,方格中的数据为对应像素单位的灰度级值。计算灰度级矩阵中各列上不同行间的灰度级梯度,可得到各列中灰度级骤然增加的像素点,该像素点连成的曲线就认为是液膜表面所对应的位置(见图1c中曲线)。最后采用统计方法得到液膜表面以下区域的像素点总个数a1。被加热表面Y方向的像素点个数Nh可通过公式(1)得到:
式中,Nhi为图1(c)中选取的任意不同位置处的像素点个数。
最终,液膜平均厚度H(见图1c)通过公式(2)计算得到:
式中λ为单位像素对应的实际长度,单位是mm/pixel。
通过图2可以看出,在本实施例的条件下,液滴喷射约0.9秒(s)后液膜波动达到动态稳定,图2中Ⅰ所示部分是液膜平均厚度随时间呈增长趋势,液膜波动尚未达到动态稳定,Ⅱ所示部分是经过约0.9s后液膜波动达到动态稳定。
采用本发明提供的方法***研究了液滴喷射启动时不同壁面温度下的液膜形成过程,揭示了液膜达到动态稳定后流量和壁面温度对平均液膜厚度的影响规律。
Claims (4)
1.一种测量固体壁面上波动液膜层平均厚度的方法,其步骤如下:
a)采用放大摄影***对液滴产生、开始撞击固体壁面、壁面上液膜厚度增加直到液膜厚度达到动态稳定的整个过程进行可视化拍摄;
b)定位液膜厚度在原始图片高度方向上的像素范围,对选定范围内的各像素点逐一读取,得到灰度级矩阵;计算灰度级矩阵中各列上不同行间的灰度级梯度,得到各列中灰度级骤然增加的像素点,该些像素点组成的曲线是液膜表面所对应的位置;
c)采用统计方法得到液膜表面以下区域的像素点总个数a1;
d)计算固体壁面Y方向的像素点个数Nh;
e)液膜层平均厚度通过下式计算得到:
式中λ为单位像素对应的实际长度,单位是mm/pixel;n为图片上行像素数。
2.如权利要求1所述的测量方法,其中,步骤b中定位液膜厚度在原始图片高度方向上的像素范围是由软件CamControl完成。
3.如权利要求1所述的测量方法,其中,步骤b中对选定范围内的各像素点逐一读取,得到灰度级矩阵,是由Matlab软件完成。
4.如权利要求1所述的测量方法,其中,步骤d中固体壁面Y方向的像素点个数Nh通过下式计算得到:
式中,Nhi为任意选取不同固体壁面高度处像素点个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210090079.2A CN102607436B (zh) | 2012-03-30 | 2012-03-30 | 一种测量固体壁面波动液膜层平均厚度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210090079.2A CN102607436B (zh) | 2012-03-30 | 2012-03-30 | 一种测量固体壁面波动液膜层平均厚度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102607436A CN102607436A (zh) | 2012-07-25 |
CN102607436B true CN102607436B (zh) | 2014-05-14 |
Family
ID=46525043
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210090079.2A Expired - Fee Related CN102607436B (zh) | 2012-03-30 | 2012-03-30 | 一种测量固体壁面波动液膜层平均厚度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102607436B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107860322A (zh) * | 2017-12-07 | 2018-03-30 | 东华大学 | 一种液膜厚度测量装置及方法 |
CN113602018B (zh) * | 2021-07-06 | 2022-08-09 | 华中科技大学 | 一种柔性电子喷印制膜漏打印缺陷在线补偿方法及*** |
CN114323659B (zh) * | 2021-12-28 | 2023-01-06 | 上海交通大学 | 一种壁面可替换的温度可控的撞壁液膜测量装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1556914A (zh) * | 2001-09-21 | 2004-12-22 | Kmac株式会社 | 利用二维反射计测量多层薄膜的厚度轮廓和折射率分布的装置及其测量方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006153770A (ja) * | 2004-11-30 | 2006-06-15 | Omron Corp | 分光計測装置 |
US20100022893A1 (en) * | 2008-07-24 | 2010-01-28 | Hart Douglas P | Self-inflating bladder |
-
2012
- 2012-03-30 CN CN201210090079.2A patent/CN102607436B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1556914A (zh) * | 2001-09-21 | 2004-12-22 | Kmac株式会社 | 利用二维反射计测量多层薄膜的厚度轮廓和折射率分布的装置及其测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102607436A (zh) | 2012-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pandit et al. | Effect of pin fin to channel height ratio and pin fin geometry on heat transfer performance for flow in rectangular channels | |
Wilson et al. | Uncertainty on PIV mean and fluctuating velocity due to bias and random errors | |
CN102607436B (zh) | 一种测量固体壁面波动液膜层平均厚度的方法 | |
Horacek et al. | Spray cooling using multiple nozzles: visualization and wall heat transfer measurements | |
CN102452081B (zh) | 机械手臂***参数的校正方法与校正装置 | |
Goth et al. | PTV/PIV measurements of turbulent flows in interior subchannels of a 61-pin wire-wrapped hexagonal fuel bundle | |
CN108875719B (zh) | 基于深度学习与红外图像识别的空冷器积灰状态感知***与计算方法 | |
CN105445319B (zh) | 一种测定钢板表面水冷换热系数的方法及装置 | |
CN104537631A (zh) | 一种基于数字图像配准技术的plif图像校正方法 | |
CN102486462B (zh) | 一种合金固件内部缺陷三维重建方法 | |
Christensen et al. | Infrared temperature measurements of the blade tip for a turbine operating at corrected engine conditions | |
Gultekin et al. | PIV measurement inside single and double droplet interaction on a solid surface | |
Lillo et al. | Plenoptic single-shot 3D imaging of in-cylinder fuel spray geometry | |
CN106468671A (zh) | 一种瞬态液晶用于测量对流换热系数的非稳态测量方法 | |
CN106407987B (zh) | 基于图像互相关匹配的电光晶体光轴出露点提取方法 | |
Jha et al. | Heat transfer from a hot moving steel plate by air-atomized spray impingement | |
Martens et al. | Dependence of the fluctuation-dissipation temperature on the choice of observable | |
CN103744269B (zh) | 光刻投影物镜波像差和成像最佳焦面的检测方法 | |
EP1178436A3 (en) | Image measurement method, image measurement apparatus and image measurement program storage medium | |
Peng et al. | A deep learning method based on patchwise training for reconstructing temperature field | |
Budziankou et al. | Modeling approach for a Wiremesh substrate in CFD simulation | |
CN103439006A (zh) | 热冲压板料温度的测量装置和测量方法 | |
Khavari et al. | Time-dependent measurements of length and area of the contact line in contact-boiling regime | |
JP2008046012A (ja) | 欠陥検出装置および欠陥検出方法 | |
Funatani et al. | High-resolution three-color PIV technique using a digital SLR camera |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140514 |
|
CF01 | Termination of patent right due to non-payment of annual fee |