An official website of the United States government
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
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
- PAR ID:
- 10416417
- 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
-
-
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
-
Abstract Viscosity is a fundamental physical property that controls lava flow dynamics, runout distance, and velocity, which are critical factors in assessing and mitigating risks associated with effusive eruptions. Natural lava viscosity is driven by a dynamic interplay among melt, crystals, and bubbles in response to the emplacement conditions. These conditions are challenging to replicate in laboratory experiments, yet this remains the most common method for quantifying lava rheology. Few in situ viscosity measurements exist, but none of those constrains the spatial evolution of viscosity along an entire active lava flow field. Here, we present the first real-time, in situ viscosity map of active lava as measured in the field at Litli-Hrútur, Iceland. We precisely measured a lava viscosity increase of over two orders of magnitude, associated with a temperature decrease, crystallinity increase, and vesicularity decrease from near-vent to distal locations, crossing the pāhoehoe–‘a‘ā transition. Our data expand the limited database of three-phase lava viscosity, which is crucial for improvements and validation of the current numerical, experimental, and petrological approaches used to estimate lava viscosity. Further, this study showcases that field viscometry provides a rapid, accurate, and precise assessment of lava viscosity that can be implemented in eruptive response modeling of lava transport.more » « less
-
Abstract Melt composition, temperature, and crystallinity are often seen as the three most important characteristics driving lava rheology, which controls eruptive behavior. Traditional methods of measuring the viscosity of crystallizing basalts often yield different mineral characteristics to natural samples and are typically bubble-free. To quantify the viscosity of basalts inclusive of bubble and crystal cargo, we developed a new technique to measure high-temperature three-phase isothermal lava viscosity and applied it to samples from the 2018 eruption of Kīlauea. This new experimental technique begins at subliquidus temperatures, preserving original phenocrysts. A short experimental duration allows for the retention of most of the original bubble population (19%–31% vs. 36% in the original lava) and accurate replication of crystal textures from field samples, as documented in quenched postexperiment samples. The observed rheological behavior in these experiments, conducted at syneruptive temperatures (1150–1105 °C) and strain rates (0.4–18 s–1), should therefore be representative of the lava flows. We measured average viscosities of 116 Pa·s at 1150 °C to 167 Pa·s at 1115 °C (i.e., only 10%–25% higher than calculated liquid viscosities at those temperatures) and a maximum of 1800 Pa·s at 1105 °C. These results are much lower than viscosity measured in traditional bubble-free experiments, which plateaued at ~14,000 Pa·s at 1115 °C. Our results suggest the effect of bubbles in three-phase magmas may be greater than predicted by models based on two-phase bubbly liquids, and this effect must be included in realistic lava flow rheology models. The method proposed here supplies a framework for providing the necessary experimental constraints.more » « less
-
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
