skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Surface bubble coalescence
We present an experimental study of bubble coalescence at an air–water interface and characterize the evolution of both the underwater neck and the surface bridge. We explore a wide range of Bond number, $Bo$ , which compares gravity and capillary forces and is a dimensionless measure of the free surface's effect on bubble geometry. The nearly spherical $$Bo\ll 1$$ bubbles exhibit the same inertial–capillary growth of the classic underwater dynamics, with limited upper surface displacement. For $Bo>1$ , the bubbles are non-spherical – residing predominantly above the free surface – and, while an inertial–capillary scaling for the underwater neck growth is still observed, the controlling length scale is defined by the curvature of the bubbles near their contact region. With it, an inertial–capillary scaling collapses the neck contours across all Bond numbers to a universal shape. Finally, we characterize the upper surface with a simple oscillatory model which balances capillary forces and the inertia of liquid trapped at the centre of the liquid-film surface.  more » « less
Award ID(s):
1849762
PAR ID:
10308881
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
915
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, Journal of Fluid Mechanics)
    null (Ed.)
    We study bubble motion in a vertical capillary tube under an external flow. Bretherton ( J. Fluid Mech. , vol. 10, issue 2, 1961, pp. 166–188) has shown that, without external flow, a bubble can spontaneously rise when the Bond number ( $${Bo} \equiv \rho g R^2 / \gamma$$ ) is above the critical value $${Bo}_{cr}=0.842$$ , where $$\rho$$ is the liquid density, $$g$$ the gravitational acceleration, $$R$$ the tube radius and $$\gamma$$ the surface tension. It was then shown by Magnini et al. ( Phys. Rev. Fluids , vol. 4, issue 2, 2019, 023601) that the presence of an imposed liquid flow, in the same (upward) direction as buoyancy, accelerates the bubble and thickens the liquid film around it. In this work we carry out a systematic study of the bubble motion under a wide range of upward and downward external flows, focusing on the inertialess regime with Bond numbers above the critical value. We show that a rich variety of bubble dynamics occurs when an external downward flow is applied, opposing the buoyancy-driven rise of the bubble. We reveal the existence of a critical capillary number of the external downward flow ( $${Ca}_l \equiv \mu U_l/\gamma$$ , where $$\mu$$ is the fluid viscosity and $$U_l$$ is the mean liquid speed) at which the bubble arrests and changes its translational direction. Depending on the relative direction of gravity and the external flow, the thickness of the film separating the bubble surface and the tube inner wall follows two distinct solution branches. The results from theory, experiments and numerical simulations confirm the existence of the two solution branches and reveal that the two branches overlap over a finite range of $${Ca}_l$$ , thus suggesting non-unique, history-dependent solutions for the steady-state film thickness under the same external flow conditions. Furthermore, inertialess symmetry-breaking shape profiles at steady state are found as the bubble transits near the tipping points of the solution branches, which are shown in both experiments and three-dimensional numerical simulations. 
    more » « less
  2. ; ; ; (, Journal of Fluid Mechanics)
    We describe the rising trajectory of bubbles in isotropic turbulence and quantify the slowdown of the mean rise velocity of bubbles with sizes within the inertial subrange. We perform direct numerical simulations of bubbles, for a wide range of turbulence intensity, bubble inertia and deformability, with systematic comparison with the corresponding quiescent case, with Reynolds number at the Taylor microscale from 38 to 77. Turbulent fluctuations randomise the rising trajectory and cause a reduction of the mean rise velocity$$\tilde {w}_b$$compared with the rise velocity in quiescent flow$$w_b$$. The decrease in mean rise velocity of bubbles$$\tilde {w}_b/w_b$$is shown to be primarily a function of the ratio of the turbulence intensity and the buoyancy forces, described by the Froude number$$Fr=u'/\sqrt {gd}$$, where$$u'$$is the root-mean-square velocity fluctuations,$$g$$is gravity and$$d$$is the bubble diameter. The bubble inertia, characterised by the ratio of inertial to viscous forces (Galileo number), and the bubble deformability, characterised by the ratio of buoyancy forces to surface tension (Bond number), modulate the rise trajectory and velocity in quiescent fluid. The slowdown of these bubbles in the inertial subrange is not due to preferential sampling, as is the case with sub-Kolmogorov bubbles. Instead, it is caused by the nonlinear drag–velocity relationship, where velocity fluctuations lead to an increased average drag. For$$Fr > 0.5$$, we confirm the scaling$$\tilde {w}_b / w_b \propto 1 / Fr$$, as proposed previously by Ruthet al.(J. Fluid Mech., vol. 924, 2021, p. A2), over a wide range of bubble inertia and deformability. 
    more » « less
  3. ; ; (, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences)
    Three-phase suspensions, of liquid that suspends dispersed solid particles and gas bubbles, are common in both natural and industrial settings. Their rheology is poorly constrained, particularly for high total suspended fractions (≳0.5). We use a dam-break consistometer to characterize the rheology of suspensions of (Newtonian) corn syrup, plastic particles and CO 2 bubbles. The study is motivated by a desire to understand the rheology of magma and lava. Our experiments are scaled to the volcanic system: they are conducted in the non-Brownian, non-inertial regime; bubble capillary number is varied across unity; and bubble and particle fractions are 0 ≤  ϕ gas  ≤ 0.82 and 0 ≤  ϕ solid  ≤ 0.37, respectively. We measure flow-front velocity and invert for a Herschel–Bulkley rheology model as a function of ϕ gas , ϕ solid , and the capillary number. We find a stronger increase in relative viscosity with increasing ϕ gas in the low to intermediate capillary number regime than predicted by existing theory, and find both shear-thinning and shear-thickening effects, depending on the capillary number. We apply our model to the existing community code for lava flow emplacement, PyFLOWGO, and predict increased viscosity and decreased velocity compared with current rheological models, suggesting existing models may not adequately account for the role of bubbles in stiffening lavas. 
    more » « less
  4. ; ; ; (, Journal of Fluid Mechanics)
    We investigate the modes of deformation of an initially spherical bubble immersed in a homogeneous and isotropic turbulent background flow. We perform direct numerical simulations of the two-phase incompressible Navier–Stokes equations, considering a low-density bubble in the high-density turbulent flow at various Weber numbers (the ratio of turbulent and surface tension forces) using the air–water density ratio. We discuss a theoretical framework for the bubble deformation in a turbulent flow using a spherical harmonic decomposition. We propose, for each mode of bubble deformation, a forcing term given by the statistics of velocity and pressure fluctuations, evaluated on a sphere of the same radius. This approach formally relates the bubble deformation and the background turbulent velocity fluctuations, in the limit of small deformations. The growth of the total surface deformation and of each individual mode is computed from the direct numerical simulations using an appropriate Voronoi decomposition of the bubble surface. We show that two successive temporal regimes occur: the first regime corresponds to deformations driven only by inertial forces, with the interface deformation growing linearly in time, in agreement with the model predictions, whereas the second regime results from a balance between inertial forces and surface tension. The transition time between the two regimes is given by the period of the first Rayleigh mode of bubble oscillation. We discuss how our approach can be used to relate the bubble lifetime to the turbulence statistics and eventually show that at high Weber numbers, bubble lifetime can be deduced from the statistics of turbulent fluctuations at the bubble scale. 
    more » « less
  5. ; ; (, Journal of Fluid Mechanics)
    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 $$r 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email