Simulation of bubbly flows with special numerical treatments of the semi-conservative and fully conservative two-fluid model

2017 ◽  
Vol 174 ◽  
pp. 25-39 ◽  
Author(s):  
Dongyue Li ◽  
Hasse Christian
Keyword(s):  
Fluid Model ◽  
Bubbly Flows ◽  
Two Fluid Model ◽  
Numerical Treatments ◽  
Two Fluid

scholarly journals On the hyperbolicity of the two-fluid model for gas–liquid bubbly flows

2018 ◽  
Vol 57 ◽  
pp. 432-447 ◽  
Author(s):  
N. Panicker ◽  
A. Passalacqua ◽  
R.O. Fox
Keyword(s):  
Fluid Model ◽  
Bubbly Flows ◽  
Two Fluid Model ◽  
Two Fluid

A One Dimensional Two-Fluid Model for Bubbly Flows in Fixed Beds

2002 ◽  
Author(s):  
M.-L. Bordas ◽  
A. Cartellier ◽  
P. Se´chet
Keyword(s):  
Pressure Drop ◽  
Void Fraction ◽  
Bubble Size ◽  
Fluid Model ◽  
Bubbly Flows ◽  
Two Phase ◽  
The Mean ◽  
Two Fluid Model ◽  
Fixed Beds ◽  
Two Fluid

Pressure drop and gas void fraction are important parameters for the design of multiphase packed bed reactors which are widely used in petrochemical industry. Several experimental studies have been devoted to the hydrodynamics of two-phase cocurrent upflow or downflow through fixed beds, and various correlations of limited range of validity are available in the literature. However, there is not yet a clear agreement on the form of the momentum equations to be used in such systems. Early attempts devoted to the pressure drop estimate were based on an extension of the Lockhart-Martinelli approach (Sweeney 1967), Rao et al. 1983). More recently, Attou at al. (1999) proposed the first serious attempt to adapt the Eulerian two-fluid model to cocurrent bubbly flows through packed beds. From an analysis of their proposal, it happens that the basic mechanical equilibrium for the gas phase needs to be reconsidered. In this scope, we derived a new model on the basis of the so-called hybrid approach initially developed for bubbly flows in ducts in absence of shear-induced turbulence (Achard and Cartellier 2000). As a first application, we considered a mean unidirectional flow of a bubbly mixture through a porous medium composed of beads uniform in size. For steady and fully established flows, and assuming a flat void fraction (α) profile, the resulting momentum equations for each phase write: Liquidphase:−dpdz=ρLg+fLS−fLG1−α(1)Gasphase:−dpdz=ρGg+fLS+fLGα(2) where fLS is the resultant of the liquid shear stress exerted on beads surface and on exterior walls, and where the quantity fLG = α F* / Vp represents the interaction force density between the gas and the liquid (F* is the mean force on bubbles and Vp = 4πa3/3 denotes the bubble volume, a being the bubble radius). The main difference with the model derived by Attou et al. is the presence of the fLS term in the gas phase equation. Without this term, the relative velocity of bubbles would be controlled by the axial pressure gradient dP/dz even in non accelerating flows which is unphysical. On the opposite, in the present model (1–2) the relative movement of bubbles is simply due to buoyancy. The set of equations (1–2) provides a mean to exploit the experimental data to derive the required closures, namely the evolution of the friction fLS with the gas content and that of the momentum exchange between phases fLG. Notably, from (1) and (2), one gets fLG=α(1−α)(ρL−ρG)g(3) In order to establish reliable closures, available experimental data of the literature are currently revisited under this framework. For the friction term, which is the principal contribution to the pressure drop, the usual closure law for fLS as given by an Ergun equation adapted to two-phase flows is under analysis. For the interfacial momentum transfer, the objective is to evaluate an “apparent” drag coefficient defined as Cd = F*/[ρL Ur2 π a2 / 2] where the mean relative velocity Ur is defined as the difference between the mean gas and liquid velocities averaged over a volume. Indeed, paralleling an approach already exploited for bubbly flows in ducts (Riviere and Cartellier 1999), it happens that the mean void fraction can be derived from equations (1) and (2) assuming a flat void fraction profile: β(1−β)−α(1−α)=(4π/3)α(1−α)[gδ2VSLνc](aδ)2fd(4) where δ is the typical size of the pores and where fd = (π/2) Rep Cd is expected to be a function of the bubble size, the porosity ε and the void fraction. To extract fd or Cd from (4), a characteristic bubble size must be specified. As shown Fig.1, the bubble size is controlled by the bed geometry and evolves between 0.2 δ and 3 δ in the dilute limit (Bordas et al. (2001)). Analysis of the existing data will be presented based on these size estimates, and comparison will be performed of this “apparent” drag with values measured for isolated bubbles in fixed beds (Fig.2).

Effect of Collision Force on Well-Posedness and Stability of the Two-Fluid Model for Vertical Bubbly Flows

2016 ◽  
Vol 184 (3) ◽  
pp. 353-362 ◽  
Author(s):  
Avinash Vaidheeswaran ◽  
William D. Fullmer ◽  
Martin Lopez de Bertodano
Keyword(s):  
Fluid Model ◽  
Bubbly Flows ◽  
Well Posedness ◽  
Two Fluid Model ◽  
Collision Force ◽  
Two Fluid

A center-averaged two-fluid model for wall-bounded bubbly flows

2006 ◽  
Vol 35 (4) ◽  
pp. 429-461 ◽  
Author(s):  
F.J. Moraga ◽  
A.E. Larreteguy ◽  
D.A. Drew ◽  
R.T. Lahey
Keyword(s):  
Fluid Model ◽  
Bubbly Flows ◽  
Two Fluid Model ◽  
Two Fluid

scholarly journals CFD Two Fluid Model for Adiabatic and Boiling Bubbly Flows in Ducts

10.5772/7097 ◽  
2010 ◽  
Author(s):  
Martin Lopez de Bertodano ◽  
Deoras Prabhudharwadkar
Keyword(s):  
Fluid Model ◽  
Bubbly Flows ◽  
Two Fluid Model ◽  
Two Fluid

New Two-Fluid Model Near-Wall Averaging and Consistent Matching for Turbulent Bubbly Flows

2016 ◽  
Vol 139 (1) ◽  
Author(s):  
Avinash Vaidheeswaran ◽  
Deoras Prabhudharwadkar ◽  
Paul Guilbert ◽  
John R. Buchanan ◽  
Martin Lopez de Bertodano
Keyword(s):  
Void Fraction ◽  
Fluid Model ◽  
Bubbly Flows ◽  
Wall Region ◽  
Two Phase ◽  
Void Fraction Distribution ◽  
Fraction Distribution ◽  
Two Fluid Model ◽  
Two Fluid ◽  
Near Wall

A new two-fluid model averaging in the near-wall region is proposed to ensure consistent matching of the two-phase k–ε turbulence model with the two-phase logarithmic law of the wall (Marie J. L., Moursali, E., and Tran-Cong, S., 1997, “Similarity Law and Turbulence Intensity Profiles in a Bubbly Boundary Layer,” Int. J. Multiphase Flow, 23(2), pp. 227–247). The void fraction distribution obtained with the averaging procedure is seen to conform to the two-phase wall function approach which is based on a double step function void fraction distribution. In particular, the proposed averaging technique is shown to achieve grid convergence in the near-wall region, which could not be obtained otherwise. Computational fluid dynamics (CFD) results with the proposed technique are in good agreement with experiments on upward bubbly flows over a flat plate, and upward and downward flows in pipes. An additional advantage of the proposed technique is that it replaces the wall force model, which has a significant degree of uncertainty in turbulent flow modeling, with a simpler geometric constraint.

Sign in / Sign up

Export Citation Format

Share Document