Direct numerical simulations of bubble-mediated gas transfer and dissolution in quiescent and turbulent flows

We perform direct numerical simulations of a gas bubble dissolving in a surrounding liquid. The bubble volume is reduced due to dissolution of the gas, with the numerical implementation of an immersed boundary method, coupling the gas diffusion and the Navier-Stokes equations. The methods are validated against planar and spherical geometries' analytical moving boundary problems, including the classic Epstein-Plesset problem. Considering a bubble rising in a quiescent liquid, we show that the mass transfer coefficient can be described by the classic Levich formula, with and the time-varying bubble size and rise velocity, and the gas diffusivity in the liquid. Next, we investigate the dissolution and gas transfer of a bubble in homogeneous and isotropic turbulence flow, extending Farsoiya et al. (J. Fluid Mech., vol. 920, 2021, A34). We show that with a bubble size initially within the turbulent inertial subrange, the mass transfer coefficient in turbulence is controlled by the smallest scales of the flow, the Kolmogorov and Batchelor microscales, and is independent of the bubble size. This leads to the non-dimensional transfer rate scaling as, where is the macroscale Reynolds number, with the velocity fluctuations, the integral length scale, the liquid viscosity, and the Schmidt number. This scaling can be expressed in terms of the turbulence dissipation rate as, in agreement with the model proposed by Lamont and Scott (AIChE J., vol. 16, issue 4, 1970, pp. 513-519) and corresponding to the high regime from Theofanous et al. (Intl J. Heat Mass Transfer, vol. 19, issue 6, 1976, pp. 613-624).