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

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

Nonlinearwaves in a liquid film-gas flow system

The instability and regular nonlinear waves in the film of a heavy viscous liquid flowing along the wall of a round tube and interacting with a gas flow are investigated. The solutions for the wave film flows are numerically obtained in the regimes from free flow-down in a counter-current gas stream to cocurrent upward flow of the film and the gas at fairly large gas velocities. Continuous transition from the counter-current to the cocurrent flow via the state with a maximum amplitude of nonlinear waves and zero values of the liquid flow rate and the phase velocity is investigated. The Kapitsa-Shkadov method is used to reduce a boundary value problem to a system of evolutionary equations for the local values of the layer thickness and the liquid flow rate.     

Experimental study on hydrodynamic characteristics of upward gas–liquid slug flow

To clarify the impacts of the hydrodynamic boundary layer and the diffusion boundary layer in the near wall zone on gas–liquid two-phase flow induced corrosion in pipelines, the hydrodynamic characteristics of fully developed gas–liquid slug flow in an upward tube are investigated with limiting diffusion current probes, conductivity probes and digital high-speed video system. The Taylor bubble and the falling liquid film characteristics are studied, the effects of various factors are examined, and the experimental results are compared with the data and models available in literature. The length of Taylor bubble, the local void fraction of the slug unit and the liquid slug, the shear stress and mass transfer coefficient in the near wall zone, are all increased with the increase of superficial gas velocity and decreased with the increase of superficial liquid velocity, whereas the length of liquid slug and the liquid slug frequency are changed contrarily. The alternate wall shear stress due to upward gas–liquid slug flow is considered to be one of the major causes for the corrosion production film fatigue cracking. A normalized formula for mass transfer coefficient is obtained based on the experimental data.     

A new approach to analyze entrance region mass transfer within a falling film

A boundary layer model is developed to analyze diffusion through a laminar falling film for incomplete penetration of the dissolved gas. Regarding the rather intractable nature of the problem, Kantorovich Integral Method is chosen. Accordingly, a mass transfer boundary layer is assumed to keep growing perpendicularly to the falling film until it hits the wall. Such an approach is superior to its preceding artworks based on Higbie’s penetration theory in terms of implementation of more realistic conditions/modeling assumptions. Furthermore, unlike penetration model, this approach gives a criterion whether the diffusion is complete. Comparing the two models, boundary layer model estimates up to 2.3 % larger mass transfer coefficients. Moreover, a sensitivity analysis of liquid velocity distribution upon Sherwood number is conducted. It is found that the local velocity at gas–liquid interface is of highest dominance. Experimental data of SO₂ absorption in water reported in literature is exploited to validate the model. It is shown that boundary layer model better fits the experimental data.     

Wavy wall influence on the hydrodynamic instability of a liquid film flowing along an inclined plane

The film dynamic of a thin liquid along an inclined and wavy wall was numerically depicted in a weighted-residual integral boundary layer equation. A qualitative and quantitative analysis was initially carried out and accurate comparisons were obtained from experimental data on film instability along a flat and inclined as well as a wavy wall. To pinpoint the effect of waviness on film instability, 20 wavy wall periods in the computational CFD domain were considered. Several waviness parameters were studied and shown to have taken on a major role in the film instability process. Finally, a wide range of main wall inclination angles was taken into account, and consequent numerical data permitted identification of a threshold angle value. For wall angles higher than the threshold angle, the film behaved as though no corrugations were present. For lower angles, the film was repeatedly altered during the acceleration and deceleration phases.     

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.     

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

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

Numerical study of condensing a small concentration of vapour inside a vertical tube

The purpose of this study is to analyse the combined heat and mass transfer of liquid film condensation from a small steam–air mixtures flowing downward along a vertical tube. Both liquid and gas stream are approached by two coupled laminar boundary layer. An implicit finite difference method is employed to solve the coupled governing equations for liquid film and gas flow together with the interfacial matching conditions. The effects of a wide range of changes of three independent variables (inlet pressure, inlet Reynolds number and wall temperature) on the concentration at exit tube, local Nusselt and Sherwood numbers, film thickness, accumulated condensate rate and temperature are carefully examined. The numerical results indicate that in the case of condensing a small concentration of vapours from a mixture, the resistance to heat and mass transfer by non-condensable gas becomes very intense. The comparisons of average Nusselt number and local condensate heat transfer coefficient with the literature results are in good agreement.     

Quadratic programming algorithm for wall slip and free boundary pressure condition

Wall slip is often observed in a highly sheared fluid film in a solid gap. This makes a difficulty in mathematical analysis for the hydrodynamic effect because fluid velocity at the liquid–solid interfaces is not known a priori. If the gap has a convergent–divergent wedge, a free boundary pressure condition, i.e. Reynolds pressure boundary condition, is usually used in the outlet zone in numerical solution. This paper, based on finite element method and parametric quadratic programming technique, gives a numerical solution technique for a coupled boundary non‐linearity of wall slip and free boundary pressure condition. It is found that the numerical error decreases with the number of elements in a negative power law having an index larger than 2. Our method does not need an iterative process and can simultaneously gives rise to fluid film pressure distribution, wall slip velocity and surface shear stress. Wall slip always decreases the hydrodynamic pressure. Large wall slip even causes a null hydrodynamic pressure in a pure sliding solid gap. Copyright © 2005 John Wiley & Sons, Ltd.     

下游泵送螺旋槽密封空化试验及性能分析

基于自主设计的可视化试验装置及膜温和膜厚测量方法,对下游泵送螺旋槽密封空化特征及性能参数进行试验研究.探讨了油压和转速对不同螺旋槽密封液膜中空穴发生位置、空穴分布及空穴边界的影响,拟合了不同螺旋槽密封空穴边界的试验关系式,并对不同空化模型的理论泄漏量和膜厚与相应试验值进行了对比分析.结果表明:内槽型和中槽型密封的空穴均发生在螺旋槽内,但两者的空穴区形状明显不同;油压的增加有助于抑制液膜中空穴的发生,而转速的增加反之;尤其是内槽型,油压对其空穴发生影响更为显著;在低转速或高油压时,Reynolds和JFO两空化模型均可用于预测泄漏量和膜厚;在高转速或低油压时,JFO空化模型预测值更准确,而Reynolds空化模型预测值偏大.     

Dynamic film thickness between bubbles and wall in a narrow channel

The present paper describes a novel technique to characterize the behavior of the liquid film between gas bubbles and the wall in a narrow channel. The method is based on the electrical conductance. Two liquid film sensors are installed on both opposite walls in a narrow rectangular channel. The liquid film thickness underneath the gas bubbles is recorded by the first sensor, while the void fraction information is obtained by measuring the conductance between the pair of opposite sensors. Both measurements are taken on a large two-dimensional domain and with a high speed. This makes it possible to obtain the two-dimensional distribution of the dynamic liquid film between the bubbles and the wall. In this study, this method was applied to an air–water flow ranging from bubbly to churn regimes in the narrow channel with a gap width of 1.5 mm.     

Pressure wave breakdown of a gas film in a liquid-filled slit

In [1–4] the results of investigating the breakdown of gas bubbles by medium-intensity pressure waves are presented and various bubble breakdown mechanisms are proposed. It is shown that breakdown may occur as a result of the formation of a cumulative jet on the boundary of the bubble or as a result of instability due to the relative motion of the bubble in the wave. In [5] experimental data on the pressure wave breakdown of a gas film in a liquid on a solid wall are reported. It is shown that at wave amplitudes p/p₀1 a liquid jet is formed at the edge of the gas film. The jet, traveling along the wall, strips off the film and carries it into the surrounding liquid. Below we investigate the pressure wave behavior of a gas film in a liquid-filled slit.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.5, pp. 175–178, September–October, 1992.     

过冷下降液膜破裂的分析

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

Nonlinear three-dimensional stretched flow of an Oldroyd-B fluid with convective condition,thermal radiation,and mixed convection

The effect of non-linear convection in a laminar three-dimensional Oldroyd-B fluid flow is addressed. The heat transfer phenomenon is explored by considering the non-linear thermal radiation and heat generation/absorption. The boundary layer assumptions are taken into account to govern the mathematical model of the flow analysis. Some suitable similarity variables are introduced to transform the partial differential equations into ordinary differential systems. The Runge-Kutta-Fehlberg fourth-and fifth-order techniques with the shooting method are used to obtain the solutions of the dimensionless velocities and temperature. The effects of various physical parameters on the fluid velocities and temperature are plotted and examined. A comparison with the exact and homotopy perturbation solutions is made for the viscous fluid case, and an excellent match is noted. The numerical values of the wall shear stresses and the heat transfer rate at the wall are tabulated and investigated. The enhancement in the values of the Deborah number shows a reverse behavior on the liquid velocities. The results show that the temperature and the thermal boundary layer are reduced when the nonlinear convection parameter increases. The values of the Nusselt number are higher in the non-linear radiation situation than those in the linear radiation situation.     

Forced convective film condensation inside vertical tubes

A theoretical study of forced convective film condensation inside vertical tubes is presented. We propose a unified procedure for predicting the pressure gradient and condensation heat transfer coefficient of a vapor flowing turbulently in the core and associated with laminar or turbulent film on the tube wall. The analysis for the vapor flows is performed under the condition that the velocity profiles are locally self-similar. The laminar and turbulent film models equate the gravity, pressure and viscous forces, and consider the effect of interfacial shear. The transition from laminar to turbulent film depends not only on the liquid Reynolds number but also on the interfacial shear stress. In this work we also proposed a new eddy viscosity model which is divided into three regions: the inner region in liquid condensate near the wall; the interface region including both liquid and vapor; and the outer region for the vapor core. Comparisons of the theory with some published experimental data showed good agreement.     

气泡碰撞亲疏水曲壁的行为特性研究

气泡碰撞固壁行为和影响因素的研究一直以来都是科学界关注的重点之一, 其在矿物浮选、气膜减阻等工业领域中的应用也极具科研价值. 论文聚焦曲壁对于气泡撞击行为特性的影响研究. 采用高速摄像技术记录气泡碰撞不同曲率半径下亲疏水曲壁的撞击过程, 分析了曲壁润湿性、曲率半径对气泡碰撞固体曲壁的影响规律. 结果表明, 气泡碰撞亲水曲壁时会发生多次弹跳直至离开曲壁; 曲率半径越大, 弹跳次数越少, 且第一次反弹的最远距离越近, 再次发生碰壁时的速度越小. 而碰撞疏水曲壁时会出现碰撞?滑移?附着的现象, 此外针对液膜挤压破裂的现象, 建立理论模型推导出液膜诱导时间的预测公式, 其主要与液膜厚度、液膜临界破裂厚度和液膜被压缩速度有关, 预测误差小于5.0%.      

Instability of a liquid layer under periodic influence: Falling film in an alternating electric field

The linear analysis of stability of a plane-parallel time-periodic flow is carried out. The numerical method which makes it possible to reduce the spectral problem for the time-dependent Orr–Sommerfeld equation to an algebraic eigenvalue problem is used. The film of viscous conducting liquid which flows down a vertical wall in the normal electric field is considered and parametric resonances are revealed.     

Oscillating heat pipe simulation considering dryout phenomena

In heat transport devices such as oscillating heat pipe (OHP), dryout phenomena is very important and avoided in order to give the optimum performance. However, from the previous studies (including our studies), the dryout phenomena in OHP and its mechanism are still unclear. In our studies of OHP (Senjaya and Inoue in Appl Thermal Eng 60:251–255, 2013; Int J Heat Mass Transfer 60:816–824, 2013; Int J Heat Mass Transfer 60:825–835, 2013), we introduced the importance and roles of liquid film in the operating principle of OHP. In our previous simulation, the thickness of liquid film was assumed to be uniform along a vapor plug. Then, dryout never occurred because there was the liquid transfer from the liquid film in the cooling section to that in the heating section. In this research, the liquid film is not treated uniformly but it is meshed similarly with the vapor plugs and liquid slugs. All governing equations are also solved in each control volume of liquid film. The simulation results show that dryout occurs in the simulation without bubble generation and growth. Dryout is started in the middle of vapor plug, because the liquid supply from the left and right liquid slugs cannot reach until the liquid film in the middle of vapor plug, and propagates to the left and right sides of a vapor plug. By inserting the bubble generation and growth phenomena, dryout does not occur because the wall of heating section is always wetted during the bubble growth and the thickness of liquid film is almost constant. The effects of meshing size of liquid film and wall temperature of heating section are also investigated. The results show that the smaller meshing size, the smaller liquid transfer rate and the faster of dryout propagation. In the OHP with higher wall temperature of heating section, dryout and its propagation also occur faster.     

Effect of wall friction and vortex generation on the radial distribution of different phases

Arguments are presented to prove the existence of rolling vortices in single-phase and two-phase flow. In the liquid phase, they appear in a boundary layer near a wall while in the continuous vapor phase they occur near the interface with a liquid film. The intensity and size of these vortices depend on the local velocity gradients normal to the wall. The interaction between the rotational field associated with such vortices and bubbles in liquid flow or droplets in vapor flow is discussed. This interaction may be called the wall-vortex effect. It appears that several, apparently unrelated, phenomena observed in two-phase flow systems may be interpreted in terms of this mechanism. Among these are: (i) radial void peaking near the walls (ii) vapor velocities less than liquid velocity observed also in vertical upward flow (iii) reduced droplet diffusion near the liquid film and (iv) reduced vapor mixing between subchannels at low steam qualities. The cause of secondary flows in non-circular channels may also be explained in terms of rolling vortices near the walls. Finally, a comparison is made with the well known Magnus effect.     

Minimum wetting rate for a decelerating liquid film

The paper presents an analytical study of a laminar decelerating liquid film falling along a vertical plate. Approximate solutions are obtained for the boundary layer within the film, film thickness, entrance length, and minimum wetting rate. It is shown that the analysis may be extended also to a film flowing over a horizontal cylinder. The theory is in reasonable agreement with the parametric trends observed in experiments on horizontal cylinders.