Scale separation and dependence of entrainment bubble-size distribution in free-surface turbulence
Xiangming Yu,Kelli Hendrickson,Dick K. P. Yue
DOI: https://doi.org/10.1017/jfm.2019.986
IF: 3.7
2019-12-27
Journal of Fluid Mechanics
Abstract:We consider the size spectrum of entrained bubbles under strong free-surface turbulence (SFST). We investigate the entrainment bubble-size spectrum per unit (mean) interface area, ${\mathcal{N}}_{a}^{E}(r)$ , with dimension length $^{-3}$ , and develop a physical/mechanistic model for ${\mathcal{N}}_{a}^{E}(r)$ through energy arguments. The model obtains two distinct regimes of ${\mathcal{N}}_{a}^{E}(r)$ , separated by bubble-size scale $r_{0}$ . For bubble radius $r>r_{0}$ , the effects of gravity $g$ dominate those of the surface tension force $\unicode[STIX]{x1D70E}/\unicode[STIX]{x1D70C}$ , and ${\mathcal{N}}_{a}^{E}(r)\propto g^{-1}\unicode[STIX]{x1D716}^{2/3}r^{-10/3}$ , where $\unicode[STIX]{x1D716}$ is the turbulence dissipation rate. For $r<r_{0}$ , surface tension is more important and ${\mathcal{N}}_{a}^{E}(r)\propto (\unicode[STIX]{x1D70E}/\unicode[STIX]{x1D70C})^{-1}\unicode[STIX]{x1D716}^{2/3}r^{-4/3}$ . From the model, we show that $r_{0}\approx r_{c}=1/2\sqrt{\unicode[STIX]{x1D70E}/\unicode[STIX]{x1D70C}g}$ , the capillary length scale, and not the generally assumed Hinze scale $r_{H}$ . For an air–water interface and Earth gravity, $r_{c}\approx$ 1.5 mm. The model provides an $\unicode[STIX]{x1D716}$ – $r$ entrainment regime map that identifies a critical dissipation rate $\unicode[STIX]{x1D716}_{cr}$ (constant for given $g$ and $\unicode[STIX]{x1D70E}/\unicode[STIX]{x1D70C}$ ) above which there is appreciable air entrainment, thus separating SFST and weak FST. We confirm the theoretical model and its predictions using two-phase, high-fidelity direct numerical simulations of a canonical FST flow using the conservative volume-of-fluid method: the respective power laws of ${\mathcal{N}}_{a}^{E}(r)\propto r^{-10/3}$ and $r^{-4/3}$ for $r>r_{0}$ and $r<r_{0}$ ; the value $r_{0}=r_{c}$ ; the scaling ${\mathcal{N}}_{a}^{E}(r)\propto \unicode[STIX]{x1D716}^{2/3}$ ; and the predictions of the $\unicode[STIX]{x1D716}$ – $r$ entrainment regime map.
mechanics,physics, fluids & plasmas