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.


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:
10504380
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
129
Issue:
5
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract We investigate the temporal accuracy of two generalized‐ schemes for the incompressible Navier‐Stokes equations. In a widely‐adopted approach, the pressure is collocated at the time steptn + 1while the remainder of the Navier‐Stokes equations is discretized following the generalized‐ scheme. That scheme has been claimed to besecond‐order accurate in time. We developed a suite of numerical code using inf‐sup stable higher‐order non‐uniform rational B‐spline (NURBS) elements for spatial discretization. In doing so, we are able to achieve high spatial accuracy and to investigate asymptotic temporal convergence behavior. Numerical evidence suggests that onlyfirst‐order accuracyis achieved, at least for the pressure, in this aforesaid temporal discretization approach. On the other hand, evaluating the pressure at the intermediate time step recovers second‐order accuracy, and the numerical implementation is simplified. We recommend this second approach as the generalized‐ scheme of choice when integrating the incompressible Navier‐Stokes equations. 
    more » « less
  3. Abstract Explosive eruptions expel volcanic gases and particles at high pressures and velocities. Within this multiphase fluid, small ash particles affect the flow dynamics, impacting mixing, entrainment, turbulence, and aggregation. To examine the role of turbulent particle behavior, we conducted an analogue experiment using a particle‐laden jet. We used compressed air as the carrier fluid, considering turbulent conditions at Reynolds numbers from approximately 5,000 to 20,000. Two different particles were examined: 14‐μm diameter solid nickel spheres and 13‐μm diameter hollow glass spheres. These resulted in Stokes numbers between 1 and 35 based on the convective scale. The particle mass percentage in the mixture is varied from 0.3% to more than 20%. Based on a 1‐D volcanic plume model, these Stokes numbers and mass loadings corresponded to millimeter‐scale particle diameters at heights of 4–8 km above the vent during large, sustained eruptions. Through particle image velocimetry, we measured the mean flow behavior and the turbulence statistics in the near‐exit region, primarily focusing on the dispersed phase. We show that the flow behavior is dominated by the particle inertia, with high Stokes numbers reducing the entrainment by more than 40%. When applied to volcanic plumes, these results suggest that high‐density particles can greatly increase the probability of column collapse. 
    more » « less
  4. 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
  5. Abstract Slowly evolving stratified flow over rough topography is subject to substantial drag due to internal motions, but often numerical simulations are carried out at resolutions where this “wave” drag must be parameterized. Here we highlight the importance of internal drag from topography with scales that cannot radiate internal waves, but may be highly nonlinear, and we propose a simple parameterization of this drag that has a minimum of fit parameters compared to existing schemes. The parameterization smoothly transitions from a quadratic drag law () for lowNh/u0(linear wave dynamics) to a linear drag law () for highNh/u0flows (nonlinear blocking and hydraulic dynamics), whereNis the stratification,his the height of the topography, andu0is the near-bottom velocity; the parameterization does not have a dependence on Coriolis frequency. Simulations carried out in a channel with synthetic bathymetry and steady body forcing indicate that this parameterization accurately predicts drag across a broad range of forcing parameters when the effect of reduced near-bottom mixing is taken into account by reducing the effective height of the topography. The parameterization is also tested in simulations of wind-driven channel flows that generate mesoscale eddy fields, a setup where the downstream transport is sensitive to the bottom drag parameterization and its effect on the eddies. In these simulations, the parameterization replicates the effect of rough bathymetry on the eddies. If extrapolated globally, the subinertial topographic scales can account for 2.7 TW of work done on the low-frequency circulation, an important sink that is redistributed to mixing in the open ocean. 
    more » « less