Highly scalable DNS solver for turbulent bubble-laden channel flow

A simulation of a downflow configuration of a turbulent channel loaded with a total of 10000 bubbles illustrates the computational capabilities. First and second order statistics of the velocity field were computed as well as the profiles of the average gas volume fraction field in the statistically steady state. In the case considered, a 47% increase of the wall shear stress was observed, brought about by turbulence modification arising from the embedded bubbles. The interaction between turbulence and bubbles at high volume fraction resulted in a strong attenuation of the rms of the velocity fluctuations in a wide region of the core of the channel.