Formation of confined free film in the window on a vertical perforated plate
-
摘要: 液膜是常见的气液接触方式之一,依据单或双自由面可分为壁面膜和自由膜。不同于传统液膜研究,板上开窗或开孔后会同时存在这2种液膜。实验采用了6种液体(Ka数52~3000)和10余种不同尺寸的矩形窗(9~1152mm2),观察了垂直板上下降液膜通过开窗区产生的自由膜流动行为。结果表明,在表面张力、惯性力、重力及黏性力的共同作用下,液体通过开窗区的方式有绕流、过流、偏流和"背流"4种主要形式。随着流量的增加,过流又可以有液滴、液柱、液膜以及它们的组合等丰富流型,同时窗内自由膜在表面张力主导下易受扰动而聚并、破碎,流型呈现多样性和动态性。当流量达到临界值时,窗口内会形成完整稳定的受限自由液膜,呈现板上壁面膜与窗口内受限自由膜交替共存的特殊液膜——孪生液膜。本文关联了临界Re数、Ka数和无量纲几何参数Nx,给出了成膜条件的经验判别式。孪生液膜具有的特殊波列结构、成膜过程的迟滞现象以及2种液膜和固体壁间相互作用引起的多重干扰等,都体现出与传统的壁面液膜和自由液膜不同的流动性质。实验结果将有助于工业过程中开窗和开孔的几何优化,强化传递过程,并进一步丰富传统的液膜研究领域。Abstract: Liquid film is a common contact way between gas/liquid and can be divided into wall-bounded film and free film according to the film formation process and the number of free surfaces. When a plate is perforated, both kinds of films exist. This paper describes some experimental observations of free-surface flows arising when a thin liquid film flows through the window on a vertical perforated plate. Dozens of rectangular windows (9~1152mm2) and six fluids (Ka from 52 to 3000) were used to investigate the flow mechanisms. Several typical flow patterns including pass-around flow, pass-through flow, bias flow and back-side flow are well defined under such a flow condition. Here a special focus is given to the window region, where various free-surface flow patterns composed of droplets, columns, sheets and their combinations were observed with increasing flow rate. Meantime, free film in the window is surface tension dominated and susceptible to disturbance, resulting in coalescence or break of liquid columns and films. At a critical flow rate, liquid film is able to full fill the window, forming a stable complete confined free film. Mutual influence between confined free film in the window region and wall-bounded film around shows special wavy trains, which is also called "twin liquid film". Based on experimental data and scaling analysis, an empirical equation which relates Reynolds number Re, Kapitza number Ka and a dimensionless length Nx is proposed to characterize the film formation conditions. It is found that critical film formation Reynolds number increases with Kapitz number and window size. And hysteresis phenomenon is manifested by obviously different flow transition Re for confined free film formation and breaking. The results can help window geometry optimization industrial processes to improve local heat and mass transfer. It can also enrich the traditional film flow investigations.
-
Keywords:
- twin-liquid film /
- flow transition /
- film formation /
- perforated plate
-
0 引 言
螺旋桨动力系统对飞机气动特性的影响(动力影响)分为直接影响和间接影响两部分。直接影响是指螺旋桨本身的气动载荷对全机的影响;间接影响则是指气流经过螺旋桨加速后,形成的滑流扫掠飞机其他部件产生的气动特性变化。螺旋桨滑流内部流动特性较为复杂,不仅其轴向速度较来流速度大、还存在一定的旋转运动。当滑流扫掠机翼时,其主要作用表现为提高机翼升力、改善飞机起降性能;当滑流扫掠到平尾时,会对飞机俯仰特性产生严重影响,在平尾相对靠近螺旋桨时影响更为显著。
对螺旋桨滑流影响的研究主要有风洞试验和数值模拟两种方法。带动力的风洞试验通过模拟螺旋桨的飞行参数,可以较为准确地得到螺旋桨滑流对飞机气动特性的影响。随着风洞试验技术的不断发展和完善,对螺旋桨滑流的研究也越来越深入,研究内容由滑流影响量获取发展到滑流特性等领域。Müller等[1]在低速风洞中研究了螺旋桨滑流对A400M气动特性的影响。李兴伟等[2]采用螺旋桨飞机动力模拟风洞试验技术及粒子图像测速技术,研究了双发常规布局涡桨飞机的螺旋桨滑流对飞机纵向气动特性的影响规律,将螺旋桨位置变化对全机气动特性的影响进行了较为深入的研究,对后续螺旋桨飞机的布局设计有重要的指导作用。
近年来,计算流体力学的发展使螺旋桨滑流数值模拟技术获得了长足进步。徐静[3]、张刘[4]、曾卓雄[5]等采用计算技术对螺旋桨滑流模拟方法、螺旋桨滑流影响进行了研究。任晓峰[6]、王伟[7]、陈波[8]等采用数值模拟方法,对滑流影响下的飞机俯仰力矩特性变化进行了分析,指出滑流对机翼、平尾的气动干扰是影响全机俯仰力矩特性变化的重要因素,并对飞机尾翼布局优化提出了建议。
当前对螺旋桨滑流的研究[9-10]已经较为深入,但国内的研究大多集中在数值模拟[11-15]方面,通过风洞试验[16-17]来解决问题的研究还相对较少。
某飞机采用大展弦比机翼,翼吊式双螺旋桨发动机。研究发现该飞机起降构型在螺旋桨大拉力情况下,中小迎角俯仰力矩特性不够理想,表现为俯仰静稳定性裕量降低,甚至出现静不稳定现象,严重影响起降安全。本文通过对飞机气动数据的深入研究,分析了问题出现的原因,给出了优化方案。
1 试验简介
试验模型为某飞机缩比全金属模型和六叶螺旋桨模型,图1给出了该模型的示意图。该飞机的系列试验在FL-13风洞第二试验段进行。FL-13风洞是一座直流式、闭口、串列双试验段的大型风洞,其第二试验段宽为8 m、高为6 m、长为15 m,有效截面积为47.4 m2,常用风速范围为20~85 m/s。
试验支撑机构为FL-13风洞特大迎角支撑机构,模型以斜腹撑方式支撑。该设备可实现模型的迎角α和偏航角β的变化,变化范围为:α=–20º~120º,β=–30º~30°。
全机载荷由布置于模型机身内部的主天平测量,另有两台动力天平分别测量螺旋桨气动载荷;螺旋桨动力通过内置电机带动桨模型实现。带动力影响试验按照间接模拟法的固定拉力系数法进行,当给定一个拉力系数后,在全部试验姿态角范围内固定不变,即在试验中改变模型姿态角时不改变螺旋桨桨叶角和转速。
2 俯仰稳定性问题
该飞机俯仰稳定性问题主要出现在起降构型大拉力状态,图2给出了起飞构型下的俯仰力矩曲线。无动力情况下,该飞机在失速前俯仰力矩曲线线性度较好,没有明显的拐折出现、俯仰静稳定导数基本保持稳定;带动力后,其拉力系数
${T_{\rm c}}$ =0.1时、在零迎角附近曲线出现了拐折,负迎角俯仰静稳定导数明显降低;随着${T_{\rm c}}$ 的增大,曲线拐点逐渐右移,拐点左侧逐渐转为俯仰静不稳定。在${T_{\rm c}}$ =0.4时,拐点已进入飞行包线范围,严重影响飞机的起降安全。3 问题分析
3.1 主要影响原因分析
飞机大襟翼构型、中小迎角时平尾(下表面)易失速,此时平尾提供的(负向)升力损失,导致全机俯仰力矩特性发生改变、俯仰力矩曲线出现拐点。在带螺旋桨动力情况下,机翼升力系数增加,平尾处下洗增大,表面更易失速。但通过对测压结果的分析表明,当飞机俯仰力矩特性出现异常时,平尾并没有出现明显的失速迹象(见图3),说明该飞机俯仰力矩特性异常并非源于平尾的分离。
图2表明,俯仰力矩曲线拐点的出现与螺旋桨动力影响有密切的关系。为梳理其主要影响要素,对动力影响进行了必要的分解。利用部件试验结果和螺旋桨天平测量结果,可以把飞机各主要部件对全机俯仰力矩的贡献表示为:
$$ {C}_{m,{\text{全机}}}\approx {C}_{m,{\text{去尾翼}}}+{C}_{m,{\text{尾翼}}}+{C}_{m,{\text{直接影响}}} $$ (1) 图4给出了式(1)各项贡献的试验结果。由图可知,螺旋桨的直接影响对全机俯仰力矩的贡献随迎角呈线性变化,且其斜率为正。这表明螺旋桨的直接影响是降低俯仰静稳定裕量,但对俯仰力矩的拐折基本无影响。尽管翼身组合体(即去尾翼)俯仰力矩曲线也存在拐折现象,但与无动力试验结果比较(见图5),在中小迎角下,俯仰静稳定裕量变化更小,因此判断翼身不是全机俯仰力矩曲线拐折原因。在中小迎角(
$ \alpha \leqslant $ 10°)下,随着$ \alpha $ 减小,尾翼对俯仰力矩的贡献呈现出明显的非线性特性,是全机俯仰力矩曲线产生拐折的主要原因。尾翼对全机俯仰力矩特性的贡献,可表示为两部分:
$$ {C}_{m,{\text{尾翼}}}\approx {C}_{m,{\text{尾翼}}({\text{无动力}})}+{C}_{m,{\text{尾翼}}({\text{滑流影响}})} $$ (2) 利用有无动力、有无尾翼试验结果,可以分解出有、无动力情况下尾翼对俯仰力矩的贡献(见图6)。图中两条曲线差量即为滑流对尾翼俯仰贡献的影响。结果表明:无动力中小迎角情况下,尾翼贡献随迎角变化呈线性变化;而带动力情况下随迎角变化呈非线性变化。该非线性变化的主原因在于,随着迎角的变化,尾翼与螺旋桨滑流扫掠区域位置关系发生了明显改变,即随迎角减小,尾翼逐渐上抬、向滑流区上部移动,并最终移出滑流区。在此过程中尾翼受滑流影响的区域逐渐减小(即滑流影响降低)。
3.2 影响滑流与尾翼相互位置关系的主要参数
滑流与尾翼相互位置关系主要受两个参数影响:螺旋桨轴线与平尾弦线的相对高差和滑流扫掠区域洗流角。
3.2.1 螺旋桨轴线与平尾弦线的相对高差影响分析
试验中通过改变平尾上反角的方法,对平尾有效高度影响进行了研究。结果如图7所示,可以看出,调整平尾有效高度能够影响平尾进出滑流的过程,从而改变俯仰力矩拐点的出现位置。不同拉力系数下,减小平尾有效高度(减小上反角),俯仰力矩曲线拐点均明显左移。在
${T_{\rm c}}$ =0.4时,曲线拐点迎角减小约4°,极大地改善了俯仰力矩特性。3.2.2 滑流扫掠区域洗流角影响分析
对比空中与近地试验结果(见图8)发现,起降构型下近地状态相对空中状态,俯仰力矩曲线拐点明显左移。这主要是受地面效应影响,近地状态机翼在滑流扫掠区域产生的下洗角明显减小,导致滑流位置相对上移,使平尾脱离滑流区的迎角减小,推迟了拐点的出现。
俯仰力矩曲线变化的另一种表现是拐点随升降舵偏度变化而规律性变化(见图9),升降舵上偏时拐点右移、升降舵下偏时拐点左移。分析认为,这也是滑流扫掠区域下洗变化对滑流影响的一种表现。升降舵上偏时,平尾负向升力(向下)增大,诱导前缘下洗增大,使平尾处滑流区下移,平尾会在较大迎角脱离滑流区,从而使俯仰力矩拐点右移,反之亦然。
3.3 螺旋桨影响分析
通过上述分析可以得出以下结论:1)螺旋桨动力直接影响俯仰静稳定性,但不影响俯仰力矩曲线拐折。2)螺旋桨轴线与平尾相对高差、滑流扫掠区域洗流角变化造成平尾在迎角变化过程中进出滑流区,是导致俯仰力矩曲线拐折迎角变化的主要原因;适当减小螺旋桨轴线与平尾相对高差有利于改善中小迎角俯仰静不稳定问题。
4 优化效果
依据上述研究成果,以减小螺旋桨轴线与平尾相对高差为主要手段,对某飞机总体布局进行了优化研究。具体措施包括:降低平尾基准高度、减小平尾上反角、提高螺旋桨轴线位置等(见图10)。
验证试验结果表明,降低平尾基准高度和减小平尾上反角等措施均可有效改善大拉力情况下起降构型中小迎角俯仰力矩曲线拐折问题,使拐点左移。最终的优化方案确定为将平尾基准高度降低约70 mm(模型)、上反角减小5°,依此方案获得的全机俯仰力矩特性(见图11)满足设计要求和起降需要。
降低平尾高度改进了全机纵向静稳定性,但对升阻特性、横航向特性也会产生一定影响。图12给出了巡航构型下平尾高度调整前后升阻特性对比。由图可知,降低平尾高度后,升阻特性变化不大。图13给出的航向特性对比表明,降低平尾高度对航向稳定性影响较为明显。降低平尾后,结合垂尾布局调整,
$ {C_n} $ 曲线在小的负侧滑出现的非线性明显减弱。5 结 论
通过以上研究,可以得出以下结论:
1)螺旋桨滑流对平尾的影响是造成某飞机初始方案在起降构型大拉力中小迎角下出现俯仰静不稳定的主要原因。
2)螺旋桨轴线与平尾相对高差、滑流扫掠区域洗流角变化,造成平尾在迎角变化过程中进出滑流区,是导致飞机俯仰力矩曲线拐折迎角变化的主要原因。
3)通过降低平尾有效高度的方法,减小螺旋桨轴线与平尾相对高差,能有效减小飞机出现俯仰力矩曲线拐折的迎角,改善中小迎角俯仰特性。
-
表 1 实验流体的物理性质
Table 1 Physical properties of the liquids used
Ka 密度
ρ/(kg·m-3)黏度
μ/(N·s·m-2)表面张力
σ/(N·m-1)68%甘油水溶液 52.7 1160 0.0175 0.058 63%甘油水溶液 127 1150 0.0090 0.057 55%甘油水溶液 295 1110 0.0046 0.055 63%酒精溶液 406 880 0.0014 0.037 无水乙醇 670 790 0.0027 0.030 去离子水 2.93×103 1000 0.00098 0.072 表中: ,g为重力加速度 -
[1] 胡军, 胡国辉, 孙德军.沿平板下落薄膜流动的研究综述[J].力学进展, 2005, 35(2):161-169. DOI: 10.6052/1000-0992-2005-2-J2004-109 Hu J, Hu G H, Sun D J. A review on thin films falling along an inclined plate[J]. Advances in Mechanics, 2005, 35(2):161-169. DOI: 10.6052/1000-0992-2005-2-J2004-109
[2] Batchelor G K, Moffatt H K, Worster M. Perspectives in fluid dynamics:a collective introduction to current research[M]. Cambridge:Cambridge University Press, 2002.
[3] Kapitza P. Dynamic stability of the pendulum with vibrating suspension point[J]. Soviet Physics-JETP, 1951, 21(5):588-597. https://www.mendeley.com/research-papers/dynamic-stability-pendulum-vibrating-suspension-point/
[4] Wasden F K, Dukler A E. A numerical study of mass transfer in free falling wavy films[J]. American Institute of Chemical Engineers Journal, 1990, 36(9):1379-1390. DOI: 10.1002/(ISSN)1547-5905
[5] Alekseenko S, Nakoryakov V E, Pokusaev B G. Wave flow of liquid films[M]. New York:Begell House, 1994.
[6] Chang H-H, Demekhin E A. Complex wave dynamics on thin films[M]. Elsevier, 2002.
[7] Kalliadasis S, Ruyer-Quil C, Scheid B, et al. Falling liquid films[M]. Springer Science & Business Media, 2011.
[8] Squire H B. Investigation of the instability of a moving liquid film[J]. British Journal of Applied Physics, 1953, 4(6):167-169. DOI: 10.1088/0508-3443/4/6/302
[9] Brown D. A study of the behaviour of a thin sheet of moving liquid[J]. Journal of Fluid Mechanics, 1961, 10(2):297-305. DOI: 10.1017/S002211206100024X
[10] Lin S, Roberts G. Waves in a viscous liquid curtain[J]. Journal of Fluid Mechanics, 1981, 112:443-458. DOI: 10.1017/S0022112081000505
[11] Lin S P. Breakup of liquid sheets and jets[M]. England:Cambridge University Press Cambridge, England, 2003.
[12] Sirignano W A, Mehring C. Review of theory of distortion and disintegration of liquid streams[J]. Progress in Energy & Combustion Science, 2000, 26(4):609-655. http://www.academia.edu/21244677/Review_of_theory_of_distortion_and_disintegration_of_liquid_streams
[13] Pritchard W. Instability and chaotic behaviour in a free-surface flow[J]. Journal of Fluid Mechanics, 1986, 165:1-60. DOI: 10.1017/S0022112086002987
[14] Mackowiak J. Fluid dynamics of packed columns[M]. London, New York:Springer, 2010.
[15] Fair J R, Seibert A F, Behrens M, et al. Structured packing performance experimental evaluation of two predictive models[J]. Industrial & Engineering Chemistry Research, 2000, 39(6):1788-1796. https://www.researchgate.net/publication/231394040_Structured_Packing_PerformanceExperimental_Evaluation_of_Two_Predictive_Models
[16] Pavlenko A, Pecherkin N, Chekhovich V, et al. Hydrodynamics in falling liquid films on surfaces with complex geometry[J]. Microgravity Science and Technology, 2009, 21(1):207-213. DOI: 10.1007/s12217-009-9157-1
[17] Yao Y, Pavlenko A, Volodin O. Effects of layers and holes on performance of wire mesh packing[J]. Journal of Engineering Thermophysics, 2015, 24(3):222-236. DOI: 10.1134/S1810232815030042
[18] Kolev N, Kralev B, Kolev D. Gas side controlled mass transfer in a new packing with stamped horizontal lamellae operating at extremely low liquid loads[J]. Chemical Engineering and Processing:Process Intensification, 2013, 63:44-49. DOI: 10.1016/j.cep.2012.07.004
[19] 戴干策, 胡剑光, 于建国, 等. 适应粘性吸收剂的有壁与无壁液膜交替的规整填料: 中国, 201210234334. 6. 2013-07-10. Dai G C, Hu J G, Yu J G, et al. Wall liquid film and wall-free liquid film alternate structured filler adapting to viscous absorbent:China, 201210234334.6. 2013-07-10.
[20] Hu J, Liu J, Yu J, et al. CO2 absorption into highly concentrated dea solution flowing over a vertical plate with rectangular windows[J]. International Journal of Greenhouse Gas Control, 2013, 19:13-18. DOI: 10.1016/j.ijggc.2013.08.007
[21] 胡剑光, 刘佳特, 袁猛, 等.新型垂直板规整填料流体力学及传质性能[J].化工学报, 2014, 65(1):116-122. http://www.cnki.com.cn/Article/CJFDTOTAL-HGSZ201401019.htm Hu J G, Liu J T, Yuan M, et al. Hydrodynamics and mass transfer characteristics of a novel vertical-sheet structured packing[J]. Journal of Chemical Industry and Engineering, 2014, 65(1):116-122. http://www.cnki.com.cn/Article/CJFDTOTAL-HGSZ201401019.htm
[22] 王良生, 戴干策.圆盘反应器成膜性的持液量研究[J].化学反应工程与工艺, 2000, 16(2):127-135. http://www.cnki.com.cn/Article/CJFDTOTAL-HXFY200002005.htm Wang L S, Dai G C. A study of film forming and hold-up in a rotating disc-ring reactor[J]. Chemical Reaction Engineering and Technology, 2000, 16(2):127-135. http://www.cnki.com.cn/Article/CJFDTOTAL-HXFY200002005.htm
[23] 邓斌, 戴干策.圆盘反应器液膜表面更新数值模拟[J].化工学报, 2015, 66(4):1407-1416. DOI: 10.11949/j.issn.0438-1157.20141688 Deng B, Dai G C. Numerical simulation of surface renewal frequency on vertically rotating disc[J]. Journal of Chemical Industry and Engineering, 2015, 66(4):1407-1416. DOI: 10.11949/j.issn.0438-1157.20141688
[24] Hu J, Yang X, Dai G. Numerical investigation on hydrodynamics of vertically confined free film[J]. The Canadian Journal of Chemical Engineering, 2016, 94(2):340-348. DOI: 10.1002/cjce.v94.2
[25] Oliver J F, Huh C, Mason S G. Resistance to spreading of liquids by sharp edges[J]. Journal of Colloid and Interface Science, 1977, 59(3):568-581. DOI: 10.1016/0021-9797(77)90052-2
[26] Brunet P, Flesselles J-M, Limat L. Dynamics of a circular array of liquid columns[J]. The European Physical Journal B, 2007, 55(3):297-322. DOI: 10.1140/epjb/e2007-00057-y
[27] Schmuki P, Laso M. On the stability of rivulet flow[J]. Journal of Fluid Mechanics, 1990, 215:125-143. DOI: 10.1017/S0022112090002580
[28] Ruan B, Jacobi A M, Li L. Effects of a countercurrent gas flow on falling-film mode transitions between horizontal tubes[J]. Experimental Thermal and Fluid Science, 2009, 33(8):1216-1225. DOI: 10.1016/j.expthermflusci.2009.07.009