Analysis of unsteady gas–liquid flows in a rectangular tank: Comparison of Euler–Eulerian and Euler–Lagrangian simulations

We study the dynamics of gas–liquid flows experimentally and computationally in a rectangular bubble column where the gas source is introduced at the corner. The flow in this reactor is complex and inherently unsteady in nature. The two-dimensional liquid phase velocity field is calculated by an Eulerian approach solving the unsteady Reynolds Averaged Navier Stokes equations. The conservation equations are closed using a two parameter turbulence model. The two-way coupling was accounted for by adding source terms in the conservation equations of the continuous phase to take into account the interaction with the dispersed phase. Bubble tracking is achieved through a Lagrangian approach. Here the equations of motion are solved taking into account the drag, pressure, buoyancy and gravity forces. The time-averaged flows along with the variables which characterize turbulence are analyzed for a wide range of gas flow-rates using Euler–Lagrangian simulations. These simulation predictions are validated with Euler–Eulerian simulations where the gas-phase distribution is captured as a void fraction and PIV experiments. The motion of bubbles induces turbulence in the flow. The applicability of two parameter models for turbulence like the standard k–ε model on time-averaged flow properties is addressed. From the results of the time averaged velocity field, turbulence intensity, turbulent viscosity and gas hold-up profiles, it is concluded that the Euler–Lagrangian model is applicable at lower gas flow-rates. The Euler–Eulerian approach was found to be valid at lower as well as higher gas flow-rates.     

Solution of the Navier–Stokes equations in velocity–vorticity form using a Eulerian–Lagrangian boundary element method

This paper describes the Eulerian–Lagrangian boundary element model for the solution of incompressible viscous flow problems using velocity–vorticity variables. A Eulerian–Lagrangian boundary element method (ELBEM) is proposed by the combination of the Eulerian–Lagrangian method and the boundary element method (BEM). ELBEM overcomes the limitation of the traditional BEM, which is incapable of dealing with the arbitrary velocity field in advection‐dominated flow problems. The present ELBEM model involves the solution of the vorticity transport equation for vorticity whose solenoidal vorticity components are obtained iteratively by solving velocity Poisson equations involving the velocity and vorticity components. The velocity Poisson equations are solved using a boundary integral scheme and the vorticity transport equation is solved using the ELBEM. Here the results of two‐dimensional Navier–Stokes problems with low–medium Reynolds numbers in a typical cavity flow are presented and compared with a series solution and other numerical models. The ELBEM model has been found to be feasible and satisfactory. Copyright © 2000 John Wiley & Sons, Ltd.     

Numerical simulation of circulation in gas-liquid column reactors: isothermal, bubbly, laminar flow

The gas-liquid flow inside a circular, isothermal column reactor with a vertical axis has been studied using numerical simulations. The flow is assumed to be in the laminar, bubbly flow regime which is characterized by a suspension of discrete air bubbles in a continuous liquid phase such as glycerol water. The mathematical formulation is based on the conservation of mass and momentum principle for the liquid phase. The gas velocity distribution is calculated via an empirically prescribed relative velocity as a function of void fraction. The interface viscous drag forces are prescribed empirically. For some cases a profile shape is assumed for the void ratio distribution. The influence of various profile shapes is investigated. The results are compared with those where the void ratio distribution is calculated from the conservation of mass equation. The mathematical model has been implemented by modifying a readily available computer code for single-phase newtonian fluid flows. The numerical discretization is based on a finite volume approach. The predictions show a good agreement with measurements. The circulation pattern seems not to be so sensitive to the actual shape of the void fraction profiles, but the inlet distribution of it is important. A significantly different flow pattern results when the void fraction distribution is calculated from the transport equation, as compared to those with a priori prescribed profiles. When the void fraction is uniformly distributed over the whole distributor plate, no circulation is observed. Calculations also show that even the two-phase systems with a few discrete bubbles can be simulated successfully by a continuum model.     

The numerical solution of the transient two‐phase flow in rigid pipelines

Consideration is given in this paper to the numerical solution of the transient two‐phase flow in rigid pipelines. The governing equations for such flows are two coupled, non‐linear, hyperbolic, partial differential equations with pressure dependent coefficients. The fluid pressure and velocity are considered as two principle dependent variables. The fluid is a homogeneous gas–liquid mixture for which the density is defined by an expression averaging the two‐component densities where a polytropic process of the gaseous phase is admitted. Instead of the void fraction, which varies with the pressure, the gas–fluid mass ratio (or the quality) is assumed to be constant, and is used in the mathematical formulation. The problem has been solved by the method of non‐linear characteristics and the finite difference conservative scheme. To verify their validity, the computed results of the two numerical techniques are compared for different values of the quality, in the case where the liquid compressibility and the pipe wall elasticity are neglected. Copyright © 1999 John Wiley & Sons, Ltd.     

On the numerical study of bubbly flow created by ventilated cavity in vertical pipe

The dispersion of bubbles into a down-liquid flow in a vertical pipe is investigated. At low flow rates, the intended design of a swarm of discrete bubbles is achieved. At high flow rates, a ventilated cavity is nonetheless formed, which is attached close to the gas sparger. Behind this ventilated cavity, three different flow regimes characterize the complex bubbly flow field downstream of the down-liquid flow: vortex region with high void fraction, transitional region and pipe flow region. In this study, a numerical model that solved the entire development of the gas–liquid flow including the extended single-phase liquid region upstream to the wall-jet and recirculating-vortex zones in order to allow a more realistic determination of the boundary conditions of the down-liquid flow was adopted. Coupling with the Eulerian–Eulerian two-fluid model to solve the respective gas and liquid phases, a population balance model was also applied to predict the bubble size distribution in the wake right below the cavity base as well as further downstream in the transitional and fully-developed pipe flow regions. The numerical model was evaluated by comparing the numerical results against the data derived from theoretical, numerical and experimental approaches. Prediction of the Sauter mean bubble diameter distributions by the population balance approach at different axial locations confirmed the dominance of breakage due to the high turbulent intensity below the ventilated cavity which led to the generation of small gas bubbles at high void fraction. Further downstream, the coalescence effect dominated leading to merging of the small bubbles to form bigger bubbles.     

NUMERICAL PREDICTION OF TWO-PHASE FLOW IN BUBBLE COLUMNS

A numerical model is described for the prediction of turbulent continuum equations for two-phase gas–liquid flows in bubble columns. The mathematical formulation is based on the solution of each phase. The two-phase model incorporates interfacial models of momentum transfer to account for the effects of virtual mass, lift, drag and pressure discontinuities at the gas–liquid interface. Turbulence is represented by means of a two-equation k–ϵ model modified to account for bubble-induced turbulence production. The numerical discretization is based on a staggered finite-volume approach, and the coupled equations are solved in a segregated manner using the IPSA method. The model is implemented generally in the multipurpose PHOENICS computer code, although the present appllications are restricted to two-dimensional flows. The model is applied to simulate two bubble column geometries and the predictions are compared with the measured circulation patterns and void fraction distributions.     

Prediction of mono-disperse gas–liquid turbulent flow in a vertical pipe

A two-fluid model in the Eulerian–Eulerian framework has been implemented for the prediction of gas volume fraction, mean phasic velocities, and the liquid phase turbulence properties for gas–liquid upward flow in a vertical pipe. The governing two-fluid transport equations are discretized using the finite volume method and a low Reynolds number $k-\varepsilon$ model is used to predict the turbulence field for the continuous liquid phase. In the present analysis, a fully developed one-dimensional flow is considered where the gas volume fraction profile is predicted using the radial force balance for the bubble phase. The current study investigates: (1) the turbulence modulation terms which represent the effect of bubbles on the liquid phase turbulence in the k–ε transport equations; (2) the role of the bubble induced turbulent viscosity compared to turbulence generated by shear; and (3) the effect of bubble size on the radial forces which results in either a center-peak or a wall-peak in the gas volume fraction profiles. The results obtained from the current simulation are generally in good agreement with the experimental data, and somewhat improved over the predictions of some previous numerical studies.     

A direct numerical simulation of axisymmetric liquid–gas two‐phase laminar flows in a pipe

An accurate finite‐volume‐based numerical method for the simulation of an isothermal two‐phase flow, consisting of a liquid slug translating in a non‐reacting gas in a circular pipe is presented. This method is built on a sharp interface concept and developed on an Eulerian Cartesian fixed‐grid with a cut‐cell scheme and marker points to track the moving interface. The unsteady, axisymmetric Navier–Stokes equations in both liquid and gas phases are solved separately. The mass continuity and momentum flux conditions are explicitly matched at the true surface phase boundary to determine the interface shape and movement. A quadratic curve fitting algorithm with marker points is used to yield smooth and accurate information of the interface curvatures. It is uniquely demonstrated for the first time with the current method that conservation of mass is strictly enforced for continuous infusion of flow into the domain of computation. The method has been used to compute the velocity and pressure fields and the deformation of the liquid core. It is also shown that the current method is capable of producing accurate results for a wide range of Reynolds number, Re, Weber number, We, and large property jumps at the interface. Copyright © 2010 John Wiley & Sons, Ltd.     

Modeling of Axisymmetric Two-phase Dilute Flows

This paper deals with the numerical treatment of Eulerian approach for dilute two-phase compressible flows (gas-particles mixtures) in axisymmetric configurations. For dilute flows, two classes of models depending on the dispersed phase volumetric fraction can be found. The volume occupied by the particles may be considered, that yields a model in which the gas phase and the dispersed phase equations are coupled through the void fraction and the source terms (Delhaye model). The void fraction effects can be neglected, that means the gas phase is a carrier phase for the particles (Ishii model). The mathematical nature of the two models is demonstrated from analysis of characteristic directions. For the Delhaye's model, a centered scheme is used to solve the system of partial differential equations, while an upwind TVD scheme is used for the Ishii's one. Then, it is shown that the problem of symmetry boundary conditions does not depend on the physical approach, as long as the flow remains dilute. However, a classical treatment for symmetry boundary conditions at the geometrical axis leads to large errors. A particular treatment for this boundary is presented: a new class of particles, described by a supplementary system of equations, is required.     

Wall shear stress and void fraction in Poiseuille bubbly flows: Part I: simple analytic predictions

Upward, co-current bubbly flows in a vertical rectangular duct are investigated at low liquid Reynolds numbers. The conditions considered are such that the pseudo-turbulent stresses remain negligible compared to the viscous stresses. The void fraction transverse distribution is idealised as step-functions and is then inserted in the conservation equations supplemented by appropriate closure laws. Analytical expressions are then obtained for the axial velocity profiles, for the lineic gas fraction and for the wall friction. The sensitivity of these quantities to the void distribution, characterised by the void fraction and the width of the three layers introduced, is discussed. It is shown that differential buoyancy effects govern the modification of the liquid velocity profiles. Notably, void peaking near walls is able to induce a wall shear stress many times higher than its single-phase flow counterpart at the same liquid flow rate. Also, the presence of a near wall region free of gas favours the onset of downward directed secondary flows. All these features correspond to experimental observations, and a few quantitative comparisons are also presented which support the validity of the model even in case of void coring. A companion paper (part II) will be devoted to systematic comparisons between predictions and experiments in the case of axisymmetric Poiseuille bubbly flows.     

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

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

NUMERICAL MODELLING OF THREE-DIMENSIONAL FLY-ASH FLOW IN POWER UTILITY BOILERS

A two-fluid Eulerian model in combination with a particle–wall collision model and generalized Eulerian boundary conditions for the particulate phase is employed to predict complex three- dimensional fly-ash flows which often cause severe erosion to boiler tubes located in power utility boilers. Mean momentum and mass conservation equations are solved for each phase using a finite volume scheme with two-way coupling and a modified renormalization group (RNG)-based k –ϵ turbulence model. Comparison of predicted particle concentration with measured data is made and excellent agreement is obtained. The detailed character of the particulate velocity field and concentration just downstream of the 180° bend shows a marked dependence on the Stokes number not previously reported. © 1997 by John Wiley & Sons, Ltd.     

A New Approach in Modeling Phase Distribution in Fully Developed Bubbly Pipe Flow

Turbulent two-phase flow equations are derived and solved for fully developed pipe flow using a composite eddy-viscosity model and a new void-fraction equation. The void fraction profile is first specified from experiments and the velocity field is calculated to validate the eddy-viscosity model. Consequently, a new equation is presented for calculation of the void fraction. This void-fraction equation incorporates the gradient of turbulent normal stresses in the radial direction, the conventional lift force, and a contribution from the unsteady drag force. The implications of this new equation, for the bubbly flow regime, are investigated by calculating the void-fraction distribution for a given velocity field. Inclusion of the normal turbulent stresses in the radial direction is shown to simulate correctly the experimentally observed trends of the phase distribution, both for upward and downward bubbly flow, without the need for a fictitious term such as the so called ``lubrication force'. This revised version was published online in July 2006 with corrections to the Cover Date.     

Prediction of turbulent gas–solids flow in curved ducts using the Eulerian–Lagrangian method

The flow of particulate two‐phase flow mixtures occur in several components of solid fuel combustion systems, such as the pressurised fluidised bed combustors (PFBC) and suspension‐fired coal boilers. A detailed understanding of the mixture characteristics in the conveying component can aid in refining and optimising its design. In this study, the flow of an isothermal, dilute two‐phase particulate mixture has been examined in a high curvature duct, which can be representative of that transporting the gas–solid mixture from the hot clean‐up section to the gas turbine combustor in a PFBC plant. The numerical study has been approached by utilising the Eulerian–Lagrangian methodology for describing the characteristics of the fluid and particulate phases. By assuming that the mixture is dilute and the particles are spherical, the governing particle momentum equations have been solved with appropriately prescribed boundary conditions. Turbulence effects on the particle dispersion were represented by a statistical model that accounts for both the turbulent eddy lifetime and the particle transit time scales. For the turbulent flow condition examined it was observed that mixtures with small particle diameters had low interphase slip velocities and low impaction probability with the pipe walls. Increasing the particle diameters (>50 μm) resulted in higher interphase slip velocities and, as expected, their impaction probability with the pipe walls was significantly increased. The particle dispersion is significant for the smaller sizes, whereas the larger particles are relatively insensitive to the gas turbulence. The main particle impaction region, and locations most prone to erosion damage, is estimated to be within an outer duct length of two to six times the duct diameter, when the duct radius of curvature to the duct diameter ratio is equal to unity. Copyright © 1999 John Wiley & Sons, Ltd.     

Bubble and liquid turbulence characteristics of bubbly flow in a large diameter vertical pipe

The bubble and liquid turbulence characteristics of air–water bubbly flow in a 200 mm diameter vertical pipe was experimentally investigated. The bubble characteristics were measured using a dual optical probe, while the liquid-phase turbulence was measured using hot-film anemometry. Measurements were performed at six liquid superficial velocities in the range of 0.2–0.68 m/s and gas superficial velocity from 0.005 to 0.18 m/s, corresponding to an area average void fraction from 1.2% to 15.4%. At low void fraction flow, the radial void fraction distribution showed a wall peak which changed to a core peak profile as the void fraction was increased. The liquid average velocity and the turbulence intensities were less uniform in the core region of the pipe as the void fraction profile changed from a wall to a core peak. In general, there is an increase in the turbulence intensities when the bubbles are introduced into the flow. However, a turbulence suppression was observed close to the wall at high liquid superficial velocities for low void fractions up to about 1.6%. The net radial interfacial force on the bubbles was estimated from the momentum equations using the measured profiles. The radial migration of the bubbles in the core region of the pipe, which determines the shape of the void profile, was related to the balance between the turbulent dispersion and the lift forces. The ratio between these forces was characterized by a dimensionless group that includes the area averaged Eötvös number, slip ratio, and the ratio between the apparent added kinetic energy to the actual kinetic energy of the liquid. A non-dimensional map based on this dimensionless group and the force ratio is proposed to distinguish the conditions under which a wall or core peak void profile occurs in bubbly flows.     

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.     

Modeling of three-phase heavy oil–water–gas bubbly flow in upward vertical pipes

A bubbly gas–bubbly oil flow pattern may occur when water, heavy oil and gas flow simultaneously in vertical pipes in such a way that water is the continuous phase. In this work, a one-dimensional, thermal, transient two-fluid mathematical model, for such flow, is presented. The model consists of mass, momentum and energy conservation equations for every phase whose numerical solution is based on the finite difference technique in the implicit scheme. The model is able to predict pressure, temperature, volumetric fraction and velocity profiles. For accurate modeling of multiphase flows, the key issue is to specify the adequate closure relationships, thus drag and virtual mass forces for the gas and oil phases were taken into account and special attention was paid on the gas–oil drag force. When this force was included into the model it was found that: (1) such force had the same order of magnitude than the oil drag force and both forces were smaller than the gas drag force, (2) the pressure, gas and oil velocities and gas and oil volume fraction profiles were affected, (3) the numerical stability was increased. The model predictions are in agreement with experimental data reported in literature.     

Studies on two-phase co-current air/non-Newtonian shear-thinning fluid flows in inclined smooth pipes

In this work, co-current flow characteristics of air/non-Newtonian liquid systems in inclined smooth pipes are studied experimentally and theoretically using transparent tubes of 20, 40 and 60 mm in diameter. Each tube includes two 10 m long pipe branches connected by a U-bend that is capable of being inclined to any angle, from a completely horizontal to a fully vertical position. The flow rate of each phase is varied over a wide range. The studied flow phenomena are bubbly flow, stratified flow, plug flow, slug flow, churn flow and annular flow. These are observed and recorded by a high-speed camera over a wide range of operating conditions. The effects of the liquid phase properties, the inclination angle and the pipe diameter on two-phase flow characteristics are systematically studied. The Heywood–Charles model for horizontal flow was modified to accommodate stratified flow in inclined pipes, taking into account the average void fraction and pressure drop of the mixture flow of a gas/non-Newtonian liquid. The pressure drop gradient model of Taitel and Barnea for a gas/Newtonian liquid slug flow was extended to include liquids possessing shear-thinning flow behaviour in inclined pipes. The comparison of the predicted values with the experimental data shows that the models presented here provide a reasonable estimate of the average void fraction and the corresponding pressure drop for the mixture flow of a gas/non-Newtonian liquid.