Efficient simulation of gas-liquid pipe flows using a generalized population balance equation coupled with the algebraic slip model

The inhomogeneous generalized population balance equation, which is discretized with the direct quadrature method of moment (DQMOM), is solved to predict the bubble size distribution (BSD) in a vertical pipe flow. The proposed model is compared with a more classical approach where bubbles are characterized with a constant mean size. The turbulent two-phase flow field, which is modeled using a Reynolds-Averaged Navier-Stokes equation approach, is assumed to be in local equilibrium, thus the relative gas and liquid (slip) velocities can be calculated with the algebraic slip model, thereby accounting for the drag, lift, and lubrication forces. The complex relationship between the bubble size distribution and the resulting forces is described accurately by the DQMOM. Each quadrature node and weight represents a class of bubbles with characteristic size and number density, which change dynamically in time and space to preserve the first moments of the BSD. The predictions obtained are validated against previously published experimental data, thereby demonstrating the advantages of this approach for large-scale systems as well as suggesting future extensions to long piping systems and more complex geometries.