Void distribution and bubble motion in bubbly flows in a 4×4 rod bundle. Part II: numerical simulation

Abstract

Numerical simulations of bubbly flows in a four by four rod bundle are carried out using a multi-fluid model to examine effects of the numerical treatment of phase distribution and drag model. The transport equations of bubble number density and void fraction are used as the continuity equation of the gas phase. Two drag models are tested: one of them accounts for the bubble deformation (aspect ratio), whereas the other does not. The rod diameter, the rod pitch and the hydraulic diameter of the rod bundle are 10, 12.5 and 9.1 mm, respectively. The gas and liquid volume fluxes are JG = 0.06 m/s and JL = 0.9 and 1.5 m/s, respectively. The bubble diameter ranges from 1 to 5 mm. Comparisons between the numerical and measured data show that (1) the restriction on bubble lateral motion due to the presence of rods can be taken into account by using the transport equation of bubble number density, whereas that of the void fraction cannot deal with the restriction and causes large errors in the distribution of void fraction and (2) the reduction in the bubble-relative velocity near the wall is predictable by using the drag model accounting for the bubble deformation effect.