We study bubble break-up in homogeneous and isotropic turbulence by direct numerical simulations of the two-phase incompressible Navier–Stokes equations. We create the turbulence by forcing in physical space and introduce the bubble once a statistically stationary state is reached. We perform a large ensemble of simulations to investigate the effect of the Weber number (the ratio of turbulent and surface tension forces) on bubble break-up dynamics and statistics, including the child bubble size distribution, and discuss the numerical requirements to obtain results independent of grid size. We characterize the critical Weber number below which no break-up occurs and the associated Hinze scale $d_h$ . At Weber number close to stable conditions (initial bubble sizes $d_0\approx d_h$ ), we observe binary and tertiary break-ups, leading to bubbles mostly between $0.5d_h$ and $d_h$ , a signature of a production process local in scale. For large Weber numbers ( $d_0> 3d_h$ ), we observe the creation of a wide range of bubble radii, with numerous child bubbles between $0.1d_h$ and $0.3d_h$ , an order of magnitude smaller than the parent bubble. The separation of scales between the parent and child bubble is a signature of a production process non-local in scale. The formation mechanism of these sub-Hinze scale bubbles relates to rapid large deformation and successive break-ups: the first break-up in a sequence leaves highly deformed bubbles which will break again, without recovering a spherical shape and creating an array of much smaller bubbles. We discuss the application of this scenario to the production of sub-Hinze bubbles under breaking waves.
more »
« less
Experimental observations and modelling of sub-Hinze bubble production by turbulent bubble break-up
We 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
- Award ID(s):
- 1844932
- NSF-PAR ID:
- 10388588
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 951
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We present high-resolution three-dimensional (3-D) direct numerical simulations of breaking waves solving for the two-phase Navier–Stokes equations. We investigate the role of the Reynolds number ( Re , wave inertia relative to viscous effects) and Bond number ( Bo , wave scale over the capillary length) on the energy, bubble and droplet statistics of strong plunging breakers. We explore the asymptotic regimes at high Re and Bo , and compare with laboratory breaking waves. Energetically, the breaking wave transitions from laminar to 3-D turbulent flow on a time scale that depends on the turbulent Re up to a limiting value $Re_\lambda \sim 100$ , consistent with the mixing transition in other canonical turbulent flows. We characterize the role of capillary effects on the impacting jet and ingested main cavity shape and subsequent fragmentation process, and extend the buoyant-energetic scaling from Deike et al. ( J. Fluid Mech. , vol. 801, 2016, pp. 91–129) to account for the cavity shape and its scale separation from the Hinze scale, $r_H$ . We confirm two regimes in the bubble size distribution, $N(r/r_H)\propto (r/r_H)^{-10/3}$ for $r>r_H$ , and $\propto (r/r_H)^{-3/2}$ for $rnull (Ed.)We present simultaneous three-dimensional measurements of deformable finite-sized bubbles and surrounding turbulent flows. The orientations of bubbles are linked to two key mechanisms that drive bubble deformation: the turbulent strain rate and slip velocity between the two phases. The strongest preferential alignment is between the bubbles and slip velocity, indicating the latter plays a dominant role. We also compared our experimental results with the deformation of ideal material elements with no slip velocity or surface tension. Without these, material elements show highly different orientations, further confirming the importance of the slip velocity in the bubble orientation. In addition to deformation, when bubbles begin to break, their relative orientations change significantly. Although the alignment of the severely deformed bubbles with the eigenvectors of the turbulent strain rate becomes much stronger, the bubble semi-major axis becomes aligned with (rather than perpendicular to) the slip velocity through an almost $90^{\circ }$ turn. This puzzling orientation change occurs because the slip velocity contains the contributions from both the bubble and the background flow. As the bubble experiences strong deformation, the rapid elongation of its semi-major axis leads to a large bubble velocity, which dominates the slip velocity and forces its alignment with the bubble's semi-major axis. The slip velocity thereby switches from a driving mechanism to a driven result as bubbles approach breakup. The results highlight the complex coupling between the bubble orientation and the surrounding flow, which should be included when modelling the bubble deformation and breakup in turbulence.more » « lessmore » « less
We conducted a spectral analysis of the turbulence kinetic energy (TKE) budget in a bubble plume using particle image velocimetry with fluorescent particles. Our findings confirmed the hypothesis of an inverse energy cascade in the bubble plume, where TKE is transferred from small to large eddies. This is attributed to direct injection of TKE by bubble passages across a wide range of scales, in contrast to canonical shear production of TKE in large scales. Turbulence dissipation was identified as the primary sink of the bubble-produced TKE and occurred at all scales. The decomposition of velocities using the critical length scale of inter-scale energy transfer allowed us to distinguish between large- and small-scale motions in the bubble plume. The large-scale turbulent fluctuations exhibited a skewed distribution and were likely associated with the return flow after bubble passage and the velocities induced by the bubble wake. The small-scale turbulent fluctuations followed a Gaussian distribution relatively well. The large-scale motions contributed to over half of the Reynolds stresses, while there were significant small-scale contributions to the normal stresses near the plume center but not to the shear stress. The large-scale motions in the vorticity field induced a street of vertically elongated vortex pairs, while the small-scale vortices exhibited similar sizes in both horizontal and vertical directions.
more » « lessAbstract Nanoparticle‐shelled bubbles, prepared with glass capillary microfluidics, are functionalized to produce catalytic micromotors that exhibit novel assembly and disassembly behaviors. Stable microbubble rafts are assembled at an air–solvent interface of nonaqueous propylene carbonate (PC) solvent by creating a meniscus using a glass capillary. Upon the addition of hydrogen peroxide fuel, catalytic microbubbles quickly break free from the bubble raft by repelling from each other and self‐propelling at the air–fuel interface (a mixture of PC and aqueous hydrogen peroxide). While most of micromotors generate oxygen bubbles on the outer catalytic shell, some micromotors contain cracks and eject bubbles from the hollow shells containing air. Nanoparticle‐shelled bubbles with a high buoyancy force are particularly attractive for studying novel propulsion modes and interactions between catalytic bubble micromotors at air–fuel interfaces.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2022.604
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.
