垂直上升气液环状流流动参数的数值计算

提出了一个新的气核-液膜耦合模型来求解垂直上升气液环状流在充分发展段的流动参数.本模型考虑了液膜、气核以及它们之间的相互影响和作用.模型中基本的气核区域和液膜区域的质量和动量方程由Fluent6.3.26进行求解,而液滴方程以及相界面上的夹带和沉积作用通过用户自定义接口函数UDF(User Defined Functi...     

Interpretation of atomization rates of the liquid film in gas-liquid annular flow

For vertical gas-liquid annular flow the fraction of the liquid in the gas is controlled by the rate of atomization of the liquid film flowing along the wall and the rate of deposition of droplets entrained in the gas. Measurements of the rate of atomization are interpreted by a Kelvin-Helmholtz mechanism. Small wavelets on the liquid film are visualized to be entrained when wave-induced variations in the gas pressure cannot be counterbalanced by surface tension effects.     

Experimental and CFD investigation of flow behavior and sand erosion pattern in a horizontal pipe bend under annular flow

The internal erosion of pipelines in oil and gas storage and transportation engineering is highly risky. High gas velocity of annular flow entrained sand will cause damage to the pipelines, and may further result in thinning of the wall. If this damage lasts for a long time, it may cause pipeline leakage and cause huge economic losses and environmental problems. In this research, an experimental device for studying multiphase flow erosion is designed, including an erosion loop and an experimental elbow that can test the erosion rate. The annular flow state and pipe wall erosion morphology can also be tested by the device. The computational fluid dynamics (CFD) method is combined with the experiment to further study the annular flow erosion mechanism in the pipeline. The relationship between gas-liquid-solid distribution and erosion profile was studied. The results show that the most eroded region occurs between 22.5° and 45° in the axial angle direction and between 90° and 135° in the circumferential angle direction of the elbow. The pits and deep scratches form on the surface of the sample after the sand collision.     

Gas—liquid annular flow under microgravity conditions: a temporal linear stability study

A temporal linear stability study was performed for a gas—liquid annular flow configuration under microgravity conditions. Data used to validate the modeling includes that generated by Texas A&M as well as all the other known data in two-phase flow under reduced gravity conditions. Following a discussion of theoretical considerations on the growth rates of different instabilities, it is shown that given the fluid properties, pipe diameter and phasic flow rates, one can predict with a high level of confidence the flow regime in the pipe. Acceptable confidence levels (80%) are achieved when one differentiates between slug, slug—annular, and annular flow. Higher confidence levels (90%) are found when one differentiates between slug and annular flow by merging the annular and slug—annular categories.     

环形通道内湍流旋流流动的数值模拟

本文应用一种考虑湍流－旋流相互作用及湍流脉动各向异性的新的代数Ｒｅｙｎｏｌｄｓ应力模型,对环形通道内的湍流旋流流动进行了数值模拟,研究了改主为旋流流数,进口轴向速度及半径比等参数对环形通道内湍流流动的影响,以及对强化环形通道内传热的作用。     

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.     

Structure of thin liquid films in gas-liquid horizontal flow

The structure of the liquid film in horizontal annular flow is studied visually using the refractive index matching technique. The liquid film is found to contain significant amount of air bubbles, which are continuously entrained, broken up and released by the rolling motion within the film. A new conceptual picture of the gas-liquid interface is presented.     

Numerical simulation for the air entrainment of aerated flow with an improved multiphase SPH model

Aerated flow is a complex hydraulic phenomenon that exists widely in the field of environmental hydraulics. It is generally characterised by large deformation and violent fragmentation of the free surface. Compared to Euler methods (volume of fluid (VOF) method or rigid-lid hypothesis method), the existing single-phase Smooth Particle Hydrodynamics (SPH) method has performed well for solving particle motion. A lack of research on interphase interaction and air concentration, however, has affected the application of SPH model. In our study, an improved multiphase SPH model is presented to simulate aeration flows. A drag force was included in the momentum equation to ensure accuracy of the air particle slip velocity. Furthermore, a calculation method for air concentration is developed to analyse the air entrainment characteristics. Two studies were used to simulate the hydraulic and air entrainment characteristics. And, compared with the experimental results, the simulation results agree with the experimental results well.     

The steady annular extrusion of a Newtonian liquid under gravity and surface tension

The steady extrusion of a Newtonian liquid through an annular die and its development outside and away from the die are studied under the influence of gravitational and surface tension forces. The finite element method (FEM) is used for the simulations. The positions of the inner and outer free surface profiles are calculated simultaneously with the other unknown fields, i.e. using the Newton–Raphson iterative scheme. The effects of three relevant parameters, i.e. the Reynolds, the Stokes and the capillary numbers, on the shape of the annular film are studied for two values of the inner to the outer diameter ratio, corresponding to a thick and a thin annular film respectively. A one‐dimensional model for the extrudate region, valid for thin annular films, is also presented, and its predictions are compared with the two‐dimensional finite element calculations. Despite the fact that it is valid away from the die exit, the one‐dimensional model predicts satisfactorily the effects of the Stokes and capillary numbers. Copyright © 2000 John Wiley & Sons, Ltd.     

Time-dependent behaviour of the liquid film in horizontal annular flow

A crucial point still to be established in the prediction of the film thickness distribution in horizontal annular two-phase flow is the mechanism(s) for transporting liquid from the bottom to the top part of the tube. To resolve this issue, the time-dependent behaviour of the liquid film is studied. Wave characteristics such as velocity and frequency are measured around the circumference. It is inferred from the autospectral density functions of film thickness variation that disturbance waves play an important, but as yet unclear, role in the formation of a liquid film in the top part of the tube. A new mechanism, based on the shape of disturbance waves is proposed.     

内管自转且公转时环管幂律流的受力分析

数值模拟了环管中内管偏心自转且公转时由轴向压力所驱动的幂律流体充分发展层流,分析了内管上的流体作用力。结果表明,内管偏心自转时流体作用力具有推动内管作和自转同向公转的效果。当只有外力矩驱动内管自转时,由于流体的作用,随内管线密度的不同,内管能达到的受力平衡态也不同:线密度较小时内管仅能在同心自转时达到受力平衡;线密度较大时内管能在作具有不变角速度和偏心率公转时达到受力平衡,且内管线密度越大,对应的受力平衡的公转的偏心率也越大。     

A numerical study of an annular liquid jet in a compressible gas medium

An annular liquid jet in a compressible gas medium has been examined using an Eulerian approach with mixed-fluid treatment. The governing equations have been solved by using highly accurate numerical methods. An adapted volume of fluid method combined with a continuum surface force model was used to capture the gas–liquid interface dynamics. The numerical simulations showed the existence of a recirculation zone adjacent to the nozzle exit and unsteady large vortical structures at downstream locations, which lead to significant velocity reversals in the flow field. It was found that the annular jet flow is highly unstable because of the existence of two adjacent shear layers in the annular configuration. The large vortical structures developed naturally in the flow field without external perturbations. Surface tension tends to promote the Kelvin–Helmholtz instability and the development of vortical structures that leads to an increased liquid dispersion. A decrease in the liquid sheet thickness resulted in a reduced liquid dispersion. It was identified that the liquid-to-gas density and viscosity ratios have opposite effects on the flow field with the reduced liquid-to-gas density ratio demoting the instability and the reduced liquid-to-gas viscosity ratio promoting the instability characteristics.     

Numerical study on the two-phase flow distribution in a T-junction

The main objective of this paper is to investigate the ability of a two-dimensional two-fluid computer code to predict the phase separation in a T-junction. A new semi-implicit numerical scheme is developed for solving the two-fluid model equations. Special attention is directed to the modelling of the constitutive for the interfacial friction term. Detailed distribution of void fraction, pressure and velocities are obtained for an air–water mixture in a vertical tee. Good agreement was obtained between the computer code results and the experimental data for the phase separation in the T-junction.     

Numerical investigation of transient nozzle flow

Abstract. The starting process of two-dimensional and axisymmetric nozzle flows has been investigated numerically. Special attention has been paid to the early phase of the starting process and to the appearance of a strong secondary shock wave. For both cases, shock intensities and velocities are obtained and discussed. The flow evolution in the axisymmetric case is proved to be more complex and the transient starting process is slower than in the plane case. Finally, the effects of changing the nozzle angle and the incident shock wave Mach number on the transient flow are addressed. It is shown that a faster start-up can be induced either by decreasing the nozzle angle or increasing the Mach number of the incident shock wave. Received 16 November 2001 / Accepted 24 September 2002 / Published online 4 December 2002 Correspondence to:A.-S. Mouronval (e-mail: mouronv@coria.fr)     

Numerical study of hysteresis in annular swirling jets with a stepped‐conical nozzle

This study investigates the experimentally observed hysteresis in the mean flow field of an annular swirling jet with a stepped‐conical nozzle. The flow is simulated using the Reynolds‐averaged Navier–Stokes (RANS) approach for incompressible flow with a k–ε and a Reynolds stress transport (RSTM) turbulence model. Four different flow structures are observed depending on the swirl number: ‘closed jet flow’, ‘open jet flow low swirl’, ‘open jet flow high swirl’ and ‘coanda jet flow’. These flow patterns change with varying swirl number and hysteresis at low and intermediate swirl numbers is revealed when increasing and subsequently decreasing the swirl. The influence of the inlet velocity profile on the transitional swirl numbers is investigated. When comparing computational fluid dynamics with experiments, the results show that both turbulence models predict the four different flow structures and the associated hysteresis and multiple solutions at low and intermediate swirl numbers. Therefore, a good agreement exists between experiments and numerics. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical study on effects of check valve closure flow conditions on pressure surges in pumping station with air entrainment

This paper proposes a numerical procedure to better compute the characteristics of pressure surges when check valves close under different flow conditions in a pumping station. Studies have shown that the effects of check valve closure on the pressure transients are predominantly dependent on the magnitude and gradient of the flow velocities immediately downstream of the check valve at the instant of valve closure. Through the present study, it was noted that the transient flow velocities near the check valve of a fluid system are also dependent on the characteristics of the air entrained into the fluid system. An improved numerical computational procedure for the fluid system with air entrainment under different transient conditions downstream of the check valve is also proposed in this paper. With a fluid system operating within the critical range of air entrainment values, the present analysis showed that there is a possibility of ‘high pressure surges’ when the check valves were closed at flow rates other than the positive flow conditions. This phenomenon was confirmed through field observations. This study thus concludes that a detailed numerical transient analysis of the fluid system, with various assumed amounts of entrained air, is necessary whenever there is the possibility of air entrainment into the fluid system, and that the flow conditions at the instant of check valve closure need to be modelled. Copyright © 2001 John Wiley & Sons, Ltd.     

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.     

Numerical simulation of the modulated flow pattern for vertical upflows by the phase separation concept

The multiphase heat transfer could be enhanced by creating thin liquid film on the wall. The phase separation concept is called due to the separated flow paths of liquid and gas over the tube cross section to yield thin liquid film. Our proposed heat transfer tube consists of an annular region close to the wall and a core region, interfaced by a suspending mesh cylinder in the tube. The heat transfer tube is a multiscale system with micron scale of mesh pores, miniature scale of annular region and macroscale of tube diameter and length. Great effort has been made to link from micron scale to macroscale. The Volume of Fluid (VOF) method simulates air/water two-phase flow for vertical upflow. The three-dimensional system was successfully converted to a two-dimensional one by using three equivalent criteria for mesh pores. The non-uniform base grid generation and dynamic grid adaption method capture the bubble interface. The numerical results successfully reproduce our experimental results. The numerical findings identify the following mechanisms for the enhanced heat transfer: (a) counter-current flow exists with upward flow in the annular region and downward flow in the core region; (b) void fractions are exact zero in the core region and higher in the annular region; (c) the liquid film thicknesses are decreased to 1/6–1/3 of those in the bare tube section; (d) the gas–liquid mixture travels much faster in the annular region than in the bare tube; (e) three-levels of liquid circulation exists: meter-scale bulk liquid circulation, moderate-scale liquid circulation around a single-elongated-ring-slug-bubble, and microliquid circulation following the ring-slug-bubble tails. These liquid circulations promote the fluid mixing over the whole tube length and within the radial direction. The modulated parameters of void fractions, velocities and liquid film thicknesses in the annular region and three-levels of liquid circulation are greatly beneficial for the multiphase heat transfer enhancement.     

环形方管道中爆轰胞格的三维研究

通过实验和三维数值模拟研究了爆轰波在环形管道内的传播。实验采用烟迹板记录了爆轰波的胞格结构。数值模拟基于带化学反应的三维Euler方程,采用五阶精度的WENO格式捕捉激波,采用具有TVD性质的三阶Runge-Kutta法处理时间项,并结合并行技术,对爆轰波的传播进行了数值研究。结果表明,环形管道外壁为收敛壁面,由于其对流场的压缩效应,外壁面及附近的胞格较小,且较均匀。而内壁为发散壁面,其对流场起稀疏效应,内壁面及附近的胞格较大,且呈周期性变化。同时, 不同壁面的胞格结构均出现了拍波(slapping wave),其形状呈弯曲的折线。