Quartz is an abundant mineral in Earth's crust whose mechanical behavior plays a significant role in the deformation of the continental lithosphere. However, the viscoplastic rheology of quartz is difficult to measure experimentally at low temperatures without high confining pressures due to the tendency of quartz (and other geologic materials) to fracture under these conditions. Instrumented nanoindentation experiments inhibit cracking even at ambient conditions, by imposing locally high mean stress, allowing for the measurement of the viscoplastic rheology of hard materials over a wide range of temperatures. Here we measure the indentation hardness of four synthetic quartz specimens and one natural quartz specimen with varying water contents over a temperature range of 23°C to 500°C. Yield stress, which is calculated from hardness but is model dependent, is fit to a constitutive flow law for low‐temperature plasticity to estimate the athermal Peierls stress of quartz. Below 500°C, the yield stresses presented here are lower than those obtained by extrapolating a flow law constrained by experiments at higher temperatures irrespective of the applied model. Indentation hardness and yield stress depend weakly on crystallographic orientation but show no dependence on water content.
Strain localization and resulting plasticity and failure play an important role in the evolution of the lithosphere. These phenomena are commonly modeled by Stokes flows with viscoplastic rheologies. The nonlinearities of these rheologies make the numerical solution of the resulting systems challenging, and iterative methods often converge slowly or not at all. Yet accurate solutions are critical for representing the physics. Moreover, for some rheology laws, aspects of solvability are still unknown. We study a basic but representative viscoplastic rheology law. The law involves a yield stress that is independent of the dynamic pressure, referred to as von Mises yield criterion. Two commonly used variants, perfect/ideal and composite viscoplasticity, are compared. We derive both variants from energy minimization principles, and we use this perspective to argue when solutions are unique. We propose a new stress‐velocity Newton solution algorithm that treats the stress as an independent variable during the Newton linearization but requires solution only of Stokes systems that are of the usual velocity‐pressure form. To study different solution algorithms, we implement 2‐D and 3‐D finite element discretizations, and we generate Stokes problems with up to 7 orders of magnitude viscosity contrasts, in which compression or tension results in significant nonlinear localization effects. Comparing the performance of the proposed Newton method with the standard Newton method and the Picard fixed‐point method, we observe a significant reduction in the number of iterations and improved stability with respect to problem nonlinearity, mesh refinement, and the polynomial order of the discretization.
more » « less- NSF-PAR ID:
- 10453922
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 21
- Issue:
- 9
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract We derive a class of exact solutions for Stokes flow in infinite and semi‐infinite channel geometries with permeable walls. These simple, explicit, series expressions for both pressure and Stokes flow are valid for all permeability values. At the channel walls, we impose a no‐slip condition for the tangential fluid velocity and a condition based on Darcy's law for the normal fluid velocity. Fluid flow across the channel boundaries is driven by the pressure drop between the channel interior and exterior; we assume the exterior pressure to be constant. We show how the ground state is an exact solution in the infinite channel case. For the semi‐infinite channel domain, the ground‐state solutions approximate well the full exact solution in the bulk and we derive a method to improve their accuracy at the transverse wall. This study is motivated by the need to quantitatively understand the detailed fluid dynamics applicable in a variety of engineering applications including membrane‐based water purification, heat and mass transfer, and fuel cells.
-
more » « less
Abstract The dynamics of geophysical dilute turbulent gas‐particles mixtures depends to a large extent on particle concentration, which in turn depends predominantly on the particle settling velocity. We experimentally investigate air‐particle mixtures contained in a vertical pipe in which the velocity of an ascending air flux matches the settling velocity of glass particles. To obtain local particle concentrations in these mixtures, we use acoustic probing and air pressure measurements and show that these independent techniques yield similar results for a range of particle sizes and particle concentrations. Moreover, we find that in suspensions of small particles (78 μm) the settling velocity increases with the local particle concentration due to the formation of particle clusters. These clusters settle with a velocity that is four times faster than the terminal settling velocity of single particles, and they double settling speeds of the suspensions. In contrast, in suspensions of larger particles (467 μm) the settling velocity decreases with increasing particle concentration. Although particle clusters are still present in this case, the settling velocity is decreased by 30%, which is captured by a hindered settling model. These results suggest an interplay between hindered settling and cluster‐induced enhanced settling, which in our experiments occur respectively at Stokes number O(100) and O(1). We discuss implications for volcanic plumes and pyroclastic currents. Our study suggests that clustering and related enhanced or hindered particle settling velocities should be considered in models of volcanic phenomena and that drag law corrections are needed for reliable predictions and hazard assessment.
-
more » « less
Microorganisms play critical roles in sea ice biogeochemical processes. However, microbes living within sea ice can be challenging to sample for scientific study. Because most techniques for microbial analysis are optimized for liquid samples, sea ice samples are typically melted first, often applying a buffering method to mitigate osmotic lysis. Here, we tested commonly used melting procedures on three different ice horizons of springtime, first year, land-fast Arctic sea ice to investigate potential methodological impacts on resulting measurements of cell abundance, photophysiology, and microbial community structure as determined by 16S and 18S rRNA gene amplicon sequencing. Specifically, we compared two buffering methods using NaCl solutions (“seawater,” melting the ice in an equal volume of 35-ppt solution, and “isohaline,” melting with a small volume of 250-ppt solution calculated to yield meltwater at estimated in situ brine salinity) to direct ice melting (no buffer addition) on both mechanically “shaved” and “non-shaved” samples. Shaving the ice shortened the melting process, with no significant impacts on the resulting measurements. The seawater buffer was best at minimizing cell lysis for this ice type, retaining the highest number of cells and chlorophyll a concentration. Comparative measurements of bacterial (16S) community structure highlighted ecologically relevant subsets of the community that were significantly more abundant in the buffered samples. The results for eukaryotic (18S) community structure were less conclusive. Taken together, our results suggest that an equivalent-volume seawater-salinity buffered melt is best at minimizing cell loss due to osmotic stress for springtime Arctic sea ice, but that either buffer will reduce bias in community composition when compared to direct melting. Overall, these findings indicate potential methodological biases that should be considered before developing a sea ice melting protocol for microbiological studies and afterwards, when interpreting biogeochemical or ecological meaning of the results.
-
Inspired by the recent realization of a two-dimensional (2-D) chiral fluid as an active monolayer droplet moving atop a 3-D Stokesian fluid, we formulate mathematically its free-boundary dynamics. The surface droplet is described as a general 2-D linear, incompressible and isotropic fluid, having a viscous shear stress, an active chiral driving stress and a Hall stress allowed by the lack of time-reversal symmetry. The droplet interacts with itself through its driven internal mechanics and by driving flows in the underlying 3-D Stokes phase. We pose the dynamics as the solution to a singular integral–differential equation, over the droplet surface, using the mapping from surface stress to surface velocity for the 3-D Stokes equations. Specializing to the case of axisymmetric droplets, exact representations for the chiral surface flow are given in terms of solutions to a singular integral equation, solved using both analytical and numerical techniques. For a disc-shaped monolayer, we additionally employ a semi-analytical solution that hinges on an orthogonal basis of Bessel functions and allows for efficient computation of the monolayer velocity field, which ranges from a nearly solid-body rotation to a unidirectional edge current, depending on the subphase depth and the Saffman–Delbrück length. Except in the near-wall limit, these solutions have divergent surface shear stresses at droplet boundaries, a signature of systems with codimension-one domains embedded in a 3-D medium. We further investigate the effect of a Hall viscosity, which couples radial and transverse surface velocity components, on the dynamics of a closing cavity. Hall stresses are seen to drive inward radial motion, even in the absence of edge tension.more » « less
