Modelling of gas-liquid bubbly flows is achieved by coupling a population balance equation with the three-dimensional, two-fluid, hydrodynamic model. For gas-liquid bubbly flows, an average bubble number density transport equation has been incorporated in the CFD code CFX 5.7 to describe the temporal and spatial evolution of the gas bubbles population. The coalescence and breakage effects of the gas bubbles are modeled. The coalescence by the random collision driven by turbulence and wake entrainment is considered, while for bubble breakage, the impact of turbulent eddies is considered. Local spatial variations of the gas volume fraction, interfacial area concentration, Sauter mean bubble diameter, and liquid velocity are compared against experimental data in a horizontal pipe, covering a range of gas (0.25 to 1.34?m/s) and liquid (3.74 to 5.1?m/s) superficial velocities and average volume fractions (4% to 21%). The predicted local variations are in good agreement with the experimental measurements reported in the literature. Furthermore, the development of the flow pattern was examined at three different axial locations of L/D = 25, 148, and 253. The first location is close to the entrance region where the flow is still developing, while the second and the third represent nearly fully developed bubbly flow patterns. 1. Introduction Gas-liquid, two-phase flow in horizontal pipes is encountered often in a number of industrial processes. Common applications include chemical plants, evaporators, oil wells and pipelines, fluidized bed combustors, and evaporators. Horizontal bubbly flows have received less attention in the literature than vertical flows, even though this flow orientation is equally important in industrial applications such as hydrotransport, an important technology in bitumen extraction. Experimental observations are also difficult in this case, as the migration of dispersed bubbles towards the top of the pipe, due to buoyancy, causes a highly nonsymmetric volume distribution in the pipe cross-section. This density stratification is not often accompanied by a strong secondary flow. Gas volume fraction, interfacial area concentration, and mean bubble diameter are the three characterizing field variables that characterize the internal flow structure of two-phase, gas-liquid flows in horizontal pipe. In various industrial processes, the gas volume fraction parameter is required for hydrodynamic and thermal design. The interfacial transport of mass, momentum, and energy is proportional to the interfacial area and the driving forces. This is an
References
[1]
G. Kocamustafaogullari and W. D. Huang, “Internal structure and interfacial velocity development for bubbly two-phase flow,” Nuclear Engineering and Design, vol. 151, no. 1, pp. 79–101, 1994.
[2]
R. Krishna and J. M. van Baten, “Scaling up bubble column reactors with the aid of CFD,” Chemical Engineering Research and Design, vol. 79, no. 3, pp. 283–309, 2001.
[3]
T. Hibiki and M. Ishii, “Development of one-group interfacial area transport equation in bubbly flow systems,” International Journal of Heat and Mass Transfer, vol. 45, no. 11, pp. 2351–2372, 2002.
[4]
A. Kitagawa, K. Sugiyama, and Y. Murai, “Experimental detection of bubble-bubble interactions in a wall-sliding bubble swarm,” International Journal of Multiphase Flow, vol. 30, no. 10, pp. 1213–1234, 2004.
[5]
Z. Xiao and R. B. H. Tan, “A model for bubble-bubble and bubble-wall interaction in bubble formation,” AIChE Journal, vol. 52, no. 1, pp. 86–98, 2006.
[6]
D. Colella, D. Vinci, R. Bagatin, M. Masi, and E. Abu Bakr, “A study on coalescence and breakage mechanisms in three different bubble columns,” Chemical Engineering Science, vol. 54, no. 21, pp. 4767–4777, 1999.
[7]
E. Olmos, C. Gentric, Ch. Vial, G. Wild, and N. Midoux, “Numerical simulation of multiphase flow in bubble column reactors. Influence of bubble coalescence and break-up,” Chemical Engineering Science, vol. 56, no. 21-22, pp. 6359–6365, 2001.
[8]
F. Lehr, M. Millies, and D. Mewes, “Bubble-size distributions and flow fields in bubble columns,” AIChE Journal, vol. 48, no. 11, pp. 2426–2443, 2002.
[9]
F. B. Campos and P. L. C. Lage, “A numerical method for solving the transient multidimensional population balance equation using an Euler-Lagrange formulation,” Chemical Engineering Science, vol. 58, no. 12, pp. 2725–2744, 2003.
[10]
P. Chen, J. Sanyal, and M. P. Dudukovic, “CFD modeling of bubble columns flows: implementation of population balance,” Chemical Engineering Science, vol. 59, no. 22-23, pp. 5201–5207, 2004.
[11]
J. Sanyal, D. L. Marchisio, O. Fox, and K. Dhanasekharan, “On the comparison between population balance models for CFD simulation of bubble columns,” Industrial and Engineering Chemistry Research, vol. 44, no. 14, pp. 5063–5072, 2005.
[12]
H. A. Jakobsen, H. Lindborg, and C. A. Dorao, “Modeling of bubble column reactors: progress and limitations,” Industrial and Engineering Chemistry Research, vol. 44, no. 14, pp. 5107–5151, 2005.
[13]
Z. Sha, A. Laari, and I. Turunen, “Multi-phase-multi-size group model for the inclusion of population balances into the CFD simulation of gas-liquid bubbly flows,” Chemical Engineering and Technology, vol. 29, no. 5, pp. 550–559, 2006.
[14]
T. Wang, J. Wang, and Y. Jin, “A CFD-PBM coupled model for gas-liquid flows,” AIChE Journal, vol. 52, no. 1, pp. 125–140, 2006.
[15]
S. Lo, “Application of the MUSIG model to bubbly flows,” AEAT -1096, AEA Technology, 1996.
[16]
G. H. Yeoh and J. Y. Tu, “Population balance modelling for bubbly flows with heat and mass transfer,” Chemical Engineering Science, vol. 59, no. 15, pp. 3125–3139, 2004.
[17]
G. Kocamustafaogullari and Z. Wang, “An experimental study on local interfacial parameters in a horizontal bubbly two-phase flow,” International Journal of Multiphase Flow, vol. 17, no. 5, pp. 553–572, 1991.
[18]
H. Luo and H. F. Svendsen, “Theoretical model for drop and bubble break-up in turbulent dispersions,” AIChE Journal, vol. 42, no. 5, pp. 1225–1233, 1996.
[19]
M. J. Prince and H. W. Blanch, “Bubble coalescence and break-up in air-sparged bubble columns,” AIChE Journal, vol. 36, no. 10, pp. 1485–1499, 1990.
[20]
T. O. Oolman and H. W. Blanch, “Bubble coalescence in stagnant liquids,” Chemical Engineering Communications, vol. 43, no. 4–6, pp. 237–261, 1986.
[21]
A. K. Chesters and G. Hofman, “Bubble coalescence in pure liquids,” Applied Scientific Research, vol. 38, no. 1, pp. 353–361, 1982.
[22]
R. D. Kirkpatrick and M. J. Lockett, “The influence of approach velocity on bubble coalescence,” Chemical Engineering Science, vol. 29, no. 12, pp. 2363–2373, 1974.
[23]
J. W. Kim and W. K. Lee, “Coalescence behavior of two bubbles in stagnant liquids,” Journal of Chemical Engineering of Japan, vol. 20, no. 5, pp. 448–453, 1987.
[24]
J. C. Rotta, Turbulente Stromungen, B. G. Teubner, Stuttgart, Germany, 1974.
[25]
H. Anglart and O. Nylund, “CFD application to prediction of void distribution in two-phase bubbly flows in rod bundles,” Nuclear Engineering and Design, vol. 163, no. 1-2, pp. 81–98, 1996.
[26]
R. T. Lahey Jr. and D. A. Drew, “The analysis of two-phase flow and heat transfer using a multidimensional, four field, two-fluid model,” Nuclear Engineering and Design, vol. 204, no. 1–3, pp. 29–44, 2001.
[27]
J. B. Joshi, “Computational flow modelling and design of bubble column reactors,” Chemical Engineering Science, vol. 56, no. 21-22, pp. 5893–5933, 2001.
[28]
R. Clift, J. R. Grace, and M. E. Weber, Bubbles, Drops, and Particles, Academic Press, New York, NY, USA, 1978.
[29]
J. B. Joshi, U. V. Parasu, Ch. V. Prasad, et al., “Gas holdup structures in bubble column reactors,” Proceedings of the Indian National Science Academy, vol. 64, pp. 441–567, 1998.
[30]
I. Zun, “The transverse migration of bubbles influenced by walls in vertical bubbly flow,” International Journal of Multiphase Flow, vol. 6, no. 6, pp. 583–588, 1980.
[31]
N. H. Thomas, T. R. Auton, K. Sene, and J. C. R. Hunt, “Entrapment and transport of bubbles by transient large eddies in turbulent shear flow,” in Proceedings of the BHRA International Conference on the Physical Modelling of Multi-Phase Flow.
[32]
D. A. Drew and S. L. Passman, Theory of Multicomponent Fluids, Springer, New York, NY, USA, 1999.
[33]
A. Tomiyama, “Drag, lift and virtual mass forces acting on a single bubble,” in Proceedings of the 3rd International Symposium on Two-Phase Flow Modelling and Experimentation, pp. 22–24, Pisa, Italy, September 2004.
[34]
S. P. Antal, R. T. Lahey, and J. E. Flaherty, “Analysis of phase distribution in fully developed laminar bubbly two-phase flow,” International Journal of Multiphase Flow, vol. 17, no. 5, pp. 635–652, 1991.
[35]
M. A. Lopez de Bertodano, Turbulent bubbly two-phase flow in a triangular duct, Ph. D. dissertation, Rensselaer Polytechnic Institute, 1992.
[36]
M. Ishii and N. Zuber, “Drag coefficient and relative velocity in bubbly, droplet or particulate flows,” AIChE Journal, vol. 25, no. 5, pp. 843–855, 1979.
[37]
T. R. Auton, “The lift force on a spherical body in a rotational flow,” Journal of Fluid Mechanics, vol. 183, pp. 199–218, 1987.
[38]
M. A. Lopez de Bertodano, R. T. Lahey Jr., and O. C. Jones, “Phase distribution in bubbly two-phase flow in vertical ducts,” International Journal of Multiphase Flow, vol. 20, no. 5, pp. 805–818, 1994.
[39]
M. Lance and M. A. Lopez de Bertodano, “Phase distribution and wall effects in bubbly two-phase flows,” Multiphase Science and Technology, vol. 8, no. 1–4, pp. 69–123, 1994.
[40]
A. Tomiyama, A. Sou, I. Zun, N. Kanami, and T. Sakaguchi, “Effects of Eotvos number and dimensionless liquid volumetric flux on lateral motion of a bubble in a laminar duct flow,” in Advances in Multiphase Flow, pp. 3–15, 1995.
[41]
A. Tomiyama, H. Tamai, I. Zun, and S. Hosokawa, “Transverse migration of single bubbles in simple shear flows,” Chemical Engineering Science, vol. 57, no. 11, pp. 1849–1858, 2002.
[42]
Y. Sato and K. Sekoguchi, “Liquid velocity distribution in two-phase bubble flow,” International Journal of Multiphase Flow, vol. 2, no. 1, pp. 79–95, 1975.
[43]
R. D. S. Cavalcanti, S. R. D. Neto, and E. O. Vilar, “A computational fluid dynamics study of hydrogen bubbles in an electrochemical reactor,” Brazilian Archives of Biology and Technology, vol. 48, pp. 219–229, 2005.
[44]
J. Y. Tu, G. H. Yeoh, G. C. Park, and M. O. Kim, “On population balance approach for subcooled boiling flow prediction,” Journal of Heat Transfer, vol. 127, no. 3, pp. 253–264, 2005.
[45]
D. Legendre and J. Magnaudet, “The lift force on a spherical bubble in a viscous linear shear flow,” Journal of Fluid Mechanics, vol. 368, pp. 81–126, 1998.
[46]
G. Kocamustafaogullari, W. D. Huang, and J. Razi, “Measurement and modeling of average void fraction, bubble size and interfacial area,” Nuclear Engineering and Design, vol. 148, no. 2-3, pp. 437–453, 1994.
[47]
A. Iskandrani and G. Kojasoy, “Local void fraction and velocity field description in horizontal bubbly flow,” Nuclear Engineering and Design, vol. 204, no. 1–3, pp. 117–128, 2001.