Search NASA⌕ Search

SEARCH · Search NASA

Results for “Slip localization”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Impact of the Earthquake Rupture on Ground-Motion Variability of the 24 August 2016 M w 6.2 Amatrice, Italy, Earthquake

Here, the devastating 24 August 2016 M w 6.2 earthquake that struck Amatrice, Italy, marked the beginning of a prolonged seismic sequence dominated by three subsequent M w ≥6.0 events in the central Apennines region. The earthquake destroyed Amatrice’s historic center, claiming the lives of 299 individuals and causing widespread damage in the neighboring villages. The severity of the ground shaking, with a recorded maximum acceleration of 850 cm/s 2 on the east–west component at the Amatrice station, was far greater than the predicted acceleration based on the Italian ground-motion model (GMM). As pointed out by several investigations, the observed ground-motion amplitude and its spatial variability during the earthquake can be linked to specific rupture characteristics, including slip distribution and rupture directivity effects revealed by the observed data (Tinti et al., 2016; Pischiutta et al., 2021). In this study, we conducted physics-based 3D numerical simulations of ground motion for the Amatrice earthquake for frequencies up to 3 Hz. We employed a series of kinematic rupture models and a well-constrained local 3D velocity model incorporating surface topography. The kinematic rupture realizations were generated using multiscale hybrid and fully stochastic models, following the technique proposed by Graves and Pitarka (2016). We focused on assessing the sensitivity of near-fault ground-motion amplitudes to earthquake rupture characteristics, in particular, the spatial slip pattern. To evaluate the quality of our simulations, we employed goodness-of-fit measurements performed in comparisons of simulated and recorded ground motions. The simulated ground motions compare well with the recorded data and predictions from GMMs for Italy, ITA18 (Lanzano et al., 2019). However, we found that the simulated interevent ground-motion variability (randomness in the source process) of peak ground velocity, σ (PGV) is higher than the constant σ (PGV) predicted by conventional GMMs. Our simulations using several rupture scenarios demonstrate that the near-fault ground-motion amplification pattern is directly related to the slip distribution pattern.

3D Ground Motions Simulations↗
DOE OSTI2025

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↗
DOE OSTI2020

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Statistics of internal stress fluctuations in dislocated crystals and relevance to density-based dislocation dynamics models

A statistical analysis of internal stress fluctuations, defined as the difference between the local mean stress and stress on dislocations, is presented for deforming crystals with 3D discrete dislocation systems. Dislocation realizations are generated using dislocation dynamics simulations and the associated stress field is computed as a superposition of a regularized stress field of dislocation lines within the domain of the solution and a complementary stress field computed via a finite-element boundary value problem. The internal stress fluctuations of interest are defined by an ensemble of the difference between the stress on dislocation lines and the local mean field stress in the crystal. The latter is established in a piecewise fashion over small voxels in the crystal thus allowing the difference between the local average stress and stress on segments to be easily estimated. The results show that the Schmid stress (resolved shear stress) and Escaig stress fluctuations on various slip systems sampled over a random set of points follow a Cauchy (Lorentz) distribution at all strain levels, with the amplitude and width of the distribution being dependent on the strain. Finally, the implications of the Schmid and Escaig internal stress fluctuations are discussed from the points of view of dislocation cross-slip and the dislocation motion in continuum dislocation dynamics.

36 MATERIALS SCIENCE↗
NASA NTRS1990

A crystallographic model for the tensile and fatigue response for Rene N4 at 982 C

An anisotropic constitutive model based on crystallographic slip theory was formulated for nickel-base single-crystal superalloys. The current equations include both drag stress and back stress state variables to model the local inelastic flow. Specially designed experiments have been conducted to evaluate the existence of back stress in single crystals. The results showed that the back stress effect of reverse inelastic flow on the unloading stress is orientation-dependent, and a back stress state variable in the inelastic flow equation is necessary for predicting inelastic behavior. Model correlations and predictions of experimental data are presented for the single crystal superalloy Rene N4 at 982 C.

Sheh, M. Y.↗
NASA NTRS1977

The 1974 Ethiopian rift geodimeter survey

The field techniques and methods of data reduction for five successive geodimeter surveys in the Ethiopian rift valley are enlarged upon, with the considered conclusion that there is progressive accumulation of upper crustal strain, consonant with on-going rift extension. The extension is restricted to the Quaternary volcanotectonic axis of the rift, namely the Wonji fault belt, and is occurring at rates of 3 to 6 mm/yr in the northern sector of the rift valley. Although this concurs with the predictions of platetectonic analysis of the Afar triple junction, it is considered premature to endorse such a concurrence on the basis of only 5 years of observations. This is underlined by the detection of local tectonic contractions and expansions associated with geothermal and gravity anomalies in the central sector of the rift valley. There is a hint of a component of dextral slip along some of the rift-floor fault zones, both from geological evidence and from the strain patterns detected in the present geodetic surveys.

Mohr, P.↗
NASA NTRS1977

Combustion of liquid sprays at high pressures

The combustion of pressure atomized fuel sprays in high pressure stagnant air was studied. Measurements were made of flame and spray boundaries at pressures in the range 0.1-9 MPa for methanol and n-pentane. At the higher test pressure levels, critical phenomena are important. The experiments are compared with theoretical predictions based on a locally homogeneous two-phase flow model. The theory correctly predicted the trends of the data, but underestimates flame and spray boundaries by 30-50 percent, indicating that slip is still important for the present experiments (Sauter mean diameters of 30 microns at atmospheric pressure under cold flow conditions). Since the sprays are shorter at high pressures, slip effects are still important even though the density ratio of the phases approach one another as the droplets heat up. The model indicates the presence of a region where condensed water is present within the spray and provides a convenient means of treating supercritical phenomena.

Shearer, A. J.↗
NASA NTRS1986

Thermofluid analysis of the SSME preburner using a gas-gas diffusion model for oxygen and hydrogen combustion at supercritical pressures

The paper discusses the thermofluid analysis of the Space Shuttle Main Engine (SSME) fuelside preburner. The governing equations have been solved numerically to predict flow, heat transfer, mixing, and combustion. A two-fluid approach is adopted in which oxygen is regarded as one fluid and hydrogen is regarded as the other fluid. The chemical kinetics is assumed to be very fast so that combustion is primarily controlled by the rate of mixing between oxygen and hydrogen. The preburner pressure is much greater than the critical pressures of oxygen and hydrogen; hence, a gas-gas diffusion model (rather than an evaporation model) has been developed to compute the rate of interphase mixing. Empirical correlations have been incorporated to account for the effect of slip on the interphase exchange. A sensitivity study has been performed with various model parameters. It is observed that the model can predict possibility of incomplete combustion and local regions of high temperatures under steady operating conditions. Some of these anomalies have been observed in actual tests, and the numerical model is useful for understanding possible causes and remedies. At least some measurements are needed for quantitative verification of the model.

Prakash, C.↗
NASA NTRS2000

A Unified Wall Boundary Treatment for Viscous and Inviscid Flows in the CE/SE Method

In the setting of the conservation element-solution element (CE/SE) method, a new and unified wall boundary treatment for the Navier-Stokes and Euler Equations is proposed. In essence, the shear stress exerted on the fluid by a wall is modeled as a source term as a part of local spacetime flux conservation in the vicinity of a wall boundary. When the fluid is inviscid, the source term vanishes and the boundary condition reduces to the usual 'slip' condition. On the other hand, when the fluid is viscous, the source-term effect is consistent with the traditional no-slip condition. Numerical results show that the new treatment is robust, efficient, and accurate for viscous and inviscid flows.

Chang, Sin-Chung↗
DOE OSTI2021

Locally commensurate charge-density wave with three-unit-cell periodicity in YBa 2 Cu 3 O y

In order to identify the mechanism responsible for the formation of charge-density waves (CDW) in cuprate superconductors, it is important to understand which aspects of the CDW’s microscopic structure are generic and which are material-dependent. Here, we show that, at the local scale probed by NMR, long-range CDW order in YBa 2 Cu 3 O y is unidirectional with a commensurate period of three unit cells ( λ = 3 b ), implying that the incommensurability found in X-ray scattering is ensured by phase slips (discommensurations). Furthermore, NMR spectra reveal a predominant oxygen character of the CDW with an out-of-phase relationship between certain lattice sites but no specific signature of a secondary CDW with λ = 6 b associated with a putative pair-density wave. These results shed light on universal aspects of the cuprate CDW. In particular, its spatial profile appears to generically result from the interplay between an incommensurate tendency at long length scales, possibly related to properties of the Fermi surface, and local commensuration effects, due to electron-electron interactions or lock-in to the lattice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Modeling the role of local crystallographic correlations in microstructures of Ti-6Al-4V using a correlated structure visco-plastic self-consistent polycrystal plasticity formulation

This paper presents a multi-level crystal plasticity-based simulation framework for modeling mechanical response and microstructure evolution of Ti-6Al-4V with α-lath/lamellar microstructures. Here, the model is a correlated structure visco-plastic self-consistent (CS-VPSC) formulation linking three scales: a single crystals micro-scale, a lath/lamellar colony meso-scale, and a lath/lamellar aggregate macro-scale. A selected hardening law for the evolution of critical resolved shear stress per slip system used in CS-VPSC is phenomenological. However, it adjusts the resistances of basal and prismatic slip systems based on the geometry of slip transfer between adjacent lamellae. Consistent with experimental evidences, the resolved shear stress on the pyramidal slip planes is dependent not only on the stress in the direction of slip but also on the two orthogonal shear stress components and the three normal stress components (non-Schmid effects). Electron backscatter diffraction (EBSD) data in conjunction with a procedure relying on α→β phase transformation is used to construct paired variants of α-laths/lamellae satisfying their local crystallographic correlations. The procedure fits volume fractions of individual laths/lamellae with the experimental EBSD data and selects a distribution of habit planes between adjacent variants with respect to the loading direction. The simulation framework is applied to interpret the deformation behavior in tension and compression along two sample directions of Ti-6Al-4V fabricated via laser powder bed fusion. Moreover, the model is used to simulate texture evolution during rolling of the material to large strains. It is demonstrated that the model is capable of predicting plastic anisotropy/asymmetry and the concomitant texture evolution. While the model reveals a significant effect of habit plane inclination with respect to the loading direction on yield stress, the comparison of the data and model predictions shows that a random distribution of habit planes fits the flow response. It is further inferred that the tension-compression asymmetry arises from the non-Schmid effects.

36 MATERIALS SCIENCE↗
NASA NTRS1984

The Corfu Landslide: Analog to Giant Landslides on Mars

In an analog to the great landslides of the Vales Marineris, Mars, a detailed study was made of the Corfu Landslide in south-central Washington. This prehistoric slide is located on the northern flank of the Saddle Mountains, southwest of Othello, Washington. The slide covers a 13 square km area centered on section 11 of T.15N., R.27E., Willamette Meridian, adjacent to the Corfu townsite. Approximately 1 cubic km of material is involved in sliding that was probably initiated by Missoula flooding through the Channeled Scabland. It is concluded that there were four primary factors involved in the initiation of the Corfu landsliding: (1) A slip surface was present at the right orientation; (2) Glacial flooding undercut the slope; (3) Wetter climatic conditions prevailed during that time period; and (4) Some seismic vibrations, known to occur locally, probably acted as a trigger. These factors show that special conditions were required in conjunction to produce landsliding. Studies in progress of the Vales Marieneris suggest that the same factors probably contributed to landsliding there.

Scott W. Lewis↗
NASA NTRS1988

Geodetic Measurement of Deformation East of the San Andreas Fault in Central California

Triangulation and trilateration data from two geodetic networks located between the western edge of the Great Valley and the San Andreas fault have been used to calculate shear strain rates in the Diablo Range and to estimate the slip rate along the Calaveras and Paicines faults in Central California. Within the Diablo Range the average shear strain rate was determined for the time period between 1962 and 1982 to be 0.15 + or - 0.08 microrad/yr, with the orientation of the most compressive strain at N 16 deg E + or - 14 deg. The orientation of the principal compressive strain predicted from the azimuth of the major structures in the region is N 25 deg E. It is inferred that the measured strain is due to compression across the folds of this area: the average shear straining corresponds to a relative shortening rate of 4.5 + or - 2.4 mm/yr. From an examination of wellbore breakout orientations and the azimuths of P-axes from earthquake focal mechanisms the inferred orientation of maximum compressive stress was found to be similar to the direction of maximum compressive strain implied by the trend of local fold structures. Results do not support the hypothesis of uniform fault-normal compression within the Coast Ranges. From trilateration measurements made between 1972 and 1987 on lines that are within 10 km of the San Andreas fault, a slip rate of 10 to 12 mm/yr was calculated for the Calaveras-Paicines fault south of Hollister. The slip rate of the Paicines fault decreases to 4 mm/yr near Bitter.

Sauber, Jeanne M.↗
DOE OSTI2022

A particle resolved simulation approach for studying shock interactions with moving, colliding solid particles

This work applies a new combination of techniques for the fully resolved simulation of compressible, gas–particle multiphase flows. Here, the adaptive wavelet collocation method is used to dynamically, and efficiently, adapt the computational grid to localized flow features and the particles. A characteristic-based volume penalization method that imposes arbitrary Dirichlet, Neumann, or Robin-type immersed boundary conditions, is used to enforce the no-slip condition at particle surfaces. A hard-sphere collision model is applied to capture the particle–particle collisions. Proof of concept test cases are presented, showcasing the dynamic grid adaptation and fully resolved two-way coupling between the phases that is possible with this approach. Results for a shock-driven single cylinder under viscous and inviscid conditions are presented along with a demonstration of a shock interacting with a cloud of randomly distributed cylinders and spheres.

42 ENGINEERING↗
NASA NTRS1992

Environmental fatigue of an Al-Li-Cu alloy. Part 2: Microscopic hydrogen cracking processes

Based on a fractographic analysis of fatigue crack propagation (FCP) in Al-Li-Cu alloy 2090 stressed in a variety of inert and embrittling environments, microscopic crack paths are identified and correlated with intrinsic da/dN-delta K kinetics. FCP rates in 2090 are accelerated by hydrogen producing environments (pure water vapor, moist air, and aqueous NaCl), as defined in Part 1. For these cases, subgrain boundary fatigue cracking (SGC) dominates for delta K values where the crack tip process zone, a significant fraction of the cyclic plastic zone, is sufficiently large to envelop 5 micron subgrains in the unrecrystallized microstructure. SGC may be due to strong hydrogen trapping at T1 precipitates concentrated at sub-boundaries. At low delta K, the plastic zone diameter is smaller than the subgrain size and FCP progresses along (100) planes due to either local lattice decohesion or aluminum-lithium hydride cracking. For inert environments (vacuum, helium, and oxygen), or at high delta K where the hydrogen effect on da/dN is small, FCP is along (111) slip planes; this mode does not transition with increasing delta K and plastic zone size. The SGC and (100) crystallographic cracking modes, and the governing influence of the crack tip process zone volume (delta K), support hydrogen embrittlement rather than a surface film rupture and anodic dissolution mechanism for environmental FCP. Multi-sloped log da/dN-log delta K behavior is produced by changes in process zone hydrogen-microstructure interactions, and not by purely micromechanical-microstructure interactions, in contradiction to microstructural distance-based fatigue models.

Robert S Piascik↗
NASA NTRS2005

Slip compensation for a Mars Rover

This paper describes the design, implementation, and experimental results of an integrated system that enables Mars rover to navigate in a high slip environment.

navigation↗
NASA NTRS2009

Simulation of RCC Crack Growth Due to Carbon Oxidation in High-Temperature Gas Environments

The carbon wall oxidation technique coupled with a CFD technique was employed to study the flow in the expanding crack channel caused by the oxidation of the channel carbon walls. The recessing 3D surface morphing procedure was developed and tested in comparison with the arcjet experimental results. The multi-block structured adaptive meshing was used to model the computational domain changes due to the wall recession. Wall regression rates for a reinforced carbon-carbon (RCC) samples, that were tested in a high enthalpy arcjet environment, were computationally obtained and used to assess the channel expansion. The test geometry and flow conditions render the flow regime as the transitional to continuum, therefore Navier-Stokes gas dynamic approach with the temperature jump and velocity slip correction to the boundary conditions was used. The modeled mechanism for wall material loss was atomic oxygen reaction with bare carbon. The predicted channel growth was found to agree with arcjet observations. Local gas flow field results were found to affect the oxidation rate in a manner that cannot be predicted by previous mass loss correlations. The method holds promise for future modeling of materials gas-dynamic interactions for hypersonic flight.

Titov, E. V.↗
DOE OSTI2023

Dislocation generation in diamond under extreme loading

Diamond is, by virtue of its high bonding forces and Peierls-Nabarro barrier stresses, among the hardest materials on earth due to the difficulty of generating and moving dislocations. We demonstrate, using molecular dynamics and analytical calculations, that the generation of defects is dependent on loading orientation. Shock-loading single-crystal diamond along [001] and [011] to a stress of 137 GPa did not reveal dislocations, whereas loading along [111] generated profuse dislocations. The introduction of a void at the latter orientation generated dislocations at 72.6 GPa. Two slip systems were identified: <011>{100} and <112>{111}. The threshold for plastic deformation in diamond is orientation dependent and significantly lower than reported in previous studies: 0.14 G (shear modulus). The void collapse generated localized amorphization for [001]. In conclusion, these results on the anisotropy of plasticity and void effects have relevance for the improvement of symmetry in the collapse of diamond capsules in inertial fusion experiments.

36 MATERIALS SCIENCE↗