Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “physics of faulting”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Defect-driven selective metal oxidation at atomic scale

Nanoscale materials modified by crystal defects exhibit significantly different behaviours upon chemical reactions such as oxidation, catalysis, lithiation and epitaxial growth. However, unveiling the exact defect-controlled reaction dynamics (e.g. oxidation) at atomic scale remains a challenge for applications. Here, using in situ high-resolution transmission electron microscopy and first-principles calculations, we reveal the dynamics of a general site-selective oxidation behaviour in nanotwinned silver and palladium driven by individual stacking-faults and twin boundaries. The coherent planar defects crossing the surface exhibit the highest oxygen binding energies, leading to preferential nucleation of oxides at these intersections. Planar-fault mediated diffusion of oxygen atoms is shown to catalyse subsequent layer-by-layer inward oxide growth via atomic steps migrating on the oxide-metal interface. These findings provide an atomistic visualization of the complex reaction dynamics controlled by planar defects in metallic nanostructures, which could enable the modification of physiochemical performances in nanomaterials through defect engineering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Desktop Audits: Overview, Case Studies, Guidance

The Energy Independence and Security Act (EISA) of 2007 requires that federal agencies audit 25% of their 3.1 billion square feet of covered facilities each year. Several newer, advanced options such as monitoring-based commissioning, fault detection and diagnostics, and energy management systems can help support data collection and analysis to be used for EISA compliance. The Federal Energy Management Program's (FEMP’s) Audit Decision Tree, Audit Definitions, and FEMP Guidance can be used by Agencies to determine the best approaches to evaluate their facilities. FEMP continues to develop additional resources to support agencies in meeting EISA compliance. This presentation gives an overview of desktop audits, including case studies and guidance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Exploiting Kubernetes to Simplify the Deployment and Management of the Multi-purpose CMS Pilot Job Factory

GlideinWMS, a widely utilized workload management system in high-energy physics (HEP) research, serves as the backbone for efficient job provisioning across distributed computing resources. It is utilized by various experiments and organizations, including CMS, OSG, Dune, and FIFE, to create HTCondor pools as large as 600k cores. In particular, a shared factory service historically deployed at UCSD has been configured to interface with more than 500 routes to compute clusters. As part of our team’s initiative to modernize infrastructure and enhance scalability, we undertook the migration of the GlideinWMS factory service into the Kubernetes environment. Leveraging the flexibility and orchestration capabilities of Kubernetes, we successfully deployed the factory service within the OSG Tiger Kubernetes cluster. The major benefits Kubernetes gives us is it streamlines the management and monitoring of the factory infrastructure, and improves fault tolerance through its resilient deployment strategies. Through this case study, we aim to share insights, challenges, and best practices encountered during the migration process. Our experience underscores the benefits of embracing containerization and Kubernetes orchestration for HEP computing infrastructure, paving the way for scalability and resilience in distributed computing environments.

Dost, Jeffrey Michael [UC, San Diego (main)]↗
DOE OSTI2020

Ambient pressure synthesis and characterization of layered honeycomb Li 2 PdO 3

Single-phase polycrystalline Li 2 PdO 3 has been synthesized at 640°C in oxygen for the first time under ambient pressure. X-ray and neutron diffraction analyses show that the sample possesses a monoclinic layered structure belonging to the C2/m space group. Rietveld refinements of neutron powder diffraction data indicate ~10% Li–Pd site exchange and DIFFaX modelling manifest ~2% stacking faults present within LiPd 2 layers. A band gap of ~2.23 eV was calculated for the golden Li 2 PdO 3 using absorbance measurements. Thermogravimetric analysis of the sample shows that Li 2 PdO 3 is stable up to 730°C under oxygen. Here a Curie tail is observed at low temperature magnetic measurements (T < 50K), yielding an effective moment of 0.038 μ B , possibly due to spin ½ impurities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Preliminary Analysis of Source Physics Experiment Explosion-Triggered Microseismicity Using the Back-Projection Method

A series of four chemical explosions were detonated in a deep borehole within the Yucca Flat Dry Alluvium Geology (DAG) at the Nevada National Security Site between 2018 and 2019. The two larger chemical explosions of 50 tons (DAG-2) and 10 tons (DAG-4) TNT equivalent yield triggered energetic aftershock sequences numbering 1392 and 347 microearthquakes, respectively, within the first 10 days. No significant aftershock activity was observed for the two smaller 1-ton explosions (DAG-1 and DAG-3). Here, we used a back-projection method based on travel-time migration and stacking of signal-to-noise ratio traces to detect, associate and locate aftershocks from a subset of 22-geophones within a larger 2 urn:x-wiley:21699313:media:jgrb54910:jgrb54910-math-0001 2 km seismic array surrounding the borehole. The aftershocks located within 300 m of the borehole and the depths were above the working points of 300 and 50 m depths of DAG-2 and DAG-4, respectively, ruling out triggering slip on geologic faults or disturbances beneath neighboring collapse craters. DAG-2 and DAG-4 aftershocks decayed at similar rates, with power-law exponents of p = 1.48 and p = 1.49, respectively. These decay rates are comparable to aftershocks sequences triggered by earthquakes and historical nuclear explosions at Yucca Flat. A smooth power-law aftershock decay within the first 10 days suggests a triggering mechanism from explosion generated stress relaxation due to the diffusion of high gas pressures in the cavity and radial fractures. A more random and episodic aftershock rate would be expected due to cavity collapse or falling rubble in chimney formation.

58 GEOSCIENCES↗
DOE OSTI2025

Explainable AI for Multivariate Time Series Pattern Exploration: Latent Space Visual Analytics With Temporal Fusion Transformer and Variational Autoencoders in Power Grid Event Diagnosis

Detecting and analyzing complex patterns in multivariate time-series data is crucial for decision-making in urban and environmental system operations. However, challenges arise from the high dimensionality, intricate complexity, and interconnected nature of complex patterns, which hinder the understanding of their underlying physical processes. Existing AI methods often face limitations in interpretability, computational efficiency, and scalability, reducing their applicability in real-world scenarios. This paper proposes a novel visual analytics framework that integrates two generative AI models, Temporal Fusion Transformer (TFT) and Variational Autoencoders (VAEs), to reduce complex patterns into lower-dimensional latent spaces and visualize them in 2D using dimensionality reduction techniques such as PCA, t-SNE, and UMAP with DBSCAN. These visualizations, presented through coordinated and interactive views and tailored glyphs, enable intuitive exploration of complex multivariate temporal patterns, identifying patterns’ similarities and uncover their potential correlations for a better interpretability of the AI outputs. The framework is demonstrated through a case study on power grid signal data, where it identifies multi-label grid event signatures, including faults and anomalies with diverse root causes. Additionally, novel metrics and visualizations are introduced to validate the models and assess the performance, efficiency, and consistency of latent maps generated by VAE, which have been utilized in prior studies for latent space cartography and used as a benchmark in this study, and the emerging TFT architecture under various configurations. These analyses provide actionable insights for model parameter tuning and reliability improvements. Comparative results highlight that TFT achieves shorter run times and superior scalability to diverse time-series data shapes compared to VAE. This work advances fault diagnosis in multivariate time series, fostering explainable AI to support critical system operations.

Explainable AI↗
DOE OSTI2022

A Sparse and Low Rank Penalized Signal Decomposition Model with Constraints: Anomaly Detection in PV Systems

Recently, robust PCA has seen its wide application in various industries for its ability to perform the task of anomaly detection. The essence of robust PCA approach is to break down the signal into a low rank component and sparse component. In many applications, a simple breakdown of the signal without accounting for the signs of low rank components and sparse components would violate the physical constraints of the decomposed signal. In addition, often times, the signals in the real world collected for a long duration has smooth changes within a day and between days. As an example, the power signals collected in a photovoltaic (PV) system are cyclostationary, exhibiting these characteristics. Neglecting the smoothness of signals would result in miss detection of anomalous signals which are smooth within a day but non-smooth between days and vice versa. In this paper, we developed a signal decomposition approach for the purpose of anomaly detection based on the idea of low rank and sparse decomposition taking into consideration the signs of the decomposed low rank and sparse components and the within-day and between-day smooth changes in the original signals. The proposed unsupervised approach for fault detection eliminates the need for faulty samples required by other machine learning methods. It does not require the full I-V characteristics to work. Furthermore, there is no need for complex modelling of PV systems as in the case of power loss analysis. Using Monte Carlo simulations, we demonstrate the ability of our proposed approach for detecting anomalies of different duration and severity in PV systems.

14 SOLAR ENERGY↗
DOE OSTI2023

ANN-based ground motion model for Turkey using stochastic simulation of earthquakes

SUMMARY Turkey is characterized by a high level of seismic activity attributed to its complex tectonic structure. The country has a dense network to record earthquake ground motions; however, to study previous earthquakes and to account for potential future ones, ground motion simulations are required. Ground motion simulation techniques offer an alternative means of generating region-specific time-series data for locations with limited seismic networks or regions with seismic data gaps, facilitating the study of potential catastrophic earthquakes. In this research, a local ground motion model (GMM) for Turkey is developed using region-specific simulated records, thus constructing a homogeneous data set. The simulations employ the stochastic finite-fault approach and utilize validated input-model parameters in distinct regions, namely Afyon, Erzincan, Duzce, Istanbul and Van. To overcome the limitations of linear regression-based models, artificial neural network is used to establish the form of equations and coefficients. The predictive input parameters encompass fault mechanism (FM), focal depth (FD), moment magnitude (Mw), Joyner and Boore distance (RJB) and average shear wave velocity in the top 30 m (Vs30). The data set comprises 7359 records with Mw ranging between 5.0 and 7.5 and RJB ranging from 0 to 272 km. The results are presented in terms of spectral ordinates within the period range of 0.03–2.0 s, as well as peak ground acceleration and peak ground velocity. The quantification of the GMM uncertainty is achieved through the analysis of residuals, enabling insights into inter- and intra-event uncertainties. The simulation results and the effectiveness of the model are verified by comparing the predicted values of ground motion parameters with the observed values recorded during previous events in the region. The results demonstrate the efficacy of the proposed model in simulating physical phenomena.

Karimzadeh, Shaghayegh (ORCID:0000000337531676)↗
DOE OSTI2020

Parametric Study to Assess Technical Prospect Feasibility for Offshore CO 2 Storage

This report is a white paper on the parametric study to assess technical prospect feasibility for CO 2 storage with and without CO 2 -EOR. The SECARB Offshore GOM team reviewed key parameters that previous modeling experience suggests have the largest impact on CO 2 plume size and CO 2 -EOR performance (SSEB, 2019) These parameters include tectonic impacts (reservoir dip and fault block size), reservoir properties (permeability) and development strategies (well field payout, waterflood, and CO 2 flooding strategies). These parameters are further evaluated in this study in order to identify upper and lower boundary conditions for the parametric values. The boundary conditions are needed to set constraints for key parameters when performing modeling assessments. Understanding the key parametric values for the GOM will facilitate locating successful offshore carbon capture utilization and storage (CCUS) projects in terms of plume size, CO 2 stored and, in the case of EOR, oil recovery. Importantly, this white paper focuses on both technical factors (geologic and engineering), but does not address economic, commercial, or stakeholder considerations. The report is organized into four sections and includes a review of the datasets used to define the parameter ranges, the results of a full-physics reservoir simulation, a summary of a reduced order modeling approach to evaluate the impact of the study parameters, and, finally, a discussion of the work that will occur during the next phase of the project.

02 PETROLEUM↗
DOE OSTI2020

Stacking faults and alternate crystal structures for the shape-memory alloy NiTi

Here, we use ab initio calculations to study the role of stacking faults in connecting the high-temperature B2 and the theoretically predicted low-temperature B33 NiTi phases. In contrast with prior work, we describe the B 2 → B 33 phase transformation in terms of alternate bilayer shifts by 1 2 [100] on the ( 011 ) B 2 plane, obtaining a viable pathway; the same mechanism could also work with the B19 parent phase. We then examine B33-like structures built from alternate stacking sequences of B19 bilayers, constructed to have monoclinic tilt angles close to the experimentally reported NiTi B19' martensite, and find four low-energy stacking-fault variants with energies 5.8–8.5 meV/atom above the calculated B19' martensite structure, suggesting that such structures might appear as a part of the NiTi martensite phase at low temperatures. Investigating further the occurrence of specific coordinated planar shifts in NiTi systems, we report a dynamically stable NiTi B27 phase and find that it is only 1.2 meV/atom above the calculated B33 ground-state structure, thus having a potential to also play a role in NiTi martensitic phase transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

WHOLESCALE - Water & Hole Observations Leverage Effective Stress Calculations And Lessen Expenses (Final Technical Report 2020 - 2024)

The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States. To reach this goal, the WHOLESCALE team has developed a methodology to incorporate and interpret data from four methods of measurement into a multi-physics model that couples thermal, hydrological, and mechanical (T H-M) processes. The WHOLESCALE team has applied this methodology at the San Emidio geothermal field, located ~100 km north of Reno, Nevada in the northwestern Basin and Range province. The WHOLESCALE team includes 30 individuals working at two universities, two national laboratories, and one industry partner. Two master-degree students and five post-doctoral researchers have gained professional experience and earned partial financial support via the WHOLESCALE project. The WHOLESCALE team has taken advantage of the perturbations created by changes in pumping operations during planned shutdowns in 2016, 2021, and 2022 to infer temporal changes in the state of stress in the geothermal system at San Emidio, Nevada, U.S. The WHOLESCALE results support the working hypothesis that increasing pore-fluid pressure reduces the effective normal stress acting across fault zones. During normal operations, pumping in deep production wells decreases fluid pressures and thus increases the effective normal stresses on faults, reducing microseismicity. During planned shutdowns, the cessation of production increases pore-fluid pressure and reduces effective normal stress. The WHOLESCALE products generated during the 4-year period between 2020 and 2024 include: three articles published in the open-access, peer-reviewed scientific literature, two master’s theses, 20 presentations or papers at scientific conferences, and 17 data sets available on public repositories. The WHOLESCALE project has been completed in two phases that included three performance periods separated by two Go/No-go Stage Gate Reviews. Tasks were classified by data type (i.e., Geologic Structure, Borehole, Geodesy, Hydrology, Seismology, and Modeling). The first phase of the project started July 31, 2020 and included ongoing project coordination (Task 1), a project kickoff (Task 2), analysis of existing data (Task 3), development of the initial stress model & deployment design (Task 4), and Go/No-go Decision Point #1 (Task 5). Phase II began with implementing the 2022 deployment (Task 6), followed by Go/No-go Decision Point #2 (Task 7) The remainder of Phase II consisted of analyzing data collected during deployment (Task 8), calibration of the stress model on all observations (Task 9), and the Final Review (August 23, 2024) & Reporting (Task 10).

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Structure of V-defects in long wavelength GaN-based light emitting diodes

The V-defect is a naturally occurring inverted hexagonal pyramid structure that has been studied in GaN and InGaN growth since the 1990s. Strategic use of V-defects in pre-quantum well superlattices or equivalent preparation layers has enabled record breaking efficiencies for green, yellow, and red InGaN light emitting diodes (LEDs) utilizing lateral injection of holes through the semi-polar sidewalls of the V-defects. In this article, we use advanced characterization techniques such as scattering contrast transmission electron microscopy, high angle annular dark field scanning transmission electron microscopy, x-ray fluorescence maps, and atom probe tomography to study the active region compositions, V-defect formation, and V-defect structure in green and red LEDs grown on (0001) patterned sapphire and (111) Si substrates. We identify two distinct types of V-defects. The “large” V-defects are those that form in the pre-well superlattice and promote hole injection, usually nucleating on mixed (Burgers vector b = ±a±c⁠) character threading dislocations. In addition, “small” V-defects often form in the multi-quantum well region and are believed to be deleterious to high-efficiency LEDs by providing non-radiative pathways. The small V-defects are often associated with basal plane stacking faults or stacking fault boxes. Furthermore, we show through scattering contrast transmission electron microscopy that during V-defect filling, the threading dislocation, which runs up the center of the V-defect, will “bend” onto one of the six {10$\overline{1}$1} semi-polar planes. In conclusion, this result is essential to understanding non-radiative recombination in V-defect engineered LEDs.

42 ENGINEERING↗
DOE OSTI2024

Single-ancilla ground state preparation via Lindbladians

We design a quantum algorithm for ground state preparation in the early fault tolerant regime. As a Monte Carlo style quantum algorithm, our method features a Lindbladian where the target state is stationary. The construction of this Lindbladian is algorithmic and should not be seen as a specific approximation to some weakly coupled system-bath dynamics in nature. Our algorithm can be implemented using just one ancilla qubit and efficiently simulated on a quantum computer. It can prepare the ground state even when the initial state has zero overlap with the ground state, bypassing the most significant limitation of methods like quantum phase estimation. As a variant, we also propose a discrete-time algorithm, demonstrating even better efficiency and providing a near-optimal simulation cost depending on the desired evolution time and precision. Numerical simulations using Ising and Hubbard models demonstrate the efficacy and applicability of our method. Published by the American Physical Society 2024

Ding, Zhiyan (ORCID:000000018863403X)↗
DOE OSTI2020

Updated Available Drawdowns for Big Hill SPR Caverns - Model Including the Caprock Fault

This report updates the estimated values of the baseline available drawdowns for the caverns at the Big Hill storage facility, and an updated table listing the available drawdowns. An updated finite element numerical analysis model, which included a fault in the caprock layers, was constructed and the daily data of actual wellhead pressures and oil-brine interfaces was used. The number of available drawdowns for each of the Big Hill SPR caverns is estimated using the new model. All caverns are predicted to have five available drawdowns remaining from a geomechanical perspective. BH-101 and 105 have a region of concern at the floor edge and/or sloping floor, where tensile and dilatant stresses are predicted to occur during each workover. The tensile state is predicted to occur because of the geometries of the edge and floor. Therefore, geomechanical examination for two caverns would be recommended after a drawdown leach. The well integrity of each cavern is not investigated in this report. The estimates for the number of baseline available drawdowns are subject to change in the future as the knowledge of physical phenomena at the sites, and the further development of the models of geomechanical behavior at the sites, evolve over time.

02 PETROLEUM↗
DOE OSTI2026

Simulating plasma wave propagation on a superconducting quantum chip

Quantum computers may one day enable the efficient simulation of strongly coupled plasmas that lie beyond the reach of classical computation in regimes where quantum effects are important and the scale separation is large. Here, in this article, we take a first step toward efficient simulation of quantum plasmas by demonstrating linear plasma wave propagation on a superconducting quantum chip. Using high-fidelity and highly expressive device-native gates, combined with an error-mitigation technique, we simulate the scattering of laser pulses from inhomogeneous plasmas. Our approach is made feasible by the identification of a suitable local spin model whose excitations mimic plasma waves, and whose circuit implementation requires a lower gate count than other proposed approaches that would require a future fault-tolerant quantum computer. This work opens avenues to study more complicated phenomena that cannot be simulated efficiently on classical computers, such as nonlinear quantum dynamics when strongly coupled plasmas are driven out of equilibrium.

general physics↗
DOE OSTI2023

Physics-based Induced Earthquake Forecasting: Process Understanding, and Hazards Mitigation

Disposal of saltwater co-produced with oil and gas is linked to elevated seismicity in the Central and Midwest US. There is a concern that these events may lead to widespread damage and an overall increase in seismicity. Thus an improved understanding of the spatially and temporally variable deformation and stress field associated with fluid injection operation is critically necessary for evaluating time-varying seismic hazards. Despite the improvements in seismic monitoring capacity and the resulting decrease in the magnitude detection threshold, estimates of induced earthquake probability remain elusive due to insufficient models incapable of accounting for the complex physics governing the process of induced seismicity. The proposed research effort will comprehensively analyze, integrate, and interpret geodetic, injection, and seismic data in the vicinity of the injection sites in Oklahoma to resolve the 4-dimensional distribution of pore pressure and stress in the shallow crust. This project, in particular, is focused on exploring the statistical relation between injection operation and increased earthquake hazard. The amplitude of and the extent to which pore pressure changes are determined by some factors, in particular, the hydrogeological properties of the rocks, such as diffusivity. Thus the available deformation data will be used to constrain hydrogeological properties of the medium, to accurately resolve the evolution of crustal stresses due to fluid injection. Having the time-varying models of stress changes, a statistical framework will be implemented to estimate the time-dependent probability of large earthquakes on the nearby fault systems. These data and models help to improve seismic hazard estimates and aid in constructing operational-induced earthquake forecast models. This information can also be integrated into the updated U.S. National Seismic Hazard Map, which local communities and authorities use in their earthquake risk estimates and mitigation efforts.

58 GEOSCIENCES↗
DOE OSTI2025

Characterization of the microstructure of yttrium hydride under proton irradiation

High moderation per unit volume solid moderator materials like yttrium hydride (YH x ) are necessary for compact nuclear microreactors. However, the phase stability and hydrogen transport processes of YH x under high-temperature irradiation are largely unknown. Proton irradiation was conducted on YH x at 300 °C and 580 °C to 0.2 dpa using 1 MeV or 2 MeV protons in a high-vacuum environment. The hydrogen concentration was determined before and after irradiation using elastic recoil detection analysis, and microstructural evolution was examined via post-irradiation scanning transmission electron microscopy and Raman spectroscopy. Dislocation loops and cavities were observed in all conditions; their distribution was correlated with the bombarding proton energy and ion irradiation temperature. This work revealed that hydrogen retention is proportional to the formation of traps for hydrogen gas atoms and identified pathways for hydrogen release. The relative contributions of bulk or fast diffusion paths, such as grain boundaries, delamination boundaries, and stacking faults are discussed; the primary mechanisms of hydrogen loss are likely based on diffusion, ruling out artefacts of the experimental design. In conclusion, the study suggests proton irradiation may be a strong surrogate to study hydrogen transport in hydride moderator materials under irradiation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

GronOR: Massively Parallel and GPU-Accelerated Non-Orthogonal Configuration Interaction for Large Molecular Systems

GronOR is a program package for non-orthogonal configuration interaction calculations for an electronic wave function built in terms of anti-symmetrized products of multi-configuration molecular fragment wave functions. The two-electron integrals that have to be processed may be expressed in terms of atomic orbitals or in terms of an orbital basis determined from the molecular orbitals of the fragments. The code has been specifically designed for execution on distributed memory massively parallel and Graphics Processing Unit (GPU)-accelerated computer architectures, using an MPI+OpenACC/OpenMP programming approach. The task-based execution model used in the implementation allows for linear scaling with the number of nodes on the largest pre-exascale architectures available, provides hardware fault resiliency, and enables effective execution on systems with distinct central processing unit-only and GPU-accelerated partitions. The code interfaces with existing multi-configuration electronic structure codes that provide optimized molecular fragment orbitals, configuration interaction coefficients, and the required integrals. Algorithm and implementation details, parallel and accelerated performance benchmarks, and an analysis of the sensitivity of the accuracy of results and computational performance to thresholds used in the calculations are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗