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

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

下降液膜在逆向流动空气作用下的换热

以下一代核电站反应堆安全壳非能动冷却为背景，对竖直管内过冷下降液膜在逆向空气．水蒸汽流动作用下的换热提出了一个分析模型，对液膜采用边界层近似理论，对空气．水蒸汽混合气采用由热质传递比拟理论引出的关联式进行分析，考察了各种条件下液膜的厚度、温度、吸收热量与蒸发热量沿流向的变化，与相关研究进行了对比，表明了模型的合理性。     

矩形液池内热毛细对流的流动不稳定性

热毛细对流及其不稳定性是微重力流体科学研究的重要内容. 对该问题的研究不仅有利于人们对微重力环境下流体行为、对流不稳定性和湍流转捩过程等基础物理现象的进一步认识,而且也将促进晶体生长、薄膜制备等空间和地面高新技术的发展. 实验研究了矩形液池中浅液层在水平温度梯度作用下产生的热毛细对流及其稳定性. 实验中,成功地利用PIV (particle image velocimetry) 技术对1mm2/s 硅油液层内的浮力热毛细对流流场结构进行了大量观测. 结果表明,液层中的流场结构经历了多种状态的转变,该过程会受到液层厚度的影响. 当液层厚度较小时,比如当d=2:8mm 时,随着液池两端温差的增大,液层中的流场结构会经历单胞对流到双胞对流再到多胞对流的转变,到达多胞对流状态之后,继续增大温差,对流涡胞的数量会有所减少,而当温差进一步增大到一定程度以后,整个液层转变为三维非定常流动;当液层厚度较大时,比如当d=4:5mm 时,随着温差的增大,流动模式的转变主要体现在水平截面流场截面上面,当温差增大到一定程度以后,在靠近高温端的附近区域会出现具有明显三维效应的"梭形结构",该梭形结构的尺寸随着温差的增大而增长,并在温差超过某个临界值时失去对称性,整个液层转变为三维非定常流动.     

微重力下液封液桥内振荡热毛细对流的数值模拟

对微重力条件下液封液桥(流体采用KF-96硅油和FC-70氟化物,桥高1.4~4毫米,内层直径2或3毫米,外层直径5毫米)和单层液桥(流体采用KF-96硅油,桥高1~1.6毫米,直径2毫米)内的热毛细对流进行了数值模拟,模拟条件与Majima等的实验条件相同。从理论上证实了大温差条件下将出现振荡热毛细对流,确定了发生振荡的临界条件并与相关实验结果进行了对比,同时还计算了振荡频率。     

圆管内大液气比下液雾在高温气流中蒸发特性的研究

本文获得了液雾在过热蒸发状态下的液气两相能量方程,给合运动学方程、粒径分布及轨道模型,组成了描述液雾的完整数学方程组。提出将时间变量隐含在空间变量中来克服由于大液气比引起的“物性漂移”,考虑了壁面捕获。计算结果表明:大滴径粒子与小滴径粒子,参数变化具有不同的特征,并对气相温度场的模型计算值与实验值进行了比较.     

大尺寸液桥热毛细对流失稳性地面实验研究

开展大尺寸液桥浮力-热毛细对流地面实验, 探究流场转捩的临界条件及临界状态附近的流动情况. 通过粒子图像测速方法(PIV) 获得流体速度场, 研究液桥内部定常和转捩后的流场结构以及流体运动规律;并用红外热像仪测量液桥自由面温度分布, 研究流体流动的时空演化和温度振荡. 实验发现大尺寸半浮区液桥浮力-热毛细对流临界值与几何参数有关, 在大普朗特(Prandtl) 数情况下, 流场存在由稳定态向不稳定态再到混沌的转捩过程, 在临界马兰哥尼(Marangoni) 数附近, 流场内会出现行波现象, 流动模式也会随高径比的变化而发生变化;当继续增大马兰哥尼数, 流动会进入混沌状态.      

过冷下降液膜破裂的分析

用简化的横向形变模型和干斑模型对过冷受热下降液膜的破裂进行分析 ,取消了膜内温度为线性分布的假设 ,考虑了壁面尺度和液固配合的影响 ,提出了形变区长度和破裂膜厚的表达式 ,讨论了分析的合理性     

双自由面溶质?热毛细液层的不稳定性

溶质?热毛细对流是流体界面的浓度和温度分布不均导致的表面张力梯度驱动的流动, 它主要存在于空间微重力环境、小尺度流动等表面张力占主导的情况中, 例如晶体生长、微流控、合金浇筑凝固、有机薄液膜生长等. 对其流动进行稳定性分析具有重要意义. 本文采用线性稳定性理论研究了双自由面溶质?热毛细液层对流的不稳定性, 得到了两种负毛细力比(η)下的临界Marangoni数与Prandtl数(Pr)的函数关系, 并分析了临界模态的流场和能量机制. 研究发现: 溶质?热毛细对流和纯热毛细对流的临界模态有较大的差别, 前者是同向流向波、逆向流向波、展向稳态模态和逆向斜波, 后者是逆向斜波和逆向流向波. 在Pr较大时, Pr增加会降低流动稳定性; 在其他参数下, Pr增加会增强流动稳定性. 在中低Pr, 溶质毛细力使流动更加不稳定; 在大Pr时, 溶质毛细力的出现可能使流动更加稳定; 在其他参数下, 溶质毛细力会减弱流动稳定性. 流动稳定性不随η单调变化. 在多数情况下, 扰动浓度场与扰动温度场都是相似的. 能量分析表明: 扰动动能的主要能量来源是表面张力做功, 但其中溶质毛细力和热毛细力做功的正负性与参数有关.      

气体速度对液膜在预膜板表面流动形态影响的LBM模拟

为了探究气流剪切作用对航空燃油在气动雾化喷嘴预膜板上流动形态的影响,首先对基于相场理论的两相流格子Boltzmann模型进行修正,并通过经典算例验证了修正后模型的准确性和可靠性.随后利用该模型模拟了同向气流驱动下液膜在水平预膜板表面上的流动,分析了气流速度对液膜流动形态的影响规律.研究表明,该模型可准确追踪具有大密度比...     

环形腔内双层薄液层热毛细对毛细对流的渐近解

为了解水平温度梯度作用下环形腔内双层薄液层热毛细对流的基本特性,采用渐近线方法 获得了热毛细对流的近似解. 环形腔外壁被加热,内壁被冷却,上、下壁面绝热. 结果表明, 当环形腔宽度与内半径比趋于零时,环形腔退化为矩形腔,所得到的主流区速度场和温度场 的表达式演化为Nepomnyashchy 等得到的矩形腔内的结果;与数值模拟结果的比较 发现,在主流区渐近解与数值解吻合较好.     

Stability of a viscous liquid film flowing down a periodic surface

The paper is devoted to a theoretical analysis of linear stability of the viscous liquid film flowing down a wavy surface. The study is based on the Navier–Stokes equations in their full statement. The developed numerical algorithm allows us to obtain pioneer results in the stability of the film flow down a corrugated surface without asymptotic approximations in a wide range over Reynolds and Kapitsa’s numbers. It is shown that in the case of moderate Reynolds numbers there is a region of the corrugation parameters (amplitude and period) where all disturbances decay in time and the wall corrugation demonstrates a stabilizing effect. At the same time, there exist corrugation parameters at which the steady-state solution is unstable with respect to perturbations of the same period as the period of corrugation. In this case the waveless solution cannot be observed in reality and the wall corrugation demonstrates a destabilizing effect.     

Wave regimes on a nonisothermal film of a viscous liquid flowing down a vertical plane

A thin film flow of a viscous liquid flowing down a vertical wall in the field of the gravity force is studied. The values of temperatures on the solid wall and on the free surface are constant. The viscosity and thermal diffusivity are functions of temperature. An equation that describes the evolution of surface disturbances is derived for small flow rates in the long-wave approximation. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 2, pp. 89–97, March–April, 2008.     

Stability of conducting liquid flowing down an inclined plane at moderate Reynolds number in the presence of constant electromagnetic field

The stability of a conducting viscous film flowing down an inclined plane at moderate Reynolds number in the presence of electromagnetic field is investigated under induction-free approximation. Using momentum integral method a non-linear evolution equation for the development of the free surface is derived. The linear stability analysis of the evolution equation shows that the magnetic field stabilizes the flow whereas the electric field stabilizes or destabilizes the flow depending on its orientation with the flow. The weakly non-linear study reveals that both the supercritical stability and subcritical instability are possible for this type of thin film flow. The influence of magnetic field on the different zones is very significant, while the impact of electric field is very feeble in comparison.     

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

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

Influence of the capillarity on a creeping film flow down an inclined plane with an edge

Summary In creeping flows of thin films, the capillarity can play a dominant role. In this paper, the creeping film flow down an inclined plane with an edge is considered. The influence of the capillarity on the velocity and the film surface is studied analytically, numerically and experimentally. Received 12 April 1999; accepted for publication 9 May 1999     

Study of gas-sheared liquid film in horizontal rectangular duct using high-speed LIF technique: Three-dimensional wavy structure and its relation to liquid entrainment

The flow of a liquid film sheared by high velocity gas stream in a horizontal rectangular duct was investigated using a high-speed laser-induced fluorescence technique. Measurements of local film thickness were resolved in both longitudinal and transverse coordinates with high spatial and temporal resolution. It was found that the generation of fast and slow ripples by the disturbance waves was qualitatively the same as it was observed earlier in completely different conditions. The transverse size and curvature of the disturbance waves and ripples were measured. A relationship between the three-dimensional structure of ripples on top of disturbance waves and the two mechanisms of liquid entrainment, known as ‘bag break-up’ and ‘ligament break-up’, is proposed.     

Free- and forced-convection film condensation from a horizontal elliptic tube with a vertical plate and horizontal tube as special cases

A simple mathematical model is developed for the study of the mixed-convection film condensation with downward flowing vapors onto a horizontal elliptic tube. Analytical analysis for both the local condensate film thickness and heat transfer characteristics under simultaneous effects of interfacial vapor shear and pressure gradient has been performed by adopting a unified geometry parameter, eccentricity e. The present results for two limit cases, e = 0 (circular tube) and e = 1.0 (vertical plate) are in an excellent agreement with the earlier works. For very slow vapor flow, the present result for dimensionless mean heat transfer coefficient reduces to the same form as in the earlier works, , whose value is 0.728 for e = 0 and 0.943 for e = 1.0. As for very fast vapor flow, the dimensionless mean heat transfer coefficient, increase with increasing eccentricity under the effects of pressure gradient caused by potential flow and surface tension.     

Effects of air-flow and evaporation on the development of thin liquid film over a rotating annular disk

Axisymmetric flow of thin pure liquid film on a spinning horizontal annular disk is studied under the action of air shear at the liquid–air interface and evaporation. The non-linear evolution equation that is obtained by singular perturbation method is solved analytically, for small Reynolds number, by using the method of characteristic and numerically by the use of Newton–Kantorovich method for any Reynolds number. Font breakdown time and its location from the center of the disk is predicted both by analytically and numerically. The result shows that the thinning of the initial film increases as air stress increase, same result is also escalated in presence of evaporation.     

The drainage and rupture of a non-foaming liquid film formed upon bubble impact with a free surface

Experiments are performed to measure the thickness of a thin liquid film formed between a free surface and the apex of a bubble which is approaching at its terminal velocity. Measurements are made using bubbles of < 1 mm dia both in distilled water and alcoholic solutions (methanol, ethanol). The results of the experiments show that the models based on the lubrication approximation fail to predict the initial stage of drainage. A better agreement is obtained by modifying the model proposed by Chesters.     

Characteristics of the critical heat flux for downward flow in a vertical tube at low flow rate and low pressure conditions

The characteristics of the critical heat flux (CHF) for downward flow were studied experimentally with an Inconel 600 circular tube test section in a water test loop at low-flow rate (0 200 kg/m²s) and low-pressure (0.1 0.7 MPa) conditions. The attention was given to the effects of upstream conditions—upper plenum and inlet throttling. Two totally different kinds of CHF behaviors were observed. It seems appropriate to interpret them as flooding-type CHF and dryout in annular flow. The CHF in downward flow may vary from extremely unstable flow CHF as low as near the flooding CHF value to stable flow CHF as high as that of upflow, depending on the upstream conditions of the test section. The CHF correlation by Mishima and that by Weber were proposed for the presentation of the lower and upper limits of the CHF for downward flow in a vertical tube at low-flow rate and low-pressure conditions.