基于数字图像投影技术的三维液膜流动测量研究

宋华振, 兆环宇, 朱程香, 王正之, 田伟, 李海星, 朱春玲

宋华振,兆环宇,朱程香,等. 基于数字图像投影技术的三维液膜流动测量研究[J]. 实验流体力学,2021,35(5):106-114. DOI: 10.11729/syltlx20200031
引用本文: 宋华振,兆环宇,朱程香,等. 基于数字图像投影技术的三维液膜流动测量研究[J]. 实验流体力学,2021,35(5):106-114. DOI: 10.11729/syltlx20200031
SONG H Z,ZHAO H Y,ZHU C X,et al. Three-dimensional liquid film flow measurement based on digital image projection technology[J]. Journal of Experiments in Fluid Mechanics, 2021,35(5):106-114.. DOI: 10.11729/syltlx20200031
Citation: SONG H Z,ZHAO H Y,ZHU C X,et al. Three-dimensional liquid film flow measurement based on digital image projection technology[J]. Journal of Experiments in Fluid Mechanics, 2021,35(5):106-114.. DOI: 10.11729/syltlx20200031

基于数字图像投影技术的三维液膜流动测量研究

基金项目: 国家重点基础研究发展(973)计划(2015CB755800);国家自然科学基金(11832012;51806105);江苏省自然科学基金(BK20180442);中国博士后科学基金(2018M642253);江苏省博士后科研资助计划(2018K294C);江苏省高校优势学科建设工程资助项目;中国空气动力研究与发展中心飞行器结冰与防除冰重点实验室基金(AIADL20180402)
详细信息
    作者简介:

    宋华振: (1995),男,山东德州人,硕士研究生。研究方向:水膜流动实验研究、非接触测量。通信地址:江苏省南京市秦淮区御道街29号南京航空航天大学航空学院(210016)。E-mail:1501141520@qq.com

    通讯作者:

    朱春玲: E-mail:clzhu@nuaa.edu.cn

  • 中图分类号: V211; O359.1

Three-dimensional liquid film flow measurement based on digital image projection technology

  • 摘要: 流体薄膜流动定量测量是分析结冰相变传热过程的必要手段。基于图像处理的数字图像投影技术(DIP)可实现对流体薄膜的非侵入式定量测量。对DIP技术的基本原理、图像互相关算法和标定方法进行了介绍,设计并搭建了DIP测量系统和平板水膜流动实验台。DIP测量系统的整体误差在5%以内,证明了系统的可靠性与准确性。在平板水膜流动实验台上开展了一系列水膜流动实验,采用DIP测量系统复原了平板水膜流动的三维全貌。通过测量结果拟合得出平均水膜高度、无量纲水膜高度与水膜雷诺数之间的关系,并与理论推导和文献实验结果进行对比,整体趋势一致。
    Abstract: Quantitative measurement of fluid film flow is a necessary means to analyze the heat transfer process of the icing phase change. Digital Image Projection (DIP) based on image processing enables non-invasive quantitative measurement of fluid films. Firstly, the basic principles, cross-correlation algorithm and calibration method of the DIP technology are introduced. Based on this, the DIP system and the flat water film flow test bench are designed and built. Subsequently, the overall error of the DIP system was tested and was found to be within 5%, indicating the reliability and accuracy of the system. A series of water film flow experiments were carried out on a flat water film flow test bench. The DIP system was used to restore the three-dimensional full picture of the flat water film flow. Finally, the relationship between the average water film thickness as well as the dimensionless water film thickness and the water film Reynolds number is obtained by fitting the measurement results, and is compared with the theoretical derivation and literature experimental results. Results show the overall trend is consistent.
  • 飞机结冰是飞机飞行中常见的一种现象,对飞行安全危害极大。当飞机穿过含有过冷水滴的云层、结冰温度较高时,过冷水滴撞击飞机表面后在气流的驱动下形成水膜。水膜表面波特性分析可以有效地研究水膜流动和传热特性,对探究飞机结冰过程中的质量分布及热量传递意义重大。随着实验技术的进步,非侵入性的液膜高度分布测量技术被广泛应用,其中包括激光诱导荧光技术(Laser Induced Fluorescence,LIF)、反射波技术和图像处理技术等。

    LIF技术被广泛应用于定常流可视化测量,其原理是在液体中添加荧光剂,根据荧光剂所激发的荧光强度来判断流体的速度、高度和温度。Li等[1]开发了一种分子标记技术,可以同时测量飞溅水滴的尺寸、速度和瞬态温度等参数。Johnson等[2]开发了一种荧光成像系统,量化了等温条件下薄膜在平坦表面上的瞬态行为。常士楠等[3]利用激光诱导荧光技术测量了一定条件下水膜流动的高度和表面波动特性。

    反射波技术中较有代表性的是基于色散共焦位移计的测量方法。该方法通过色散镜头将光源射出的宽光谱复色光进行光谱色散,形成不同波长的单色光,每一个波长都对应一个光源到被测物体之间的距离,通过计算被感测到的波长换算距离。Rishikesan等[4]利用色散共焦传感器评估不同加工表面的轮廓参数,以此来测量微观水平表面的光洁度。冷梦尧等[5-6]利用色散共焦位移计对高速气流驱动下的平板水膜流动进行测量,得到了高速气流剪切作用下的气液波动界面剪切因子。

    基于图像处理的结构化光技术也是实现非侵入式高度测量的一种方式,近年来也受到流体动力学界的关注。该技术将已知的光图案主动投影到物体上,并从一个或多个视点捕获的光图案中提取三维表面形状[7]。Cazabat等[8]通过投射等距条纹来重构扩散薄膜的剖面高度,研究不同温度下的薄膜爬升特性。Cobelli等[9]使用数字条纹投影系统进行实验,以测量水流的表面变形。Hu等[10]开发出数字条纹投影系统,并使用该系统定量测量风驱动下的水膜和溪流在平板上的瞬态行为。

    在上述测量技术中,LIF技术凭借着其非侵入式测量和低信号−噪声比特性得到了广泛应用。但LIF技术的荧光强度通常依赖于温度,在不可控变温过程(如结冰)中很难使用LIF技术对流体薄膜高度进行定量测量。反射波技术虽然可以对几乎任何材料进行高精度无侵入式测量[11],但是无法实时测量三维流体流动的全貌。基于图像处理的结构化光技术虽能实现非侵入式流体全貌的复现,但也存在图像处理过程相位变化复杂、流体三维全貌预测精度较低等问题[12]

    数字图像投影技术(Digital Image Projection, DIP)具有时间分辨功能,是一种全区域水膜高度的测量方法。相比于其他基于图像处理的结构光技术,DIP技术可以通过更简单的图像处理方法在采集的图像中提取被测物体的三维表面全貌,且测量精度较好。Zhang等[13]应用该技术测量NACA0012翼型表面水膜和溪流随风速变化的流动过程,复原了水膜在不同风速下的三维全貌,但未建立水膜流动与外界条件间的关系式。国内在DIP技术方面的相关研究较为少见。

    为验证DIP技术在流体测量中的可靠性与准确性,本文自主设计搭建了DIP测量系统和平板水膜流动实验台,进行了一系列水膜流动实验,并利用标准件对DIP系统测量精度进行了测试。

    DIP测量技术的基本原理如图1所示。摄像机、投影仪和实验基板形成三角形布局。在目标流体未流过实验基板之前,投影仪投射网格图像于实验基板,此时摄像机采集实验基板上的网格图片作为参考图像。当目标流体流过实验基板时,网格图片将投射于流体表面并发生变形,摄像机记录此时的变形网格作为测试图像。通过分析测试图像和参考图像对应点的位移,就可以重构被测物体的三维形状。

    图  1  DIP技术原理图
    Fig.  1  DIP technology diagram

    为了将测试图像和参考图像对应点的位移转换为目标流体高度分布,DIP系统需要标定,以建立高度与位移之间的关系。

    建立如图1所示的坐标系,y方向高度为0的平面为后续测量的参考平面。实验基板上无任何目标流体时,任意网格点投影于点E;当流体流经实验基板时,原投影于点E的网格点将投影于关注点C。若此时摄像机的采集方向垂直于实验基板,根据相似三角形原理可得:

    $$ \frac{{\overline {DE} }}{d} = \frac{{\overline {CD} }}{{H - \overline {CD} }} $$ (1)

    则此时关注点C处的流体高度为:

    $$ \overline {CD} = \frac{{{H} \ \overline {DE} }}{{d + \overline {DE} }} $$ (2)

    距离d通常远大于$\overline {DE} $,上式可简化为:

    $$ \overline {CD} \approx \frac{H}{d}\ \overline {DE} = K \ \overline {DE} $$ (3)

    对于给定的DIP系统,H为摄像机到实验基板的高度,Hd是固定的,因此参数K为常数,可由标定方法确定;$\overline {DE} $可通过测量测试图像和参考图像之间关注点C处的位移来确定;通过式(3)就可以得出流体流经实验基板时关注点C处的高度信息,从而复原流体的整个三维全貌。

    DIP测量系统使用图像互相关算法[14-16]来确定测试图像与参考图像之间的关注点C处的位移量。在参考图像和测试图像的特定区域分别选择一个小窗口作为关联窗口,计算两关联窗口的相关系数。相关系数R以下式表示:

    $$ R = \frac{{\displaystyle\sum\limits_i {\displaystyle\sum\limits_j {({I_{\rm{r}}}(i,j) - {{\overline I }_{\rm{r}}})(I(i + r,j + s) - \overline I )} } }}{{\sqrt {\displaystyle\sum\limits_i {\displaystyle\sum\limits_j {{{({I_{\rm{r}}}(i,j) - \overline {{I_{\rm{r}}}} )}^2}\displaystyle\sum\limits_i {\displaystyle\sum\limits_j {{{(I(i + r,j + s) - \overline I )}^2}} } } } } }} $$ (4)

    其中,IrI分别表示从参考图像和测试图像中获得的关联窗口的图像强度,${\overline I _{\rm{r}}}$$ \overline I $为相应的平均强度值;ij分别表示水平方向和垂直方向,(rs)为图像位移量。通过找到相关系数的峰值位置就可确定两关联窗口的位移量,即关注点的位移量。图2展示了一个规则网格的参考图像和由测试板上的长方块产生的关注点移动的测试图像,使用上述算法即可得到红框内的网格点的位移量。

    图  2  网格点位移图
    Fig.  2  Grid intersection displacement diagram

    对于一个给定的DIP系统,测试图像和参考图像中对应点的位移随着目标物体的高度变化而呈线性变化。DIP标定可用于确定测试图像关注点的位移量与实际高度的转换系数K。如图3所示,通过数字投影仪将网格图像投影到实验基板上。实验基板水平安装于升降台顶部,升降台高度(即实验基板在y方向上的位置)用微米级驱动器调整。摄像机用于捕获投影在实验基板上的图像,所获取的图像被存储在计算机上用于图像处理。

    图  3  DIP标定示意图
    Fig.  3  DIP calibration diagram

    在DIP标定过程中,通过调整升降台的高度,将测试板以某一间隔移动到y方向上的多个位置,获取测试板上投影的网格图像,再利用图像互相关算法对图像进行处理。以y=0处的图像作为参考图像,计算测试板在其他位置时的关注点位移量,从而得到图像网格点的位移值。图4为测试板在不同高度时的测试图像网格点的像素位移量。从图中可以看出,高度与网格图像点的位移量之间的关系可以用线性函数拟合,符合式(3)的线性关系。标定过程得出的线性比例因子即为DIP测量系统的K值,且K值是一个常数。

    图  4  DIP标定结果图
    Fig.  4  DIP calibration result diagram

    图5为平板水膜流动实验台。风道气流经收缩段整流后进入实验段,为保证气流均匀,在收缩段上游设置蜂窝器和阻力网。实验段横截面为300 mm×75 mm。为方便更换实验基板,实验段顶部采用可拆卸的玻璃盖板。实验基板为抛光铝材,其下装有升降台。实验基板的一端开有出水槽,出水槽连接供水箱,供水流量通过高精度齿轮泵进行调节。

    图  5  平板水膜流动实验台
    Fig.  5  Flat water film flow test setup

    为了测试DIP系统的测量精度,使用DIP系统测量了长宽高为1200.0 mm×1000.0 mm×0.5 mm的实验基板。首先对DIP测量系统进行标定,得出系统的K值。随后获取参考图像和测试图像,使用图像互相关算法计算参考图像与测试图像之间关注点的位移量。最后,根据K值和位移量计算得到实验基板的高度分布。图6(a)给出了实验基板高度h'分布的DIP测量结果。实验基板的高度分布已知,DIP系统的测量精度可以通过测得的高度分布与实际高度分布进行定量比较评估。DIP测量误差分布如图6(b)所示。测量结果显示,实验基板中间区域测量误差较小(2%以内),上下两侧的测量误差在2%~5%之间。这是由于投影仪的投影方向与摄像机采集方向需呈一定角度,而DIP系统配置的摄像机垂直于实验基板,因此,投影仪倾斜投射出的网格在实验基板中部区域聚焦较好,网格点清晰,而两侧区域聚焦较差,网格点模糊,增大了测量误差。此外,DIP系统部件的搭建角度和摄像机分辨率的限制也是造成DIP测量误差的主要原因。由DIP测量原理可知,投影部件与摄像机部件之间的夹角越大,DIP测量系统单位距离的网格点像素位移量越大,测量精度越高。受实验条件的限制,本文DIP测量系统的部件夹角为81°。在图像互相关算法原理中,摄像机的分辨率越高,图像互相关算法计算出的最终位移量精度越高。本文实验摄像机的分辨率为4×106 pxiel,拍摄区域面积为80 mm×120 mm。

    图  6  实验台精度
    Fig.  6  Test setup precision

    为说明DIP测量系统对毫米级高度测量的精度,本文在已搭建好的水膜流动实验台上使用DIP测量技术进行3 mm高度范围内的梯度测量。引入色散共焦位移计同时测量,将其测量结果作为参照。在图3所示的实验台基础上将色散共焦位移计垂直于实验基板固定。由低到高等间隔升高升降台,每次升高500 μm,共6次。每次升高升降台后,采集网格测试图像并同时记录色散共焦位移计测量结果。该实验重复3次,最终结果取平均值,如图7所示。从图中可以看出:两种测量方式得到的结果均呈现出良好的线性趋势;随着测量高度增大,共焦位移计与DIP测量系统的测量误差均增大。相较于DIP测量,共焦位移计的测量误差随着实测距离越来越大。这是由于共焦位移计有特定的测量距离范围,当超出测量范围时,共焦位移计不能很好地聚焦待测平面,导致测量距离越大、测量误差越大;这也可以反映出在3 mm高度范围内,当测量高度较大时,DIP测量技术的精度优于色散共焦位移计。

    图  7  测量结果对比图
    Fig.  7  Comparison of measurement results

    在水膜流动中,水膜雷诺数定义为:

    $$ {{\mathop{Re}\nolimits} _\text{f}} = Q/(W \ {\mu _\text{w}}) $$ (5)

    其中,Q为水膜质量流量,W为水膜润湿宽度,μw为水的动力黏度系数。

    在实验中,流场风速Ua=17 m/s,水膜雷诺数Ref范围为16.5~188.6。为增大水表面的反射性,在水中加入1%体积比的白色乳胶漆。通过落球法测试配置液体在温度为20 ℃时的动力黏度μw=0.985×10–3 Pa·s。在实验段上部空间,使用Texas instrument DLP450型投影仪投影方形网格图像,同时使用PCO HS4型高速摄像机采集图像,采集频率为100 Hz,曝光时间为1 ms。摄像机匹配镜头型号为定焦尼克尔105 mm。图8显示了实验基板水膜的参考图像和变形的测试图像,h为水膜高度。投影的单个方形网格尺寸为2.8 mm×2.8 mm,网格条纹宽度为0.8 mm。与参考图像(图8(a))相比,测量图像(图8(b))由于存在水膜而显著变形。将参考图片与测量图片使用图像互相关算法进行计算,计算过程中的关联窗口设置为2.4 mm×2.4 mm的方形窗口。图8(c)展示了经过三维重构后获得的水膜流过铝板表面的高度分布。

    图  8  实验采集图像与计算结果对比图
    Fig.  8  Comparison of test photo and calculation

    本文开展了不同雷诺数下的水膜流动实验,水膜高度的DIP计算结果如图9所示,L为水膜流动长度。整体上,水膜的高度随着水膜雷诺数的增大而增大。水膜雷诺数在16.5~23.6范围时(图9(a)(b)),润湿区域的水量较少,风速对表面波的传导影响较大,此时水膜整体较薄,无法形成稳定的表面波传导。随着水膜雷诺数的增大(图9(c)(f)),润湿区域内的水量逐渐增多,水膜在气液界面剪切力的作用下形成表面波。随着水膜雷诺数的增大,水膜表面波越来越向润湿区域中心集中,中心润湿区域的水膜高度越来越大。当水膜雷诺数继续增大到165.0及以上时(图9(g)(h)),实验设定的风速不能够以稳定的表面波向后发展水膜。

    图  9  水膜流动DIP计算结果(Ua=17 m/s、Ref =16.5~188.6)
    Fig.  9  DIP calculation results of water film flow (Ua=17 m/s and Ref =16.5-188.6)

    为了更好地研究水膜雷诺数对水膜表面波形的影响,分别取水膜雷诺数为16.5、70.7、188.6的DIP计算结果的中间截面剖面曲线,如图10所示。Ref =16.5时,水膜表面波呈现出振幅小,周期短且不稳定的特性。Ref =70.7时,润湿区域内水膜的高度整体增大,水膜表面波振幅和周期均较稳定。Ref =188.6时,在水膜平均高度增大的同时,出现了表面波振幅较大、周期较长的现象。因此,在特定风速下随着水膜雷诺数的增大,水膜表面波呈现出振幅逐渐增大、振动频率逐渐降低的波动特性。

    图  10  水膜剖面表面波形图
    Fig.  10  Water film surface wave profile

    在水膜高度较薄、水流速度较低的条件下,水膜流动可简化为二维不可压缩层流流动。对Navier-Stokes方程进行简化得出水膜高度的关系式[3]为:

    $$ h = \sqrt {\frac{{2{\mu _{\rm{w}}}{Q_{\rm{a}}}}}{\tau }} $$ (6)

    其中Qa为上游流入的单位润湿长度水的体积流量,τ为气−液界面剪切力。借鉴气−液两相分层流研究的经验[17],界面剪切力可以表示为:

    $$ \tau = \frac{1}{2}{f_i}{\rho _{\rm{a}}}{\left( {{U_{\rm{a}}} - C} \right)^2} $$ (7)

    其中,fi为界面剪切系数,通过实验结果测得;ρa为空气密度;Ua为空气主流速度;C为气−液界面波的相位速度,由于C远小于气流速度,故在计算中忽略[18-19]。对于平板边界层模型[20],界面剪切系数${f_i} = 0.074{\left( {{{{\mathop{ Re}\nolimits} }_{\rm{a}}}} \right)^{ - 0.2}}$,其中${{\mathop{ Re}\nolimits} _{\rm{a}}} = {U_{\rm{a}}}H/{\upsilon _{\rm{a}}}$为气流雷诺数,υa为空气运动黏度系数。通过式(6)和(7)可得水膜高度与水膜雷诺数之间的关系为:

    $$ h = 7.35\ {\mu _{\rm{w}}} \cdot U_{\rm{a}}^{ - 1} \cdot {\left( {{\rho _{\rm{a}}}\ {\rho _{\rm{w}}}} \right)^{ - 0.5}} \cdot {\mathop{ Re}\nolimits} _{\rm{a}}^{0.1} \cdot {R}{e}^{0.5} $$ (8)

    本文测量了8个水膜雷诺数状态下10个时刻的水膜高度分布,并将结果平均化处理得到平均水膜高度$ \bar h $。平均水膜高度$ \bar h $与雷诺数的关系如图11所示,用幂函数拟合得到其关系式为:

    $$ \bar h = 0.038\,3{{\mathop{ Re}\nolimits} _{\rm{f}}}^{0.568\,9} $$ (9)
    图  11  平均水膜厚度对比图
    Fig.  11  Comparison of average film thickness

    图11所示,实验得出的拟合曲线与理论推导曲线整体趋势相近。但是水膜雷诺数的实测值略大于理论推导值,并且该趋势随水膜雷诺数的增大逐渐明显。这是由于随着水膜雷诺数的增大,水膜的流速增大,水膜流速增大到一定程度时形成的表面波为滚动不稳定状态。此时水膜的表面波将影响气−液界面剪切力τ的计算,随水膜雷诺数增大而增大,理论推导曲线与实验曲线差别逐渐增大。因此,当水膜雷诺数增大到一定程度时,需考虑水膜表面相位速度对水膜高度的影响。此外,本文的理论曲线基于一系列假设推导得出,这对比较结果也带来一定影响。随着水膜雷诺数的增大,水膜高度增大,DIP投影网格清晰度降低,加大了图像互相关算法的误差,这也使得图像处理结果误差增大。

    为进一步验证DIP测量技术应用于平板水膜流动的可靠性与准确性,特引入无量纲水膜高度。将平均水膜高度与剪切应力相结合进行无量纲化,定义水膜无量纲高度h+为:

    $$ {h^ + } = \frac{{{u^*}\bar h}}{{{\upsilon _{\rm{w}}}}} $$ (10)

    其中,${u^*} = \sqrt {\tau /\rho } $为水膜内的速度,υw为水的运动黏度系数,$ \bar h $取水膜平均高度实测值。由式(5)、(6)和(8)得出无量纲水膜高度$ h^ + $与水膜雷诺数Ref之间的关系为:

    $$ {h^ + } = 1.397\, 3 \ {\mathop{ Re}\nolimits} _{\rm{f}}^{0.568\, 9} $$ (11)

    表1列出了以往文献及本研究的无量纲水膜高度的计算模型,其中Kosky & Staub模型[21]为水平冷凝环状流模型,Rishikesan模型[4]为基于色散共焦位移计测量的平板水膜流动模型,Hughmark模型[22]为等温条件下的上升流模型。与前两者比较,Hugh-mark模型将气−液两相的质量流速、管直径和流体的物理性质作为研究变量来预测水膜高度,本文预测模型与之类似。各模型计算值结果如图12所示。可以看出,虽然本文研究对象与Rishikesan模型的研究对象均为平板水膜流动,但是本文使用的DIP测量方法与Rishikesan模型的单点共焦位移计测量方法最终测量误差存在差别(如图7所示),从而导致这两种模型存在偏差。本文测量水膜高度方法虽与其他文献的测量方法不同,但应用到水膜无量纲高度的模型中,与其他研究模型计算结果的趋势具有良好的一致性,更为接近Hughmark模型,误差在5%以内。

    表  1  无量纲水膜厚度模型
    Table  1  Dimensionless Film Thickness Models
    AuthorAppliedRefRelation
    Kosky & Staub[21]Horizontal
    condense flow
    13.0~370.0 $h^{+}=1.316\,0 \ { {Re} }_{ {\rm{f} } }^{0.529\,9}$
    Rishikesan[4]Plate flow26.0~128.0 $\,\, h^{+}=1.292\,0 \ {{Re} }_{ {\rm{f} } }^{0.528\,9}$
    Hughmark[22]Upward flow25.0~250.0$h^{+}=0.874\,0 \ { {Re} }_{ {\rm{f} } }^{2 / 3}$
    Present workPlate flow16.0~190.0 $h^{+}=1.397\,3\ {{Re} }_{ {\rm{f} } }^{0.568\,9}$
    下载: 导出CSV 
    | 显示表格
    图  12  无量纲水膜厚度对比图
    Fig.  12  Comparison of Dimensionless Water Film Thickness

    本文介绍了可实现非侵入定量测量的DIP技术,说明了其测量原理、标定方法,并将其应用于平板表面水膜流动的定量测量,最后将测量结果与理论推导和文献结果进行对比得到以下结论:

    1) DIP测量技术可对流体流动的三维全貌进行复原。本文复原了平板水膜表面波随雷诺数的形态变化。

    2) DIP测量技术拥有良好的精度。在本文的实验台参数设置下,DIP系统的测量误差小于5%。高度在3 mm以内的测量,DIP测量技术精度优于色散共焦位移计。

    3)平板水膜流动实验结果显示:在风速17 m/s状态下,随着水膜雷诺数的增大,平均水膜高度整体增大;当水膜雷诺数增大到一定程度(Ref >165.0),水膜向中间聚拢,呈现出中间区域高于两侧区域的现象。随着水膜雷诺数的增大,水膜表面波呈现出振幅逐渐增大、振动频率逐渐降低的波动特性。

    4)测量与理论推导发现:水膜表面波相位速度的定量测量是影响本文平均水膜高度测量的关键因素,会导致实测值与理论推导值的差别随着水膜雷诺数的增大而增大。

    5)无量纲水膜高度与水膜雷诺数的关系模型与文献中的模型吻合良好,再次证明了DIP技术在流体薄膜测量中的可靠性与准确性。

  • 图  6   实验台精度

    Fig.  6   Test setup precision

    图  1   DIP技术原理图

    Fig.  1   DIP technology diagram

    图  2   网格点位移图

    Fig.  2   Grid intersection displacement diagram

    图  3   DIP标定示意图

    Fig.  3   DIP calibration diagram

    图  4   DIP标定结果图

    Fig.  4   DIP calibration result diagram

    图  5   平板水膜流动实验台

    Fig.  5   Flat water film flow test setup

    图  7   测量结果对比图

    Fig.  7   Comparison of measurement results

    图  8   实验采集图像与计算结果对比图

    Fig.  8   Comparison of test photo and calculation

    图  9   水膜流动DIP计算结果(Ua=17 m/s、Ref =16.5~188.6)

    Fig.  9   DIP calculation results of water film flow (Ua=17 m/s and Ref =16.5-188.6)

    图  10   水膜剖面表面波形图

    Fig.  10   Water film surface wave profile

    图  11   平均水膜厚度对比图

    Fig.  11   Comparison of average film thickness

    图  12   无量纲水膜厚度对比图

    Fig.  12   Comparison of Dimensionless Water Film Thickness

    表  1   无量纲水膜厚度模型

    Table  1   Dimensionless Film Thickness Models

    AuthorAppliedRefRelation
    Kosky & Staub[21]Horizontal
    condense flow
    13.0~370.0 $h^{+}=1.316\,0 \ { {Re} }_{ {\rm{f} } }^{0.529\,9}$
    Rishikesan[4]Plate flow26.0~128.0 $\,\, h^{+}=1.292\,0 \ {{Re} }_{ {\rm{f} } }^{0.528\,9}$
    Hughmark[22]Upward flow25.0~250.0$h^{+}=0.874\,0 \ { {Re} }_{ {\rm{f} } }^{2 / 3}$
    Present workPlate flow16.0~190.0 $h^{+}=1.397\,3\ {{Re} }_{ {\rm{f} } }^{0.568\,9}$
    下载: 导出CSV
  • [1]

    LI H X,CHEN F,HU H. Simultaneous measurements of droplet size, flying velocity and transient temperature of in-flight droplets by using a molecular tagging technique[J]. Experiments in Fluids,2015,56(10):1-14. doi: 10.1007/s00348-015-2063-y

    [2]

    JOHNSON M F G,SCHLUTER R A,BANKOFF S G. Fluorescent imaging system for global measurement of liquid film thickness and dynamic contact angle in free surface flows[J]. Review of Scientific Instruments,1997,68(11):4097-4102. doi: 10.1063/1.1148352

    [3]

    CHANG S N,YU W D,SONG M J,et al. Investigation on wavy cha-racteristics of shear-driven water film using the planar laser induced fluorescence method[J]. International Journal of Multiphase Flow,2019,118:242-253. doi: 10.1016/j.ijmultiphaseflow.2019.04.016

    [4]

    RISHIKESAN V,SAMUEL G L. Evaluation of surface profile para-meters of a machined surface using confocal displacement sensor[J]. Procedia Materials Science,2014,5:1385-1391. doi: 10.1016/j.mspro.2014.07.456

    [5] 冷梦尧,常士楠,丁亮,等. 水平表面气流剪切作用下的水膜厚度[J]. 航空学报,2017,38(2):65-73.

    LENG M Y,CHANG S N,DING L,et al. Thickness of water film driven by gas stream on horizontal plane[J]. Acta Aeronautica et Astronautica Sinica,2017,38(2):65-73.

    [6]

    LENG M Y,CHANG S N,WU H. Experimental investigation of shear-driven water film flows on horizontal metal plate[J]. Experi-mental Thermal and Fluid Science,2018,94:134-147. doi: 10.1016/j.expthermflusci.2018.02.004

    [7]

    SALVI J,FERNANDEZ S,PRIBANIC T,et al. A state of the art in structured light patterns for surface profilometry[J]. Pattern Reco-gnition,2010,43(8):2666-2680. doi: 10.1016/j.patcog.2010.03.004

    [8]

    CAZABAT A M,HESLOT F,TROIAN S M,et al. Fingering in-stability of thin spreading films driven by temperature gradients[J]. Nature,1990,346(6287):824-826. doi: 10.1038/346824a0

    [9]

    COBELLI P J,MAUREL A,PAGNEUX V,et al. Global measure-ment of water waves by Fourier transform profilometry[J]. Experi-ments in Fluids,2009,46(6):1037-1047. doi: 10.1007/s00348-009-0611-z

    [10]

    HU H,WANG B,ZHANG K,et al. Quantification of transient behavior of wind-driven surface droplet/rivulet flows using a digital fringe projection technique[J]. Journal of Visualization,2015,18(4):705-718. doi: 10.1007/s12650-014-0264-8

    [11] 王津楠. 光谱共焦位移传感器研究与设计[D]. 哈尔滨: 哈尔滨工业大学, 2016: 8-10.

    WANG J N. Research and design of chromatic confocal displacement sensor[D]. Harbin: Harbin Institute of Technology, 2016: 8-10.

    [12]

    ZHANG K, ZHANG S, ROTHMAYER A, et al. Development of a digital image projection technique to measure wind-driven water film flows[R]. AIAA 2013-0247, 2013.

    [13]

    ZHANG K,WEI T,HU H. An experimental investigation on the surface water transport process over an airfoil by using a digital image projection technique[J]. Experiments in Fluids,2015,56(9):1-16. doi: 10.1007/s00348-015-2046-z

    [14]

    HU H,KOOCHESFAHANI M M. Molecular tagging velocimetry and thermometry and its application to the wake of a heated circular cylinder[J]. Measurement Science and Technology,2006,17(6):1269-1281. doi: 10.1088/0957-0233/17/6/s06

    [15]

    HU H,KOOCHESFAHANI M M. Thermal effects on the wake of a heated circular cylinder operating in mixed convection regime[J]. Journal of Fluid Mechanics,2011,685:235-270. doi: 10.1017/jfm.2011.313

    [16]

    GENDRICH C P,KOOCHESFAHANI M M. A spatial correlation technique for estimating velocity fields using molecular tagging velocimetry (MTV)[J]. Experiments in Fluids,1996,22(1):67-77. doi: 10.1007/BF01893307

    [17]

    KAREV A R,FARZANEH M,LOZOWSKI E P. Character and stability of a wind-driven supercooled water film on an icing sur-face—I. Laminar heat transfer[J]. International Journal of Thermal Sciences,2003,42(5):481-498. doi: 10.1016/S1290-0729(02)00049-2

    [18]

    ANDRITSOS N,HANRATTY T J. Influence of interfacial waves in stratified gas-liquid flows[J]. AIChE Journal,1987,33(3):444-454. doi: 10.1002/aic.690330310

    [19]

    TZOTZI C,ANDRITSOS N. Interfacial shear stress in wavy stratified gas-liquid flow in horizontal pipes[J]. International Journal of Multiphase Flow,2013,54:43-54. doi: 10.1016/j.ijmultiphaseflow.2013.03.003

    [20] 吴望一. 流体力学(下册)[M]. 北京: 北京大学出版社, 2004.
    [21]

    KOSKY P G,STAUB F W. Local condensing heat transfer coeffi-cients in the annular flow regime[J]. AIChE Journal,1971,17(5):1037-1043. doi: 10.1002/aic.690170505

    [22]

    HUGHMARK G A. Film thichness, entrainment, and pressure drop in upward annular and dispersed flow[J]. AIChE Journal,1973,19(5):1062-1065. doi: 10.1002/aic.690190533

  • 期刊类型引用(1)

    1. 陈清华,余健军,王建刚,季家东,牛天宇. 基于图像处理的管翅式换热器翅片尺寸测量方法研究. 环境技术. 2023(08): 87-93+117 . 百度学术

    其他类型引用(1)

图(12)  /  表(1)
计量
  • 文章访问数:  656
  • HTML全文浏览量:  264
  • PDF下载量:  55
  • 被引次数: 2
出版历程
  • 收稿日期:  2020-03-08
  • 修回日期:  2020-04-02
  • 网络出版日期:  2021-11-03
  • 刊出日期:  2021-11-04

目录

/

返回文章
返回
x 关闭 永久关闭