SEARCH · Search NASA

Results for “anisotropic stresses”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Anisotropy evolution of elastomeric foams during uniaxial compression measured via in-situ X-ray computed tomography

We have characterized the three-dimensional evolution of microstructural anisotropy of a family of elastomeric foams during uniaxial compression via in-situ X-ray computed tomography. Flexible polyurethane foam specimens with densities of 136, 160 and 240 kg/m 3 were compressed in uniaxial stress tests both parallel and perpendicular to the foam rise direction, to engineering strains exceeding 70%. The uncompressed microstructures show slightly elongated ellipsoidal pores, with elongation aligned parallel to the foam rise direction. The evolution of this microstructural anisotropy during deformation is quantified based on the autocorrelation of the image intensity, and verified via the mean intercept length as well as the shape of individual pores. Trends are consistent across all three methods. In the rise direction, the material remains transversely anisotropic throughout compression. Anisotropy initially decreases with compression, reaches a minimum, then increases up to large strains, followed by a small decrease in anisotropy at the largest strains as pores collapse. Compression perpendicular to the foam rise direction induces secondary anisotropy with respect to the compression axis, in addition to primary anisotropy associated with the foam rise direction. In contrast to compression in the rise direction, primary anisotropy initially increases with compression, and shows a slight decrease at large strains. These surprising non-monotonic trends and qualitative differences in rise and transverse loading are explained based on the compression of initially ellipsoidal pores. Microstructural anisotropy trends reflect macroscopic stress-strain and lateral strain response. These findings provide novel quantitative connections between three-dimensional microstructure and anisotropy in moderate density polymer foams up to large deformation, with important implications for understanding complex three-dimensional states of deformation.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Anisotropic structure and optical engineering of strontium titanate and zirconia responding to sequential hydrogen and helium irradiation

Unraveling the correlation between strain engineering with anisotropic optical properties via nanoscale defects gradually evolved into a strategy for fundamental studies and technological applications, yet it remains understudied in functional complex oxides compared to metals and semiconductors. Here, the methodology of strain engineering for the determination of the lattice parameters, parallel/normal to the sample surface, in the individual layers of single-crystalline superlattices is derived, which is based on analysis of high-angle X-ray diffraction measurements in combination with diffraction reciprocal space mapping. With modeling that takes into account the effect of elastic properties, two elastically anisotropic materials, SrTiO 3 and ZrO 2 , have been compared in terms of defect-induced elastic strain caused by individual and sequential H + and He 2+ irradiation. The anisotropic lattice swelling with corresponding refractive index, and obstructive behavior of elastic strain recovery are demonstrated in SrTiO 3 , while approaching strain behaviors accompanied by isotropic refractive index distribution in both orientations are confirmed in ZrO 2 . Under sequential He 2+ /H + irradiation, pre-existing He-vacancy complexes acted as vacancy traps, preferentially capturing H + clusters to induce lattice distortion and enhance absorption. Nanohardness increments (ΔH) calculated by the DBH model matched nanoindentation results, confirming that sequential irradiation generated higher indentation yield stress than individual irradiations.

Anisotropic expansion↗

DOE OSTI2022

An Eulerian crystal plasticity framework for modeling large anisotropic deformations in energetic materials under shocks

Here, this paper demonstrates a novel Eulerian computational framework for modeling anisotropic elastoplastic deformations of organic crystalline energetic materials (EM) under shocks. While Eulerian formulations are advantageous for handling large deformations, constitutive laws in such formulations have been limited to isotropic elastoplastic models, which may not fully capture the shock response of crystalline EM. The present Eulerian framework for high-strain rates, large deformation material dynamics of EM incorporates anisotropic isochoric elasticity via a hypo-elastic constitutive law and visco-plastic single-crystal models. The calculations are validated against atomistic calculations and experimental data and benchmarked against Lagrangian (finite element) crystal plasticity computations for shock-propagation in a monoclinic organic crystal, octahydro-1,3,5,7-tetranitro-1,3,5,7 tetrazocine (β-HMX). The Cauchy stress components and the resolved shear stresses calculated using the present Eulerian approach are shown to be in good agreement with the Lagrangian computations for different crystal orientations. The Eulerian framework is then used for computations of shock-induced inert void collapse in β-HMX to study the effects of crystal orientations on hotspot formation under different loading intensities. The computations show that the hotspot temperature distributions and the collapse profiles are sensitive to the crystal orientations at lower impact velocities (viz., 500 m/s); when the impact velocity is increased to 1000 m/s, the collapse is predominantly hydrodynamic and the role of anisotropy is modest. The present methodology will be useful to simulate energy localization in shocked porous energetic material microstructures and other situations where large deformations of single and polycrystals govern the thermomechanical response.

42 ENGINEERING↗

DOE OSTI2024

A crystal-plasticity-informed Gaussian Process Regression model to capture anisotropy in single crystal shape memory alloys

This work presents a machine learning (ML) framework that model the anisotropic actuation responses in a shape memory alloy. A Gaussian Process Regression (GPR) based ML model is trained on a set of different crystal orientations subjected to different actuation conditions. The training employed thermo-mechanical responses from a crystal-plasticity model that captures phase-transformation, stress-induced plasticity, and transformation-induced plasticity. Further, on training the GPR-ML model at fixed stress level for different orientations, it captured the thermo-mechanical responses accounting for the anisotropy, and predicted responses for new orientations with good accuracy. The GPR-ML model is able to capture the transformation temperature variations even when trained using multiple stress levels, and the transformation strain showed significant deviations. The developed GPR-ML model gave reasonable predictions for an unexplored sample set of orientations and loading conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Synchrotron self-Compton radiation from magnetically dominated turbulent plasmas in relativistic jets

ABSTRACT Relativistic jets launched by rotating black holes are powerful emitters of non-thermal radiation. Extraction of the rotational energy via electromagnetic stresses produces magnetically dominated jets, which may become turbulent. Studies of magnetically dominated plasma turbulence from first principles show that most of the accelerated particles have small pitch angles, i.e. the particle velocity is nearly aligned with the local magnetic field. We examine synchrotron self-Compton radiation from anisotropic particles in the fast cooling regime. The small pitch angles reduce the synchrotron cooling rate and promote the role of inverse Compton (IC) cooling, which can occur in two different regimes. In the Thomson regime, both synchrotron and IC components have soft spectra, νFν ∝ ν1/2. In the Klein–Nishina regime, synchrotron radiation has a hard spectrum, typically νFν ∝ ν, over a broad range of frequencies. Our results have implications for the modelling of BL Lacertae objects (BL Lacs) and gamma-ray bursts (GRBs). BL Lacs produce soft synchrotron and IC spectra, as expected when Klein–Nishina effects are minor. The observed synchrotron and IC luminosities are typically comparable, which indicates a moderate anisotropy with pitch angles θ ≳ 0.1. Rare orphan gamma-ray flares may be produced when θ ≪ 0.1. The hard spectra of GRBs may be consistent with synchrotron radiation when the emitting particles are IC cooling in the Klein–Nishina regime, as expected for pitch angles θ ∼ 0.1. Blazar and GRB spectra can be explained by turbulent jets with a similar electron plasma magnetization parameter, σe ∼ 104, which for electron–proton plasmas corresponds to an overall magnetization σ = (me/mp)σe ∼ 10.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Resonant ultrasound spectroscopy for crystalline samples containing initial strain

Resonant ultrasound spectroscopy (RUS) is a mature and robust technique for the nondestructive characterization of the elastic properties of solids capable of providing the elastic constants of anisotropic crystalline solids. The traditional method is based on assuming that the solid is linear elastic and, therefore, obeys a linearized Hookean constitutive relationship (Hooke’s law). In this work, an alternative constitutive law is provided that allows for an initial strain or prestrain to be present stemming from residual stress. Then, the constitutive relationship is integrated into the RUS framework. The model is demonstrated using a realistic prestrain field obtained by simulating shot-peening processing of a polycrystalline Cu sample. The sensitivity of the resonances to the developed prestrain is established and discussed. This work allows researchers to consider the influence of initial strain or residual stress in their samples and the potential influence on accurate elastic constant estimates. The model also supports the potential of RUS for the nondestructive characterization of prestrain in materials.

42 ENGINEERING↗

DOE OSTI2026

Oil-impregnated densified wood veneer with high electrical insulation enabled by nanosized oil channels

Growing energy demands and renewable integration are stressing the aging power grid infrastructure. Lignocellulosic oil-impregnated paper is widely used in power transformers but suffers from critical limitations, such as low dielectric strength, mechanical strength, and thermal conductivity, causing premature transformer failures. Here, we demonstrate a superior electrically insulating oil-impregnated paper design using the naturally anisotropic structure of densified wood veneer to achieve nanosized channels of oil that efficiently disrupt electrical breakdown pathways. The developed oil-impregnated densified wood (ODW) creates aligned cellulose fibers with 166 ± 87–nanometer oil nanochannels, achieving record dielectric strength of 105 kilovolts per millimeter. The structure also delivers a mechanical strength of up to 384 megapascals and a thermal conductivity of 0.33 watts per meter per kelvin, enabling enhanced longevity upon thermal aging tests. The ODW could replace conventional transformer insulation to enhance power transformer performance and improve lifetime. Moreover, its anisotropic oil-filled nanochannel design offers a general strategy for hybrid dielectrics in medium- and high-voltage applications.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Quantifying precursors to void nucleation and coalescence in aluminum

Ductile rupture is a common failure mode for engineering alloys. It is generally comprised of three mechanisms: void nucleation at secondary particles, void growth, and coalescence of voids. Conventional models for ductile rupture have limited precision for predicting macroscopic failure. This is partially because they have been corroborated using ex-situ observations of these three mechanisms and calibrated by macroscopic strain or stress metrics. In this study, in-situ high-energy X-ray characterization was conducted to correlate sites of void nucleation via particle cracking and sites of void growth with grain-scale metrics. Here it is shown that particle cracking is not predicated by elevated stress metrics, which disagrees with classical models. Instead, particle cracking tended to occur for the largest, least spherical particles. This trend persisted on the aggregate-scale and within individual neighborhoods around particles. Furthermore, an Eshelby analysis illuminated a statistically significant increase in maximum principal stress within a cracked particle, compared to an uncracked particle. In-situ observations also revealed two separate examples of intragranular void growth comprised of flat, crack-like features that grew showing alignment with slip systems of the highest resolved shear stress. Post-mortem fractography revealed a secondary population of dimples on these crack-like features, implying that void sheeting may be occurring during this crack-like void growth. Furthermore, these in-situ observations imply that classical models for void growth, that assume a homogenous medium, should be extended to account for the anisotropic grain-level behavior to correctly capture distinct features of void nucleation and growth.

Al-2219↗

DOE OSTI2025

The revised elastic field of an edge dislocation

We modify the traditional linear-elastic field of an edge dislocation. The modifications stem from symmetry and energy requirements imposed in terms of embedded (deformed) coordinates. These requirements are satisfied if a line force is added to the dislocation field. The modification to the stress field is expressed by coefficients that are a function only of Poisson’s ratio. Qualitatively, the field of the dislocation is increased in the glide direction and decreased in the climb direction. There are small changes to the field of a screw dislocation that do not entail the addition of a line force. Effects of nonlinearity, anisotropic elasticity, and image forces are briefly considered.

Hirth, J. P.↗

DOE OSTI2020

Sapphire advanced mitigation process: wet etch to expose sub-surface damage and increase laser damage resistance and mechanical strength

In this study, a novel, to the best of our knowledge, method of wet chemical etching of sapphire workpieces (such as optics, wafers, windows, and cones), called the sapphire advanced mitigation process (or sapphire AMP), has been developed that exposes sub-surface mechanical damage created during the optical fabrication process and significantly enhances the surface laser damage resistance (> 2×) and mechanical strength (up to ~ 2.6×). Sapphire AMP involves first treating the workpiece with a mixture of sulfuric and phosphoric acid ([H 2 SO 4 ]:[H 3 PO 4 ]= 1 : 3) at 220°C, followed with phosphoric acid at 160°C, then with sodium hydroxide base (NaOH) and surfactant at 40°C, and finally with a high-pressure deionized water spray rinse. Sapphire AMP has been demonstrated on both A- and C-plane sapphire workpieces. The mechanism of this etch process involves the reaction of the sapphire (Al 2 O 3 ) surface with sulfuric acid (H 2 SO 4 ) forming aluminum sulfate [Al2(SO4)3], which has low solubility. The high phosphoric acid content in the first and second steps of sapphire AMP results in the efficient conversion of Al 2 (SO 4 ) 3 to aluminum phosphate (AlPO 4 ), which is very soluble, greatly reducing reaction product redeposition on the workpiece surface. Sapphire AMP is shown to expose sub-surface mechanical damage on the sapphire surface created during the grinding and polishing processes, whose etched morphology has either isotropic or anisotropic evolution depending on the nature of the initial surface damage. Sapphire AMP was also designed to remove the key known surface, laser absorbing precursors (namely, foreign chemical impurities, the fracture surface layer of preexisting sub-surface damage, and reaction product or foreign species redeposition or precipitation). Static and sliding indention induced surface microfractures on sapphire are shown after sapphire AMP to have a significant decrease in the fast photoluminescence intensity (a known metric for measuring the degree of laser damaging absorbing precursors). In addition, the onset of laser damage (at 351 nm 3 ns) on sapphire AMP treated workpieces was shown to increase in fluence from ~ 4 to >9.5J/cm 2 . Finally, biaxial ball-on-ring mechanical tests on sapphire disks showed an increase in the failure stress from 340 MPa (with pre-existing 28 µm flaws) to ~ 900MPa after sapphire AMP, which is attributed to the blunting of the surface microfractures.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Reaction Heterogeneity in LiFePO 4 Agglomerates and the Role of Intercalation-Induced Stress

As an important battery cathode material, reaction distribution in lithium iron phosphate (LiFePO4) has been extensively studied in dispersed particle systems, but remains poorly understood for mesoscopic agglomerates (or secondary particles) that are used in most commercial batteries. Herein, we apply three-dimensional X-ray spectroscopic imaging to characterize the two-phase structure in LiFePO4 secondary particles during electrochemical cycling. (De)lithiated domains are found to not form the commonly assumed core-shell structure but develop highly anisotropic filamentary morphology that is rate independent and symmetric between charging and discharging. Phase-field simulations elucidate that the observed 1D phase growth behavior is not caused by the 1D lithium diffusivity of LiFePO4 but the elastic interaction between primary particles, which gives rise to stronger reaction heterogeneity than dispersed nanoparticles. As a result, uniform lithium (de)intercalation does not occur on the secondary particle surface even at high cycling rates.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Microstructure Anisotropy Effects on Fracture and Fatigue Mechanisms in Shape Memory Alloy Martensites

The chief scientific object was to investigate the fatigue and fracture mechanisms in NiTi shape memory alloys, in the presence of stress concentrating factors such as micro-pores and notches. In particular, the interactions between structure (factors related to the geometry of the specimen, such as holes, pores and notches) and microstructure (related to properties of the material such as grains, secondary phases, twins) in influencing the inelastic deformation, leading up to fracture in NiTi was examined. Primarily, there were two major thrusts. The first thrust focused on fracture behavior in Shape Memory Alloys where the high crystal symmetry cubic phase transformed into the low symmetry phase during deformation and propagation of cracks. This thrust revealed the interplay between size of micro-pores, size of grains and grain orientations in determining the resistance to crack growth, during cyclic tensile deformation of porous-like NiTi specimens. The role of secondary phases (impurities of Ti, often formed during fabrication) were also shown to induce local inelastic deformation gradients during tensile deformation and reduce overall extent of deformation. The second thrust focused on the fatigue and fracture behavior in the low crystal symmetry martensite phase. Unlike high symmetry phases, characterization of deformation in low symmetry phases such as martensitic NiTi has been challenging, due to a hierarchical microstructure spanning from nanoscale twins to the macroscale, a highly anisotropic microstructure stemming from the low crystal symmetry and a multitude of inelastic mechanisms for deformation. Thus, advanced X-ray techniques such as 3D high-energy diffraction microscopy and x-ray micro-Laue diffraction and spectroscopy were used in conjunction with traditional techniques such as SEM and TEM to understand the inelastic deformation in the martensite phase in NiTi, in the presence of notches. This thrust revealed a stark difference in deformation mechanisms based on the nature of stress: regions near the notch experience a highly triaxial stress while those far away experience a uniaxial stress. Additionally, while regions close to the notch were shown to deform locally by amorphization/nanocrystallization, ultimately leading to fracture, regions far away from the notch deformed in a more traditional manner, with detwinning and twin reorientation the primary mechanisms of deformation. More importantly, this thrust revealed that in the presence of structural factors such as notches, the traditional picture of reorientation (coalescence) of nanoscale martensite twins into increasingly larger domains with increased loading does not hold true, even up to fracture of the specimen. These conclusions will inform the modeling and understanding of deformation behavior of specimens, where structure and microstructure coexist and interact.

36 MATERIALS SCIENCE↗

DOE OSTI2022

The local wavenumber model for computation of turbulent mixing

We present an overview of the current status in the development of a two-point spectral closure model for turbulent flows, known as the local wavenumber (LWN) model. The model is envisioned as a practical option for applications requiring multi-physics simulations in which statistical hydrodynamics quantities such as Reynolds stresses, turbulent kinetic energy, and measures of mixing such as density-correlations and mix-width evolution, need to be captured with relatively high fidelity. In this review, we present the capabilities of the LWN model since it was first formulated in the early 1990s, for computations of increasing levels of complexity ranging from homogeneous isotropic turbulence, inhomogeneous and anisotropic single-fluid turbulence, to two-species mixing driven by buoyancy forces. The review concludes with a discussion of some of the more theoretical considerations that remain in the development of this model.

42 ENGINEERING↗

DOE OSTI2020

Solid state grain alignment of permanent magnets in near-final shape

Magnet microstructure manipulation in the solid state by controlled application of a sufficient stress in a direction during high temperature annealing in a single-phase region of heat-treatable magnet alloys, e.g., alnico-type magnets is followed by magnetic annealing and draw annealing to improve coercivity and saturation magnetization properties. The solid-state process can be termed highly controlled abnormal grain growth (hereafter AGG) and will make aligned sintered anisotropic magnets that meet or exceed the magnetic properties of cast versions of the same alloy types.

Anderson, Iver E.↗

DOE OSTI2022

Solid state grain alignment of permanent magnets in near-final shape

Anderson, Iver E.↗

DOE OSTI2020

Vulnerable window of yield strength for swelling-driven fracture of phase-transforming battery materials

Abstract Despite numerous experimental and theoretical investigations of the mechanical behavior of high-capacity Si and Ge Li-ion battery anodes, our basic understanding of swelling-driven fracture in these materials remains limited. Existing theoretical studies have provided insights into elasto-plastic deformations caused by large volume change phase transformations, but have not modeled fracture explicitly beyond Griffith’s criterion. Here, we use a multi-physics phase-field approach to model self-consistently anisotropic phase transformation, elasto-plastic deformation, and crack initiation and propagation during lithiation of Si nanopillars. Our computational results reveal that fracture occurs within a “vulnerable window” inside the two-dimensional parameter space of yield strength and fracture energy and highlight the importance of taking into account the surface localization of plastic deformation to accurately predict the magnitude of tensile stresses at the onset of fracture. They further demonstrate how the increased robustness of hollow nanopillars can be understood as a direct effect of anode geometry on the size of this vulnerable window. Those insights provide an improved theoretical basis for designing next-generation mechanically stable phase-transforming battery materials undergoing large volume changes.

25 ENERGY STORAGE↗

DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Pore-scale visualization of natural hydrate-bearing sediments

Accurate modeling of gas hydrate reservoir productivity and geomechanical risks associated with subsurface dissociation of natural gas hydrates (NGH) requires the determination of model parameters through physical testing on natural hydrate-bearing sediments (HBS). This involves investigating the hydro-mechanical behavior of undisturbed hydrate samples from nature under in situ conditions using pressure core characterization and analysis, which provides a unique opportunity for research. By employing state-of-the-art micro computed tomography imagery on cryogenically preserved, hydrate-bearing sediment samples, we can determine hydrate saturation as well as permeability with and without the presence of hydrates in the sediment. Furthermore, utilizing a machine learning based image segmentation technique, it is possible to extract pore space and grain information. Subsections of the entire image volume were used to determine anisotropic permeabilities using a finite-difference method Stokes solver (FDMSS). Additionally, permeability measurements on whole pressure and temperature preserved hydrate-bearing core were analyzed by utilizing the National Energy Technology Laboratory’s (NETL) Pressure Core Characterization and X-ray CT Visualization Tool (PCXT) to manipulate, cut, and analyze pressure preserved sediment. Permeabilities were measured under a broad range of vertical stress states to simulate expected pressure changes during production scenarios, and the results show that permeabilities derived from images are in agreement with those from traditional core derived experiments. The collected stress-dependent permeability, permeability anisotropy, and corresponding gas hydrate saturations provide valuable input into numerical simulations of reservoir productivity. These properties have been proven to be key parameters determining a long-term reservoir response under depressurization.

Liu, Mengwei [Oak Ridge Institute for Science and ↗