液桥中不同Prandtl数流体的热毛细对流

本文研究从上面加热的液桥中不同 pr 数流体的热表面张力驱动对流。在 Ma 数相同的条件下,不同 pr 数流体液桥中的温度分布和流场结构定性相似,但定量结果不同。小pr数(pr<<1)流体液桥中的粘性边界层远小于热边界层,最大流函数所在位置向冷端偏移,有较大的流动速度。结果表明,Ma 数是描述这种流动的敏感参数。     

Deformation of the Free Surface in a Moving Locally-Heated Thin Liquid Layer

Liquid film flow on a vertical surface is studied experimentally and theoretically under the determining influence of the thermocapillary forces. In the two-dimensional steady-state case the shape of the film surface is calculated numerically within the thin layer approximation with allowance for the temperature dependence of the viscosity of the liquid and redistribution of the heat flux in the heating element. A local heat source was used in the experiments to produce temperature gradients up to 10 K/mm on the liquid surface. The film thickness was determined by means of the schlieren method with reflection. The relative thickness of the roller in the upper heater edge zone, characteristic of the formation of regular structures, is measured. The thickness is h/h ₀=1.32 ±0.07, which agrees satisfactorily with the results of numerical calculations.     

Thermocapillary Vortices Induced by a Light Beam near a Bubble Surface in a Hele-Shaw Cell

This paper studies thermocapillary vortices induced by local heating of a bubble surface in a Hele-Shaw cell by a light beam. It is found that the vortex rotation frequency and its depth depend on the distance from the light-beam projection onto the layer to the bubble boundary. The surface velocity of the thermocapillary flow is calculated using the balance of the near-surface and return flows of the thermocapillary vortex and the equality of capillary and dynamic pressures. It is shown that a decrease in the surface velocity and the vortex rotation frequency with increase in the distance from the light beam to the bubble surface is due to a decrease in the temperature gradient between the illuminated and cold poles of the bubble.__________Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 5, pp. 93–99, September–October, 2005.     

半浮区液桥热毛细振荡流

采用非定常、三维直接数值模拟方法研究大Pr数半浮区液桥热毛细对流从定常流向振荡流的过渡过程．文中详细描述了热毛细振荡流的起振和振荡特征,给出了液桥横截面上振荡流的流场和温度分布．在地面引力场条件下计算的结果与地面实验的结果进行比较,得出液桥水平截面上的流场和温度分布图样以一定的速度旋转,自由表面固定点处流体的环向流速正、负交替变化的一致结论．     

半浮区液桥热毛细对流速度场特征

液桥流场的速度分布是液桥表面张力梯度驱动对流的一个重要研究内容。通过地面模拟微重力环境的液桥实验以及数值计算,研究了外加温差、液桥腰径等对流场速度分布的影响,实验结果与数值计算结果相当符合。结果还表明,当液桥腰径减小时,流场的结构发生了变化,涡心由一个转化成两个。     

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

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

A model of film deformation and rupture under the action of thermocapillary forces

A model of the deformation of a free weightless liquid film with rims fixed at a plane contour and subject to the action of thermocapillary forces is investigated. The film free-surface temperature is assumed to be a known function of the coordinates and time. The equation describing the film thickness evolution in the long-wave approximation is of the second order with respect to time and of the fourth order with respect to the longitudinal coordinates. A model plane nonstationary problem is calculated, making it possible to estimate the lifetime of the film as a function of the rate of variation of the temperature profile curvature on the free boundary. If the characteristic time of the temperature variation is large as compared with the film’s natural oscillation period, the problem can be considered in the quasi-steady approximation, which is studied in detail for the plane and axisymmetric cases. The critical values of the temperature profile curvature, for which the film thickness on the symmetry axis vanishes, resulting in the rupture of the film, are calculated.     

Experimental Investigation of Thermocapillary Convection Induced by a Local Temperature Inhomogeneity near the Liquid Surface. 2. Radiation-Induced Source of Heat

The structure and stability of a thermocapillary flow from a concentrated source of heat induced by radiation and located near the free surface of the liquid filling a deep reservoir are experimentally studied. Absolute stability of this flow for all depths and powers of the heat source is found. The profile of the liquid surface near the heat spot is considered experimentally.     

热毛细对流温度场ESPI检测研究

在地面上用二维液体盒研究了热壁下气泡周围液体中的热毛细对流现象,并应用电子散斑干涉技术（ESP）对热毛细对流温度场进行了实时检测研究,给出了部分典型的实验结果。     

Instabilities in two-liquid layers subject to a horizontal temperature gradient

We study theoretically the stability of two superposed fluid layers heated laterally. The fluids are supposed to be immiscible, the interface undeformable and of infinite horizontal extension. Combined thermocapillary and buoyancy forces give rise to a basic flow when a temperature difference is applied. The calculations are performed for a melt of GaAs under a layer of molten B₂O₃, a configuration of considerable technological importance. Four different flow patterns and five temperature configurations are found for the basic state in this system. A linear stability analysis shows that the basic state may be destabilized by oscillatory motions leading to the so-called hydrothermal waves. Depending on the relative height of the two layers these hydrothermal waves propagate parallel or perpendicular to the temperature gradient. This analysis reveals that these perturbations can alter significantly the liquid flow in the liquid-encapsulated crystal growth techniques. PACS 47.20.Dz, 47.20.Bp, 47.54.+r, 47.27.Te, 44.25.+f, 47.20.Ma     

热毛细对流温度场全息干涉检测研究

本文采用Ｈｅｌｅ－Ｓｈａｗ盒体,在地面上模拟二维热毛细对流,应用全息干涉技术对热毛细对流温度场进行了测试研究,为了获得感兴趣的热毛细对流温度增量分布,文中提出了将背景温度分布条纹与热毛细对流温度增量分布条纹相分离的方法。     

环形浅液池内热流体波的本质特征

为了了解环形浅液池内热毛细对流的特征,利用有限容积法对外壁受热、内壁冷却、厚度为1 mm的环形浅液池内硅油的热毛细对流进行了三维数值模拟.结果表明,当Marangoni(Ma)数小于临界值时,随着Ma数的增加,内、外壁附近的温度梯庹上升,稳定的二维轴对称流动的径向速率增加;超过临界Ma数后,漉动转化为三维振荡对流,形成热流体波.沿径向运动的同时,伴随着热流体波的传播流体质点成对地绕顺时针和逆时针方向旋转.热流体波的周向传播速度较快,而流体质点的周向速度很小.分析发现,热流体波为对数螺线形波纹,其传播角为常数;随着Ma数的增大,传播角增大.     

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 prediction of turbulent thermocapillary convection in superposed fluid layers with a free interface

The present paper introduces a new numerical method for predicting the characteristics of thermocapillary turbulent convection in a differentially-heated rectangular cavity with two superposed and immiscible fluid layers. The unsteady Reynolds form of the Navier–Stokes equations and energy equation are solved by using the control volume approach on a staggered grid system using SIMPLE algorithm. The turbulence quantities are predicted by applying the standard k–ε turbulence model. The level set formulation is applied for predicting the topological changes of the interface separating the two fluid layers and to provide an accurate and robust modeling of the interfacial normal and tangential stresses. The computational results obtained showed good agreement when compared with the previous experimental, numerical and analytical benchmark data for different validation cases in both laminar and turbulent regimes. The present numerical method is then applied to predict the velocity and temperature distribution in two immiscible liquid layers with undeformable interface for a wide range of Marangoni numbers. The laminar-turbulent transition is demonstrated by obtaining the turbulence features at high interfacial temperature gradient which is characterized by high Marangoni number. The effect of increasing Marangoni number on the interface dynamics in turbulent regime is also investigated.     

Direct numerical simulation of thermocapillary flow based on the Volume of Fluid method

A numerical method for direct simulation of thermal Marangoni effects at dynamically deformable interface of two-phase incompressible fluids is developed. The approach is based on the Volume of Fluid (VOF) method with special focus on the numerical treatment of the temperature surface gradient because of its decisive role as the driving force of the flow. The surface gradient calculation is split into computing its length and direction in order to satisfy the correct thermal boundary condition at the interface without losing mobility of the interface. The method is applied to three different types of thermocapillary flow, namely thermocapillary migration of a droplet in an ambient fluid with linear temperature gradient, thermocapillary convection in a liquid layer under linear temperature gradient along the interface, and Marangoni convection due to Bénard–Marangoni instability. In the first case, different aspects of the dynamics of the migration are considered for validation such as the terminal migration velocity, the initial acceleration and quantification of the wall effects. The simulation also considers high convective heat transfer and covers a wide range of Marangoni numbers up to 5000, where good agreement with both theoretical and experimental results is achieved. In the second case, the convection velocity in the liquid layer is compared with an analytical result. In the final application, pattern formation due to the Bénard–Marangoni instability in a liquid layer in square geometry of small aspect ratio is investigated for realistic Biot number and dynamically deformable fluid interface. The results show good agreement with experiments from literature, where our numerical simulation also predicts cell pattern for a particular aspect ratio which is outside the limitation of the above cited experimental work.     

磁场对双扩散液层热毛细对流的影响

通过数值模拟的方法对磁场作用下的双扩散液层热毛细对流进行了研究, 模型中同时考虑了热毛细效应和溶质毛细效应的存在. 研究结果显示, 外部磁场能够有效削弱液层内热毛细对流的强度, 改变热毛细对流的对流结构; 随着磁场强度的增大, 液层内热毛细对流的对流强度逐渐减小, 热质传递过程中扩散效应逐渐得到增强; 最终, 溶质浓度沿水平方向呈梯度分布. 因此, 当磁场强度足够大时能够实现晶体生长中所需的纯扩散条件.     

矩形液池热毛细对流转捩途径研究

主要研究矩形液池热毛细对流的分岔转捩. 通过测量流体内部温度振荡情况, 详细研究了热毛细对流的转捩过程和转捩途径. 实验发现, 矩形液池热毛细对流的转捩过程依次经历了定常、规则振荡、不规则振荡的阶段. 对于不同普朗特数的硅油在不同长高比情况下, 通向混沌的途径不同. 在转捩过程中, 随着温差的增加, 普朗特数在16 (1cSt) 以下和普朗特数为25 (1.5cSt)、长高比为26 的硅油热毛细对流主要以准周期分岔的转捩方式为主;而普朗特数为25 以上的则以倍周期分岔的转捩方式为主;两种分岔有时还会伴随有切分岔形式的出现.实验中还观察到了表面波动和对流涡胞振荡等现象.      

Numerical study on the effect of initial flow velocity on liquid film thickness of accelerated slug flow in a micro tube

Numerical simulation of air–water slug flows accelerated from steady states with different initial velocities in a micro tube is conducted. It is shown that the liquid film formed between the gas bubble and the wall in an accelerated flow is significantly thinner than that in a steady flow at the same instantaneous capillary number. Specifically, the liquid film thickness is kept almost unchanged just after the onset of acceleration, and then gradually increases and eventually converges to that of an accelerated flow from zero initial velocity. Due to the flow acceleration, the Stokes layer is generated from the wall, and the instant velocity profile can be given by superposition of the Stokes layer and the initial parabolic velocity profile of a steady flow. It is found that the velocity profile inside a liquid slug away from the bubble can be well predicted by the analytical solution of a single-phase flow with acceleration. The change of the velocity profile in an accelerated flow changes the balance between the inertia, surface tension and viscous forces around the meniscus region, and thus the resultant liquid film thickness. By introducing the displacement thickness, the existing correlation for liquid film thickness in a steady flow (Han and Shikazono, 2009) is extended so that it can be applied to a flow with acceleration from an arbitrary initial velocity. It is demonstrated that the proposed correlation can predict liquid film thickness at Re < 4600 within the range of ±10% accuracy.     

Heat transfer in a thin layer of a viscous heavy noncompressible liquid under wavy flow conditions

The article describes a method for calculating the flow of heat through a wavy boundary separating a layer of liquid from a layer of gas, under the assumption that the viscosity and heat-transfer coefficients are constant, and that a constant temperature of the fixed wall and a constant temperature of the gas flow are given. A study is made of the equations of motion and thermal conductivity (without taking the dissipation energy into account) in the approximations of the theory of the boundary layer; the left-hand sides of these equations are replaced by their averaged values over the layer. These equations, after linearization, are used to determine the velocity and temperature distributions. The qualitative aspect of heat transfer in a thin layer of viscous liquid, under regular-wavy flow conditions, is examined. Particular attention is paid to the effect of the surface tension coefficient on the flow of heat through the interface.Notation x, y coordinates of a liquid particle - t time - v and u coordinates of the velocity vector of the liquid - p pressure in the liquid - c_v, , T,, andv heat capacity, thermal conductivity coefficient, temperature, density, and viscosity of the liquid, respectively - g acceleration due to gravity - surface-tension coefficient - c phase velocity of the waves at the interface - T_w wall temperature - h₀ thickness of the liquid layer - u₀ velocity of the liquid over the layer Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 147–151, July–August, 1970.     

Pure steam condensation model with laminar film in a vertical tube

A new physical model for calculating the liquid film thickness and condensation heat transfer coefficient in a vertical condenser tube is proposed by considering the effects of gravity, liquid viscosity, and vapor flow in the core region of the flow. To estimate the velocity profile in the liquid film, the liquid film was assumed to be in Couette flow forced by the interfacial velocity at the liquid–vapor interface. For simplifying the calculation procedures, the interfacial velocity was estimated by introducing an empirical power-law velocity profile. The resulting film thickness and heat transfer coefficient from the model were compared with the experimental data and the results obtained from the other condensation models. The results demonstrated that the proposed model described the liquid film thinning effect by the vapor shear flow and predicted the condensation heat transfer coefficient from experiments reasonably well.