skip to main content


Title: Rheology of three-phase suspensions determined via dam-break experiments
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
Award ID(s):
1929008
NSF-PAR ID:
10373933
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Volume:
477
Issue:
2254
ISSN:
1364-5021
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. We present a numerical study of non-colloidal spherical and rigid particles suspended in Newtonian, shear thinning and shear thickening fluids employing an immersed boundary method. We consider a linear Couette configuration to explore a wide range of solid volume fractions ( $0.1\leqslant \unicode[STIX]{x1D6F7}\leqslant 0.4$ ) and particle Reynolds numbers ( $0.1\leqslant Re_{p}\leqslant 10$ ). We report the distribution of solid and fluid phase velocity and solid volume fraction and show that close to the boundaries inertial effects result in a significant slip velocity between the solid and fluid phase. The local solid volume fraction profiles indicate particle layering close to the walls, which increases with the nominal $\unicode[STIX]{x1D6F7}$ . This feature is associated with the confinement effects. We calculate the probability density function of local strain rates and compare the latter’s mean value with the values estimated from the homogenisation theory of Chateau et al. ( J. Rheol. , vol. 52, 2008, pp. 489–506), indicating a reasonable agreement in the Stokesian regime. Both the mean value and standard deviation of the local strain rates increase primarily with the solid volume fraction and secondarily with the $Re_{p}$ . The wide spectrum of the local shear rate and its dependency on $\unicode[STIX]{x1D6F7}$ and $Re_{p}$ point to the deficiencies of the mean value of the local shear rates in estimating the rheology of these non-colloidal complex suspensions. Finally, we show that in the presence of inertia, the effective viscosity of these non-colloidal suspensions deviates from that of Stokesian suspensions. We discuss how inertia affects the microstructure and provide a scaling argument to give a closure for the suspension shear stress for both Newtonian and power-law suspending fluids. The stress closure is valid for moderate particle Reynolds numbers, $O(Re_{p})\sim 10$ . 
    more » « less
  4. Aqueous foams are ubiquitous; they appear in products and processes that span the cosmetics, food, and energy industries. The versatile applicability of foams comes as a result of their intrinsic viscous and elastic properties; for example, foams are exploited as drilling fluids in enhanced oil recovery for their high viscosity. Recently, so-called capillary foams were discovered: a class of foams that have excellent stability under static conditions and whose flow properties have so far remained unexplored. The unique architecture of these foams, containing oil-coated bubbles and a gelled network of oil-bridged particles, is expected to affect foam rheology. In this work, we report the first set of rheological data on capillary foams. We study the viscoelastic properties of capillary foams by conducting oscillatory and steady shear tests. We compare our results on the rheological properties of capillary foams to those reported for other aqueous foams. We find that capillary foams, which have low gas volume fractions, exhibit long lasting rheological stability as well as a yielding behavior that is reminiscent of surfactant foams with high gas volume fractions. 
    more » « less
  5. The rheology of suspensions of rings (tori) rotating in an unbounded low Reynolds number simple shear flow is calculated using numerical simulations at dilute particle number densities ( n ≪ 1 ). Suspensions of non-Brownian rings are studied by computing pair interactions that include hydrodynamic interactions modeled using slender body theory and particle collisions modeled using a short-range repulsive force. Particle contact and hydrodynamic interactions were found to have comparable influences on the steady-state Jeffery orbit distribution. The average tilt of the ring away from the flow-vorticity plane increased during pairwise interactions compared to the tilt associated with Jeffery rotation and the steady-state orbit distribution. Particle stresses associated with the increased tilt during the interaction were found to be comparable to the stresses induced directly by particle contact forces and the hydrodynamic velocity disturbances of other particles. The hydrodynamic diffusivity coefficients in the gradient and vorticity directions were also obtained and were found to be two orders of magnitude larger than the corresponding values in fiber suspensions at the same particle concentrations. Rotary Brownian dynamics simulations of isolated Brownian rings were used to understand the shear rate dependence of suspension rheology. The orbit distribution observed in the regime of weak Brownian motion, P e ≫ ϕ T − 3, was surprisingly similar to that obtained from pairwise interaction calculations of non-Brownian rings. Here, the Peclet number P e is the ratio of the shear rate and the rotary diffusivity of the particle and ϕ T is the effective inverse-aspect ratio of the particle (approximately equal to 2 π times the inverse of its non-dimensional Jeffery time period). Thus, the rheology results obtained from pairwise interactions should retain accuracy even for weakly Brownian rings ( n ≪ 1 and ϕ T − 3 ≪ P e ). 
    more » « less