Direct numerical simulation of bubble rising in turbulence

We describe the rising trajectory of bubbles in isotropic turbulence and quantify the slowdown of the mean rise velocity of bubbles with sizes within the inertial subrange. We perform direct numerical simulations of bubbles, for a wide range of turbulence intensity, bubble inertia and deformability, with systematic comparison with the corresponding quiescent case, with Reynolds number at the Taylor microscale from 38 to 77. Turbulent fluctuations randomise the rising trajectory and cause a reduction of the mean rise velocity w∼ _b compared with the rise velocity in quiescent flow w_b. The decrease in mean rise velocity of bubbles w∼ _b/w_b is shown to be primarily a function of the ratio of the turbulence intensity and the buoyancy forces, described by the Froude number Fr = u^'/√gd, where u^' is the root-mean-square velocity fluctuations, g is gravity and d is the bubble diameter. The bubble inertia, characterised by the ratio of inertial to viscous forces (Galileo number), and the bubble deformability, characterised by the ratio of buoyancy forces to surface tension (Bond number), modulate the rise trajectory and velocity in quiescent fluid. The slowdown of these bubbles in the inertial subrange is not due to preferential sampling, as is the case with sub-Kolmogorov bubbles. Instead, it is caused by the nonlinear drag–velocity relationship, where velocity fluctuations lead to an increased average drag. For Fr > 0.5, we confirm the scaling w∼ _b/w_b ∝ 1/Fr, as proposed previously by Ruth et al. (J. Fluid Mech., vol. 924, 2021, p. A2), over a wide range of bubble inertia and deformability.