微扰动对液体螺旋槽机械密封性能的影响

针对非接触式液体螺旋槽机械密封,分析了周期性轴向微扰动和角微扰动对密封环端面间液膜厚度的影响规律,建立了单周期螺旋槽液膜模型,采用数值法求解时变雷诺方程,研究了微扰动对液膜承载力和泄漏率等密封特性参数的影响.结果表明,膜厚对液膜承载力和泄漏率的影响显著,膜厚增加,液膜承载力减小,泄漏率增大.液膜承载力、泄漏率的变化幅度和频率主要受到轴向扰动的影响.角扰动造成沿圆周方向分布的液膜承载力不均匀而使液膜稳定性变差.     

高超声速边界层液膜演化过程和冷却机理研究

高超声速液膜冷却技术是通过一系列狭缝或孔洞压出冷却工质,在飞行器表面边界层形成一层低温冷却膜,阻止高超声速气流对飞行器的气动加热.其作为一种主动冷却方式在高超声速飞行器表面热防护有着巨大的应用潜力.文章采用数值方法,结合VOF模型,研究25 km飞行高度和Ma=5气流条件下的液膜铺展情况,并通过不同冷却工质的入射速度、角度、表面张力和黏性系数条件,讨论了液膜在平板上的演化过程和冷却机理.结果表明,在气流作用下,液膜向壁面下游发展,液膜的存在导致边界层分离,连续液膜会在一定位置断裂为液块,然后进一步破碎为液滴.入射条件和液体性质的改变,会影响液膜沿流向的发展,具体表现在连续液膜断裂点的位置和连续液膜的厚度.在所设定的计算域内,壁面热流降低了80%～95%,液膜对壁面的冷却效率随着液膜形态的变化而变化.     

STABILITY OF ULTRATHIN LIQUID FILM EVOLUTION WITH SURFACTANT

针对固体基底上厚度小于100 nm的含活性剂超薄液膜演化过程, 基于润滑理论推导出包含分离压影响的液膜厚度和活性剂浓度的演化方程, 采用正则模态法导出了描述液膜线性稳定性的特征方程, 分析了多个特征参数对线性稳定性的影响, 数值模拟了液膜厚度和活性剂浓度演化历程, 对比了模拟所得非线性结果与线性分析预测结果的一致性.结果表明:范德华力具有促进扰动增长的作用, 较强的玻恩斥力促使扰动衰减, 使液膜趋于稳定;较小的毛细力数易使液膜凹陷处发生二次失稳, 并最终导致去润湿现象发生;液膜厚度和溶于液膜内部的活性剂浓度初值越大, 液膜稳定性越强, 液膜表面活性剂浓度影响则相反;增大吸附系数不利于液膜稳定性.     

考虑空化效应的螺旋槽机械密封液膜动力学特性研究

考虑液膜空化效应的影响，研究螺旋槽液体润滑机械密封的动力学特性. 基于液体润滑理论和小扰动法，建立了考虑液膜空化的密封微扰膜压控制方程，采用有限单元法对端面液膜三自由度微扰下的液膜刚度和阻尼系数进行了数值求解，分析了不同参数对液膜密封动力系数的影响. 螺旋槽深度在10 μm左右、槽坝比在0.75左右、槽宽比在0.4左右，螺旋角在9°左右时液膜具有最大的轴向和角向刚度系数. 螺旋槽深度在5 μm左右、槽宽比在0.6左右、螺旋角在20°左右时，两角向交叉阻尼绝对值最大. 初始偏角的存在使密封压力呈现非对称性，从而使两角向动力系数绝对值不再相等. 液膜轴向刚度k_zz在数量级上远大于其余液膜刚度值，液膜轴向阻尼d_zz、角向阻尼d_αα和d_ββ远大于液膜其余阻尼值且随着转速和间隙的增大而减小.      

高超声速溢流冷却实验研究

高超声速溢流冷却是一种新型的飞行器热防护方法,基本思想为:在高热流区布置溢流孔,控制冷却液以溢流方式流出,之后通过飞行器表面摩阻作用展布为液膜,形成热缓冲层以降低飞行器表面热流.目前,溢流冷却技术还处于探索阶段,实现工程应用前还需开展大量的实验验证和机理研究工作.本文首次开展溢流冷却的实验研究工作,采用热流测量、液膜厚度测量及液膜流动特性观测技术,搭建了完善的溢流冷却风洞实验平台,对溢流冷却热防护性能和高超声速条件下液膜流动规律进行了初步研究.研究表明:(1)高超声速流场中通过溢流能够在飞行器表面形成液膜并有效隔离外部高温气流,可降低飞行器表面热流率;(2)楔面上的液膜前缘流动是一个逐渐减速的过程,增加冷却液流量液膜厚度变化不明显,但液膜前缘运动速度增大;(3)液膜层存在表面波,在时间和空间方向发生演化,导致液膜厚度的微弱扰动;(4)液膜层存在横向展宽现象,即液膜层宽度大于溢流缝宽度.原因是液膜层与流场边界层条件不匹配,存在压力梯度,迫使冷却液向低压区流动,从而展宽液膜层,并且流量越高,横向展宽现象越明显.     

高超声速溢流冷却实验研究

高超声速溢流冷却是一种新型的飞行器热防护方法,基本思想为:在高热流区布置溢流孔,控制冷却液以溢流方式流出,之后通过飞行器表面摩阻作用展布为液膜,形成热缓冲层以降低飞行器表面热流. 目前,溢流冷却技术还处于探索阶段,实现工程应用前还需开展大量的实验验证和机理研究工作. 本文首次开展溢流冷却的实验研究工作,采用热流测量、液膜厚度测量及液膜流动特性观测技术,搭建了完善的溢流冷却风洞实验平台,对溢流冷却热防护性能和高超声速条件下液膜流动规律进行了初步研究. 研究表明:(1) 高超声速流场中通过溢流能够在飞行器表面形成液膜并有效隔离外部高温气流,可降低飞行器表面热流率;(2) 楔面上的液膜前缘流动是一个逐渐减速的过程,增加冷却液流量液膜厚度变化不明显,但液膜前缘运动速度增大;(3) 液膜层存在表面波,在时间和空间方向发生演化,导致液膜厚度的微弱扰动;(4) 液膜层存在横向展宽现象,即液膜层宽度大于溢流缝宽度. 原因是液膜层与流场边界层条件不匹配,存在压力梯度,迫使冷却液向低压区流动,从而展宽液膜层,并且流量越高,横向展宽现象越明显.      

管内液膜流动失稳现象的直接数值模拟

为分析液膜流动时垂直液膜的重力加速度分量和表面张力对扰动波振幅的影响,通过FTM模拟了液膜在重力作用下流动失稳过程。结果表明:重力分量和表面张力系数越大,扰动振幅衰减得越快;并且表面张力系数越小,扰动振幅波动的频率越高。此外还分析了促进液膜流动失稳的两个因素,即密度比和雷诺数Re。当增大密度比时,液膜会出现两种极端的变化:一种是液膜不再完全覆盖底面,会形成锲状液膜,随着速度变大和振幅脉动,当锲状液膜前锋出现新的扰动波时,新的扰动会与初始扰动发生干涉产生二次波动;另一种是增加密度比时,液膜将会在剪切力的作用下产生飞脱的现象,在Re数较大时,液膜流动失稳会在剪切力的作用下产生K-H(开尔文-赫姆霍兹)不稳定性,且初始扰动对形成开尔文-赫姆霍兹波浪具有诱导作用。     

旋流喷嘴中空旋转射流近区域流动的研究

从理论上分析了旋流喷嘴喷出的中空射流近区域的液膜的运动,在只考虑液膜表面张力的作用下,应用质量守恒和动量定理,建立了描述液膜运动的非线性常微分（积分）方程组,该方程组可以用数值方法方便地求解,结果表明,液体离开旋流喷嘴后在自由空间形成的液膜呈葫芦形状,其速度和液膜厚度等都周期性地变化。本结果是在液厝受拓动失称碎成液滴前的最基本运动状态,可以在射流的近区域内实验观察到,也是进一步从理论液膜破碎雾化过     

饱和蒸发和热毛细力作用下低雷诺数液膜在波纹竖壁上的流动

考虑表面蒸发压力和热毛细力作用情况下,对饱和蒸发状态下低雷诺数自由降落液膜在小波幅正弦型波纹壁面上的流动进行理论分析。对控制微分方程及边界条件进行量纲一化并引入流函数,对微分方程及边界条件进行摄动展开,得到了这种情况下液膜流动的简化分析模型,求出了近似解析解。讨论了壁面波纹、表面张力、蒸发压力、热毛细力对液膜流动的影响。研究表明:液膜的波动幅度随蒸发强度和热毛细力的增大而增大;液膜波动与壁面波纹的相位差随蒸发强度增大而增大,随热毛细力增大而减小。     

单个上升气泡与异质液膜相互作用的数值模拟

为了探讨在流体区域内气泡上升然后与异质液膜接触的过程中,气泡速度、形状随时间的变化,以及区域内液膜的运动特性。本文采用FTM方法模拟研究了在2D×D的空间范围内单个气泡从静止开始上升到与异质液膜接触过程中的运动特性。研究结果表明:气泡上升到与液膜一定距离后,会使液膜中间凸起变薄,并将液膜挤向两侧,气泡的竖直速度先增大后减小;液膜厚度不同时,气泡上升初期相同时刻的纵横比?基本相同;气泡上升后期相同时刻的纵横比随着液膜厚度的增厚而增大;不同表面张力下,相同时刻气泡纵横比随着表面张力的减小而减小;表面张力越小,气泡顶部和底端的压力差会越小,上升速度越小。     

Temperature pulsations and thermocapillary effects on the surface of a wavy heated liquid film

The temperature pulsations and wave characteristics in water film flow along a vertical plate with a heater are investigated. Using an infrared scanner, the temperature field on the film surface is measured for various heat flux densities on the heater. Experimental data on the variation of the temperature with time on a local segment of the liquid film surface during wave transmission are obtained. In the absence of a heat flux the data obtained are in good agreement with the results of other researchers for an isothermal liquid film. When the down-flowing liquid is heated, the thermocapillary forces lead to the formation of rivulets and a thin film between them. It is shown that in the inter-rivulet zone the relative wave amplitude increases due to the action of the thermocapillary forces.     

波纹竖板层流降膜蒸发传热的分析

对高粘度液体在等温正弦形波纹壁面上的自由降落与蒸发建立了摄动分析模型。得到了流动的分析解和蒸发传热的数值解。考察了壁面波纹的波幅和波数、液膜表面张力及贝克利数对流动与传热的影响,结果表明,加大波纹的波幅、适当选择波数、减小贝克利数可增强传热,而表面张力对蒸发传热的影响较小。     

Numerical analysis on flow dynamics and heat transfer of falling liquid films with interfacial waves

Flow dynamics and heat transfer of falling liquid films with interfacial waves flowing on a vertical plate have been studied with originally proposed numerical simulation method. To discretize basic equations a staggered grid fixed on a physical space is employed. A small amplitude disturbance generated at inflow boundary develops to a solitary wave which consists of a large amplitude roll wave and small amplitude capillary waves. Instantaneous streamwise velocity profiles at the wave crest and trough are very different from a laminar flow. A circulation flow occurs in the roll wave and it affects temperature distributions, especially the strong effect is observed for high Prandtl number liquids. The interfacial wave enhances the heat transfer by two kinds of effects which are a film thinning effect and a convection effect. The dominating effect depends on the Prandtl number. Received on 23 December 1998     

Measurement of water falling film thickness to flat plate using confocal chromatic sensoring technique

An experimental investigation of wavy water film falling down a flat plate has been carried out using confocal chromatic sensoring technique to determine the instantaneous and statistical characteristics of the film. The experiments involved three parameters: liquid feed mode, Reynolds number and plate inclination angle. The present time–average film thickness data is compared with the previous experimental and theoretical results showing a good agreement. A new correlation for the average film thickness is suggested. Our results show that the liquid feed mode has a vital influence on the film thickness and that the film thickness increases with Reynolds number and decreased plate inclination angle. The root–mean–square value of the film thickness fluctuations depends non-monotonically on the film Reynolds number. The corresponding mechanisms are analyzed.     

Heat and mass transfer analogy for condensation of humid air in a vertical channel

This study examines energy transport associated with liquid film condensation in natural convection flows driven by differences in density due to temperature and concentration gradients. The condensation problem is based on the thin-film assumptions. The most common compositional gradient, which is encountered in humid air at ambient temperature is considered. A steady laminar Boussinesq flow of an ideal gas–vapor mixture is studied for the case of a vertical parallel plate channel. New correlations for the latent and sensible Nusselt numbers are established, and the heat and mass transfer analogy between the sensible Nusselt number and Sherwood number is demonstrated. Received on 15 November 1999     

An experimental test of equations predicting the load bearing capacity of squeeze films of elastico-viscous liquid

A form of squeeze film apparatus was recently described in which the movement of one plate towards the other was simulated by the continuous volume generation of liquid over the plate area. The liquid exuded from a large number of holes in the lower plate surface and formed a “continous flow” version of squeeze film apparatus with no moving parts [1]. A later paper gave derivations of equations from which squeeze film load bearing capacity could be evaluated, taking into account viscous, inertial and normal stress effects in the liquid film [2].In order to find the total load in a squeeze film system, it was necessary to obtain the relationship between the first normal stress difference and shear rate for the liquid in use, using an experimental method. At high shear rates, the jet thrust method provided these data [3,4] and from them the load bearing capacity of squeeze films of hot, polymer-thickened oil were predicted [2].A more complete test of the method is possible with a highly elastic liquid because considerable load enhancement due to extra stress is present at moderate deformation rates in squeeze film systems [1,5,6,7]. Thus a 0.1 per cent aqueous polyacrylamide solution gives well-defined load enhancement and (quite independently) the jet thrust method gives the relationship between normal stress and shear rate from which predictions of load enhancement may be made. Furthermore, convergent nozzles may be used in the jet thrust apparatus [3] to measure the stress development in an elastic liquid which is being simulateneously sheared and stretched, a situation which more closely resembles the squeeze film case than that of steady shear.     

下降液膜层流换热发展段中的积分分析

对等温竖直壁面上层流液膜发展段的流动换热进行了理论分析 ,用积分法得到并求解了速度和温度边界层的积分方程式 ,求得了发展段的流动和换热特性的变化规律 ,讨论了进口流态和普朗特数对换热的影响。     

等热流密度加热垂直下降液膜的研究

利用等热流密度加热条件下降膜流动的三维模型方程进行线性稳定性分析和数值模拟。线性稳定性分析表明，模型方程在小到中等Reynolds数下都适用，并且流向不稳定性增长率随着Reynolds数和Marangoni数增加而增加，展向不稳定性增长率则随着Marangoni数增加而增加，随着Reynolds数增加而减小，流向和展向对扰动波数都存在一个不稳定区间。三维数值模拟表明，在等热流密度加热条件下，液膜在随机扰动的情况下最终会形成带孤立波的三维溪流状结构，液膜与气体的换热也因溪流状结构的出现而加强；在随机扰动的基础上引入占优势地位的展向最不稳定扰动会使得换热增强，液膜会提前破裂；在随机扰动的基础上引入占优势地位的流向最不稳定扰动时，液膜的换热会增强，但不会提前破裂；在随机扰动的基础上同时引入占优势地位的流向和展向最不稳定扰动时，换热会加强且液膜会提前破裂。     

Study of the thermocapillary effect on a wavy falling film using a fiber optical thickness probe

The wave flow of a water film down a vertical plate with a 150×150 mm heater has been experimentally studied. The effect of the heat flux on the film flow leads to the formation of periodically flowing rivulets and thin film between them due to the action of thermocapillary forces in spanwise direction. The local film thickness between rivulets is measured by means of a noncontact fiber optical probe. As the heat flux grows, the average film thickness continuously decreases but upon reaching about 50% of the initial thickness, the film spontaneously breaks down. It is found that the decrease of the wave amplitude between rivulets is caused by the reduction of the local Reynolds number and is in a qualitative agreement with the laws of the hydrodynamics for the isothermal case. That is, no appreciable effect of streamwise thermocapillary forces on the wave amplitudes is detected. The experimental results are in good agreement with recently published data obtained by the capacitance method.