Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture modelling”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

3D seismic imaging of a fracture damage zone controlling reservoir compartmentalization at the Raft River EGS using multi-azimuth walkaway VSP

Accurate imaging of steeply dipping fracture zones in crystalline enhanced geothermal systems (EGS) is critical for constraining permeability architecture and guiding stimulation design. However, such structures remain poorly resolved by conventional surface seismic methods. We present a fully three-dimensional (3D) elastic-waveform inversion-migration workflow applied to multi-azimuth walkaway vertical seismic profiling (VSP) data acquired in a deviated borehole at the Raft River EGS. The workflow integrates first-arrival traveltime tomography, multi-scale elastic waveform inversion (EWI), and elastic least-squares reverse-time migration (ELSRTM) to recover high-resolution compressional-and shear-wave velocity models and to image structural discontinuities in the crystalline basement. The results reveal a laterally continuous low-velocity anomaly, with shear-wave velocity reductions of 25-30%, consistent with fractureinduced mechanical weakening. Two steeply dipping discontinuities bound a 50-80 m wide fracture damage zone. Independent constraints from microseismic clustering and geochemical compartmentalization corroborate the geometry and structural significance of this feature. Synthetic modeling further confirms that structures of this scale are resolvable (~30 m). These findings indicate that the Narrows structure is a distributed fracture damage zone rather than a discrete fault plane. By resolving fracture-zone geometry at the tens-of-meters scale using a single borehole, this workflow provides a practical and transferable approach for improving structural characterization, reducing uncertainty in permeability architecture, and supporting reservoir modeling and stimulation design in fractured crystalline EGS reservoirs.

58 GEOSCIENCES↗
DOE OSTI2020

The Sealing Mechanisms of a Fracture in Opalinus Clay as Revealed by in situ Synchrotron X-Ray Micro-Tomography

The detailed mechanisms of the sealing of a single fracture, from hydration to almost complete closure by increase of confining pressure, as monitored from in situ synchrotron X-ray microtomography during the flow of carbonated water, is here shown for the first time. Different mechanisms play the key role at different stages in the evolution of the fracture. Hydration mechanically weakens the surfaces of the fracture and induces a first closure due to microcracking at the asperity contacts, increasing their size and creating choke points. Increase in confining stress promptly hydraulically seals the fracture by closing the main choke point, with a relative small deformation of the sample. Finally, the more pervasive mechanical deformation observed at higher stresses almost completely seals the whole fracture. The evolution of the sample has been described and quantified using 4D image processing, focusing on the evolution of aperture and digital volume correlation. Hydraulic properties of the sample at different stages have been modeled via Stokes flow simulation, and the results compared to the morphometric analysis, finding positive correlations with the average fracture aperture variation along the flowpath in function of confining pressure. Opalinus Clay is found to be a rock markedly prone to sealing in case of flow with carbonated water; this behavior is the result of the large fraction of clays and of its microstructure, lacking both cementing phases and large stiff particles. CO2 in this sample has no evident role in the evolution of the fracture; chemically induced weathering on the surface has not been detected, in contrast with the behavior observed in samples with carbonates as cementing phase.

58 GEOSCIENCES↗
DOE OSTI2023

Identification of Faults Susceptible to Induced Seismicity (Final Report)

Central to the work documented in this report is the capability of geocellular models to represent the geologic conceptual model updated with fault identification from machine learning and joint inversion modeling of microseismic data measured and recorded as a consequence of CO 2 injection at a field demonstration site: the Illinois Basin - Decatur Project (IBDP). This work required seven unique geocellular models with 100s of simulated variations to gain a very high degree of confidence in the identification of geologic features present that contributed to induced microseismicity at IBDP. All forward modeling: pressure modeling, stress modeling, and seismic modeling used the same geologic conceptual model and representations of that model at different scales. The pressure modeling and poroelastic modeling created “snapshots” of pore pressure and stress field changes at different times during CO 2 injection, in which microseismic events were clustered (in time). These pressure and stress snapshots, within the framework and architecture of the geologic conceptual model via the geocellular model, informed the single fault and fault network models to ascertain the likelihood of fault movement (seismic or aseismic). The outcomes of the pressure, stress, and fault/fault network (seismic) modeling confirmed that the faults in the geologic conceptual model in Task 2 were likely the source of microseismic events measured at IBDP and acted as conduits for pressure to be transmitted from the injection interval into the Precambrian crystalline basement rock. This closely coordinated and integrated unique modeling approach was conducted to prove the viability of our proposed workflow 1) to better resolve crystalline basement faults, 2) detect subseismic faults that could be activated by injection, 3) increase the certainty in fault detection and their susceptibility to release seismic energy, and 4) understand transmission of pressure vertically from the well to the underlying fractured crystalline basement. The proposed methodology was effective in guiding an iterative process of calibrating forward modeling results based on similar geocellular models while honoring the geologic conceptual model (i.e., characterization data and knowledge of regional geology); this led to higher level of certainty in the identification of fault/faults zones to control seismicity and transmission of pressure to the regions of recorded and located injection induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2020

Mechanistic Approach to Analyzing and Improving Unconventional Hydrocarbon Production [Slides]

DOE research is developing the physical basis and tools needed to manage pressure effectively to increase recovery efficiency. By coupling fast, accurate physics with machine learning, DOE is producing science-based platforms any operator can use. DOE’s research portfolio is targeting hydrocarbon transport at multiple scales, with the goal of increasing recovery efficiency. DOE’s research has led to new, fast & accurate platforms for predicting gas production from fractured shales. Using data from the MSEEL-I site to calibrate our physics-based model, we have early results on pressure management. We have shown that both mechanical and chemical processes in the matrix can negatively impact production.

08 HYDROGEN↗
DOE OSTI2024

Development of Graphite Thermal and Mechanical Modeling Capabilities in Grizzly

Nuclear-grade graphites are used extensively in the core designs of multiple types of advanced nuclear reactors. In the reactor environment, graphite is exposed over long durations to extreme conditions, including high temperatures, radiation and potentially molten salt and oxygen. Exposure to these conditions can cause several degradation mechanisms in graphite, including nonuniform volumetric strains induced by irradiation and thermal expansion, which lead to stresses that can compromise the performance of graphite components. Evaluating component integrity, predicting component performance over the reactor lifetime, and developing design standards all require robust tools for predicting fracture initiation and propagation in graphite structural components in nuclear reactors. This report documents progress in an ongoing effort to develop modeling and simulation tools in the Grizzly code for predicting the performance of graphite exposed to reactor conditions. Recent developments include a set of thermal and mechanical models that now include the IG-110, NBG-18, and H-451 graphite grades. Improvements have also been made to a nonlinear damaged plasticity model applicable to predicting damage under tension and compression to quasibrittle materials, including graphite. In addition, enhancements have been made to the extended finite element method implementation targeted at simulating graphite fracture. These include new capabilities for crack nucleation in the interior of a solid body, improved treatment of crack nucleation on free surfaces, and more robust modeling approaches for crack growth approaching free surfaces or other cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the Feasibility of Deep Geothermal Wells Using Numerical Reservoir Simulation

This study examines the geothermal energy extraction potential from the basement rock within the Denver–Julesburg Basin, focusing on the flow performance and heat extraction efficiency of different geothermal well configurations. It specifically compares U-shaped, V-shaped, inclined V-shaped, and pipe-in-pipe configurations against enhanced geothermal system setups. Through numerical modeling, we evaluated the thermal behavior of these systems under various operational scenarios and fracture conditions. The results suggest that while closed-loop systems offer moderate temperature increases, Enhanced geothermal system configurations show substantial potential for high-temperature extraction. This underscores the importance of evaluating well configurations in complex geological settings. The insights from this study aid in strategic geothermal energy planning and development, marking significant advancements in geothermal technology and setting a foundation for future explorations and optimizations.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
DOE OSTI2020

Fluid-driven fractures in granular media: Insights from numerical investigations

We investigate the mechanisms of opening-mode fracture initiation in granular media. The study is based on a simulation of grain-scale fluid-grain interactions through a coupled numerical approach in which the discrete element method is used to solve for the mechanics of a solid granular medium, and computational fluid dynamics is used to model fluid flow and drag forces. In this work, we present benchmark problems with analytical solutions and validate this numerical model against experiments on a viscous-drag-driven cavity in the literature. Additional simulation results show fracture initiation mechanisms in a random granular packing subjected to constant boundary stresses and to fluid injection with a localized source. The dimensionless variable $F_s/F_{sk}$ (ratio of seepage force $F_s$ and skeletal force $F_{sk}$) incorporates the impacts of physical properties and injection parameters including fluid viscosity, injection velocity, grain size, and effective stresses, and it has been used as a criterion separating regimes of fluid invasion and drag-driven fracture opening. Our simulation results show that $F_s/F_{sk}$ in combination with $τ_1$ (ratio of diffusion time from hydromechanical coupling and injection time) serves as a prediction of fracture opening within granular packing. We suggest a simple criterion ($F_s/F_{sk}$ > 1 or $τ_1$ > 0.17) that is valid for various types of granular media and injection conditions to determine if fracture opening will occur. Among other applications, this study is useful to predict the initiation and propagation of fractures in natural sediments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Quantum-well states in fractured crystals of the heavy-fermion material CeCoIn5

Quantum-well states appear in metallic thin films due to the confinement of the wave function by the film interfaces. Using angle-resolved photoemission spectroscopy, we unexpectedly observe quantum-well states in fractured single crystals of CeCoIn5. We confirm that confinement occurs by showing that these states’ binding energies are photon-energy independent and are well described with a phase accumulation model, commonly applied to quantum-well states in thin films. This indicates that atomically flat thin films can be formed by fracturing hard single crystals. For the two samples studied, our observations are explained by free-standing flakes with thicknesses of 206 and 101 Å. We extend our analysis to extract bulk properties of CeCoIn5. Specifically, we obtain the dispersion of a three-dimensional band near the zone center along in-plane and out-of-plane momenta. We establish part of its Fermi surface, which corresponds to a hole pocket centered at G. We also reveal a change of its dispersion with temperature, a signature that may be caused by the Kondo hybridization.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Detecting and Characterizing Fracture Zones Using a Convolutional Neural Network

This project directly supports the Geothermal Technologies Office (GTO) objectives outlined in the Multi-Year Program Plan (MYPP) by advancing two key research areas: “Exploration and Characterization” and “Data, Modeling, and Analysis.” This project has successfully demonstrated a pre-drilling ability to image and characterize the distribution and connectivity of subsurface faults and fractures, key parameters for identifying permeable pathways that enable geothermal fluids to circulate and produce energy. Specifically, we developed and implemented innovative machine learning methodologies to enhance geothermal exploration. Large-scale faults were detected using a Convolutional Neural Network (CNN), while small-scale fractures were characterized using a novel Double-Beam Neural Network (DBNN). These tools have proven both technically effective and cost-efficient by reducing reliance on expensive exploratory drilling. Through collaboration with our geothermal industry partner, this research has significantly advanced techniques for identifying hidden geothermal systems and extending the productive lifespan of existing geothermal fields. We applied our methods to two geothermal fields—Soda Lake (Nevada) and Lightning Dock (New Mexico)—to identify shallow steam-charged fracture zones and characterize deep faults at depths of 1.5-2 km. The steam zone identified at the Soda Lake geothermal field showed excellent agreement with prior drilling data, validating the effectiveness of our approaches. In addition, the analysis revealed three new prospective drilling targets for further development and verification. The outcomes of this project improve our scientific understanding of geothermal reservoir behavior, enhance exploration efficiency, extend the economic life of existing geothermal plants. Ultimately, these advancements contribute to GTO’s goal of achieving more sustainable, affordable, and data-driven geothermal energy development across the United States.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Obstruction and Interference in Low-Energy Models for Twisted Bilayer Graphene

The electronic bands of twisted bilayer graphene (TBLG) with a large-period moiré superlattice fracture to form narrow Bloch minibands that are spectrally isolated by forbidden energy gaps from remote dispersive bands. When these gaps are sufficiently large, one can study a band-projected Hamiltonian that correctly represents the dynamics within the minibands. This inevitably introduces nontrivial geometrical constraints that arise from the assumed form of the projection. Here we show that this choice has a profound consequence in a low-energy experimentally observable signature that therefore can be used to tightly constrain the analytic form of the appropriate low-energy theory. Additionally, we find that this can be accomplished by a careful analysis of the electron density produced by backscattering of Bloch waves from an impurity potential localized on the moiré superlattice scale. We provide numerical estimates of the effect that can guide experimental work to clearly discriminate between competing models for the low-energy band structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Development and Refinement of Bison Models for Higher Burnup Application

Electric Power Research Institute (EPRI) is an international research organization funded by organizations in the energy and government sectors. The company conducts research to tackle the energy challenges of its member organizations. In support of the U.S. nuclear industry, EPRI and Idaho National Laboratory were awarded a U.S. Department of Energy Technology Commercialization Fund to develop a loss of coolant accident (LOCA) analysis tool. This tool was developed to support burnup extension activities by supplementing the existing capabilities of operators of existing large light-water reactor plants. The tool couples INL developed codes Bison and RELAP5-3D for advanced simulation of LOCA behavior. The project consists of seven major tasks, one of which is the development and refinement of Bison models for higher burnup applications. This task is broken down into six targeted areas of development and refinement including: formation of high burnup structure, fuel fracture and pulverization, transient fission gas release, axial gas communication within a fuel rod, cladding ballooning and burst, and axial fuel relocation. This report summarizes the models for both the fuel and cladding in Bison for high burnup application, identifies the areas of development and refinement necessary to enhance the Bison capabilities to be used in the coupled simulation tool with RELAP5-3D. Integration and regression testing to demonstrate proper implementation of the new features in Bison is highlighted.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Microstructure based modeling of friction stir welded joint of dissimilar metals using crystal plasticity

We report the friction stir welding (FSW) process shows promising results in joining dissimilar metals which are otherwise almost impossible to join using traditional welding techniques. Being a new technique, the deformation and the failure mechanism of the joints made by the FSW process needs to be investigated. In this work, a joint between AZ31 Mg alloy and DP590 steel is modeled using phenomenological crystal plasticity formulation on the mesoscale in the form of a representative volume element (RVE). The interface of the two materials is modeled using a cohesive zone model. A parametric study has been performed to understand the effect of grain size and interface fracture toughness as well as strength on the mechanical performance of the joint. It was found that the grain size of AZ31 Mg alloy, as well as DP590 steel, has little effect on the overall joint performance. On the other hand, interface fracture toughness and strength have a significant impact on the mechanical properties of the joint

36 MATERIALS SCIENCE↗
DOE OSTI2023

Real-Time Monitoring of Fracture Dynamics with a Contrast Agent-Assisted Electromagnetic Method

In collaboration with the Advanced Energy Consortium, our team has previously demonstrated that the placement of electrically active proppants (EAPs) in a hydraulic fracture surveyed by electromagnetic (EM) methods can enhance the imaging of the stimulated reservoir volumes during hydraulic fracturing. That work culminated in constructing a well-characterized EAP-filled fracture anomaly at the Devine field pilot site (DFPS). In subsequent laboratory studies, we observed that the electrical conductivity of our EAP correlates with changes in pressure, salinity, and flow. Thus, we postulated that the EAP could be used as an in-situ sensor for the remote monitoring of these changes in previously EAP-filled fractures. This paper presents our latest field data from the DFPS to demonstrate such correlations at an intermediate pilot scale. We conducted surface-based EM surveys during freshwater (200 ppm) and saltwater (2,500 ppm) slug injections while running surfaced-based EM surveys. Simultaneously, we measured the following: 1) bottomhole pressure and salinity in five monitoring wells; 2) injection rate using high-precision data loggers; 3) distributed acoustic sensors in four monitoring wells; and 4) tiltmeter data on the survey area. We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, by comparing the electrical field traces with the bottomhole pressure, flow rate, and salinity, we concluded that the observed electric field in our study is influenced by fracture dilation and flow rate. Salinity effect was observed when saltwater was injected. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce all of the measured electric field magnitudes. Preliminary estimates showed that including streaming potential in our geophysical model may be needed to reduce the simulation mismatch. The methods developed and demonstrated during this study will lead to a better understanding of the extent of fracture networks, formation stress states, fluid leakoff and invasion, characterizations of engineered fracture systems, and other applications where monitoring subsurface flow tracking is deemed important.

02 PETROLEUM↗
DOE OSTI2025

Numerical simulation of asteroid geometry variance on airburst threat

For an atmospheric airburst the primary source of concern when assessing uncertainty is the size and velocity. Determining these properties provides the basis for threat assessment, as the total energy of the asteroid may then be estimated, and the threat investigated thoroughly. Even with clarity as to how much energy an asteroid may deposit, a great deal of uncertainty still exists for the actual energy deposition process. One such source of uncertainty is the geometry of the incoming asteroid. The geometry of an asteroid will alter the stress distribution during entry, which adds uncertainty to when fracture will occur. Here, in this study, we use Smoothed Particle Hydrodynamics to model the atmospheric airburst of Tunguska-scale asteroids with varying geometric profiles, including a sphere, ellipsoid, binary and superellipsoid. Each asteroid is modeled as a homogenous structure with strength. We assess uncertainty through a series of planar 2D simulation cases for each geometry, comparing the source of stochasticity across geometries. A single 3D airburst simulation for each geometry is also analyzed. Additionally, the 3D cases are compared to the highly uncertain Tunguska event, predicting variance in burst height across geometries, but all bounded by theoretical burst heights proposed for Tunguska.

Airburst↗
DOE OSTI2025

Graph neural networks for mechanical property prediction of 2D fiber composites

This work investigates the ability of graph neural networks (GNNs) to homogenize 2D fiber composite microstructures. We use different inhomogeneity and anisotropy indices to motivate and show that the Volume Elements (VEs) used in ML methods should ideally be far from their Representative Volume Element (RVE) size limit and, consequently, are notably anisotropic. Hence, training only the isotropic limit properties may not be acceptable. Another aspect is the need to normalize elastic stiffness values for ML, especially when high elastic contrast ratios are encountered between composite phases or in the material set. We introduce a normalization technique based on the mean-field method (MFM) to handle such high contrast ratios and train for the entire stiffness tensor. We show that the proposed GNN approaches exhibit high accuracy and efficiency compared to traditional methods and convolutional neural networks, utilizing unstructured graphs constructed from microstructure topology. Our model successfully predicts the stiffness tensor, peak strength under bulk damage, and brittle fracture initiation strength across diverse microstructure configurations while maintaining high accuracy even for extreme material contrasts and volume fractions. We also present a method to improve prediction accuracy for small dataset sizes using Voronoi partitioning.

Brittle strength↗
DOE OSTI2023

Topology-enhanced mechanical stability of swelling nanoporous electrodes

Abstract Materials like silicon and germanium offer a 10-fold improvement in charge capacity over conventional graphite anodes in lithium-ion batteries but experience a roughly threefold volume increase during lithiation, which challenges ensuring battery integrity. Nanoporous silicon, created by liquid-metal-dealloying, is a potentially attractive anode design to mitigate this challenge, exhibiting both higher capacity and extended cycle lifetimes. However, how nanoporous structures accommodate the large volume change is unknown. Here, we address this question by using phase-field modeling to produce nanoporous particles and to investigate their elastoplastic swelling behavior and fracture. Our simulations show that enhanced mechanical stability results from the network topology consisting of ligaments connected by bulbous, sphere-like nodes. The ligaments forcefully resist elongation while the nodes, behaving like isolated spherical particles, experience large stresses driving fracture. However, being smaller compared to a sphere of the same volume as the entire nanoporous particle, the nodes are more protected against fracture.

36 MATERIALS SCIENCE↗