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.


Title: Bubble pinch-off in turbulence
Although bubble pinch-off is an archetype of a dynamical system evolving toward a singularity, it has always been described in idealized theoretical and experimental conditions. Here, we consider bubble pinch-off in a turbulent flow representative of natural conditions in the presence of strong and random perturbations, combining laboratory experiments, numerical simulations, and theoretical modeling. We show that the turbulence sets the initial conditions for pinch-off, namely the initial bubble shape and flow field, but after the pinch-off starts, the turbulent time at the neck scale becomes much slower than the pinching dynamics: The turbulence freezes. We show that the average neck size, d ¯ , can be described by d ¯ ∼ ( t − t 0 ) α , where t 0 is the pinch-off or singularity time and α ≈ 0.5 , in close agreement with the axisymmetric theory with no initial flow. While frozen, the turbulence can influence the pinch-off through the initial conditions. Neck shape oscillations described by a quasi–2-dimensional (quasi-2D) linear perturbation model are observed as are persistent eccentricities of the neck, which are related to the complex flow field induced by the deformed bubble shape. When turbulent stresses are less able to be counteracted by surface tension, a 3-dimensional (3D) kink-like structure develops in the neck, causing d ¯ to escape its self-similar decrease. We identify the geometric controlling parameter that governs the appearance of these kink-like interfacial structures, which drive the collapse out of the self-similar route, governing both the likelihood of escaping the self-similar process and the time and length scale at which it occurs.  more » « less
Award ID(s):
1844932
PAR ID:
10145520
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
116
Issue:
51
ISSN:
0027-8424
Page Range / eLocation ID:
25412 to 25417
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract In this work, we experimentally measured the pinch‐off of a gas bubble on a biphilic surface, which consisted of an inner circular superhydrophobic region and an outer hydrophilic region. The superhydrophobic region had a radius ofRSHvarying from 2.8 to 19.0 mm, where the largeRSHmodeled an infinitely large superhydrophobic surface. We found that during the pinch‐off, the contact line had two different behaviors: for smallRSH, the contact line was fixed at the boundary of superhydrophobic and hydrophilic regions, and the contact angle gradually increased; in contrast, for largeRSH, the contact angle was fixed, and the contact line shrank toward the bubble center. Furthermore, we found that regardless of bubble size and contact line behavior, the minimum neck radius collapsed onto a single curve after proper normalizations and followed a power–law relation where the exponent was close to that for bubble pinch‐off from a nozzle. The local surface shapes near the neck were self‐similar. Our results suggest that the surface wettability has a negligible impact on the dynamics of pinch‐off, which is primarily driven by liquid inertia. Our findings improve the fundamental understanding of bubble pinch‐off on complex surfaces. 
    more » « less
  2. In this work, we experimentally studied bubble formation on the superhydrophobic surface (SHS) under a constant gas flow rate and at quasi-static regime. SHS with a radius RSHS ranging from 4.2 to 19.0 mm was used. We observed two bubbling modes A and B, depending on RSHS. In mode A for small RSHS, contact line fixed at the rim of SHS, and contact angle (θ) initially reduced, then maintained as a constant, and finally increased. In mode B for large RSHS, contact line continuously expanded, and θ slowly reduced. For both modes, during necking, contact line retracts, and θ was close to the equilibrium contact angle. Moreover, the pinch-off of bubble at the early stage was similar to the pinch-off of bubble from a nozzle and followed a power-law relation Rneck ∼ τ0.54, where Rneck is the minimum neck radius and τ is the time to detaching. Furthermore, we calculated the forces acting on the bubble and found a balance between one lifting force (pressure force) and two retaining forces (surface tension force and buoyancy force). Last, we found a waiting time for a finite volume to be detected for large RSHS. The detached volume was well predicted by Tate volume, which was derived based on balance between buoyancy and surface tension and was a function of bubble base radius. 
    more » « less
  3. 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
  4. Bubble-mediated gas exchange in turbulent flow is critical in bubble column chemical reactors as well as for ocean–atmosphere gas exchange related to air entrained by breaking waves. Understanding the transfer rate from a single bubble in turbulence at large Péclet numbers (defined as the ratio between the rate of advection and diffusion of gas) is important as it can be used for improving models on a larger scale. We characterize the mass transfer of dilute gases from a single bubble in a homogeneous isotropic turbulent flow in the limit of negligible bubble volume variations. We show that the mass transfer occurs within a thin diffusive boundary layer at the bubble–liquid interface, whose thickness decreases with an increase in turbulent Péclet number, $$\widetilde {{Pe}}$$ . We propose a suitable time scale $$\theta$$ for Higbie ( Trans. AIChE , vol. 31, 1935, pp. 365–389) penetration theory, $$\theta = d_0/\tilde {u}$$ , based on $$d_0$$ the bubble diameter and $$\tilde {u}$$ a characteristic turbulent velocity, here $$\tilde {u}=\sqrt {3}\,u_{{rms}}$$ , where $$u_{{rms}}$$ is the large-scale turbulence fluctuations. This leads to a non-dimensional transfer rate $${Sh} = 2(3)^{1/4}\sqrt {\widetilde {{Pe}}/{\rm \pi} }$$ from the bubble in the isotropic turbulent flow. The theoretical prediction is verified by direct numerical simulations of mass transfer of dilute gas from a bubble in homogeneous and isotropic turbulence, and very good agreement is observed as long as the thin boundary layer is properly resolved. 
    more » « less
  5. Brehm, Christop; Pandya, Shishir (Ed.)
    Paraffin wax is a prominent solid fuel for hybrid rockets. The atomization process of the paraffin wax fuel into he hybrid rocket combustion involves the droplets pinching off from the fuel surface. Therefore, droplet formation and pinch-off dynam- ics is analyzed using a one-dimensional axisymmetric approximation to understand droplet size distribution and pinch-off time. A mixed finite element formulation is used to solve the numerical problem. The computational algorithm uses adaptive mesh refinement to capture singularity and runs self-consistently to calculate droplet elongation. The code is verified using the Method of Manufactured Solution (MMS) and validated against laboratory experiments. Moreover, paraffin wax simulations are explored for varying inlet radius and it is found that the droplet size increases very slightly with the increasing inlet radius. Also, the pinch-off time increases up to a point where it starts to decrease as we increase the inlet radius. This behavior leads to a conjecture for the theoretical maximum radius that the droplet approaches as the inlet radius increases, which is a motivation for the future work. 
    more » « less