skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Friday, December 13 until 2:00 AM ET on Saturday, December 14 due to maintenance. We apologize for the inconvenience.


Title: Magnetic Resonance Imaging of Multi‐Phase Lava Flow Analogs: Velocity and Rheology
Abstract

The rheology of lavas and magmas exerts a strong control on the dynamics and hazards posed by volcanic eruptions. Magmas and lavas are complex mixtures of silicate melt, suspended crystals, and gas bubbles. To improve the understanding of the dynamics and effective rheology of magmas and lavas, we performed dam‐break flow experiments using suspensions of silicone oil, sesame seeds, and N2O bubbles. Experiments were run inside a magnetic resonance imaging (MRI) scanner to provide imaging of the flow interior. We varied the volume fraction of sesame seeds between 0 and 0.48, and of bubbles between 0 and 0.21. MRI phase‐contrast velocimetry was used to measure liquid velocity. We fit an effective viscosity to the velocity data by approximating the stress using lubrication theory and the imaged shape of the free surface. In experiments with both particles and bubbles (three‐phase suspensions), we observed shear banding in which particle‐poor regions deform with a lower effective viscosity and dominate flow propagation speed. Our observations demonstrate the importance of considering variations in phase distributions within magmatic fluids and their implications on the dynamics of volcanic eruptions.

 
more » « less
Award ID(s):
1654588
PAR ID:
10416417
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
128
Issue:
5
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract

    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.

     
    more » « less
  3. Abstract

    Caldera-forming eruptions of mushy silicic magma are among the most catastrophic natural events on Earth. In such magmas, crystals form an interlocking framework when their content reaches critical thresholds, resulting in the dramatic increase in viscous resistance to flow. Here, we propose a new mechanism for the ascent of mushy magma based on microstructural observations of crystal-rich silicic pumices and lavas from the Central Andes and decompression experiments. Microstructural data include spherical vesicles and jigsaw-puzzle association of broken crystals in pumices, whereas there is limited breakage of crystals in lavas. These observations insinuate that shearing of magma during ascent was limited. Decompression experiments reveal contrasting interaction between growing gas bubbles and the crystal framework in crystal-rich magma. Under slow decompression typical of effusive eruptions, gas extraction is promoted, whereas under rapid decompression, bubbles are retained and the crystal framework collapses. This feedback between decompression rate, retention of gas bubbles, and integrity of the crystal framework leads to strong non-linearity between magma decompression rate and eruption explosivity. We extend these findings to caldera-forming eruptions of crystal-rich magma where large overpressures are induced by caldera-collapse, resulting in magma plug-flow, rapid decompression facilitated by shear-localization at conduit margins, and explosive eruption.

     
    more » « less
  4. Understanding the orientation dynamics of anisotropic colloidal particles is important for suspension rheology and particle self-assembly. However, even for the simplest case of dilute suspensions in shear flow, the orientation dynamics of non-spherical Brownian particles are poorly understood. Here we analytically calculate the time-dependent orientation distributions for non-spherical axisymmetric particles confined to rotate in the flow–gradient plane, in the limit of small but non-zero Brownian diffusivity. For continuous shear, despite the complicated dynamics arising from the particle rotations, we find a coordinate change that maps the orientation dynamics to a diffusion equation with a remarkably simple ratio of the enhanced rotary diffusivity to the zero shear diffusion: $D_{eff}^{r}/D_{0}^{r}=(3/8)(p-1/p)^{2}+1$ , where $p$ is the particle aspect ratio. For oscillatory shear, the enhanced diffusion becomes orientation dependent and drastically alters the long-time orientation distributions. We describe a general method for solving the time-dependent oscillatory shear distributions and finding the effective diffusion constant. As an illustration, we use this method to solve for the diffusion and distributions in the case of triangle-wave oscillatory shear and find that they depend strongly on the strain amplitude and particle aspect ratio. These results provide new insight into the time-dependent rheology of suspensions of anisotropic particles. For continuous shear, we find two distinct diffusive time scales in the rheology that scale separately with aspect ratio $p$ , as $1/D_{0}^{r}p^{4}$ and as $1/D_{0}^{r}p^{2}$ for $p\gg 1$ . For oscillatory shear flows, the intrinsic viscosity oscillates with the strain amplitude. Finally, we show the relevance of our results to real suspensions in which particles can rotate freely. Collectively, the interplay between shear-induced rotations and diffusion has rich structure and strong effects: for a particle with aspect ratio 10, the oscillatory shear intrinsic viscosity varies by a factor of ${\approx}2$ and the rotational diffusion by a factor of ${\approx}40$ . 
    more » « less
  5. As a paradigmatic model of active fluids, bacterial suspensions show intriguing rheological responses drastically different from their counterpart colloidal suspensions. Although the flow of bulk bacterial suspensions has been extensively studied, the rheology of bacterial suspensions under confinement has not been experimentally explored. Here, using a microfluidic viscometer, we systematically measure the rheology of dilute Escherichia coli suspensions under different degrees of confinement. Our study reveals a strong confinement effect: the viscosity of bacterial suspensions decreases substantially when the confinement scale is comparable or smaller than the run length of bacteria. Moreover, we also investigate the microscopic dynamics of bacterial suspensions including velocity profiles, bacterial density distributions, and single bacterial dynamics in shear flows. These measurements allow us to construct a simple heuristic model based on the boundary layer of upstream swimming bacteria near confining walls, which qualitatively explains our experimental observations. Our study sheds light on the influence of the boundary layer of collective bacterial motions on the flow of confined bacterial suspensions. Our results provide a benchmark for testing different rheological models of active fluids and are useful for understanding the transport of microorganisms in confined geometries. 
    more » « less