Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fracture prediction”

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

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Influence of extreme fracture flow channels on the thermal performance of open-loop geothermal systems at commercial scale

Abstract Adequate stewardship of geothermal resources requires accurate forecasting of long-term thermal performance. In enhanced geothermal systems and other fracture-dominated reservoirs, predictive models commonly assume constant-aperture fractures, although spatial variations in aperture can greatly affect reservoir permeability, fluid flow distribution, and heat transport. Whereas previous authors have investigated the effects of theoretical random aperture distributions on thermal performance, here we further explore the influence of permeability heterogeneity considering field-constrained aperture distributions from a meso-scale field site in northern New York, USA. Using numerical models of coupled fluid flow and heat transport, we conduct thermal–hydraulic simulations for a hypothetical reservoir consisting of a relatively impervious porous matrix and a single, horizontal fracture. Our results indicate that in highly channelized fields, most well design configurations and operating conditions result in extreme rates of thermal drawdown (e.g., 50% drop in production well temperatures in under 2 years). However, some other scenarios that account for the risks of short-circuiting can potentially enhance heat extraction when mass flow rate is not excessively high, and the direction of geothermal extraction is not aligned with the most permeable features in the reservoir. Through a parametric approach, we illustrate that well separation distance and relative positioning play a major role in the long-term performance of highly channelized fields, and both can be used to help mitigate premature thermal breakthrough.

Rangel-Jurado, Nicolás↗
DOE OSTI2022

From Fluid Flow to Coupled Processes in Fractured Rock: Recent Advances and New Frontiers

Abstract Quantitative predictions of natural and induced phenomena in fractured rock is one of the great challenges in the Earth and Energy Sciences with far‐reaching economic and environmental impacts. Fractures occupy a very small volume of a subsurface formation but often dominate fluid flow, solute transport and mechanical deformation behavior. They play a central role in CO 2 sequestration, nuclear waste disposal, hydrogen storage, geothermal energy production, nuclear nonproliferation, and hydrocarbon extraction. These applications require predictions of fracture‐dependent quantities of interest such as CO 2 leakage rate, hydrocarbon production, radionuclide plume migration, and seismicity; to be useful, these predictions must account for uncertainty inherent in subsurface systems. Here, we review recent advances in fractured rock research covering field‐ and laboratory‐scale experimentation, numerical simulations, and uncertainty quantification. We discuss how these have greatly improved the fundamental understanding of fractures and one's ability to predict flow and transport in fractured systems. Dedicated field sites provide quantitative measurements of fracture flow that can be used to identify dominant coupled processes and to validate models. Laboratory‐scale experiments fill critical knowledge gaps by providing direct observations and measurements of fracture geometry and flow under controlled conditions that cannot be obtained in the field. Physics‐based simulation of flow and transport provide a bridge in understanding between controlled simple laboratory experiments and the massively complex field‐scale fracture systems. Finally, we review the use of machine learning‐based emulators to rapidly investigate different fracture property scenarios and accelerate physics‐based models by orders of magnitude to enable uncertainty quantification and near real‐time analysis.

58 GEOSCIENCES↗
DOE OSTI2019

Prediction of the Shear Failure of Opened Rock Fractures and Implications for Rock Slope Stability Evaluation

Rock slope commonly fails due to the shear failure of rock fractures. Shear strength of rock fractures are reduced substantially once the fracture surfaces are mismatched or opened. We propose a new criterion to predict the shear strength of rock fractures in different opening states. The degree of interlocking representing the true asperity contact area is incorporated into the modified model of Saeb and Amadei. The effect of fracture opening on asperity dilation and degradation is separately considered. The transitional stress that is a critical parameter involved in the model is analytically determined based on energy consideration. The new model is validated with experimental results from direct shear tests on synthetic fractures with regular‐shaped asperities. Good agreement between the analytical solution and the experimental data confirms the capacity of the proposed model. Therefore, the model has great potential for assessing the stability of rock slopes where fractures are often opened due to stress relief and engineering disturbances.

Li, Yingchun (ORCID:0000000201083704)↗
DOE OSTI2021

Uncertainty bounds for multivariate machine learning predictions on high-strain brittle fracture

Simulation of the crack network evolution on high strain rate impact experiments performed in brittle materials is very compute-intensive. The cost increases even more if multiple simulations are needed to account for the randomness in crack length, location, and orientation, which is inherently found in real-world materials. Constructing a machine learning emulator can make the process faster by orders of magnitude. There has been little work, however, on assessing the error associated with their predictions. Estimating these errors is imperative for meaningful overall uncertainty quantification. In this work, we extend the heteroscedastic uncertainty estimates to bound a multiple output machine learning emulator. Overall, we find that the response prediction is accurate within its predicted errors, but with a somewhat conservative estimate of uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Phase-field modeling and experiments of dynamic fracture in single crystal quartz

Predicting the onset and characteristics of brittle fracture is important for a wide range of engineering and geological material applications. In this paper, we study important aspects of brittle fracture in α-quartz by phase-field modeling and experiments using a top-down approach. In the modeling framework, the work term in the Griffith energy balance is replaced with internal energy contributions that represent surface energy, thermal energy, and elastic strain energy stored in defects. This allows parametrization of individual energy contributions in terms of internal state variables and keeps track of energy partitioning after the onset of fracture. The path and history dependence of fracture is included in evolution laws for internal state variables, e.g., entropy evolution, while the energy remains a true potential. In the experimental part, dynamic compression experiments coupled with X-ray phase contrast imaging are performed on cube-like samples with a hole. In the top-down analysis, dynamic compression and three point bending experiments from the literature are simulated with the developed phase-field damage model. In conclusion, the fitted model highlights the strain rate, size, and stress state dependence of damage nucleation and evolution in single crystal α-quartz.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Radiation-induced Swelling and Microcracking in SiC Cladding for LWRs

After the Fukushima accident in 2011, SiC has attracted considerable attention as a potential material for fuel cladding that can offer increased accident tolerance in light water reactors (LWRs) compared to zircaloy cladding. Existing measurements on nuclear grade SiC/SiC composites made of high-purity, stoichiometric SiC have already shown very good radiation tolerance in various harsh environments. However, there are a number of challenges that need to be addressed before SiC composites can be successfully used in LWR technologies. One of them is irradiation-induced swelling, which can lead to development of significant stresses in the matrix, followed by microcracking and subsequent release of fission products. Models of irradiation-induced swelling had been largely empirical and therefore they cannot be used to predict the microstructural dependence of microcracking. Existing models are also limited by the lack of understanding of defects that form due to irradiation in the range of temperatures relevant to fuel cladding in LWRs (<1000°C). Many of the defects in this regime of temperature are too small to be detected with traditional transmission electron microscopy (TEM) techniques. In fact the defects observed in traditional TEM account only for 10-45% of the swelling measured in irradiated SiC. Here, we developed a multi-scale simulation methodology combined with state of the art experimental imaging techniques based on aberration corrected scanning transmission electron microscopy (STEM) to provide a scientific basis for predictive models of radiation swelling and microcracking in SiC. Specifically, we have developed a phase field fracture model that predicts microcrack initiation and growth in polycrystalline SiC. The model accounts for single crystal anisotropy, including anisotropy in the elastic constants and anisotropic surface energy. It also considers the reduced fracture strength of grain boundaries. The model was implemented using the MARMOT mesoscale nuclear materials tool. We have parameterized the model for SiC using molecular dynamics simulation results. The model has been verified using single crystal, bicrystal, and polycrystal examples. Swelling model has also been developed based on cluster dynamics simulations of defect evolution in SiC. We found that point defects and small defect clusters have a significant contribution to swelling. This finding is important because point defects and very small defect clusters are not visible in experiments (e.g., in TEM). Combined modeling and experiments can bring critical insights into how swelling depends on irradiation conditions and the microstructure. In situ nanoindentation tests and ex situ nanoindentation tests were performed on assynthesized 3C-SiC and 3C-SiC irradiated under different radiation conditions. Without radiation, bend contour movement and residual bend contours after indentation demonstrate some room temperature plastic flow and residual plastic deformation. Plastic deformation recovery is also observed on unirradiated samples during unloading. Radiation-induced embrittlement is demonstrated by differences in crack geometry in in situ tests, and is consistent with ex situ test results. Failure by single, straight, cracks indicates decreased fracture toughness due to irradiation in 3C-SiC. Future work includes coupling of the swelling and fracture models and development of a master equation that can be used in the BISON code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗
DOE OSTI2023

Fabric controls on fracture surface roughness of an architected rock material

Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.

3D printed rock↗
DOE OSTI2026

Oil and Gas Recovery Capability

Oil and gas recovery increasingly relies on stimulation and production strategies that involve complex interactions between rock deformation, fractures, and fluid flow underground. Conventional modeling approaches often simplify or separate these effects, which can limit their ability to accurately predict performance in fractured or geologically complex reservoirs.

02 PETROLEUM↗
DOE OSTI2023

Multiscale spatial analysis of fracture nodes in two dimensions

Spatial arrangement of fractures as a function of scale is an important component of fracture quantification for inferential and predictive modeling. Available methods that analyze fracture spatial arrangement are based on one-dimensional spacing data; therefore, they are limited to semi-parallel fractures. Such methods cannot be applied to fracture networks in higher dimensions, particularly when fractures have different orientations. Here to characterize fracture arrangements in two dimensions, we propose using Ripley’s K-function, as a method of point pattern analysis, to quantify spatial arrangement of fracture nodes. Fracture nodes, such as barycenters, intersection points, and tips, are point-based representations of fracture locations and connectivity within the fracture network. We introduce formulations for isotropic as well as directional analyses of spatial arrangement. In addition, we derive formulations for edge correction in circular and rectangular study domains. Finally, we demonstrate applications of Ripley’s K-function on two natural fracture datasets. Our proposed method supports quantification and characterization of fracture spatial arrangements that allow practitioners to build representative models of fractures in the subsurface.

02 PETROLEUM↗
DOE OSTI2025

Learning the factors controlling mineral dissolution in three-dimensional fracture networks: applications in geologic carbon sequestration

We perform a set of high-fidelity simulations of geochemical reactions within three-dimensional discrete fracture networks (DFN) and use various machine learning techniques to determine the primary factors controlling mineral dissolution. The DFN are partially filled with quartz that gradually dissolves until quasi-steady state conditions are reached. At this point, we measure the quartz remaining in each fracture within the domain as our primary quantity of interest. We observe that a primary sub-network of fractures exists, where the quartz has been fully dissolved out. This reduction in resistance to flow leads to increased flow channelization and reduced solute travel times. However, depending on the DFN topology and the rate of dissolution, we observe substantial variability in the volume of quartz remaining within fractures outside of the primary subnetwork. This variability indicates an interplay between the fracture network structure and geochemical reactions. We characterize the features controlling these processes by developing a machine learning framework to extract their relevant impact. Specifically, we use a combination of high-fidelity simulations with a graph-based approach to study geochemical reactive transport in a complex fracture network to determine the key features that control dissolution. We consider topological, geometric and hydrological features of the fracture network to predict the remaining quartz in quasi-steady state. We found that the dissolution reaction rate constant of quartz and the distance to the primary sub-network in the fracture network are the two most important features controlling the amount of quartz remaining. This study is a first step towards characterizing the parameters that control carbon mineralization using an approach with integrates computational physics and machine learning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modeling Chemo-Mechanics with Electrolyte Infiltration to Quantify Degradation of Cathode Particles

One of the main goals in modeling lithium-ion batteries is to improve/predict longevity and resilience of new chemistries. To that end, this talk investigates the formation of stress-induced fracture within polycrystalline cathode particles and the impact on capacity loss. Physically based cathode aging dynamics is simulated in a 3D, continuum-level chemo-mechanical model. The model captures anisotropic Li diffusion within a single polycrystalline particle comprised of hundreds to thousands of randomly oriented grains. A recent addition to this model includes electrolyte infiltration, which occurs when the electrolyte seeps into surface cracks within the particle. The model predicts that particle fracture is primarily due to non-ideal grain interactions with slight dependence on high-rate charge demands. Essentially, when neighboring grains are misaligned, they expand a different rates relative to one another leading to high stresses and ultimately the formation of intraparticle cracks. The model predicts that small particles with large grains develop significantly less damage than larger particles with small grains. Finally, the model predicts most of the chemo-mechanical damage accumulates in the first charge after formation. This chemo-mechanical "damage saturation" effect indicates that initial particle fracture occurs within the first few cycles, while long-term cathode degradation is not solely chemo-mechanically induced. The principle contribution of this research is the use of an anisotropic chemo-mechanical model to test how particle geometry affect capacity fade, which predicts that particle size has a stronger effect on capacity fade than grain size and ultimately that small particles with large grains have the least capacity fade.

cathode↗
DOE OSTI2023

Enhancing Lithium-Ion Battery Aging Simulations: Coupling a High-Resolution, 3D, Grain-Scale Electromechanical Model to a Single-Particle Model

One of the main goals in modeling lithium-ion batteries is to improve/predict longevity and resilience of new chemistries. Unfortunately, this requires simulation of thousands of charge/discharge cycles, which can be rather time consuming depending on the fidelity of the simulation. The purpose of this talk is to discuss a new modeling framework that couples a high-resolution, continuous damage model (CDM) to a single particle model (SPM) resulting in a good combination of speed and fidelity. In previous work, a 3D, continuum-level chemo-mechanical model was developed to investigate cracking within a single cathode particle comprised of hundreds to thousands of randomly oriented grains. The CDM predicted that particle fracture is primarily due to non-ideal grain interactions with slight dependence on high-rate charge demands. Essentially, when neighboring grains were misaligned, they expanded different rates relative to one another leading to high stresses and ultimately the formation of intraparticle cracks. The model predicted that small particles with large grains develop significantly less damage than larger particles with small grains. Finally, the model predicted most of the chemo-mechanical damage accumulates in the first charge after formation. This chemo-mechanical "damage saturation" effect indicated that initial particle fracture occurs within the first few cycles, while long-term cathode degradation is not solely chemo-mechanically induced. This led to a need for simulating fatigue-like mechanism that degrade the battery over longer time scales. In order to reach the time scales, need to resolve fatigue-like degradation, the CDM needs to be complemented by a faster model. Therefore, recent efforts have been focused on using results from the CDM to inform parameters within the SPM. These parameters are homogenization factors that are associated with diffusion, particle radius, and/or exchange current density. By coupling the CDM to the SPM, the aging simulation is broken up into two domains: short-term and long-term degradation. The short-term degradation occurs over a single cycle and is handled by the CDM due of its high fidelity, but relatively expensive computational cost. Such mechanism include break-in crack caused by mismatches in grain orientation. The long-term degradation occurs over tens of cycles and is handled by the SPM due it's computational efficiency. Most of the fatigue-like mechanism fall into this category. The eventual goal of this modeling framework is to upscale to a psudo-2D model allowing for full electrode simulation, which are informed by high-fidelity grain-scale simulations.

electrochemistry↗
DOE OSTI2024

Learning nonlinear operators in latent spaces for real-time predictions of complex dynamics in physical systems

Abstract Predicting complex dynamics in physical applications governed by partial differential equations in real-time is nearly impossible with traditional numerical simulations due to high computational cost. Neural operators offer a solution by approximating mappings between infinite-dimensional Banach spaces, yet their performance degrades with system size and complexity. We propose an approach for learning neural operators in latent spaces, facilitating real-time predictions for highly nonlinear and multiscale systems on high-dimensional domains. Our method utilizes the deep operator network architecture on a low-dimensional latent space to efficiently approximate underlying operators. Demonstrations on material fracture, fluid flow prediction, and climate modeling highlight superior prediction accuracy and computational efficiency compared to existing methods. Notably, our approach enables approximating large-scale atmospheric flows with millions of degrees, enhancing weather and climate forecasts. Here we show that the proposed approach enables real-time predictions that can facilitate decision-making for a wide range of applications in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Conditions for stable operation of solid oxide electrolysis cells: oxygen electrode effects

Solid oxide electrolysis cells (SOECs) convert renewable electricity to fuels with efficiency substantially higher than other electrolysis technologies. However, questions remain regarding degradation mechanisms that limit SOEC long-term stability. One of the key degradation mechanisms is oxygen electrode delamination; although prior studies have improved the understanding of this mechanism, it is still difficult to predict how degradation depends on SOEC materials and operating conditions, i.e., temperature, voltage, and current density. Here we present a study aimed at developing a quantitative understanding of oxygen electrode delamination. Experimentally, a life test study of symmetric and full cells with yttria-stabilized zirconia (YSZ) electrolytes and Gd-doped ceria (GDC) barrier layers was done with three different perovskite oxygen electrode materials. Fracture was observed at the perovskite–GDC interface above a critical current density and below a critical operating temperature. A theory is presented that combines a calculation of the effective oxygen pressure across the electrolyte with an estimation of the pressure required for fracture. Fracture is correctly predicted for a critical oxygen partial pressure of ~7200 atm and an associated electrode overpotential of ~0.2 V, occurring at the electrode/GDC interface because of the relatively low perovskite fracture toughness. As a result, damage at the GDC/YSZ interface was also observed in some cases and explained by a peak in the oxygen pressure at this interface.

30 DIRECT ENERGY CONVERSION↗