Multi-fluid simulation of turbulent bubbly pipe flows

Radial distributions of void fraction αGαG, bubble aspect ratio E, phasic velocities VG and VL and turbulent kinetic energy k in bubbly pipe flows are measured using an image processing method and a laser Doppler velocimetry. Multi-fluid simulations are conducted to examine applicability of state-of-the-art closure relations to the turbulent bubbly pipe flows. The experimental results indicate that aspect ratio of bubbles in the near wall region takes a higher value than that of free rising bubbles due to the presence of wall, and that the change in the aspect ratio induces decrease in relative velocity between bubbles and liquid in the near wall region. Drag coefficient CD of a bubble in a bubbly pipe flow tends to increase with magnitude of shear flow, and the effect of shear flow on CD is estimated by the correlation proposed by Legendre and Magnaudet (1998). Comparison between the simulated and the measured results indicate that the effects of bubble shape and shear flow on drag force acting on bubbles should be taken into account for accurate predictions of bubbly pipe flows. The turbulence models proposed by Lopez de Bertodano et al. (1994) and by Hosokawa and Tomiyama (2004a) give good predictions for turbulence modification caused by bubbles.