A mathematical model for turbulent gas-liquid flows
with mass transfer and chemical reactions is presented
and a robust solution strategy based on nested iterations
is proposed for the numerical treatment of the intricately
coupled PDEs. In particular, the incompressible
Navier-Stokes equations are solved by a discrete projection
scheme from the family of Pressure Schur Complement
methods. Novel high-resolution finite element schemes
of FCT and TVD type are employed for the discretization
of unstable convective terms. The implementation of a
modified k-epsilon turbulence model is described in
detail. The performance of the developed simulation
tools is illustrated by a number of three-dimensional
numerical examples.