Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHEAR”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

Effect of Shear Displacement and Stress Changes on Fracture Hydraulic Aperture and Flow Anisotropy

Fluid flow in fractures has the potential to drastically change the economic, environmental, and safety risks associated with a subsurface operation. This work focuses on the analysis of experimental shear fracturing and subsequent shear displacement of a shale specimen that contains multiple pre-existing natural fractures. A triaxial direct shear apparatus with concurrent µCT imaging allowed continuous monitoring of sample permeability and changes in fracture geometry during fracture creation, displacement, and changes in stress. Steady-state, low Reynolds number fluid flow simulations were used to determine the changes of each individual fracture’s hydraulic aperture due to displacement and effective stress. The development of anisotropy in the hydraulic aperture of a generated shear fracture was determined for shear displacements of 0.38, 0.74, and 2.04 mm. At high confining stress (30 MPa), shear displacement significantly reduced the hydraulic aperture of the shear fracture to below flow detectable limits and X-ray imaging resolution (23 × 23 × 23 µm 3 ). The shear generated fracture was found to be more sensitive to changes in effective stress compared to the pre-existing natural fractures. Existing models derived from fluid flow equations accounting for stochastic fracture roughness were able to partially fit the simulated flow properties of the pre-existing fractures, but no model was found to fit the flow properties of the newly formed shear fracture. These results are the first report of simultaneous changes in the geometry and flow properties of surrounding, parallel natural fractures as a shear fracture is generated with varying stress and shear displacement conditions.

58 GEOSCIENCES↗

Bending and shear improvements in 3D-printed core sandwich composites through modification of resin uptake in the skin/core interphase region

Fused deposition modeling (FDM) printed polymers are rarely used as a structural material due to anisotropic and low mechanical properties compared with conventional composites. In recent years, greater need has been expressed for recycling of materials, such as recyclable FDM, at the end of service life to reduce environmental pollution and manufacture cost. However, how the amount of resin uptake in the skin and skin/core interphase affects the bending and shear performance of the sandwich composites when replacing the low strength and ductile core (conventional core) with a high strength and brittle core (FDM printed PLA (polylactic acid) core) still remains unclear. A new manufacturing routine is needed to improve the incorporation of FDM printed polymers in composite structures. In this work, FDM printed PLA was used as core material and sandwiched between two unidirectional glass fiber reinforced polymer (GFRP) skins to form a sandwich composite by compression-molding (CM) process, which provides a good manufacturing strategy for skin/core interphase modification. The significance of the CM process is proved by investigating the effect of resin uptake on bending and in-plane/out-of-plane shear performances. Current first order shear deformation (FSDT) theory lacks a direct connection between the in-plane shear stress and out-of-shear stress in the core region of sandwich composites. With the help of DIC, a connection between the in-plane shear and the out-of-plane shear strain was built and in-plane shear properties can acquire through out-of-plane shear properties, hence reducing the redundancy of sample preparation or the need for simulation. A significant improvement was found compared with the optimized resin uptake (Optimized resin uptake range: 20.43%–22.86 wt%) 3D-printed PLA core sandwich composite and lowest performance sandwich composite (Improvement: in-plane shear strength (~34%)/modulus (~29%), out-of-plane shear strength (~25%)/modulus (~31%), specific peak bending load (~19%)). Finally, compared with balsa core sandwich composites, the 3D-printed cores are suitable for use in composite sandwich structures in many applications with a satisfactory strength-to-weight ratio.

3D-printed core sandwich composite↗

How does vertical wind shear influence entrainment in squall lines?

The influence of vertical wind shear on updraft entrainment in squall lines is not well understood. To address this knowledge gap, a suite of high-resolution idealized numerical model simulations of squall lines were run in various vertical wind shear (hereafter “shear”) environments to study the effects of shear on entrainment in deep convective updrafts. Low-level horizontal mass flux into the leading edge of the cold pool was strongest in the simulations with the strongest low-level shear. These simulations consequently displayed wider updrafts, less entrainment-driven dilution, and larger buoyancy than the simulations with comparatively weak low-level shear. An analysis of vertical accelerations along trajectories that passed through updrafts showed larger net accelerations from buoyancy in the simulations with stronger low-level shear, which demonstrates how less entrainment-driven dilution equated to stronger updrafts. Additionally, the effects of upper-level shear on entrainment and updraft vertical velocities were generally less pronounced than the effects of low-level shear. We argue that in addition to the outflow boundary-shear interactions and their effect on updraft tilt established by previous authors, decreased entrainment-driven dilution is yet another beneficial effect of strong low-level shear on squall line updraft intensity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of sheared E × B flow in self-organized, improved confinement states in magnetized plasmas

A great scientific success story of magnetic fusion research in the past several decades is the theoretical development and experimental testing of the process of turbulence decorrelation and stabilization by sheared E×B flow which shows that E×B shear effects are ubiquitous in magnetized plasmas. This concept of turbulence decorrelation and stabilization has the universality needed to explain the H-mode edge transport barriers seen in limiter and divertor tokamaks, stellarators, and mirror machines; the broader edge transport barrier seen in VHmode plasmas; and the core transport barriers formed in tokamaks. Similar effects are seen in linear devices. These examples of confinement improvement are of considerable physical interest; it is not often that a system self-organizes to reduce transport when an additional source of free energy is applied to it. The transport decrease associated with E×B velocity shear is also of great practical benefit to fusion research, since it contributed to substantially increased fusion yield in all DT magnetic fusion experiments conducted to date. The fundamental physics involved in transport reduction is the effect of E×B shear on the growth, radial extent, and phase correlation of turbulent eddies in the plasma. The same basic transport reduction process can be operational in various portions of the plasma because there are a number of ways to change the radial electric field Er. An important secondary theme in this area is the synergistic effect of E×B velocity shear and magnetic shear. Although the E×B velocity shear appears to have an effect on broader classes of microturbulence, magnetic shear can mitigate some potentially harmful effects of E×B velocity shear and facilitate turbulence stabilization. Our present understanding in this area is the result of a multi-decade, intertwined effort in theory, modeling and diagnostic development combined with continuing experimental investigations. These experiments have clearly demonstrated that increased E×B shear causes reductions in turbulence and transport. Finally, the experimental findings are generally consistent with the basic theoretical models although considerable work remains to be done before we have a fully predictive theory of transport in magnetized plasmas including E×B shear effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mesoscale simulations of pressure-shear loading of granular tungsten carbide

Numerical simulations of pressure-shear loading of a granular material are performed using the shock physics code CTH. A simple mesoscale model for the granular material is used that consists of a randomly packed arrangement of solid circular or spherical grains of uniform size separated by vacuum. The grain material is described by a simple shock equation of state, elastic perfectly plastic strength model, and fracture model with baseline parameters for WC taken from previous mesoscale modeling work. Simulations using the baseline material parameters are performed at the same initial conditions of pressure-shear experiments on dry WC powders. Except for some localized flow regions appearing in simulations with an approximate treatment of sliding interfaces among grains, the samples respond elastically during shear, which is in contrast to experimental observations. By extending the simulations to higher shear wave amplitudes, macroscopic shear failure of the simulated samples is observed with the shear strength increasing with increasing stress confinement. The shear strength is also found to be strongly dependent on the grain interface treatment and on the fracture stress of the grains, though the variation in shear strength due to fracture stress decreases with increasing stress confinement. At partial compactions, the transverse velocity histories show strain-hardening behavior followed by formation of a shear interface that extends through the transverse dimensions of the sample. Near full compaction, no strain hardening is observed and, instead, the sample transitions sharply from an elastic response to formation of an internal shear interface. Agreement with experiment is shown to worsen with increasing confinement stress with simulations overpredicting the shear strengths measured in experiment. The source of the disagreement can be ultimately attributed to the Eulerian nature of the simulations, which do not treat contact and fracture realistically.

36 MATERIALS SCIENCE↗

Mitigating Shear-dependent Object Detection Biases with Metacalibration

M ETACALIBRATION is a new technique for measuring weak gravitational lensing shear that is unbiased for isolated galaxy images. Here we test M ETACALIBRATION with overlapping, or "blended" galaxy images. Using standard M ETACALIBRATION , we find a few percent shear measurement bias for galaxy densities relevant for current surveys, and that this bias increases with increasing galaxy number density. We show that this bias is not due to blending itself, but rather to shear-dependent object detection. If object detection is shear-independent, no deblending of images is needed, in principle. We demonstrate that detection biases are accurately removed when including object detection in the M ETACALIBRATION process, a technique we call M ETADETECTION . This process involves applying an artificial shear to images of small regions of sky and performing detection on the sheared images, as well as measurements that are used to calculate a shear response. We demonstrate that the method can accurately recover weak shear signals even in highly blended scenes. In the M ETACALIBRATION process, the space between objects is sheared coherently, which does not perfectly match the real universe in which some, but not all, galaxy images are sheared coherently. We find that even for the worst case scenario, in which the space between objects is completely unsheared, the resulting shear bias is at most a few tenths of a percent for future surveys. We discuss additional technical challenges that must be met in order to implement M ETADETECTION for real surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Influence of alkali-silica reaction on the shear capacity of reinforced concrete beams with minimum transverse reinforcement

A growing number of critical concrete infrastructure are affected by alkali-silica reaction (ASR) damage such as the Seabrook Nuclear Power Plant in New Hampshire, Parker arch-gravity dam in Arizona, and several highway bridges in California and Texas. ASR causes expansion and cracking and degrades the concrete mechanical properties. Despite a wealth of material level studies there is still limited large-scale experimental data regarding the effects of ASR on reinforced concrete (RC) members. Due to the brittle nature of the shear failure in RC structures, this study focuses on the shear response of full-scale ASR damaged RC beams with minimum shear reinforcement. Six RC beams were built with different levels of ASR susceptibility and conditioned in different environments during which continuous expansion monitoring was performed. The beams all contained reactive fine aggregate (sand) and two of them had additional alkali, 1.25% by weight of cement, to accelerate ASR. The highest expansion rate of the beams happened during the first 150 days and the expansions stayed constant after 240 days and 330 days for beams conditioned in outdoor and laboratory conditions, respectively, until the last measurement at 575 days. The highest expansion, 0.4%, was seen in the beams with additional alkali and conditioned outside with regular water spray. Out of the six beams, three were selected at different levels of ASR damage and two shear tests were performed on the minimally reinforced spans close to the ends of each beam. Results indicated that beams gain shear strength from ongoing cement hydration in the presence of moisture for ASR expansions less than 0.2%. Compared to one of the beams with 0.2% ASR expansion, one of the other samples with 0.4% expansion lost 6% of its shear strength, 25% of its shear stiffness, and showed about two times larger shear cracks and shear deformations at peak load. Finally, the shear reinforcement yielded at 20% less load in the beam with 0.4% expansion compared to the beam with 0.2% expansion.

42 ENGINEERING↗

Microstructural evolution in shear-punch tests: A comparative study of pure Cu and Cu-Cr alloy

Understanding the mechanisms behind microstructural evolution during shear deformation has been a long-standing area of interest. However, establishing a connection between microstructure, mechanical properties, and extent of shear deformation is challenging and requires refined experimental approaches. Shear-punch testing (SPT) provides a controlled method to introduce shear into small volumes of material that later can be subjected to detailed microstructural characterization. In this study, we utilize an SPT device to induce shear deformation to pure copper and a binary copper-chromium alloy. Electron backscatter diffraction and transmission electron microscopy were used to study the mechanisms of plastic deformation after SPT. Our results indicate that shear deformation of pure Cu produces a dense network of intercepting microshear bands upon sustained deformation. Twin boundaries undergo degradation into high angle grain boundaries due to simultaneous deviation from the axis-angle pair condition of 60° misorientation on [111] direction. The presence of 50% volume Cr particles in the soft Cu matrix fundamentally altered the shear deformation mechanism. Preferential deformation of the Cu matrix led to accelerated shear-induced formation of low and high angle grain boundaries, and subsequent grain refinement. Comparatively, no grain refinement occurred in pure Cu at a strain ~10 times larger (ϵ = 4.73) than that of the copper-chromium case (ϵ = 0.42). Overall, our study sheds light on the microstructural evolution during shear deformation and highlights the significant influence of particle reinforcements on the shear deformation mechanisms of metals.

36 MATERIALS SCIENCE↗

Influence of Shear Strength Assumptions on BISON Debonding Simulations

Accurately predicting the thermomechanical response of buffer–IPyC debonding in TRISO fuel particles requires reliable mechanical property inputs for each coating layer, particularly the normal and shear strengths that influence interlayer delamination and stress concentrations. Micro tensile testing of AGR-2 fuel particles provided experimentally measured normal strengths for the buffer, IPyC, and buffer–IPyC interface; however, shear strength was not measured. As a result, BISON simulations of interface debonding must rely on assumed shear strength values, typically estimated as 20–40% of the measured ultimate tensile strength. This study evaluates how these assumed shear strength values influence cohesive zone model (CZM) predictions of buffer–IPyC separation in AGR 2 TRISO particles. Using micro tensile data from three AGR 2 compacts (2 1 3, 5 1 3, and 6 3 3), BISON simulations were performed with multiple shear strength assumptions to quantify their effect on radial and tangential stress evolution, debonding, and gap propagation. The results show that shear strength is a high sensitivity parameter: increasing the assumed shear strength significantly alters the stress distribution at the buffer–IPyC junction, shifts the predicted debonding location, and changes the extent of partial gap formation. While normal strength controls the initiation of interface separation, shear strength strongly influences the mode mixity of the failure process and the resulting stress concentrations transmitted to the IPyC and SiC layers. These findings highlight a critical gap in current TRISO mechanical characterization. Without experimentally measured shear strength, BISON simulations must rely on approximations that introduce uncertainty into predictions of coating layer integrity and fission product barrier performance. Future fuel qualification campaigns should therefore consider measurement of shear strength at the interlayer interfaces to reduce model uncertainty and improve the fidelity of TRISO fuel performance simulations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A High-Speed Rotational Diamond Anvil Cell for In Situ Analysis of Hierarchical Microstructural Evolution of Metallic Alloys during Extreme Shear Deformation

High speed shear deformation is ubiquitous in engineering applications, ranging from material processing methods such as friction stir processing/extrusion and in tribological contacts. However, analyzing the microstructural evolution of materials while they are undergoing high speed shear deformation have been a long-standing challenge. This led to predominant reliance on ex situ microscopy before and after shear deformation. But ex situ microscopy lacks the ability to analyze dynamic and transient hierarchical microstructural evolution mechanisms that could occur during shear deformation of materials. Therefore, to better understand the dynamic mechanisms of mass and energy transfer in materials under shear deformation, we developed a first of its kind high-speed rotational diamond anvil cell (HS-RDAC) for synchrotron-based in situ high-energy x-ray diffraction (XRD). We studied the time resolved lattice strain evolution, XRD peak broadening and changes in spatial variation of shear deformation induced alloying in pure metal and metal alloy sheets and powder mixture using the HS-RDAC. These in situ results were combined with detailed ex situ microstructural characterization before and after the shear deformation using transmission electron microscopy and atom probe tomography, which revealed the different stages of evolution of a shear deformation induced hierarchical nanostructure. Multiscale computational simulations including computational fluid dynamics, crystal plasticity, molecular dynamic simulation and density functional theory uncovered the mechanisms behind morphological changes, evolution of defect structures and changes in driving force for shear deformation induced intermixing. In conclusion, this in situ HS-RDAC capability, in combination with ex situ microstructural characterization and computational simulations, can provide new insights into the hierarchical microstructural evolution pathway during shear deformation.

36 MATERIALS SCIENCE↗

Shallow Cumulus Entrainment Dynamics in a Sheared Environment

Abstract Vertical wind shear has long been known to tilt convective towers and reduce thermal ascent rates. The purpose of this study is to better understand the physical mechanisms responsible for reduced ascent rates in shallow convection. In particular, the study focuses on cloud-edge mass flux to assess how shear impacts mass-flux profiles of both the ensemble and individual clouds of various depths. A compositing algorithm is used to distill large-eddy simulation (LES) output to focus on up- and down-shear cloud edges that are not influenced by complex cloud geometry or nearby clouds. A direct entrainment algorithm is used to estimate the mass flux through the cloud surface. We find that the dynamics on the up- and down-shear sides are fundamentally different, with the entrainment of environmental momentum and dilution of buoyancy being primarily responsible for the reduced down-shear ascent rates. Direct estimates of fluid flow through the cloud interface indicate a counter-shear organized flow pattern that entrains on the down-shear side and detrains on the up-shear side, resulting from the subcloud shear being lifted into the cloud layer by the updraft. In spite of organized regions of entrainment and detrainment, the overall net lateral mass flux remains unchanged with respect to the no shear run, with weak detrainment present throughout cloud depth.

Meteorology & Atmospheric Sciences↗

A New Simplified Discrete Fracture Model for Shearing of Intersecting Fractures and Faults

Abstract Shearing of fractures and faults is important because it can result in permeability change or even induce seismicity—both are keys for efficient and safe energy recovery and storage in Earth systems. Quantitative analysis of shearing of intersecting fractures and faults is challenging because it can involve dynamic frictional contacts that are complicated by deformation of the rock matrix. To predict the shearing of intersecting fractures/faults, we attempt to answer the question of how intersections impact the shearing of a fracture network and whether we can simplify the description as compared to classical discrete fracture network (DFN) models. To answer these questions, we conducted a series of numerical simulations on scenarios for variable numbers of intersecting fractures. All these examples yield consistent results: the results of using DFNs are consistent with those of using hypothetical major paths. This leads to a new model, which we name simplified discrete fracture network model , to analyze shearing of intersecting fractures/faults using major path(s). We found that the intersections of fractures do not fundamentally change the shearing of two intersecting fractures if the intersecting angles are small. Furthermore, increasing the number of fractures/faults may relax the stress as more fractures/faults become available for shearing and distributing the stress. The simplified DFN model, which can capture efficiently the shearing behavior of each major paths from a large number of intersecting fractures/faults, will be a promising conceptual model that is complementary to existing equivalent continuum and discrete fracture models to analyze shearing of intersecting fractures/faults.

58 GEOSCIENCES↗

An atomistic survey of shear coupling in asymmetric tilt grain boundaries and interpretation using the disconnections framework

Grain Boundaries (GB) play an important role in determining the behavior of polycrystalline materials. While the mechanisms of motion and associated shear response for symmetric tilt grain boundaries (STGBs) are well studied, the same is not true for asymmetric tilt grain boundaries (ATGBs) despite their greater prevalence in polycrystals. Here, this study aims to investigate the shear response of a large collection of asymmetric tilt grain boundaries (ATGBs) using molecular dynamics (MD) simulations and interpret the data using a discrete disconnections model that works remarkably well for STGBs. MD simulations of shear-driven ATGBs show that the plastic shear (shear coupling factor) in the region swept by a GB exhibits a complex dependence on the inclination angle, and this dependence changes with the misorientation of the GB. In addition, the shear response was observed to be highly sensitive to the applied shear rate and temperature. Recognizing the spatial and temporal scale limitations of MD simulations, we extended the discrete disconnections mesoscale model of Khateret al. (2012) to calculate the nucleation barriers of disconnection modes and predict the effective shear coupling of an ATGB. We observed that the mesoscale model’s predictions of the shear coupling factor of ATGBs do not agree with those observed in MD simulations. Finally, we examine the hypotheses of our mesoscale model that contribute to disagreements between MD simulations and the mesoscale model and propose improvements to the mesoscale model for future work.

36 MATERIALS SCIENCE↗

Insight into ideal shear strength of Ni-based dilute alloys using first-principles calculations and correlational analysis

Here the present work examines the effect of alloying elements (denoted X) on the ideal shear strength for 26 dilute Ni-based alloys, Ni11X, as determined by first-principles calculations of pure alias shear deformations. The variations in ideal shear strength are quantitatively explored with correlational analysis techniques, showing the importance of atomic properties such as size and electronegativity. The shear moduli of the alloys are affirmed to show a strong linear relationship with their ideal shear strengths, while the shear moduli of the individual alloying elements were not indicative of alloy shear strength. Through combination with available ideal shear strength data on Mg alloys, a potential application of the Ni alloy data is demonstrated in the search for a set of atomic features suitable for machine learning applications to mechanical properties. As another illustration, the calculated Ni ideal shear strengths play a key role in a predictive multiscale framework for deformation behavior of single crystal alloys at large strains, as shown by simulated stress–strain curves.

36 MATERIALS SCIENCE↗

Shear Response of Ionizable Polymer Melts at the Crossover from Ionomers to Polyelectrolytes

Nonlinear shear response of polymers is affected by inherent barriers for diffusive motion, including entanglements and topology. In ionizable polymers, ionic clusters further constrain the intrinsic dynamics of the polymers, significantly enhancing their viscosity. Here, using fully atomistic molecular dynamics simulations, the nonlinear shear response of ionizable polymers is presented, across the transition from the ionomer regime where distinctive clusters dominate the structure to the polyelectrolyte regime where clusters percolate, in polystyrene randomly sulfonated with fractions of SO 3 − groups of f = 0.20 and 0.35, in pristine and tetrahydrofuran (THF) swollen polyelectrolyte melts. For f = 0.20, the ionic clusters first fracture into smaller clusters followed by splitting into individual ionic groups and eventually reform. At higher f, the clusters morph in shape but do not break under high shear. At very high shear rates, all of the chains stretch and recoil rapidly. As the shear rate is reduced, some chains stretch and recoil, while others remain largely unaffected by the shear. Macroscopically, for all systems, the shear viscosity displays initially an elastic response, followed by nonlinear shear stress overshoot and, eventually, a steady state. The evolution of viscosity with time and shear reflects that of the ionizable domains.

cluster chemistry↗

Hydro-Mechanical Measurements of Sheared Crystalline Rock Fractures With Applications for EGS Collab Experiments 1 and 2

We present hydro-mechanical measurements that characterize shear on natural fractures in schist, amphibolite, and rhyolite specimens from the enhanced geothermal system (EGS) Collab Project's Experiment 1 and 2 sites (E1 and E2) at the Sanford Underground Research Facility. We employed a triaxial direct shear method augmented with X-ray imaging to perform hydroshearing (injection-induced shearing) and mechanical shearing on naturally fractured specimens at in situ stress conditions. Measurements included fracture permeability, strength, stress-dependent aperture, shear dilation, and frictional strength. Results reveal that in situ natural fractures must be permeable, weak, and shear-oriented to be hydrosheared, and only a subset of the observable in situ fractures were suitable. When sheared, the fracture permeability typically increased by a factor of 10 or more and this increase was retained over time. However, shear slip did not always result in permeability increase. High phyllosilicate content associated with exceptionally weak fractures exhibited poor or even decreased permeability after stimulation. These measurements in combination with site data were used to conduct a slip-tendency analysis for different fracture sets, and we selected the top candidate natural fractures for hydroshearing at the EGS Collab sites. We also found that the lower in situ shear stress and stronger fractures at the E2 site make hydroshearing more challenging than at the E1 site. Overall, the methods and analysis used in our work can be applied to any geothermal project to identify in situ joint sets that are best suited for hydroshearing, which in turn can help to optimize well placement and energy production.

58 GEOSCIENCES↗

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Medium-density amorphous ice unveils shear rate as a new dimension in water’s phase diagram

Recent experiments revealed a new amorphous ice phase, medium-density amorphous ice (MDA), formed by ball-milling ice I h at 77 K [Rosu-Finsen et al., Science 379, 474–478 (2023)]. MDA has density between that of low-density amorphous (LDA) and high-density amorphous (HDA) ices, adding to the complexity of water’s phase diagram, known for its glass polyamorphism and two-state thermodynamics. The nature of MDA and its relation to other amorphous ices and liquid water remain unsolved. Here, we use molecular simulations under controlled pressure and shear rate at 77 K to produce and investigate MDA. Here. we find that MDA formed at constant shear rate is a steady-state nonequilibrium shear-driven amorphous ice (SDA), that can be produced by shearing ice I h , LDA, or HDA. Our results suggest that MDA could be obtained by ball-milling water glasses without crystallization interference. Increasing the shear rate at ambient pressure produces SDAs with densities ranging from LDA to HDA, revealing shear rate as a new thermodynamic variable in the nonequilibrium phase diagram of water. Indeed, shearing provides access to amorphous states inaccessible by controlling pressure and temperature alone. SDAs produced with shearing rates as high as 10 6 s −1 sample the same region of the potential energy landscape than hyperquenched glasses with identical density, pressure, and temperature. Intriguingly, SDAs obtained by shearing at ~10 8 s −1 have density, enthalpy, and structure indistinguishable from those of water “instantaneously” quenched from room temperature to 77 K over 10 ps, making them good approximants for the “true glass” of ambient liquid water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗