Introducing PyFaultSlip: A free and open-source tool for the assessment of induced fault slip hazards from deep fluid injection
Not Available
Polun, Sean G. (ORCID:0000000261873687)↗
SEARCH · Search NASA
Results for “Slip”
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.
New molecular dynamics studies of hydrogen effects on cross-slip energy barriers in austenitic stainless steels
Explore the source record for details and available documents.
Zhou, Xiaowang [Sandia National Laboratories (SNL-↗
Heterogeneous fatigue damage in a nickel-based single-crystal superalloy unraveled using correlative 3D X-ray technology
Nickel-based single-crystal (Ni-SX) superalloys under cyclic stress are susceptible to cracking at stress-concentration sites, eventually leading to low-cycle fatigue (LCF) failure. LCF cracks typically originate from intrinsic defects (e.g., voids and carbides) within solidified dendrites. However, systematic quantitative experimental analyses of defect-mediated local damage remain limited. To thoroughly understand the microscopic origins and evolution of LCF damage, correlated 3D mapping of dendrites across various regions is essential. Here, in this study, macroscale micro-computed tomography (μ-CT) was initially used to capture internal interdendritic secondary cracks within bulk DD413 superalloy after LCF testing at 760 °C. Subsequently, a multimodal methodology combining synchrotron 3D microdiffraction (3D-μXRD), high-resolution μ-CT, and electron microscopy was established. This approach allowed precise localization of internal damage zones near interdendritic secondary cracks and detailed mapping of the 3D correlated distributions of dendrites, defects, and residual stress/strain fields within these zones at submicron spatial resolution. Finally, the same approach was applied to specimens subjected to interrupted loading at approximately 40 % of the fatigue life to uncover the early damage states of dendrites. The dendrite cores (DCs) and interdendritic regions (IDs) exhibit microscale heterogeneous mechanical responses: nearly defect-free DCs accumulate local irreversible slip along specific slip systems to generate slip bands, while the IDs containing various defects accommodate local microplasticity through the activation of multiple slip systems around these defects. The local tensile stress near defects in the IDs exceeds that in the DC slip band regions by more than threefold, leading to the generation of local damage zones within the IDs. Chain-like defect distributions facilitate the interconnection of these local zones into a continuous damage region, further elevating the overall tensile stress in the IDs. Additionally, geometrically necessary dislocations alone are insufficient as indicators of LCF damage; both the internal stress state and its magnitude must be considered. These experimental results provide critical data and insights for the development of multi-physics fatigue models.
Localized deformation↗
Evidence for a Single Holocene Paleoseismic Event on the Pajarito Fault, Northern New Mexico
Low-slip rate fault systems tend to be less studied than their high-slip rate counterparts, and paleoseismic techniques used to study them may pose challenges in interpretation that differ from high-slip rate systems. A good example of this is the Pajarito fault system (PFS), a normal fault complex within the Rio Grande rift. Despite numerous previous paleoseismic trenching studies conducted on the PFS between 1990 and 2003, considerable uncertainty remains regarding its Holocene paleoseismic history, particularly for the primary Pajarito fault (PF). To further clarify the PF paleoseismic history, we present data from paleoseismic investigations of 6 trenches at 3 distinct locations along the PF. Though the totality of the age and structural data obtained in this study is complex and not entirely consistent with any one interpretation, a single Holocene paleoearthquake occurring younger than ∼1,600 to 2,300 kcal yr BP is the simplest interpretation. It is possible that the PF records two Holocene events, with a penultimate event 6.9–2.4 kcal yr BP event and the aforementioned most recent event (MRE) between 2.3 and 1.6 kcal yr BP. However, only a single wall of one trench, out of a total of 12 walls in our 6 trenches, provides evidence supporting that interpretation. This study finds evidence of a single late Holocene paleoseismic event on the PF and sparse evidence for 2 Holocene paleoseismic events on the PF and highlights the benefits of logging multiple trench walls to better understand the complexity that results from this low-slip rate, low-deposition-rate fault system.
58 GEOSCIENCES↗
Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models
In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.
Crystal plasticity↗
In-situ observations of cyclic deformation in an extruded Mg-2Nd-1Y-0.1Zr-0.1Ca alloy
In this study, the evolution of deformation mechanisms during cyclic loading in an extruded, solution-treated Mg–2Nd–1Y–0.1Zr–0.1Ca alloy was investigated using a combination of in-situ loading, scanning electron microscopy (SEM), electron backscatter diffraction (EBSD), and focused ion beam (FIB) nanofabrication. The initial microstructure exhibited a random crystallographic texture with no preferred grain orientation. Flat, rectangular dog-bone specimens were subjected to load-controlled, fully reversed fatigue for 50 cycles, during which the same region was sequentially mapped to track microstructural changes. After 10 cycles of loading deformation twins were observed. During tensile reloading detwinning or narrowing of those twinned regions occurred. After 20 cycles, detwinning ceased and residual twins remained in the material. SEM imaging revealed numerous surface slip traces after cyclic loading. EBSD-assisted slip trace analysis identified the activation of prismatic and pyramidal < c+a> slip systems during low-cycle fatigue. Site-specific scanning transmission electron microscopy (STEM) further revealed that deformation was also accommodated by basal < a> slip and the dissociation of < c+a> dislocations. Center-of-symmetry (COS) analysis confirmed that the dissociation of < c+a> dislocations resulted in the formation of I₁ intrinsic stacking faults after cyclic loading. These findings provide new insights into the complex interplay of dislocation mechanisms governing fatigue deformation in rare-earth-containing Mg alloys.
Cyclic deformation↗
Link statistics of dislocation network during strain hardening
Dislocations are line defects in crystals that multiply and self-organize into a complex network during strain hardening. The length of dislocation links, connecting neighboring nodes within this network, contains crucial information about the evolving dislocation microstructure. By analyzing data from Discrete Dislocation Dynamics (DDD) simulations in face-centered cubic (fcc) Cu, we characterize the statistical distribution of link lengths of dislocation networks during strain hardening on individual slip systems. Here, our analysis reveals that link lengths on active slip systems follow a double-exponential distribution, while those on inactive slip systems conform to a single-exponential distribution. The distinctive long tail observed in the double-exponential distribution is attributed to the stress-induced bowing out of long links on active slip systems, a feature that disappears upon removal of the applied stress. We further demonstrate that both observed link length distributions can be explained by extending a one-dimensional Poisson process to include different growth functions. Specifically, the double-exponential distribution emerges when the growth rate for links exceeding a critical length becomes super-linear, which aligns with the physical phenomenon of long links bowing out under stress. This work advances our understanding of dislocation microstructure evolution during strain hardening and elucidates the underlying physical mechanisms governing its formation.
Crystal plasticity↗
Slow Quasiparticle Dynamics and Anyonic Statistics in a Fractional Quantum Hall Fabry-Pérot Interferometer
Anyons are two-dimensional particles with fractional exchange statistics that emerge as elementary excitations of fractional quantum Hall phases. Experimentally, their exchange statistics can be measured in the edge-state Fabry-Pérot interferometer, wherein the presence of 𝑁 𝑞𝑝 localized anyons contributes a phase 𝑁 𝑞𝑝 𝜃 𝑎 to the interference pattern where 𝜃 𝑎 is twice the exchange phase. Here we report the observation of large, hysteretic phase jumps in a monolayer graphene Fabry-Pérot interferometer at 𝜈 = 1/3. When the filling factor is increased from 𝜈 < 1/3 toward the center of the plateau, we observe phase slips with magnitude Δ𝜃 ≈ 2𝜋/3, consistent with the addition of individual quasiparticles to the interferometer bulk. These phase slips occur as instantaneous jumps in the interference signal, with intervals between the jumps indicating quasiparticle equilibration times exceeding 20 min. We use this long timescale to investigate the effect of changes in interferometer area 𝐴 𝐼 and 𝑁 𝑞𝑝 independently at fixed magnetic field, revealing a striking memory effect in the phase slip magnitude. In particular, as the 𝜈 =1/3 plateau is approached from higher filling, we observed phase slips with Δ𝜃 significantly larger than 2𝜋/3 over the same range of gate voltage where quantized jumps are seen for increasing 𝜈. We discuss this asymmetry in terms of bulk-edge coupling of quasiparticles localized near the edge or in the bulk, and argue that this effect can be qualitatively reconciled with theoretical expectations for strongly interacting quasiparticles in the presence of weak disorder and strongly nonequilibrium charge dynamics. Besides providing a replication of interferometric measurements sensitive to 𝜃 𝑎 , our results highlight the key role played by charge dynamics on signatures of the anyon phase, and demonstrate that fractional quasiparticles can be indefinitely localized in nonequilibrium configurations.
anyons↗
AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES
Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.
Fraze, Matthew↗
Assessing Seismic Risk for CO2 Geologic Storage: Comparative Analysis of the Delaware Basin and Basin and Range Province Projects
ABSTRACT: Effective management of induced seismicity is critical for safe and sustainable CO2 storage. This study evaluates fault slippage risks in the Delaware Basin (Texas) and Basin and Range Province (Utah), integrating geological, operational, and geomechanical parameters to assess fault stability and seismic hazard mitigation. In the Delaware Basin, two sites were analyzed under an injection rate of 20,000 bbl/day over 25 years. One site showed low fault slip risk, while the other exhibited higher reactivation potential due to proximity to critically stressed faults. Sensitivity analysis revealed that increased pore pressure significantly heightened slip potential, highlighting the necessity of precise pressure control and real-time monitoring. In the Basin and Range Province, fault stability was evaluated at Neck of the Desert, Escalante Desert, Parowan, and Beaver sites under injection rates of 8,750 bbl/day per site over 30 years. Minimal fault slip risk was observed at Neck of the Desert and Escalante Desert sites, whereas Parowan and Beaver sites exhibited elevated slip potential due to semi-critically stressed faults sensitive to modest pore pressure increases. The findings demonstrate that fault slippage analysis, combined with sensitivity analysis of pore pressure and friction coefficients, is essential for understanding seismic risks. Continuous monitoring, adaptive injection management, and rigorous geomechanical analysis are key strategies for minimizing induced seismicity in CO2 sequestration projects.
58 GEOSCIENCES↗
Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2
This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.
Fault permeability↗
Comparative studies of cross-phase dynamics in turbulent momentum flux and particle flux at the tokamak edge
Turbulent transport events, including turbulent transport flux of momentum (i.e., turbulent momentum flux or Reynolds stress) and turbulent transport flux of particle (i.e., turbulent particle flux), have important effects on the confinement performance of magnetic confinement fusion devices. Poloidal Reynolds stress is the ensemble average of the product of radial velocity fluctuations and poloidal velocity fluctuations, i.e., $\langle {\widetilde{v}}_{r}{\widetilde{v}}_{\theta }\rangle$. Turbulent particle flux is the ensemble average of the product of radial velocity fluctuations and density fluctuations, i.e., $\langle \widetilde{n}{\widetilde{v}}_{r}\rangle$. Changes in either amplitude of fluctuations or cross phase between fluctuations can cause changes in turbulent transport. In this paper, cross-phase dynamics in the Reynolds stress and turbulent particle flux at the tokamak edge are studied in detail. Reynolds stress and turbulent particle flux are, respectively, written as the product of fluctuation amplitudes and an average cross-phase factor. The mathematical expressions of the average cross-phase factors are derived. The average cross-phase factors and the power spectra of cross phase are obtained using experimental measurement data. It is found that the cross-phase dynamics in Reynolds stress and particle flux are very different. Reynolds stress is found to be more sensitive to cross phase than particle flux is. In the strong $E\times B$ shear layer, spatial slips of cross phase lead to the obvious radial gradient of Reynolds stress. In the no/weak $E\times B$ shear region, the cross phase in Reynolds stress tends to lock. Here, phase locking refers to that the power spectra of phase tend to distribute around a fixed phase which does not change with radial position, while phase slip means that the power spectra of cross phase tend to distribute around a phase that varies with radial position. Phase slip or locking mainly describes the central phase weighted by the power spectra, while the phase scattering mainly describes the dispersion of the power spectrum distribution of the phase. The increased scattering of cross phase, which indicates the power spectra distribution of the phase is more dispersed, contributes to the decreased Reynolds stress for higher collisionality. The cross phase in particle flux tends to lock in both strong and no/weak shear regions. The degree of scattering of cross phase in the particle flux does not change obviously as collisionality increases. For higher collisionality, it is the increased density fluctuation amplitude rather than cross-phase dynamics that leads to the increased particle flux. The underlying physical mechanism that causes Reynolds stress and particle flux to exhibit different phase dynamics is discussed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Coupled geomechanical investigation of depletion-induced fault reactivation
Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.
Coupled geomechanical modeling↗
Hydrogen influences thermal activation parameters for dislocation glide during low cycle fatigue of 316L stainless steel
Measurements of activation areas are used to investigate the effect of hydrogen on the kinetics of dislocation glide during cyclic deformation in cold-worked 316L stainless steel. Non-charged and hydrogen-precharged (H-precharged) specimens were tested in low cycle fatigue (LCF) under plastic strain control. A series of plastic strain rate changes was performed periodically at the peak true plastic strain from the first cycle to half-life, and at various plastic strain values around stable hysteresis loops near half-life to determine the operational activation area, Δ a ∗. Both material conditions demonstrate a rapid increase in Δ a ∗ during the initial rapid softening followed by a region of approximately constant values coinciding with a reduced rate of softening. Near half-life, hydrogen reduces Δ a ∗ at a given true stress due to its effect on the activation distance and obstacle spacing. The magnitudes of Δ a ∗ reveal that bypassing solutes, cutting forest dislocations, and initiating cross slip are important mechanisms of thermally activated dislocation glide at all amplitudes, except hydrogen suppresses cross slip at the lowest plastic strain amplitudes. These results are supported by electron microscopy characterization of deformed microstructures. A Haasen plot analysis indicates that forest dislocations control the kinetics of deformation in both material conditions. It also reveals the presence of athermal obstacles in both non-charged and H-precharged conditions, likely dense dislocation tangles and cell walls. Additionally, the effect of hydrogen on microstructure evolution (by reducing the propensity for cross slip) leads to a dependence of athermal stress on plastic strain amplitude.
Activation area↗
Understanding the deformation behavior of the γ rich transformative Fe 38.5 Mn 20 Co 20 Cr 15 Si 5 Cu 1.5 complex concentrated alloy using in situ synchrotron diffraction
In situ tensile testing coupled with synchrotron x-ray diffraction was used to study the deformation behavior of metastability-engineered Fe 38.5 Mn 20 Co 20 Cr 15 Si 5 Cu 1.5 complex concentrated alloy. Monitoring the evolution of phase fraction and strain hardening response allowed the determination of true critical stress for the onset of the transformation-induced plasticity (TRIP) to be ∼375 MPa, preceded by slip starting at ∼255 MPa. In situ EBSD was used to validate the critical stress for transformation at the microstructure level and observe slip traces to confirm prior slip activity before the transformation. Further, a modeling framework based on stacking fault energy (SFE) was developed to predict the critical stress for transformation. Modeling suggested the SFE of the alloy to fall nearly 15 mJ/m 2 , which agrees well with SFE values calculated using synchrotron peak shifting (12 mJ/m 2 ) and thermodynamic calculation (11 mJ/m 2 ). As a result of γ-fcc to ε-hcp phase transformation, new {0002} ε planes emerged parallel to unaligned {111} γ planes with the tensile loading following S-N orientation relationship. Such selective emergence of new diffraction rings corresponding to ε phase is understood based on the reorientation of γ crystals with reference to tensile axis. In conclusion, this approach can be extended to effectively design alloys based on critical stress required for activating different deformation mechanisms to further push the limits of the strength-ductility envelope.
Complex concentrated alloy↗
Deformation of Hydrous Phases Egg [AlSiO 3 (OH)], δ [AlO(OH)] and Stishovite [Si 1‐n H 4n O 2 ] Relevant to Anisotropy of the Earth's Mantle
Seismic anisotropy of the Earth's mantle has been mostly attributed to crystallographic preferred orientation (CPO) generated during subduction and convection of an anhydrous mantle. But some hydrous phases are also stable at mantle conditions. Here we present results from diamond-anvil cell deformation experiments at high pressure and temperature on hydrous phases Egg [AlSiO 3 (OH)], δ [AlO(OH)] and hydrous stishovite [Si 1-n H 4n O 2 ], transformed from the clay mineral kaolinite. They develop strong CPO during axial compression, suggesting that they likely contribute to seismic anisotropy and heterogeneity in the mantle. Comparing experimental results with viscoplastic polycrystal plasticity models suggest that phase Egg deforms dominantly by (001) slip, δ by (010) slip and stishovite by {100} slip which could be incorporated in future models of mantle geodynamics.
58 GEOSCIENCES↗
Spherical and sessile droplet dynamics by fluctuating hydrodynamics
We simulate the mesoscopic dynamics of droplets formed by phase-separated fluids at nanometer scales where thermal fluctuations are significant. Both spherical droplets fully immersed in a second fluid and sessile droplets which are also in contact with a solid surface are studied. Our model combines a Cahn–Hilliard formulation with incompressible fluctuating hydrodynamics; for sessile droplets, the fluid–solid contact angle is specified as a boundary condition. Deterministic simulations with an applied body force are used to measure the droplets' mobility from which a diffusion coefficient is obtained using the Einstein relation. Stochastic simulations are independently used to obtain a diffusion coefficient from a linear fit of the variance of a droplet's position with time. In some scenarios, these two measurements give the same value but not in the case of a spherical droplet initialized near a slip wall or in the case of sessile droplets with large contact angles (≥90°) on both slip and no-slip surfaces.
Bell, John B↗