Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Fault structure, stress, or pressure control of the seismicity in shale? Insights from a controlled experiment of fluid-induced fault reactivation

It is clear that clay formations are present in reservoirs and earthquake faults, but questions remain on their mechanical behavior, as they can vary from ductile (aseismic) to brittle (seismic). An experiment, at a scale of 10 m, aims to reactivate a natural fault by fluid pressure in shale materials. The injection area was surrounded by a dense monitoring network comprising pressure, deformation, and seismicity sensors, in a well-characterized geological setting. Thirty-two microseismic events were recorded during several injection phases in five different locations within the fault zone. Their computed magnitude ranged between -4.3 and -3.7. Their spatiotemporal distribution, compared with the measured displacement at the injection points, shows that most of the deformation induced by the injection is aseismic. Whether the seismicity is controlled by the fault architecture, mineralogy of fracture filling, fluid, and/or stress state is then discussed. The fault damage zone architecture and mineralogy are of crucial importance, as seismic slip mainly localizes on the sealed-with-calcite fractures which predominate in the fault damage zone. As no seismicity is observed in the close vicinity of the injection areas, the presence of fluid seems to prevent seismic slips. The fault core acts as an impermeable hydraulic barrier that favors fluid confinement and pressurization. Therefore, the seismic behavior seems to be strongly sensitive to the structural heterogeneity (including permeability) of the fault zone, which leads to a heterogeneous stress response to the pressurized volume.

58 GEOSCIENCES↗
DOE OSTI2017

Toward full exploitation of coherent and incoherent information in Sentinel-1 TOPS data for retrieving surface displacement: Application to the 2016 Kumamoto (Japan) earthquake

Sentinel-1's continuous observation program over all major plate boundary regions makes it well suited for earthquake studies. However, decorrelation due to large displacement gradients and limited azimuth resolution of the Terrain Observation by Progressive Scan (TOPS) data challenge acquiring measurements in the near field of many earthquake ruptures and prevent measurements of displacements in the along-track direction. Here we propose to fully exploit the coherent and incoherent information of TOPS data by using standard interferometric synthetic aperture radar (InSAR), split-bandwidth interferometry in range and azimuth, swath/burst-overlap interferometry, and amplitude cross correlation to map displacements in both the line-of-sight and the along-track directions. Application to the 2016 Kumamoto earthquake sequence reveals the coseismic displacements from the far field to the near field. By adding near-field constraints, the derived slip model reveals more shallow slip than obtained when only using far-field data from InSAR, highlighting the importance of exploiting all coherent and incoherent information in TOPS data.

58 GEOSCIENCES↗
DOE OSTI2017

Strength anisotropy of shales deformed under uppermost crustal conditions

Conventional triaxial tests were performed on three sets of samples of Tournemire shale along different orientations relative to bedding (0°, 45°, and 90°). Experiments were carried out up to failure at increasing confining pressures ranging from 2.5 to 160 MPa, at strain rates ranging between 3 × 10 -7 s -1 and 3 × 10 -5 s -1 . This allowed us to determine the entire anisotropic elastic compliance matrix as a function of confining pressure. Results show that the orientation of principal stress relative to bedding plays an important role on the brittle strength, with 45° orientation being the weakest. We fit our results with a wing crack micromechanical model and an anisotropic fracture toughness. We found low values of internal friction coefficient and apparent friction coefficient in agreement with friction coefficient of clay minerals (between 0.2 and 0.3) and values of K /c comparable to that already published in the literature. We also showed that strain rate has a strong impact on peak stress and that dilatancy appears right before failure and hence highlighting the importance of plasticity mechanisms. Although brittle failure was systematically observed, stress drops and associated slips were slow and deformation always remained aseismic (no acoustic emission were detected). Lastly, this confirms that shales are good lithological candidates for shallow crust aseismic creep and slow slip events.

58 GEOSCIENCES↗
DOE OSTI2017

Temporal variation of intermediate-depth earthquakes around the time of the M 9.0 Tohoku-oki earthquake

Here, the temporal evolution of intermediate-depth seismicity before and after the 2011 M9.0 Tohoku-oki earthquake reveals interactions between plate interface slip and deformation in the subducting slab. We investigate seismicity rate changes in the upper and lower planes of the double seismic zone beneath northeast Japan using both a β statistic approach and a temporal epidemic type aftershock sequence model. We do not observe an anomalous precursory increase in intermediate-depth earthquake activity preceding the mainshock; however, following the mainshock, we observe a rate increase for the intermediate-depth earthquakes in the upper plane. The average ratio of upper plane to lower plane activity and the mean deep aseismic slip rate both increased by factor of 2. An increase of downdip compression in the slab resulting from coseismic and postseismic deformation enhanced seismicity in the upper plane, which is dominated by events accommodating downdip shortening from plate unbending.

58 GEOSCIENCES↗
DOE OSTI2017

Bounding the moment deficit rate on crustal faults using geodetic data: Methods

Here, the geodetically derived interseismic moment deficit rate (MDR) provides a first-order constraint on earthquake potential and can play an important role in seismic hazard assessment, but quantifying uncertainty in MDR is a challenging problem that has not been fully addressed. We establish criteria for reliable MDR estimators, evaluate existing methods for determining the probability density of MDR, and propose and evaluate new methods. Geodetic measurements moderately far from the fault provide tighter constraints on MDR than those nearby. Previously used methods can fail catastrophically under predictable circumstances. The bootstrap method works well with strong data constraints on MDR, but can be strongly biased when network geometry is poor. We propose two new methods: the Constrained Optimization Bounding Estimator (COBE) assumes uniform priors on slip rate (from geologic information) and MDR, and can be shown through synthetic tests to be a useful, albeit conservative estimator; the Constrained Optimization Bounding Linear Estimator (COBLE) is the corresponding linear estimator with Gaussian priors rather than point-wise bounds on slip rates. COBE matches COBLE with strong data constraints on MDR. We compare results from COBE and COBLE to previously published results for the interseismic MDR at Parkfield, on the San Andreas Fault, and find similar results; thus, the apparent discrepancy between MDR and the total moment release (seismic and afterslip) in the 2004 Parkfield earthquake remains.

58 GEOSCIENCES↗
DOE OSTI2017

Dominant Controls of Downdip Afterslip and Viscous Relaxation on the Postseismic Displacements Following the M w 7.9 Gorkha, Nepal, Earthquake

Here, we analyze three-dimensional GPS coordinate time series from continuously operating stations in Nepal and South Tibet and calculate the initial 1 year postseismic displacements. We first investigate models of poroelastic rebound, afterslip, and viscoelastic relaxation individually and then attempt to resolve the trade-offs between their contributions by evaluating the misfit between observed and simulated displacements. We compare kinematic inversions for distributed afterslip with stress-driven afterslip models. The modeling results show that no single mechanism satisfactorily explains near- and far-field postseismic deformation following the Gorkha earthquake. When considering contributions from all three mechanisms, we favor a combination of viscoelastic relaxation and afterslip alone, as poroelastic rebound always worsens the misfit. The combined model does not improve the data misfit significantly, but the inverted afterslip distribution is more physically plausible. The inverted afterslip favors slip within the brittle-ductile transition zone downdip of the coseismic rupture and fills the small gap between the mainshock and largest aftershock slip zone, releasing only 7% of the coseismic moment. Our preferred model also illuminates the laterally heterogeneous rheological structure between India and the South Tibet. The transient and steady state viscosities of the upper mantle beneath Tibet are constrained to be greater than 10 18 Pa s and 10 19 Pa s, whereas the Indian upper mantle has a high viscosity ≥10 20 Pa s. The viscosity in the lower crust of southern Tibet shows a clear trade-off with its southward extent and thickness, suggesting an upper bound value of ~8 × 10 19 Pa s for its steady state viscosity.

58 GEOSCIENCES↗
DOE OSTI2018

Spatial and Temporal Variations in the Moment Tensor Solutions of the 2008 Wenchuan Earthquake Aftershocks and Their Tectonic Implications

Following the mainshock of the 2008 M8 Wenchuan Earthquake, there were more than 300 M L ≥ 4.0 aftershocks that occurred between 12 May 2008 and 8 September 2010. We analyzed here the broadband waveforms for these events and found 160 events with sufficient signal-to-noise levels to invert for seismic moment tensors. Considering the length of the activated fault and the distances to the recording stations, four velocity models were employed to account for variability in crustal structure. The moment tensor solutions show considerable variations with a mixture of mainly reverse and strike-slip mechanisms and a small number of normal events and ambiguous events. We analyzed the spatial and temporal distribution of the aftershocks and their mechanism types to characterize the structure and the deformation occurring in the Longmen Shan fold and thrust belt. Our results suggest that the stress is very complex at the Longmen Shan fault zone. The moment tensors have both a spatial segmentation with two major categories of the moment tensor of thrust and strike slip; and a temporal pattern that the majority of the aftershocks gradually migrated to thrust-type events. The variability of aftershock mechanisms is a strong indication of significant tectonic release and stress reorganization that activated numerous small faults in the system.

58 GEOSCIENCES↗
DOE OSTI2025

Dislocation‐Driven Formation of Oriented Macroperiodic Metastructures of Curved Single Crystal Lattices in Glass

Abstract Single crystals fabricated in glass by localized heating can develop uniquely deformed lattices stabilized by the surrounding amorphous medium. The development of lattice curvature appears to be intrinsic to the crystal growth process in some systems, while the result of the locally changing crystallography in others. In this work, a model laser‐fabricated rotating lattice Sb 2 S 3 crystal grown in stoichiometric glass is used to demonstrate fabrication of novel macroperiodic metastructures that utilize intrinsic lattice curvature superimposed with subtle crystallographic influences. The limited availability of slip systems drives the lattice curvature magnitude to vary with crystal growth direction, maximizing for lattices aligned with the predominant Burgers vector along with corresponding increases in dislocation density. Misaligned lattice orientations form smaller secondary lattice curvatures arising from misaligned Burgers vectors with further elastic contributions. Over extended crystal growth, these secondary components align the lattice to rotate about either the <001> or <010> crystal axes forming repeating metastructures of lattice orientation with periodicity 20–160 microns in length. The mechanistic approach used in this work may be expanded to other systems with known slip systems to better understand and design macroperiodic metastructures.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Revisiting Johnson and Jackson boundary conditions for granular flows

In this article, we revisit Johnson and Jackson boundary conditions for granular flows. The oblique collision between a particle and a flat wall is analyzed by adopting the classic rigid-body theory and a more realistic semianalytical model. Based on the kinetic granular theory, the input parameter for the partial-slip boundary conditions, specularity coefficient, which is not measurable in experiments, is then interpreted as a function of the particle-wall restitution coefficient, the frictional coefficient, and the normalized slip velocity at the wall. An analytical expression for the specularity coefficient is suggested for a flat, frictional surface with a low frictional coefficient. The procedure for determining the specularity coefficient for a more general problem is outlined, and a working approximation is provided.

boundary condition↗
DOE OSTI2023

Simulation-based characterization of the variability of earthquake risk to buildings in the near-field

Recent advancements in high performance computing platforms and computational workflow for regional-scale simulations are enabling unprecedented modeling of fault-to-structure earthquake processes. Regional simulations resolving ground motions at frequencies relevant to engineered systems are becoming computationally viable and provide a new capability to improve understanding of the geographical distribution and intensity of risk to buildings and critical infrastructure. As computational capabilities advance, it is essential to move beyond illustrative single rupture realizations for scenario earthquake events towards the development of a full suite of rupture realizations that appropriately characterize the range of risk to building systems. The work described in this article investigates the application of a suite of fault rupture realizations with the objective of assessing near-fault, site-specific seismic demand variability for building structures. A representative high-performance regional-scale computational model is utilized to execute ground motion and building response simulations based on 18 kinematic rupture realizations of an M7 strike-slip scenario earthquake. The fault rupture models for the scenario earthquake are created by systematically perturbing the hypocenter location and stochastically generating rupture parameters (slip, rise time, rake angle) to represent a breadth of ground motion intensities resulting from the spatial and temporal variabilities of an earthquake rupture process. The resulting seismic demand variability for three-story (short period) and forty-story (long period) steel moment-resisting frame buildings is characterized in terms of the median and distribution of peak inter-story drift ratio for a range of near-fault sites. The full suite of 18 fault rupture realizations and approximately 280,000 nonlinear dynamic building simulations indicate that the three-story building undergoes higher median seismic demand and significantly greater variability of demand at a given site than the forty-story building, which has important implications for the level of certainty in predicting building performance during an earthquake. The simulations performed provide deeper insight into the relationship between fault rupture parameterization and building response, which is essential information for developing a representative suite of rupture realizations for specific earthquake scenarios.

58 GEOSCIENCES↗
DOE OSTI2015

Thermal and capillary effects on the caprock mechanical stability at In Salah, Algeria

Thermo-mechanical effects are important in geologic carbon storage because CO 2 will generally reach the storage formation colder than the rock, inducing thermal stresses. Capillary functions, i.e., retention and relative permeability curves, control the CO 2 plume shape, which may affect overpressure and thus, caprock stability. To analyze these thermal and capillary effects, we numerically solve non-isothermal injection of CO 2 in deformable porous media considering the In Salah, Algeria, CO 2 storage site. Here, we find that changes in the capillary functions have a negligible effect on overpressure and thus, caprock stability is not affected by capillary effects. But, we show that for the strike slip stress regime prevalent at In Salah, stability decreases in the lowest parts of the caprock during injection due to cooling-induced thermal stresses. Simulations show that shear slip along pre-existing fractures may take place in the cooled region, whereas tensile failure is less likely to occur. Indeed, only the injection zone and the lowest tens of meters of the 900-m-thick caprock at In Salah might be affected by cooling effects, which would thus not jeopardize the overall sealing capacity of the caprock. Furthermore, faults are likely to remain stable far away from the injection well because outside the cooled region the injection-induced stress changes are not sufficient to exceed the anticipated shear strength of minor faults. Nonetheless, we recommend that thermal effects should be considered in the site characterization and injection design of future CO 2 injection sites to assess caprock stability and guarantee a permanent CO 2 storage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Tailoring the mechanical properties of nanoparticle networks that encompass biomimetic catch bonds

Biological “catch” bonds display the distinctive attribute that the bond lifetime can increase under an applied force. Insertion of biomimetic catch bonds into hybrid materials could lead to composites that exhibit remarkable mechanical properties. We model the tensile behavior of polymer-grafted nanoparticle (PGN) networks interconnected by a mixture of catch bonds and conventional “slip” bonds, whose lifetimes decrease with force. We formulate a kinetic master equation that provides the complete probabilistic description of a system of two PGNs. The master equation is used to analyze the parameter space that determines the rupture behavior of the catch bonds in the two-particle system. We then utilize two exemplary sets of catch bond parameters in three-dimensional computer simulations of larger PGN networks under strain-controlled tensile deformation. We demonstrate that the strain at break and toughness of the networks can be altered by “tuning” the attributes of the catch bonds and varying the fraction of catch bonds in the network. We show that networks encompassing the catch bonds could exhibit a several-fold increase in the strain-at-break and toughness relative to those interconnected solely by the slip bonds. Our studies can provide valuable guidelines for tailoring the mechanical properties of novel, bio-inspired nanocomposites.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Numerical Simulation of M9 Megathrust Earthquakes in the Cascadia Subduction Zone

We estimate ground motions in the Pacific Northwest urban areas during M9 subduction scenario earthquakes on the Cascadia megathrust by simulating wave propagation from an ensemble of kinematic source descriptions. Velocities and densities in our computational mesh are defined by integrating the regional Cascadia Community Velocity Model (CVM) v1.6 (Stephenson et al. P-and S-wave velocity models incorporating the Cascadia subduction zone for 3D earthquake ground motion simulations—update for open-file report 2007–1348, US Geological Survey, 2017) including the ocean water layer with a local velocity model of the Georgia basin (Molnar, Predicting earthquake ground shaking due to 1D soil layering and 3D basin structure in SW British Columbia, Canada, 2011), including additional near-surface velocity information. We generate six source realizations, each consisting of a background slip distribution with correlation lengths, rise times and rupture velocities consistent with data from previous megathrust earthquakes (e.g., 2011 M 9 Tohoku or 2010 M 8.8 Maule). We then superimpose M ~ 8 subevents, characterized by short rise times and high stress drops on the background slip model to mimic high-frequency strong ground motion generation areas in the deeper portion of the rupture (Frankel, Bull Seismol Soc Am 107(1):372–386, 2017). The wave propagation is simulated using the discontinuous mesh (DM) version of the AWP finite difference code. We simulate frequencies up to 1.25 Hz, using a spatial discretization of 100 m in the fine grid, resulting in surface grid dimensions of 6540 × 10,728 mesh points. At depths below 8 km, the grid step increases to 300 m. We obtain stable and accurate results for the DM method throughout the simulation time of 7.5 min as verified against a solution obtained with a uniform 100 m grid spacing. Peak ground velocities (PGVs) range between 0.57 and 1.0 m/s in downtown Seattle and between 0.25 and 0.54 m/s in downtown Vancouver, while spectral accelerations at 2 s range between 1.7 and 3.6 m/s 2 and 1.0 and 1.3 m/s 2 , respectively. These long-period ground motions are not significantly reduced if plastic Drucker-Prager yielding in shallow cohesionless sediments is taken into account. Effects of rupture directivity are significant at periods of ~ 10 s, but almost absent at shorter periods. We find that increasing the depth extent of the subducting slab from the truncation at 60 km in the Cascadia CVM version 1.6 to ~ 100 km increases the PGVs by 15% in Seattle and by 40% in Vancouver.

finite difference method↗
DOE OSTI2018

The Ideal Strengths of Superconducting MgCNi{sub 3} and CdCNi{sub 3}

The stress-strain curves under tensile deformation in the 〈100〉, 〈110〉, and 〈111〉 directions and under shear deformation in the (001)〈110〉, (110)〈1̄10〉, (111)〈11̄0〉, and (111)〈112̄〉 slip systems have been systematically calculated by first-principles method to study the ideal strengths of superconducting MgCNi{sub 3} and CdCNi{sub 3}. The ideal strengths in the three tensile directions are found to be reduced in the order of 〈100〉 → 〈110〉 → 〈111〉 and those for the four shear slip systems in the order of (110)〈1̄10〉→(111)〈112̄〉→(111)〈11̄0〉→(001)〈110〉 for both superconductors. Their lowest ideal tensile strengths are found to be larger than the corresponding highest ideal shear strengths, which could explain why both superconductors have the ductility. The obtained lattice constants and elastic properties coincide well with the the available experimental and theoretical values.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Predicting Fluid Flow Regime, Permeability, and Diffusivity in Mudrocks from Multiscale Pore Characterisation

In geoenergy applications, mudrocks prevent fluids to leak from temporary (H 2 , CH 4 ) or permanent (CO 2 , radioactive waste) storage/disposal sites and serve as a source and reservoir for unconventional oil and gas. Understanding transport properties integrated with dominant fluid flow mechanisms in mudrocks is essential to better predict the performance of mudrocks within these applications. In this study, small-angle neutron scattering (SANS) experiments were conducted on 71 samples from 13 different sets of mudrocks across the globe to capture the pore structure of nearly the full pore size spectrum (2 nm–5 μm). We develop fractal models to predict transport properties (permeability and diffusivity) based on the SANS-derived pore size distributions. The results indicate that transport phenomena in mudrocks are intrinsically pore size-dependent. Depending on hydrostatic pore pressures, transition flow develops in micropores, slip flow in meso- and macropores, and continuum flow in larger macropores. Fluid flow regimes progress towards larger pore sizes during reservoir depletion or smaller pore sizes during fluid storage, so when pressure is decreased or increased, respectively. Capturing the heterogeneity of mudrocks by considering fractal dimension and tortuosity fractal dimension for defined pore size ranges, fractal models integrate apparent permeability with slip flow, Darcy permeability with continuum flow, and gas diffusivity with diffusion flow in the matrix. This new model of pore size-dependent transport and integrated transport properties using fractal models yields a systematic approach that can also inform multiscale multi-physics models to better understand fluid flow and transport phenomena in mudrocks on the reservoir and basin scale.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microstructure-Based Estimation of Strength and Ductility Distributions for \(\alpha +\beta \) Titanium Alloys

Titanium alloys are processed to develop a wide range of microstructure configurations and therefore material properties. While these properties are typically measured experimentally, a framework for property prediction could greatly enhance alloy design and manufacturing. In this work, a microstructure-sensitive framework is presented for the prediction of strength and ductility as well as estimates of the bounds in variability for these properties. The framework explicitly considers distributions of microstructure via new approaches for instantiation of structure in synthetic samples. The parametric evaluation strategy, including the finite element simulation package FEpX, is used to create and test virtual polycrystalline samples to evaluate the variability bounds of mechanical properties in Ti-6Al-4V. Critical parameters for the property evaluation framework are provided by measurements of single crystal properties and advanced characterization of microstructure and slip system strengths in 2D and 3D. Property distributions for yield strength and ductility are presented, along with the validation and verification steps undertaken. Comparisons between strain localization and slip activity in virtual samples and in experimental grain-scale strain measurements are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accumulative Roll Bonding of Alloy 2205 Duplex Steel and the Accompanying Impacts on Microstructure, Texture, and Mechanical Properties

The mechanical and microstructural evolution of Alloy 2205 during severe plastic deformation is examined in this study. A combination of accumulative roll bonding (ARB) and cold rolling results in the successful formation of a nanograined dual-phase microstructure of austenite and ferrite with some transformed martensite. Severe deformation to cumulative reductions of 80.5, 92.5, 95, and 97 pct were performed. Microscopy indicates that grain dimensions in the sheet normal direction is less than 100 nm for reductions ≥ 92.5 pct. Shear banding is observed at reductions ≥ 95 pct while twinning is only observed at reductions < 92.5 pct. Neutron diffraction measurements indicated the presence of martensite for reductions ≥ 95 pct at ~ 8 pct volume fraction. Taken in conjunction, it appears that during initial ARB processing, both slip and twinning are active plastic mechanisms. As twinning becomes exhausted, martensitic transformation, slip, and intermittent shear banding account for the active plasticity mechanisms. Material hardness saturates at 92.5 pct reduction, with a maximum hardness of 45 HRC. Sub-sized tensile testing confirms this approximate hardness with measurements indicating a UTS of ~ 1440 MPa. Texture analysis of crystal orientation distributions in the plate normal direction suggest an approximate Kurdjumov–Sachs orientation relationship at all reductions above 80 pct indicating stability of the orientation relationship at high strains. The intragranular structure develops a fine scale sub-grain content with increasing deformation, resulting in a continual evolution of texture up to and including 97 pct reduction. The final structure presents strong components of Goss and rotated cube texture in both the austenite and ferrite. In this body of work we aim to compare ARB of an industrially relevant FCC/BCC system (Alloy 2205) to historical model FCC/BCC systems such as Cu/Nb.

36 MATERIALS SCIENCE↗