Numerical simulation of multiphase flow in bubble column reactors. Influence of bubble coalescence and break-up

Abstract

Population balance equations have been combined to a classical hydrodynamic Euler/Euler simulation to investigate the operation of a cylindrical bubble column. The MUSIG (mutiple-size-group) model implemented in the CFX 4.3 commercial software has been used. Hydrodynamic experimental variables, i.e. local axial liquid velocity and local gas hold-up, have been compared to the corresponding calculated values, showing a quite good agreement, except for the gas hold-up when the column is no more operating in the homogeneous regime. Bubble sizes have been investigated, showing that two domains of superficial gas velocities can be distinguished. In the first domain, coalescence occurs predominantly, Sauter diameter increases with the superficial gas velocity, bubble size distribution is narrow and Sauter diameter is continuously evolving along the column axis. In the second domain, break-up becomes more intensive and compensates coalescence, bubble size distribution becomes wider, since more small bubbles are formed, an equilibrium Sauter diameter appears when the superficial gas velocity increases. Furthermore an equilibrium Sauter diameter appears along the column axis, and it can be noticed that this phenomenon appears lower in the column when the gas flow rate is increased. In these two domains the characteristics of the bubbles are typical of those of the homogeneous and transition regimes.