We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
We investigate air-entraining flows where degassing, rather than fragmentation, plays a significant role. Of interest is the power-law slope 
$\beta$ of the bulk bubble size distribution 
$N(a)$ during the air-generating period, when the total volume of bubbles is increasing. We study a canonical air-entraining flow created by strong underlying free-surface turbulence. We perform analysis using the population balance equation (PBE) and computations using direct numerical simulations (DNS) with bubble tracking. We quantify the importance of degassing by the ratio of degassing flux (
$Q_D$) to entrainment flux (
$Q_I$), 
$\mathcal {D}=Q_D/Q_I$, and the ratio of degassing rate (
$\varLambda (a)$) to fragmentation rate (
$\varOmega (a)$) for a bubble of radius 
$a$, 
$\varLambda (a)/\varOmega (a)$. For a broad range of large Froude numbers 
${{Fr}}=U/\sqrt {L g}$, DNS give 
$\mathcal {D}=\operatorname {O}(1)$ (independent of 
${{Fr}}$), showing that degassing is relevant, and 
$\varLambda (a) \gg \varOmega (a)$, showing that the bubble population is degassing-dominated. In contrast to fragmentation-dominated populations, such as those due to wave breaking, where 
$\beta =-10/3$, degassing-dominated populations have qualitatively different 
$N(a)$ during air entrainment. Analysis using the PBE shows that degassing-dominated 
$\beta$ is a function of 
$\varLambda (a)$, which has a turbulence-driven regime (
$a< a_\varLambda$) and a buoyancy-driven regime (
$a>a_\varLambda$). Here, 
$a_\varLambda$ is the bubble radius where terminal buoyant rise velocity equals 
$u_{rms}$. Consequently, 
$N(a)$ exhibits a split power with 
$\beta (a< a_\varLambda )=-4.\bar {3}$ and 
$\beta (a>a_\varLambda )=-5.8\bar {3}$ for moderate bubble Reynolds numbers 
${{Re}}_b$. For large 
${{Re}}_b$, 
$\beta (a>a_\varLambda )=-4.8\bar {3}$. The DNS strongly confirm these findings for moderate 
${{Re}}_b$. By identifying and describing degassing-dominated bubble populations, this work contributes to the understanding and interpretation of broad types of air-entraining problems where degassing plays a relevant role.
We investigate the fundamental time scales that characterise the statistics of fragmentation under homogeneous isotropic turbulence for air–water bubbly flows at moderate to large bubble Weber numbers, 
$We$. We elucidate three time scales: 
$\tau _r$, the characteristic age of bubbles when their subsequent statistics become stationary; 
$\tau _\ell$, the expected lifetime of a bubble before further fragmentation; and 
$\tau _c$, the expected time for the air within a bubble to reach the Hinze scale, radius 
$a_H$, through the fragmentation cascade. The time scale 
$\tau _\ell$ is important to the population balance equation (PBE), 
$\tau _r$ is critical to evaluating the applicability of the PBE no-hysteresis assumption, and 
$\tau _c$ provides the characteristic time for fragmentation cascades to equilibrate. By identifying a non-dimensionalised average speed 
$\bar {s}$ at which air moves through the cascade, we derive 
$\tau _c=C_\tau \varepsilon ^{-1/3} a^{2/3} (1-(a_{max}/a_H)^{-2/3})$, where 
$C_\tau =1/\bar {s}$ and 
$a_{max}$ is the largest bubble radius in the cascade. While 
$\bar {s}$ is a function of PBE fragmentation statistics, which depend on the measurement interval 
$T$, 
$\bar {s}$ itself is independent of 
$T$ for 
$\tau _r \ll T \ll \tau _c$. We verify the 
$T$-independence of 
$\bar {s}$ and its direct relationship to 
$\tau _c$ using Monte Carlo simulations. We perform direct numerical simulations (DNS) at moderate to large bubble Weber numbers, 
$We$, to measure fragmentation statistics over a range of 
$T$. We establish that non-stationary effects decay exponentially with 
$T$, independent of 
$We$, and provide 
$\tau _r=C_{r} \varepsilon ^{-1/3} a^{2/3}$ with 
$C_{r}\approx 0.11$. This gives 
$\tau _r\ll \tau _\ell$, validating the PBE no-hysteresis assumption. From DNS, we measure 
$\bar {s}$ and find that for large Weber numbers (
$We>30$), 
$C_{\tau }\approx 9$. In addition to providing 
$\tau _c$, this obtains a new constraint on fragmentation models for PBE.
