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.
Attention:The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 7:00 AM ET to 7:30 AM ET on Friday, April 24 due to maintenance. We apologize for the inconvenience.


Title: Viscosity Measurements at High Pressures: A Critical Appraisal of Corrections to Stokes' Law
Abstract Fluids and melts in planetary interiors significantly influence geodynamic processes from volcanism to global‐scale differentiation. The roles of these geofluids depend on their viscosities (η). Constraining geofluidηat relevant pressures and temperatures relies on laboratory‐based measurements and is most widely done using Stokes' Law viscometry with falling spheres. Yet small sample chambers required by high‐pressure experiments introduce significant drag on the spheres. Several correction schemes are available for Stokes' Law but there is no consensus on the best scheme(s) for high‐pressure experiments. We completed high‐pressure experiments to test the effects of (a) the relative size of the sphere diameter to the chamber diameter and (b) the top and bottom of the chamber, that is, the ends, on the sphere velocities. We examined the influence of current correction schemes on the estimated viscosity using Monte Carlo simulations. We also compared previous viscometry work on various geofluids in different experimental setups/geometries. We find the common schemes for Stokes' Law produce statistically distinct values ofη. When inertia of the sphere is negligible, the most appropriate scheme may be the Faxén correction for the chamber walls. Correction for drag due to the chamber ends depends on the precision in the sinking distance and may be ineffective with decreasing sphere size. Combining the wall and end corrections may overcorrectη. We also suggest the uncertainty inηis best captured by the correction rather than propagated errors from experimental parameters. We develop an overlying view of Stokes' Law viscometry at high pressures.  more » « less
Award ID(s):
2246803 2246802 1753125
PAR ID:
10507580
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
AGU
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Edition / Version:
1.0
Volume:
129
Issue:
5
ISSN:
2169-9313
Page Range / eLocation ID:
1-30
Subject(s) / Keyword(s):
Silicate melt viscosity high pressure Stoke's law synchrotron
Format(s):
Medium: X Size: 5.1MB Other: pdf
Size(s):
5.1MB
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The earliest form of continental crust was produced by tonalite‐trondhjemite‐granodiorite (TTG) magmas. Molten albite (NaAlSi3O8) is representative of TTGs and also a major component of modern crust‐forming magma. The viscosity of the melt controls the magma ascent rate and hence influences the production of new continental crust. It is well known that the viscosity (η) of albitic melt exhibits an anomalous pressure (P) dependence. However, prior results on the meltηat high‐Pdiffer significantly which limits our ability to predict the movement of crust‐forming magma at depth. In this study, we more tightly constrained theP‐effect onηin anhydrous albitic melt via high‐Pand high‐temperature (T) falling sphere experiments. We limited undesirable drag effects by using small sphere‐to‐capsule diameter ratios (d/D) such thatd/D ≤ 0.12, and evaluated uncertainties due to such drag using a Monte Carlo approach. Our results show that meltηfirst decreases withP(i.e., ∂η/∂P < 0) and then increases with continued compression (∂η/∂P > 0) with a well‐definedηminimum (ηmin) at ∼6 GPa along a ∼2,000 K isotherm. We find that the viscosity of the melt can be described by an Arrhenius formalism with an activation volume that varies withPandT. The results indicate thatηof aluminosilicate magmas decrease with depth and temperature in the crust, thereby mobilizing the magmas to promote rapid volcanic eruptions. The results also suggest that TTG magmas relevant for the early Earth could pond during ascent due to the anomalousP‐effect onη. 
    more » « less
  2. Abstract This paper concerns the existence of global weak solutionsá la Lerayfor compressible Navier–Stokes–Fourier systems with periodic boundary conditions and the truncated virial pressure law which is assumed to be thermodynamically unstable. More precisely, the main novelty is that the pressure law is not assumed to be monotone with respect to the density. This provides the first global weak solutions result for the compressible Navier-Stokes-Fourier system with such kind of pressure law which is strongly used as a generalization of the perfect gas law. The paper is based on a new construction of approximate solutions through an iterative scheme and fixed point procedure which could be very helpful to design efficient numerical schemes. Note that our method involves the recent paper by the authors published in Nonlinearity (2021) for the compactness of the density when the temperature is given. 
    more » « less
  3. Many viscous liquids behave effectively as incompressible under high pressures but display a pronounced dependence of viscosity on pressure. The classical incompressible Navier-Stokes model cannot account for both features, and a simple pressure-dependent modification introduces questions about the well-posedness of the resulting equations. This paper presents the first study of a second-gradient extension of the incompressible Navier-Stokes model, recently introduced by the authors, which includes higher-order spatial derivatives, pressure-sensitive viscosities, and complementary boundary conditions. Focusing on steady flow down an inclined plane, we adopt Barus' exponential law and impose weak adherence at the lower boundary and a prescribed ambient pressure at the free surface. Through numerical simulations, we examine how the flow profile varies with the angle of inclination, ambient pressure, viscosity sensitivity to pressure, and internal length scale. 
    more » « less
  4. We revisit Brenner's seminal work on the Stokes resistance of a slightly deformed sphere (Chem. Engng Sci., vol. 19, 1964, p. 519), evaluate its range of validity and extend its applicability to higher deformations for axisymmetric particles, using hydrodynamic radius as the measure of Stokes resistance. Brenner's method solves the flow around a slightly deformed sphere through two mapping steps: the first mapping translates the surface velocity on the deformed sphere to that over a reference sphere of arbitrary radius using an asymptotic expansion of the flow field in terms of deformation amplitude and a Taylor expansion of the velocity field around the surface of the reference sphere. Subsequently, the second mapping extrapolates the velocity field from the surface of the reference sphere to any point in the fluid using Lamb's general solution for Stokes flow. While the original work addresses slightly deformed spheres to a linear order in deformation amplitude, we demonstrate that the first mapping, in combination with axisymmetric spectral modes (J. Fluid Mech., vol. 936, 2022, R1), can accommodate significant deformations to arbitrary orders of perturbation, and thus is not limited to slightly deformed spheres. Also, while first-order analysis is suitable for nearly spherical particles, second-order terms can provide a reasonable range for significantly higher deformations. 
    more » « less
  5. We report experiments on flow of wormlike micellar solutions past a falling sphere. By increasing the salt-to-surfactant concentration ratio, and beyond a viscosity peak, wormlike micelles experience a transition from linear to branched microstructure. Two viscoelastic wormlike micelles with salt to surfactant concentrations on each side of the viscosity peak are considered. Our results indicate three significant differences in flows of branched and linear micelles. First, while the sphere drag correction factor rapidly decreases upon increasing Weissenberg number in linear micelles, it shows an apparent local maximum at Wi ≈ 3 in branched micelles. Second, despite its high viscoelasticity, the time-averaged flow of branched micelles around the falling sphere exhibits a fore-and-aft symmetry, while a strong negative wake is observed in linear micelles at relatively weaker flows. Third, branched micelles exhibit a stronger flow-induced birefringence than linear micelles in an otherwise identical condition. Our hypothesis is that subject to strong flows around the falling sphere, branched micelles can relax much more efficiently than linear wormlike micelles through sliding of the branched junctions. This additional stress relaxation mechanism may facilitate micellar orientation, produce a marginal sphere drag reduction and a Newtonian-like flow profile around the falling sphere. Finally, unsteady flow is observed in both linear and branched micellar solutions beyond some critical thresholds of the extensional Weissenber number. Our results corroborate a recently proposed criterion for onset of instability in flow of wormlike micelles past a falling sphere, thereby, suggesting that micellar branching does not affect the mechanism of flow instability. 
    more » « less