In this paper results of direct numerical simulation (DNS) of bubbles rising in viscous Newtonian liquids with high-density ratio are presented. The simulations are carried out with the highly parallelized code FS3D, which employs the Volume-of-Fluid (VOF) method. The high degree of parallization of the code allows high resolution of the computational domain, such that the Kolmogorov length scale inside the liquid phase is resolved for the simulations. For validation of the numerical results the terminal rise velocities, bubble shapes and flow fields are compared to experimental data as well as to approximate analytical solutions. For high Morton numbers terminal rise velocities and aspect ratios agree very well with experimental values. For lower Morton numbers there is an increasing difference between experimental and numerical rise velocities. The aspect ratios of ellipsoidal bubbles match both with experimental measurements and with theoretical values of Taylor and Acrivos. At very low Reynolds numbers (ReB < 1) the velocity fields in and outside of the bubble show good semi-quantitative agreement with the analytical creeping flow solution of Hadamard and Rybczynski.

This content is only available via PDF.
You do not currently have access to this content.