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

[+] Author Affiliations
M.-L. Bordas, A. Cartellier, P. Séchet

Laboratoire des Ecoulements Géophysiques et Industriels, Grenoble, France

Paper No. FEDSM2002-31382, pp. 669-670; 2 pages
  • ASME 2002 Joint U.S.-European Fluids Engineering Division Conference
  • Volume 1: Fora, Parts A and B
  • Montreal, Quebec, Canada, July 14–18, 2002
  • Conference Sponsors: Fluids Engineering Division
  • ISBN: 0-7918-3615-0 | eISBN: 0-7918-3600-2
  • Copyright © 2002 by ASME


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:

Liquid phase:
  −dpdz = ρLg + fLS − fLG1 − α    (1)
Gas phase:
  −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 Ur 2 π 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).

Copyright © 2002 by ASME
Topics: Fluids , Bubbly flow
This article is only available in the PDF format.



Interactive Graphics


Country-Specific Mortality and Growth Failure in Infancy and Yound Children and Association With Material Stature

Use interactive graphics and maps to view and sort country-specific infant and early dhildhood mortality and growth failure data and their association with maternal

Citing articles are presented as examples only. In non-demo SCM6 implementation, integration with CrossRef’s "Cited By" API will populate this tab (http://www.crossref.org/citedby.html).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In