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.     

Lattice Boltzmann study of the interaction between a single solid particle and a thin liquid film

The isothermal single-component multi-phase lattice Boltzmann method(LBM) combined with the particle motion model is used to simulate the detailed process of liquid film rupture induced by a single spherical particle.The entire process of the liquid film rupture can be divided into two stages.In Stage 1,the particle contacts with the liquid film and moves into it due to the interfacial force and finally penetrates the liquid film.Then in Stage 2,the upper and lower liquid surfaces of the thin film are driven by the capillary force and approach to each other along the surface of the particle,resulting in a complete rupture.It is found that a hydrophobic particle with a contact angle of 106.7° shows the shortest rupture duration when the liquid film thickness is less than the particle radius.When the thickness of the liquid film is greater than the immersed depth of the particle at equilibrium,the time of liquid film rupture caused by a hydrophobic particle will be increased.On the other hand,a moderately hydrophilic particle can form a bridge in the middle of the liquid film to enhance the stability of the thin liquid film.     

Simultaneous heat and mass transfer inside a vertical tube in evaporating a heated falling alcohols liquid film into a stream of dry air

A numerical study of the evaporation in mixed convection of a pure alcohol liquid film: ethanol and methanol was investigated. It is a turbulent liquid film falling on the internal face of a vertical tube. A laminar flow of dry air enters the vertical tube at constant temperature in the downward direction. The wall of the tube is subjected to a constant and uniform heat flux. The model solves the coupled parabolic governing equations in both phases including turbulent liquid film together with the boundary and interfacial conditions. The systems of equations obtained by using an implicit finite difference method are solved by TDMA method. A Van Driest model is adopted to simulate the turbulent liquid film flow. The influence of the inlet liquid flow, Reynolds number in the gas flow and the wall heat flux on the intensity of heat and mass transfers are examined. A comparison between the results obtained for studied alcohols and water in the same conditions is made.     

Formation of a secondary droplet over the crown upon the liquid film rupture under the action of a laser beam

The free surface of a liquid film exposed to a laser beam is deformed and suffers a rupture. Depending on the thermal load intensity and the thermal properties of the liquid the rupture can be accompanied by the formation of secondary droplets over the film crown. This process is investigated using a mathematical model describing the motion of the thin layer of a viscous nonisothermal liquid. The model is based on the two-dimensional Navier–Stokes equations. The boundary conditions at the film-gas and film-liquid interfaces necessary for the solution of these equations are derived in the explicit form. The results of the solution of model problems are presented.     

Simulation and measurement of 3D shear-driven thin liquid film flow in a duct

Three-dimensional flow behavior of thin liquid film that is shear-driven by turbulent air flow in a duct is measured and simulated. Its film thickness and width are reported as a function of air velocity, liquid flow rate, surface tension coefficient, and wall contact angle. The numerical component of this study is aimed at exploring and assessing the suitability of utilizing the FLUENT-CFD code and its existing components, i.e. Volume of Fluid model (VOF) along with selected turbulence model, for simulating the behavior of 3D shear-driven liquid film flow, through a comparison with measured results. The thickness and width of the shear-driven liquid film are measured using an interferometric technique that makes use of the phase shift between the reflections of incident light from the top and bottom surfaces of the thin liquid film. Such measurements are quite challenging due to the dynamic interfacial instabilities that develop in this flow. The results reveal that higher air flow velocity decreases the liquid film thickness but increases its width, while higher liquid flow rate increases both its thickness and width. Simulated results provide good estimates of the measured values, and reveal the need for considering a dynamic rather than a static wall contact angle in the model for improving the comparison with measured values.     

Gas entrainment by a liquid film falling around a stationary Taylor bubble in a vertical tube

Gas entrainment by a liquid film falling around a stationary Taylor bubble in a 0.1 m diameter vertical tube is studied experimentally with the purpose of validating a model formulated in an earlier phase of our research. According to this model for a fixed liquid velocity the gas entrainment should be proportional to the waviness of the film (its intermittency) and the wave height and inversely proportional to the film thickness. For Taylor bubble lengths ranging from 1D to 15D these film parameters have been measured with a Laser Induced Fluorescence technique. The gas entrainment has been determined from the net gas flux into the liquid column underneath the Taylor bubble by using data on gas re-coalescence into the rear of the Taylor bubble. These data are available for lengths ranging from 4.5D to 9D. The model results with the measured film characteristics compare well with the observed gas entrainment. The fact that the net gas flux becomes constant for long Taylor bubbles, whereas the wave height still increases, warrants further study.     

Liquid film break-up in a model of a prefilming airblast nozzle

 The paper describes the atomisation process of a liquid in an axissymmetric shear layer formed through the interaction of turbulent coaxial jets (respectively, inner and outer jets), with and without swirl, in a model airblast prefilming atomiser. The atomisation process and spray quality was studied using different visualisation techniques, namely laser shadowgraphy and digital image acquisition. The experiments were conducted for different liquid flow rates, Reynolds numbers ranging from 6600 to 66000 and 27300 to 92900 for the inner and outer air flows, respectively, for different outer flow swirl levels, and two liquid film thicknesses −0.2 and 0.7 mm. All the tests were carried out at atmospheric pressure and using water. The results include the analysis of the film structure at break-up and of the break-up length, and suggest that the deterioration of the liquid film close to the atomising edge exhibits a periodic behaviour and is mainly dependent on the inner air velocity. Film thickness strongly affects the time and length scales of the break-up process for the lower range of air velocities. For higher inner air velocities, the break-up length and time become less dependent on liquid flow rate and initial film thickness. Received: 14 March 1997/Accepted: 27 October 1997     

Model of a wavy flow in a falling film of a viscous liquid

A new model is developed for describing long-wave perturbations in a falling film of a viscous liquid. The model is based on an integral approach and an expansion of the velocity profile into a series in linearly independent basis functions of a boundary-value problem. A linear analysis of film flow stability is performed, and dispersion dependences are obtained. Results predicted by the new model are demonstrated to be in good agreement with available experimental data on the film flow over a gently sloping surface.     

Droplet formation in the rupture of a liquid film under the action of a thermal load

The deformation and rupture of a liquid film, suspended between two solid walls, under the action of a localized thermal load is considered. For studying this process, a two-dimensional model is used, which describes the motion of a thin layer of a viscous non-isothermal liquid under microgravity conditions. For modeling the dynamics of the liquid, the Navier-Stokes equations in the “vorticity—stream function” variables are used. A numerical analysis of the influence of thermal loads on the deformation and the mechanism of the rupture of freely suspended films is performed. Is is shown that for a certain width of the thermal beam acting on the film free surface the rupture of the film may occur with the formation of a droplet. The results of the solution of model problems are presented.     

A two-phase boundary-layer model for laminar mixed-convection condensation with a noncondensable gas on inclined plates

A finite-volume-based numerical model for mixed-convection laminar film condensation from a flowing mixture of a vapor and a heavier noncondensable gas on inclined isothermal flat plates is presented. The full boundary layer equations for the liquid film and the vapor-gas mixtures (including liquid inertia and energy convection terms) are solved implicitly with appropriate liquid-mixture interface conditions. Results were obtained for three mixtures, covering wide ranges of liquid Prandtl number and free-stream gas concentration in the forced-convection, mixed-convection and free-convection flow regimes. The effects of liquid inertia were found to be significant only for low-Prandtl-number fluids and lower gas concentrations. The effects of liquid energy convection were found to be significant only for high-Prandtl-number fluids and to be most significant for mixed-convection condensation. Received on 3 March 1998     

基于湍流模型的高速螺旋槽机械密封稳态性能研究

针对高速工况下的液膜润滑螺旋槽端面机械密封,建立了其湍流润滑模型,采用有限单元法结合松弛迭代技术实现了润滑方程和液膜湍流模型的数值求解,对比分析了层流模型和湍流模型下不同螺旋槽几何参数和工况参数对密封性能的影响. 结果表明:液膜湍流效应显著提升了螺旋槽机械密封端面液膜的动力润滑效应,密封的开启力、泄漏率和刚度明显大于层流模型预测值. 在不同条件下,比较而言螺旋槽内产生更加明显的湍流效应,其内液膜流动行为远不同于层流模型. 以开启力为优化目标,湍流模型获得的优化螺旋槽几何参数在螺旋角、槽深明显不同于层流模型. 在高速和低黏度介质下,机械密封的湍流效应不可忽略.      

Development of a PLIC-VOF method for the dynamic simulation of entry region flow in a laminar falling film

The present work describes a numerical procedure to simulate the development of hydrodynamic entry region in a gravity-driven laminar liquid film flow over an inclined plane. It provides a better insight into the physics of developing film in entry region. A novel numerical approach is proposed which has the potential to provide solutions for the complex physics of liquid film spreading on solid walls. The method employs an incompressible flow algorithm to solve the governing equations, a PLIC-VOF method to capture the free surface evolution and a continuum surface force (CSF) model to include the effect of surface tension. To account for the moving contact line on the solid substrate, a precursor film model based wall treatment is implemented. Liquid film flow has been simulated for the Reynolds number range of 5 ≤ Re ≤ 37.5, and the predicted results are found to agree well with the available analytical and experimental data.     

Numerical simulation of ligament-growth on a spinning wheel

Liquid ligaments can grow from perturbations in liquid film spread on a spinning wheel due to the centrifugal force acting on the film. Typically, the growth is strongly influenced by the surface tension on the evolving liquid–air interface. This phenomenon, frequently exploited in industry for the production of fibers, was investigated numerically and the Volume-of-Fluid (VOF) methodology was used to model the interface. The freely available Gerris code with adaptive mesh refinement (AMR) was used to achieve the fine resolution of the computational grid required at the evolving liquid–air interface. The results were compared with the experimental data captured by a high-speed camera. The influence of the process operating variables on the ligament growth is also presented.     

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

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

Convective heat and mass transfer from a falling film to a laminar gas stream

A numerical study has been made of convective heat and mass transfer from a falling film to a laminar gas stream between vertical parallel plates. The effects of gas-liquid phase coupling, variable thermophysical properties, and film vaporization have been considered. Simultaneous mass, momentum and heat transfer between liquid film and gas stream is numerically studied by solving the respective governing equations for the liquid film and gas stream together. The influences of the inlet liquid temperature and liquid flowrate on the cooling of liquid film are examined for air-water and air-ethanol systems. Results show that the heat transfer from the gas-liquid interface to the gas stream is predominantly determined by the latent heat transfer connected with film evaporation. Additionally, better liquid film cooling is noticed for the system having a higher inlet liquid temperature or a lower liquid flowrate.     

Bifurcation analysis of the travelling waves on a falling power-law fluid film

Non-linear waves on the surface of a falling film of power-law fluid on a vertical plane are investigated. The waves are described by evolution equations previously derived as a generalization of the model for the Newtonian liquid. This paper presents a numerical bifurcation analysis of the steady travelling waves on a falling film described by an equation containing three parameters. It is shown that the wave regimes are sensitive to the power-law index as well as to the film parameter and the wavenumber that is typical for the Newtonian liquid.     

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

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

Mass transfer with gas desorption from the surface of a liquid film in the presence of a cocurrent flow

Results of a thcoretical and experimental study of dynamics and mass transfer during desorption of a gas from a liquid film in the presence of a cocurrent air flow are presented. The calculation model is based on solving integral momentum and diffusion relations for the gaseous and liquid phascs. Both laminar and turbulent regimes of the film flow are analyzed. The experimental study of mass transfer was conducted for carbon dioxide desorption from a water film. Criterial relations for mass transfer in the gaseous and liquid phases are obtained. The experiments showed that the heat-transfer coefficients for the case under study are one order of magnitude grcater than those for the flow of a smooth film. Possible mechanisms of such an appreciable intensification of the liquid-film mass transfer in a cocurrent gas flow are discussed. Kutateladze Institute of Thermal Physics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 4, pp. 131–138, July–August, 2000.     

过冷下降液膜破裂的分析

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

Solid dissolution in a thin liquid film on a horizontal rotating disk

A model for the rate of dissolution in liquid film on horizontal rotating disk is obtained by the method of Leveque. It as well as models found in the literature are subjected to experimental verification by dissolving disk cast of gypsum in two liquids. Satisfactory agreement with the model predictions is found. The rate with rotation is compared to that in gravitational film. Enhancements up to 2.5 times are established.