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: Universal non-monotonic drainage in large bare viscous bubbles
Abstract Bubbles will rest at the surface of a liquid bath until their spherical cap drains sufficiently to spontaneously rupture. For large film caps, the memory of initial conditions is believed to be erased due to a visco-gravitational flow, whose velocity increases from the top of the bubble to its base. Consequently, the film thickness has been calculated to be relatively uniform as it thins, regardless of whether the drainage is regulated by shear or elongation. Here, we demonstrate that for large bare bubbles, the film thickness is highly nonuniform throughout drainage, spanning orders of magnitude from top to base. We link the film thickness profile to a universal non-monotonic drainage flow that depends on the bubble thinning rate. These results highlight an unexpected coupling between drainage velocity and bubble thickness profiles and provide critical insight needed to understand the retraction and breakup dynamics of these bubbles upon rupture.  more » « less
Award ID(s):
1351466
PAR ID:
10442780
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Nature Communications
Volume:
14
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Journal of Fluid Mechanics)
    null (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 » « less
  2. ; ; (, Journal of Fluid Mechanics)
    Marangoni flow is the motion induced by a surface tension gradient along a fluid–fluid interface. In this study, we report a Marangoni flow generated when a bath of surfactant contacts a pre-wetted film of deionized water on a vertical substrate. The thickness profile of the pre-wetted film is set by gravitational drainage and so varies with the drainage time. The surface tension is lower in the bath due to the surfactant, and thus a liquid film climbs upwards along the vertical substrate due to the surface tension difference. Particle tracking velocimetry is performed to measure the dynamics in the film, where the mean fluid velocity reverses direction as the draining film encounters the front of the climbing film. The effect of the surfactant concentration and the pre-wetted film thickness on the film climbing is then studied. High-speed interferometry is used to measure the front position of the climbing film and the film thickness profile. As a result, higher surfactant concentration induces a faster and thicker climbing film. Also, for high surfactant concentrations, where Marangoni driving dominates, increasing the film thickness increases the rise speed of the climbing front, since viscous resistance is less important. In contrast, for low surfactant concentrations, where Marangoni driving balances gravitational drainage, increasing the film thickness decreases the rise speed of the climbing front while enhancing gravitational drainage. We rationalize these observations by utilizing a dimensionless parameter that compares the magnitudes of the Marangoni stress and gravitational drainage. A model is established to analyse the climbing front, either in the Marangoni-driving-dominated region or in the Marangoni-balanced drainage region. Our work highlights the effects of the gravitational drainage on the Marangoni flow, both by setting the thickness of a pre-wetted film and by resisting the film climbing. 
    more » « less
  3. ; ; (, Journal of Fluid Mechanics)
    Thin-film equations are utilised in many different areas of fluid dynamics when there exists a direction in which the aspect ratio can be considered small. We consider thin free films with Marangoni effects in the extensional flow regime, where velocity gradients occur predominantly along the film. In practice, because of the local deposition of surfactants or input of energy, asymmetric distributions of surfactants or surface tension more generally, are possible. Such examples include the surface of bubbles and the rupture of thin films. In this study, we consider the asymmetric thin-film equations for extensional flow with Marangoni effects. Concentrating on the case of small Reynolds number$$ Re $$, we study the deposition of insoluble surfactants on one side of a liquid sheet otherwise at rest and the resulting thinning and rupture of the sheet. The analogous problem with a uniformly thinning liquid sheet is also considered. In addition, the centreline deformation is discussed. In particular, we show analytically that if the surface tension isotherm$$\sigma = \sigma (\varGamma )$$is nonlinear (surface tension$$\sigma$$varies with surfactant concentration$$\varGamma$$), then accounting for top–bottom asymmetry leads to slower (faster) thinning and pinching if$$\sigma = \sigma (\varGamma )$$is convex (concave). The analytical progress reported in this paper allows us to discuss the production of satellite drops from rupture via Marangoni effects, which, if relevant to surface bubbles, would be an aerosol production mechanism that is distinct from jet drops and film drops. 
    more » « less
  4. ; (, 7th Thermal and Fluids Engineering Conference)
    Extensive research has been conducted to resolve small-scale microlayer and bubble nucleation and departure processes in flow boiling, building on controlled pool boiling studies. Large-scale two-phase flow structures, such as Taylor bubbles, are known to locally modify transport due to their wakes and varying surrounding liquid film thickness. However, the effect of interaction of such large-scale flow processes with bubble nucleation is not yet well characterized. Wakes may drive premature nucleating bubble departure, or conversely, suppress boiling due to boundary layer quenching, significantly affecting overall heat transfer. To explore such phenomena, a two-phase flow boiling visualization facility is developed to collect simultaneous high-speed visualization and infrared (IR) thermal imaging temperature distribution data. The test cell channel is 420 mm long with a 10 mm × 10 mm internal square-cross section. A transparent conductive indium tin oxide (ITO) coated sapphire window serves as a heater and IR interface for measuring the internal wall temperature. The facility is charged with a low boiling point fluid (HFE7000) to reduce uncertainties from heat loss to the laboratory environment. Vertical saturated flow boiling wake-nucleation interaction experiments are performed for varying liquid volume flow rates (0.5 − 1.5 L min-1, laminar-to-turbulent Re) and heat fluxes (0 − 100 kW m-2). Discrete vapor slugs are injected to explore interactions with nucleate boiling processes. By measuring film heater power, surface temperature distributions, and pressures, local instantaneous heat transfer coefficients (HTC) can be obtained. Results will be applied to assess simulations at matched conditions for void fraction, and size statistics of flow structures. 
    more » « less
  5. ; ; ; (, Journal of Fluid Mechanics)
    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
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Text Encoded Conversion
  • XML
  • Save / Share this Record
  • Send to Email