Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fault 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 271 records · Page 15

Mixing effects on geothermometric calculations of the Newdale geothermal area in the Eastern Snake River Plain, Idaho

The Newdale geothermal area in Madison and Fremont Counties in Idaho is a known geothermal resource area whose thermal anomaly is expressed by high thermal gradients and numerous wells producing warm water (up to 51 °C). Geologically, the Newdale geothermal area is located within the Eastern Snake River Plain (ESRP) that has a time-transgressive history of sustained volcanic activities associated with the passage of Yellowstone Hotspot from the southwestern part of Idaho to its current position underneath Yellowstone National Park in Wyoming. Locally, the Newdale geothermal area is located within an area that was subjected to several overlapping and nested caldera complexes. The Tertiary caldera forming volcanic activities and associated rocks have been buried underneath Quaternary flood basalts and felsic volcanic rocks. Two southeast dipping young faults (Teton dam fault and an unnamed fault) in the area provide the structural control for this localized thermal anomaly zone. Geochemically, water samples from numerous wells in the area can be divided into two broad groups – Na-HCO3 and Ca-(Mg)-HCO3 type waters and are considered to be the product of water-rhyolite and water-basalt interactions, respectively. Each type of water can further be subdivided into two groups depending on their degree of mixing with other water types or interaction with other rocks. For example, some bivariate plots indicate that some Ca-(Mg)-HCO3 water samples have interacted only with basalts whereas some samples of this water type also show limited interaction with rhyolite or mixing with Na-HCO3 type water. Traditional geothermometers [e.g., silica variants, Na-K-Ca (Mg-corrected)] indicate lower temperatures for this area; however, a traditional silica-enthalpy mixing model results in higher reservoir temperatures. We applied a new multicomponent equilibrium geothermometry tool (e.g., Reservoir Temperature Estimator, RTEst) that is based on inverse geochemical modeling which explicitly accounts for boiling, mixing, and CO2 degassing. RTEst modeling results indicate that the well water samples are mixed with up to 75% of the near surface groundwater. Relatively, Ca-(Mg)-HCO3 type water samples are more diluted than the Na-HCO3 type water samples. However, both water types result in similar reservoir temperatures, up to 150 °C. Samples in the vicinity of faults produced higher reservoir temperatures than samples away from the faults. Although both the silica-enthalpy mixing and RTEst models indicated promising geothermal reservoir temperatures, evaluation of the subsurface permeability and extent of the thermal anomaly is needed to define the hydrothermal potential of the Newdale geothermal resource.

Eastern Snake River Plain, Newdale, RTEst, geother↗

Cenozoic thermal, mechanical and tectonic evolution of the Rio Grande rift

Two areas of New Mexico which exhibit complex but similar Cenozoic histories of extensional tectonism are analyzed. The first study area is the Basin and Range province and southern Rio Gande rift in southern New Mexico; the second study area is the central Rio Grande rift in central and northern New Mexico, the southern San Luis basin, the Espanola basin, and the Albuquerque basin. Two phases of extension were identified: the first phase which began in mid-Oligocene was characterized by local high-strain extension events, low-angle faulting, and the development of broad, shallow basins, all indicating an approximately NE-SW + or - 25 deg extension direction; the later phase which occurred primarily in the late Miocene, was characterized by synchronous, high-angle faulting, resulting in large vertical strains which produced the modern Rio Grande rift morphology. Extension direction was approximately E-W. Geotherms were estimated and lithospheric strength curves were calculated for these two phases of extension. A high geotherm was deduced for the early phase resulting in a shallow crustal brittle transition, and insignificant mantle strength. The lithosphere subsequently cooled, resulting in a significant zone of mantle strength beneath the Moho. It is concluded that the interrelationship among regional and local prerifting, synrifting, and postrifting events in the Rio Grande rift attests to the fact that the rifting (in the region studied) should be considered in the context of other geologic events.

Morgan, P.↗

Self-stabilizing byzantine-fault-tolerant clock synchronization system and method

Systems and methods for rapid Byzantine-fault-tolerant self-stabilizing clock synchronization are provided. The systems and methods are based on a protocol comprising a state machine and a set of monitors that execute once every local oscillator tick. The protocol is independent of specific application specific requirements. The faults are assumed to be arbitrary and/or malicious. All timing measures of variables are based on the node's local clock and thus no central clock or externally generated pulse is used. Instances of the protocol are shown to tolerate bursts of transient failures and deterministically converge with a linear convergence time with respect to the synchronization period as predicted.

Malekpour, Mahyar R.↗

Detection of High-impedance Arcing Faults in Radial Distribution DC Systems

High voltage, low current arcing faults in DC power systems have been researched at the NASA Glenn Research Center in order to develop a method for detecting these 'hidden faults', in-situ, before damage to cables and components from localized heating can occur. A simple arc generator was built and high-speed and low-speed monitoring of the voltage and current waveforms, respectively, has shown that these high impedance faults produce a significant increase in high frequency content in the DC bus voltage and low frequency content in the DC system current. Based on these observations, an algorithm was developed using a high-speed data acquisition system that was able to accurately detect high impedance arcing events induced in a single-line system based on the frequency content of the DC bus voltage or the system current. Next, a multi-line, radial distribution system was researched to see if the arc location could be determined through the voltage information when multiple 'detectors' are present in the system. It was shown that a small, passive LC filter was sufficient to reliably isolate the fault to a single line in a multi-line distribution system. Of course, no modification is necessary if only the current information is used to locate the arc. However, data shows that it might be necessary to monitor both the system current and bus voltage to improve the chances of detecting and locating high impedance arcing faults

Gonzalez, Marcelo C.↗

Compartmentalization of the Coso East Flank geothermal field imaged by 3-D full-tensor MT inversion

Previous magnetotelluric (MT) studies of the high-temperature Coso geothermal system in California identified a subvertical feature of lowresistivity (2-5Ohmm) and appreciable lateral extent (>1 km) in the producing zone of the East Flank field. However, these models could not reproduce gross 3-D effects in the recorded data. We perform 3-D full-tensor inversion and retrieve a resistivity model that out-performs previous 2-D and 3-D off-diagonal models in terms of its fit to the complete 3-D MT data set as well as the degree of modelling bias. Inclusion of secondary Zxxand Zyydata components leads to a robust east-dip (60°) to the previously identified conductive East Flank reservoir feature, which correlates strongly with recently mapped surface faults, downhole well temperatures, 3-D seismic reflection data, and local microseismicity. We perform synthetic forward modelling to test the best-fit dip of this conductor using the response at a nearby MT station. We interpret the dipping conductor as a fractured and fluidized compartment, which is structurally controlled by an unmapped blind East Flank fault zone.

58 GEOSCIENCES↗

Numerical modeling of heterogeneous asperity distributions controlling the growth of shear rupture on a frictional fault

Copyright 2016 ARMA, American Rock Mechanics Association. A better understanding of the strength distributions and the seismic response of asperities failing within a cohesive (or breakdown) region of the growing shear rupture is important in many engineering disciplines, such as, hydraulic fracture and induced seismicity. In the laboratory, a fault is experimentally modeled by a 400 mm × 80 mm × 12 mm PMMA slider sheared across a large PMMA base plate. A pressure-sensitive film was employed to localize and size, and measure normal stress upon contact junctions. Moments prior to rapid sliding, multiple discrete and localized dynamic events (or foreshocks) were observed at a section of the fault that exhibited a sparse to dense transition in asperity distribution. Experiments showed that foreshocks occurred as the crack (rupture front) moved into the resistive patch of densely distributed asperities. We use the results of pressure sensitive film asperity measurements to develop a heterogeneous finite element (FE) model that studies the effect of asperities within the cohesive region have on the expanding shear crack. Slip accumulated heterogeneously during the loading cycle due to the non-uniform distribution of asperities. We study the constitutive behavior of a frictional interface where local slip is determined by the asperities and their ability to communicate elastically. A foreshock (Mw -7.4) was simulated in a quasi-static manner by instantaneoulsy removing an asperity contact during the loading cycle. The foreshock increased shear stress perturbations to neighboring regions up to distances of ~ 10 to 15 source radii.

Selvadurai, PA↗

Fault Scarp Detection Beneath Dense Vegetation Cover: Airborne Lidar Mapping of the Seattle Fault Zone, Bainbridge Island, Washington State

The emergence of a commercial airborne laser mapping industry is paying major dividends in an assessment of earthquake hazards in the Puget Lowland of Washington State. Geophysical observations and historical seismicity indicate the presence of active upper-crustal faults in the Puget Lowland, placing the major population centers of Seattle and Tacoma at significant risk. However, until recently the surface trace of these faults had never been identified, neither on the ground nor from remote sensing, due to cover by the dense vegetation of the Pacific Northwest temperate rainforests and extremely thick Pleistocene glacial deposits. A pilot lidar mapping project of Bainbridge Island in the Puget Sound, contracted by the Kitsap Public Utility District (KPUD) and conducted by Airborne Laser Mapping in late 1996, spectacularly revealed geomorphic features associated with fault strands within the Seattle fault zone. The features include a previously unrecognized fault scarp, an uplifted marine wave-cut platform, and tilted sedimentary strata. The United States Geologic Survey (USGS) is now conducting trenching studies across the fault scarp to establish ages, displacements, and recurrence intervals of recent earthquakes on this active fault. The success of this pilot study has inspired the formation of a consortium of federal and local organizations to extend this work to a 2350 square kilometer (580,000 acre) region of the Puget Lowland, covering nearly the entire extent (approx. 85 km) of the Seattle fault. The consortium includes NASA, the USGS, and four local groups consisting of KPUD, Kitsap County, the City of Seattle, and the Puget Sound Regional Council (PSRC). The consortium has selected Terrapoint, a commercial lidar mapping vendor, to acquire the data.

Harding, David J.↗

Structural review of the Vredefort dome

The structure of the older-than-3.2-Ga Archean basement and Archean-to-Precambrian sedimentary/volcanic rocks (3.07 to ca. 2.2 Ga) in the center of the Witwatersrand Basin to the southwest of Johannesburg (South Africa) is dominated by the ca. 2.0-Ga megascopic Vredefort 'Dome' structure. The effect of the 'Vredefort event' is demonstrably large and is evident within a northerly arc of about 100 km radius around the granitic core of the structure. Northerly asymmetric overturning of the strata is observed within the first 17 km (strata is horizontal in the south), followed by a 40-km-wide rim synclinorium. Fold and fault structures (normal, reverse, and strike-slip) are locally as well as regionally concentrically arranged with respect to the northern and western sides of the structure. The unusual category of brittle deformation, the so-called 'shock deformation', observed in the collar strata has attracted worldwide attention over the past two decades. These deformation phenomena include the presence of coesite and stishovite, mylonites, and pseudotachylites, cataclasis at a microscopic scale, and the ubiquitous development of multiply striated joint surfaces (which include shatter cones, orthogonal, curviplanar, and conjugate fractures). The macroscopic to microscopic deformation features have led to the formulation of various hypotheses to account for the origin of the Vredefort structure: (1) tectonic hypotheses--deep crustal shear model, doming and N-directed thrust fault model, fold interference model, and diapir model; (2) the exogenous bolide impact hypothesis; and (3) the endogenous cryptoexplosion model.

Colliston, W. P.↗

Fault geometries and extension in the Valles Marineris, Mars

The central troughs of the Valles Marineris system on Mars are generally regarded to be of tectonic origin because they are linear and bounded by faults. Their radial orientation to the Tharsis rise suggests that at least their incipient formation was influenced by the formation of the rise, but the exact role of Tharsis is not yet well understood. Collapse and erosion may have significantly modified the original tectonic troughs. To better understand the magnitude and mechanism of crustal extension across the Valles Marineris, the attitudes of trough-wall faults that parallel the long axes of the troughs were studied. Our results will help to determine if the extension is related to the Tharsis rise, to a local arch along which the troughs lie, or to other mechanisms. Fault-plane attitudes were obtained by performing a photogrammetric study of fault traces in the walls of Candor and Melas Chasmata. If the exposed wall faults are not major trough-bounding faults, it is assumed that these major faults are buried at the base of the walls and have similar geometries to the exposed faults. Three dimensional coordinates of several points along each fault trace were obtained and three-point solutions to their geometries were constructed. As expected, the results show that the faults dip toward the interior of the troughs, suggesting that they are indeed normal faults associated with trough formation.

Chadwick, D. J.↗

Topographic Slant Range Modeling and Fault Detection for Precision Planetary Landing

This work presents a novel landing site relative topographic measurement model for aslant range sensor being utilized for precision planetary landing operations. The measurement model accounts for the local terrain the slant range sensor captures and leverages knowledge of the estimated landing site provided by the navigation filter. Notably, in contrast to previous works, the new model does not rely on surface normal approximation, reducing the model’s sensitivity to noisy digital elevation maps which represent the local topography. In addition to the measurement model, this work introduces a novel fault detection method, denoted the probabilistic inspection of topographic filter altitude likelihood (PITFAL) algorithm, that implements a statistical outlier rejection algorithm. PITFAL is designed for multi-beam slant range sensors, such as the Navigation Doppler LIDAR (NDL), and identifies statistically inconsistent range estimates through a consensus check on the set of apparent altitudes computed for each individual beam. These models are numerically validated by the Safe and Precise Landing Capability Evolution (SPLICE) project’s high-fidelity terrestrial and lunar lander simulations.

Davis W Adams↗

Topographic Slant Range Modeling and Fault Detection for Precision Planetary Landing

This work presents a novel landing site relative topographic measurement model for aslant range sensor being utilized for precision planetary landing operations. The measurement model accounts for the local terrain the slant range sensor captures and leverages knowledge of the estimated landing site provided by the navigation filter. Notably, in contrast to previous works, the new model does not rely on surface normal approximation, reducing the model’s sensitivity to noisy digital elevation maps which represent the local topography. In addition to the measurement model, this work introduces a novel fault detection method, denoted the probabilistic inspection of topographic filter altitude likelihood (PITFAL) algorithm, that implements a statistical outlier rejection algorithm. PITFAL is designed for multi-beam slant range sensors, such as the Navigation Doppler LIDAR (NDL), and identifies statistically inconsistent range estimates through a consensus check on the set of apparent altitudes computed for each individual beam. These models are numerically validated by the Safe and Precise Landing Capability Evolution (SPLICE) project’s high-fidelity terrestrial and lunar lander simulations.

Davis W Adams↗

A novel vibration-based fault diagnostic algorithm for gearboxes under speed fluctuations without rotational speed measurement

The localized failures of gears introduce cyclic-transient impulses in the measured gearbox vibration signals. These impulses are usually identified from the sidebands around gear-mesh harmonics through the spectral analysis of cyclo-stationary signals. However, in practice, several high-powered applications of gearboxes like wind turbines are intrinsically characterized by nonstationary processes that blur the measured vibration spectra of a gearbox and deteriorate the efficacy of spectral diagnostic methods. Although order-tracking techniques have been proposed to improve the performance of spectral diagnosis for nonstationary signals measured in such applications, the required hardware for the measurement of rotational speed of these machines is often unavailable in industrial settings. Moreover, existing tacho-less order-tracking approaches are usually limited by the high time-frequency resolution requirement, which is a prerequisite for the precise estimation of the instantaneous frequency. To address such issues, a novel fault-signature enhancement algorithm is proposed that can alleviate the spectral smearing without the need of rotational speed measurement. This proposed tacho-less diagnostic technique resamples the measured acceleration signal of the gearbox based on the optimal warping path evaluated from the fast dynamic time-warping algorithm, which aligns a filtered shaft rotational harmonic signal with respect to a reference signal assuming a constant shaft rotational speed estimated from the approximation of operational speed. The effectiveness of this method is validated using both simulated signals from a fixed-axis gear pair under nonstationary conditions and experimental measurements from a 750-kW planetary wind turbine gearbox on a dynamometer test rig. Lastly, the results demonstrate that the proposed algorithm can identify fault information from typical gearbox vibration measurements carried out in a resource-constrained industrial environment.

17 WIND ENERGY↗

Temperature uncertainty modelling with proxy structural data as geostatistical constraints for well siting: an example applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural geology information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aims of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations. The data used to support the modelling are measured temperatures and structural proxies from 3D geologic modelling (distance to fault, distance to fault intersections and terminations, Coulomb stress change and dilation tendency), the latter considered ‘secondary’ data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally-likely temperature models (i.e. realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: V prior and entropy. V prior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because of the dipping ellipsoid shape created and that the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e. linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g. dilation v. distance to fault termination) should be used given the specific geothermal system. In general, V prior from local varying mean results identify locations that are close to high values for the structural proxies: areas with higher probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comparison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

15 GEOTHERMAL ENERGY↗

Internal structure of the San Jacinto fault zone at Jackass Flat from data recorded by a dense linear array

The internal structure of the Clark fault in the trifurcation area of the San Jacinto fault zone is imaged using seismograms recorded by a dense linear array (Jackass Flat, JF) crossing the surface trace of the fault and an adjacent array (TR) to the SW. Delay times between phase arrivals associated with ~3500 local earthquakes and nine teleseismic events are used to estimate velocity variations within the arrays. The teleseismic P waves travel faster beneath the TR than the JF array, in contrast to larger scale tomographic results. Statistical analysis of local P-wave delay times indicates that the entire JF array, with an aperture of ~400 m, is inside a low-velocity damage zone. This low-velocity zone is bounded on the NE side by a shallow bimaterial interface generating fault zone head waves, and it contains an inner zone of more intense damage generating fault zone trapped waves. The P-wave velocity contrast across the local bounding bimaterial interface is 10–15 per cent. The trapping structure is associated with a width of ~200 m, S-wave velocity reduction of ~35 per cent with respect to the surrounding rock, Q-value of ~20 and depth of ~3.5 km. The imaging results suggest that the main seismogenic fault is near the SW end of the JF array, in agreement with a prominent geomorphologic feature. The existence of intense local damage on the crustal block with faster larger scale velocity at depth is consistent with common propagation of earthquake ruptures in the area to the NW.

58 GEOSCIENCES↗

MINERALOGICAL AND TEXTURAL CONTROLS ON SHEAR STRENGTH, SLIP STABILITY AND PERMEABILITY OF FAULTS

Induced seismicity resulting from fluid injection into the subsurface related to water and CO2 disposal, hydraulic fracturing and the stimulation of geothermal reservoirs present an important societal concern. These human activities involve the injection of large volumes of pressurized fluid into the subsurface, potentially at high rates, raising local pore pressures and disturbing the pristine local stress regime by lowering effective normal stress on pre-existing faults and fractures. The reduction of effective normal stress may trigger fault/fracture reactivation and in some cases result in hazardous seismic ruptures. Effective management and engineering of anthropogenic seismic events requires substantial understanding in the mechanisms, especially for the controlling factors on coupled rheological and transport response, including fault shear strength, slip stability, and permeability evolution during such events. In this study, we explore the coupled rheological and transport response of faults and fractures during reactivation as controlled by two fundamental controlling properties, viz., mineralogy and textural features. We approach this problem through shear experiments on analog faults and fractures via laboratory and numerical experiments. Specifically, we investigate: (1) the influence of frictionally weak minerals (talc) in mixtures of mineral analogs featuring contrasting frictional properties, (2) the influence of iron oxide grain coatings on quartz aggregates, and (3) the influence of fracture roughness in mated fractures on the ensemble shear strength, slip stability, and permeability evolution during reactivation events. We address the following questions in this study: (1) how much and what distribution of frictionally weak minerals is required to induce significant weakening in faults consisting of a matrix of frictionally strong minerals, (2) how does a pre-imposed weak mineral layer influence the rheological and transport behavior of faults, (3) what is the influence of a trace amount of grain coating materials introduce on the coupled behavior of faults, and finally (4) how do asperity height and wave length control the ensemble behavior of faults. These questions are explicitly answered in the following.

58 GEOSCIENCES↗

Dynamic modeling of injection-induced fault reactivation and ground motion and impact on surface structures and human perception

We summarize recent modeling studies of injection-induced fault reactivation, seismicity, and its potential impact on surface structures and nuisance to the local human population. We used coupled multiphase fluid flow and geomechanical numerical modeling, dynamic wave propagation modeling, seismology theories, and empirical vibration criteria from mining and construction industries. We first simulated injection-induced fault reactivation, including dynamic fault slip, seismic source, wave propagation, and ground vibrations. From co-seismic average shear displacement and rupture area, we determined the moment magnitude to about M w = 3 for an injection-induced fault reactivation at a depth of about 1000 m. We then analyzed the ground vibration results in terms of peak ground acceleration (PGA), peak ground velocity (PGV), and frequency content, with comparison to the U.S. Bureau of Mines’ vibration criteria for cosmetic damage to buildings, as well as human-perception vibration limits. For the considered synthetic M w = 3 event, our analysis showed that the short duration, high frequency ground motion may not cause any significant damage to surface structures, and would not cause, in this particular case, upward CO 2 leakage, but would certainly be felt by the local population.

58 GEOSCIENCES↗

Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS

New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Dislocation nano-hydrides in nickel: Nucleation, evolution and effects on dislocation behaviors

Nano-hydrides have been predicted to precipitate at the core of edge dislocations in the Ni-H system, a mechanism that may promote hydrogen embrittlement. However, nano-hydride nucleation, growth, and effects on dislocation behavior have seldom been explored. This work combines molecular dynamics grand canonical Monte Carlo (MD-GCMC) simulations and continuum modeling to uncover a wide range of phenomena linked to dislocation nano-hydrides. Simulations reveal that nano-hydrides can be stabilized at dislocation cores with all character angles, including screw segments, due to the hydrostatic stresses around the cores of the Shockley partials. Nano-hydride nucleation takes place in these regions, and growth is dictated by the character angle θ of the perfect dislocation. The equilibrium stacking fault width d eq varies dynamically to increase the local hydrostatic stress field and facilitate the formation of the nano-hydride, forming a constriction-like feature and leading to three distinct behaviors: d eq decreases for θ > 30°, d eq remains unchanged for θ = 30°, and d eq increases for θ < 30°. Remote hydrostatic and Escaig stresses are also shown to influence the nucleation stage, implying stress concentrations such as those ahead of crack tips may facilitate nano-hydride precipitation. Moreover, we identify a new hydrogen-induced 60° dislocation reaction that emits a Shockley partial on a conjugate plane, with potential implications for twin nucleation. Testable predictions from this study are then used to reinterpret previous results from the literature. In conclusion, these findings provide a comprehensive framework to assess nano-hydride formation and evolution at dislocations in nickel and other face-centered cubic metals, with important implications to hydrogen embrittlement.

Dislocations↗