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 $k_L$ can be described by the classic Levich formula $k_L = (2/\sqrt {{\rm \pi} })\sqrt {\mathscr {D}_l\,U(t)/d(t)}$ , with $d(t)$ and $U(t)$ the time-varying bubble size and rise velocity, and $\mathscr {D}_l$ 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 $k_L$ is controlled by the smallest scales of the flow, the Kolmogorov $\eta$ and Batchelor $\eta _B$ microscales, and is independent of the bubble size. This leads to the non-dimensional transfer rate ${Sh}=k_L L^\star /\mathscr {D}_l$ scaling as ${Sh}/{Sc}^{1/2} \propto {Re}^{3/4}$ , where ${Re}$ is the macroscale Reynolds number ${Re} = u_{rms}L^\star /\nu _l$ , with $u_{rms}$ the velocity fluctuations, $L^*$ the integral length scale, $\nu _l$ the liquid viscosity, and ${Sc}=\nu _l/\mathscr {D}_l$ the Schmidt number. This scaling can be expressed in terms of the turbulence dissipation rate $\epsilon$ as ${k_L}\propto {Sc}^{-1/2} (\epsilon \nu _l)^{1/4}$ , in agreement with the model proposed by Lamont & Scott ( AIChE J. , vol. 16, issue 4, 1970, pp. 513–519) and corresponding to the high $Re$ regime from Theofanous et al. ( Intl J. Heat Mass Transfer , vol. 19, issue 6, 1976, pp. 613–624).
more »
« less
The effect of nonlinear drag on the rise velocity of bubbles in turbulence
We investigate how turbulence in liquid affects the rising speed of gas bubbles within the inertial range. Experimentally, we employ stereoscopic tracking of bubbles rising through water turbulence created by the convergence of turbulent jets and characterized with particle image velocimetry performed throughout the measurement volume. We use the spatially varying, time-averaged mean water velocity field to consider the physically relevant bubble slip velocity relative to the mean flow. Over a range of bubble sizes within the inertial range, we find that the bubble mean rise velocity $\left \langle v_z \right \rangle$ decreases with the intensity of the turbulence as characterized by its root-mean-square fluctuation velocity, $u'$ . Non-dimensionalized by the quiescent rise velocity $v_{q}$ , the average rise speed follows $\left \langle v_z \right \rangle /v_{q}\propto 1/{\textit {Fr}}$ at high ${\textit {Fr}}$ , where ${\textit {Fr}}=u'/\sqrt {dg}$ is a Froude number comparing the intensity of the turbulence to the bubble buoyancy, with $d$ the bubble diameter and $g$ the acceleration due to gravity. We complement these results by performing numerical integration of the Maxey–Riley equation for a point bubble experiencing nonlinear drag in three-dimensional, homogeneous and isotropic turbulence. These simulations reproduce the slowdown observed experimentally, and show that the mean magnitude of the slip velocity is proportional to the large-scale fluctuations of the flow velocity. Combining the numerical estimate of the slip velocity magnitude with a simple theoretical model, we show that the scaling $\left \langle v_z \right \rangle /v_{q}\propto 1/{\textit {Fr}}$ originates from a combination of the nonlinear drag and the nearly isotropic behaviour of the slip velocity at large ${\textit {Fr}}$ that drastically reduces the mean rise speed.
more »
« less
- Award ID(s):
- 1844932
- NSF-PAR ID:
- 10308644
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 924
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We experimentally investigate the rise velocity of finite-sized bubbles in turbulence with a high energy dissipation rate of $\unicode[STIX]{x1D716}\gtrsim 0.5~\text{m}^{2}~\text{s}^{-3}$ . In contrast to a 30–40 % reduction in rise velocity previously reported in weak turbulence (the Weber number ( $We$ ) is much smaller than the Eötvös number ( $Eo$ ); $We\ll 1null (Ed.)We experimentally investigate the breakup mechanisms and probability of Hinze-scale bubbles in turbulence. The Hinze scale is defined as the critical bubble size based on the critical mean Weber number, across which the bubble breakup probability was believed to have an abrupt transition from being dominated by turbulence stresses to being suppressed completely by the surface tension. In this work, to quantify the breakup probability of bubbles with sizes close to the Hinze scale and to examine different breakup mechanisms, both bubbles and their surrounding tracer particles were simultaneously tracked. From the experimental results, two Weber numbers, one calculated from the slip velocity between the two phases and the other acquired from local velocity gradients, are separated and fitted with models that can be linked back to turbulence characteristics. Moreover, we also provide an empirical model to link bubble deformation to the two Weber numbers by extending the relationship obtained from potential flow theory. The proposed relationship between bubble aspect ratio and the Weber numbers seems to work consistently well for a range of bubble sizes. Furthermore, the time traces of bubble aspect ratio and the two Weber numbers are connected using the linear forced oscillator model. Finally, having access to the distributions of these two Weber numbers provides a unique way to extract the breakup probability of bubbles with sizes close to the Hinze scale.more » « lessAbstract The transverse-momentum $$(p_{\textrm{T}})$$ ( p T ) spectra of K $$^{*}(892)^{0}~$$ ∗ ( 892 ) 0 and $$\mathrm {\phi (1020)}~$$ ϕ ( 1020 ) measured with the ALICE detector up to $$p_{\textrm{T}} $$ p T = 16 GeV/ c in the rapidity range $$-1.2< y < 0.3,$$ - 1.2 < y < 0.3 , in p–Pb collisions at the center-of-mass energy per nucleon–nucleon collision $$\sqrt{s_{\textrm{NN}}} = 5.02$$ s NN = 5.02 TeV are presented as a function of charged particle multiplicity and rapidity. The measured $$p_{\textrm{T}} $$ p T distributions show a dependence on both multiplicity and rapidity at low $$p_{\textrm{T}} $$ p T whereas no significant dependence is observed at high $$p_{\textrm{T}} $$ p T . A rapidity dependence is observed in the $$p_{\textrm{T}} $$ p T -integrated yield (d N /d y ), whereas the mean transverse momentum $$\left( \langle p_{\textrm{T}} \rangle \right) $$ ⟨ p T ⟩ shows a flat behavior as a function of rapidity. The rapidity asymmetry ( $$Y_{\textrm{asym}}$$ Y asym ) at low $$p_{\textrm{T}} $$ p T (< 5 GeV/ c ) is more significant for higher multiplicity classes. At high $$p_{\textrm{T}} $$ p T , no significant rapidity asymmetry is observed in any of the multiplicity classes. Both K $$^{*}(892)^{0}~$$ ∗ ( 892 ) 0 and $$\mathrm {\phi (1020)}~$$ ϕ ( 1020 ) show similar $$Y_{\textrm{asym}}$$ Y asym . The nuclear modification factor $$(Q_{\textrm{CP}})$$ ( Q CP ) as a function of $$p_{\textrm{T}} $$ p T shows a Cronin-like enhancement at intermediate $$p_{\textrm{T}} $$ p T , which is more prominent at higher rapidities (Pb-going direction) and in higher multiplicity classes. At high $$p_{\textrm{T}}$$ p T (> 5 GeV/ $$c$$ c ), the $$Q_{\textrm{CP}}$$ Q CP values are greater than unity and no significant rapidity dependence is observed.more » « lessWe present experiments on large air cavities spanning a wide range of sizes relative to the Hinze scale $d_{H}$ , the scale at which turbulent stresses are balanced by surface tension, disintegrating in turbulence. For cavities with initial sizes $d_0$ much larger than $d_{H}$ (probing up to $d_0/d_{H} = 8.3$ ), the size distribution of bubbles smaller than $d_{H}$ follows $N(d) \propto d^{-3/2}$ , with $d$ the bubble diameter. The capillary instability of ligaments involved in the deformation of the large bubbles is shown visually to be responsible for the creation of the small bubbles. Turning to dynamical, three-dimensional measurements of individual break-up events, we describe the break-up child size distribution and the number of child bubbles formed as a function of $d_0/d_{H}$ . Then, to model the evolution of a population of bubbles produced by turbulent bubble break-up, we propose a population balance framework in which break-up involves two physical processes: an inertial deformation to the parent bubble that sets the size of large child bubbles, and a capillary instability that sets the size of small child bubbles. A Monte Carlo approach is used to construct the child size distribution, with simulated stochastic break-ups constrained by our experimental measurements and the understanding of the role of capillarity in small bubble production. This approach reproduces the experimental time evolution of the bubble size distribution during the disintegration of large air cavities in turbulence.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2021.556
Warning: Leaving National Science Foundation Website
You are now leaving the National Science Foundation website to go to a non-government website.
Website:
NSF takes no responsibility for and exercises no control over the views expressed or the accuracy of the information contained on this site. Also be aware that NSF's privacy policy does not apply to this site.
