Most lava flows carry bubbles and crystals in suspension. From earlier works, it is known that spherical bubbles increase the effective viscosity while bubbles deformed by rapid flow decrease it. Changes in the spatial distribution of bubbles can lead to variable rheology and flow localization and thus modify the resulting lava flow structure and morphology. To understand the roles of bubble and solid phase crystal distributions, we conducted a series of analog experiments of high bubble fraction suspensions. We poured the analog lava on an inclined slope, observed its shape, calculated the velocity field, and monitored its local thickness. A region of localized rapid flow and low vesicularity, whose thickness is thinner than the surrounding area, develops at the center of the bubbly flows. These features suggest that the locally higher liquid fraction decreases the effective viscosity, increases the fluid density, and accelerates the flow. We also found that a halted particle‐bearing bubbly flow can resume flowing. We interpret this to result from the upward vertical separation of bubbles, which generates a liquid‐rich layer at the bottom of the flow. In our experiment, bubbles are basically spherical and decrease the flow velocity, while our estimate suggests that bubbles in natural lava flows could increase or decrease flow velocity. Downstream decreases in flow velocity stops the bubble deformation and can cause a sudden increase of effective viscosity. The vertical segregation of the liquid phase at the slowed flow front may be a way to generate a cavernous shelly paho’eho’e.
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
- NSF-PAR ID:
- 10397467
- Publisher / Repository:
- Nature Publishing Group
- 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
-
more » « less
Abstract -
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
-
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
-
more » « less
Abstract Squeezing bubbles in a tapered microgap has proved to be effective for improving flow stability in flow boiling. A previous study from our research group has successfully demonstrated using tapered microgap for transforming pool boiling into a self-sustained flow boiling-like system for cooling CPU through thermosiphon. To overcome the imaging challenges with nucleating vapor bubbles, the present work investigates the squeezing behaviour of air-injected bubbles between a tapered microgap with taper angles of 5°, 10°, and 15°. The air bubbles are injected at a rate of 3 ml/min, 15ml/min, and 30 ml/min in a pool of water. The bubble squeezing is recorded at 2000fps using a Photron high-speed camera. The experimental analysis compares the displacement and velocity of the advancing and receding bubble interfaces. The analysis found that in certain test cases, multiple bubbles coalesced while exiting the tapered microgap. In all the test cases, the receding interface of the bubble slingshots after detaching pushes the bubble out of the tapered microgap. The result from the current study provides an insight into the bubble flow and squeezing behavior that can be used for optimizing taper microgap geometries to enhance critical heat flux and heat transfer coefficient of two-phase, and air-injected single-phase heat transfer systems.
-
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
