An analytical relation for the void fraction distribution in a fully developed bubbly flow in a vertical pipe. (2nd October 2016)

• Record Type: Journal Article
• Authors: Marfaing, O.; Guingo, M.; Laviéville, J.; Bois, G.; Méchitoua, N.; Mérigoux, N.; Mimouni, S.
• Is Part Of: Chemical engineering science. Volume 152(2016)
• Page Start: 579
• Page End: 585
• Publication Date: 2016-10-02
• Subjects: Two-fluid model -- Analytical relation -- Verification -- Wall force -- Lift force -- Dispersion force -- Drag force -- NEPTUNE_CFD; Chemical engineering -- Periodicals; Génie chimique -- Périodiques; Chemical engineering; Periodicals; Electronic journals; 660
• DOI: 10.1016/j.ces.2016.06.041
• Languages: English
• ISSNs: 0009-2509

Abstract

Abstract: The problem of a steady, axisymmetric, fully developed adiabatic bubbly flow in a vertical pipe is studied analytically with the two-fluid model. The exchange of momentum between the phases is described as the sum of drag, lift, wall and dispersion contributions, with constant coefficients. Under these conditions, we are able to express analytically the void fraction profile as a function of the liquid velocity and pressure profiles. This relation is valid independently of the Reynolds stress model in the liquid phase – and can serve as a verification case for multiphase flow codes. The analytical void fraction profile vanishes at the wall, as a result of the balance between dispersion and wall forces. It presents a peak near the wall for upward flows, whereas its maximum is reached in the center of the pipe for downward flows. This is illustrated by calculations performed for upward and downward bubbly flows with the NEPTUNE_CFD code. Highlights: We perform an analytical study of an upward bubbly flow in a pipe. Drag, dispersion, lift, and wall forces are considered with constant coefficients. The void fraction profile is expressed analytically as a function of the liquid velocity and pressure. The relation is valid for high Reynolds number flows. At the wall, the void fraction is zero, due to the singular wall force. Numerical calculations are performed with a CFD code.