Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic fracture”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Scaling of Imbibition Front Dynamics in Heterogeneous Porous Media

Here we report a quantitative study on the dynamics of countercurrent imbibition in fractured rock samples with the application of CO 2 geosequestration, as the rock matrix is saturated with CO 2 and brine exists in the adjacent microfracture. Direct pore-scale simulations of spontaneous and forced imbibition in a three-dimensional heterogeneous porous matrix enable us to capture the effect of microscopic properties on the imbibition dynamics as these pore-level descriptions are ignored in the existing coarse-grained models. We have established new scaling classes for the propagation of the center of mass of the imbibition front and fluid-fluid interface broadening. It has become evident that the wettability is a major factor which leads to anomalies in these scaling laws. While the governing capillarity under the water-wet condition induces fast roughening dynamics of the interface, the adverse effect of capillary pressure under the intermediate-wet condition leads to a damping effect on the imbibition front fluctuations.

58 GEOSCIENCES↗
DOE OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗
DOE OSTI2020

The emergence of small-scale self-affine surface roughness from deformation

Most natural and man-made surfaces appear to be rough on many length scales. There is presently no unifying theory of the origin of roughness or the self-affine nature of surface topography. One likely contributor to the formation of roughness is deformation, which underlies many processes that shape surfaces such as machining, fracture, and wear. Using molecular dynamics, we simulate the biaxial compression of single-crystal Au, the high-entropy alloy Ni 36.67 Co 30 Fe 16.67 Ti 16.67 , and amorphous Cu 50 Zr 50 and show that even surfaces of homogeneous materials develop a self-affine structure. By characterizing subsurface deformation, we connect the self-affinity of the surface to the spatial correlation of deformation events occurring within the bulk and present scaling relations for the evolution of roughness with strain. These results open routes toward interpreting and engineering roughness profiles.

42 ENGINEERING↗
DOE OSTI2024

SMART – A Comprehensive Research and Development Program to Demonstrate Application of Machine Learning for Supporting CCS Deployment

Presentation material for a paper presented at the GHGT-17 conference, Calgary, Canada, October 20-24, 2024. The objective of the US Department of Energy’s SMART Initiative, i.e., Science-informed Machine Learning (ML) for Accelerating Real-Time Decisions in Subsurface Applications, is to showcase how the utilization of ML can significantly improve efficiency and effectiveness of field-scale commercial carbon storage operations. This paper will present the results from the current phase of SMART (field deployment) for demonstrating the applicability of ML-based tools and workflows for: (a) virtual learning during the pre-injection permitting phase, (b) advanced storage reservoir imaging to better characterize fractures and faults, and (c) dynamic storage reservoir modelling and optimization to inform operational decision making and visualization of system evolution.

CO2 geologic storage↗
DOE OSTI2024

SMART – A Comprehensive Research and Development Program to Demonstrate Application of Machine Learning for Supporting CCS Deployment

The objective of the US Department of Energy’s SMART Initiative, i.e., Science-informed Machine Learning (ML) for Accelerating Real-Time Decisions in Subsurface Applications, is to showcase how the utilization of ML can significantly improve efficiency and effectiveness of field-scale commercial carbon storage operations. This paper will present the results from the current phase of SMART (field deployment) for demonstrating the applicability of ML-based tools and workflows for: (a) virtual learning during the pre-injection permitting phase, (b) advanced storage reservoir imaging to better characterize fractures and faults, and (c) dynamic storage reservoir modelling and optimization to inform operational decision making and visualization of system evolution.

Siriwardane, Hema↗
DOE OSTI2024

Vacancy-Mediated Increases in Brine–Salt Surface Energies

Salt formations have been explored for the permanent isolation of spent nuclear fuel based on their high thermal conductivity, self-healing nature, and low hydraulic permeability to brine flow. Vacancy defect concentrations in salt complicate fracture mechanics not driven by dislocation dynamics and can influence the resulting surface structure. Classical molecular dynamic simulations were used to simulate tensile testing of salt crystals (halite) with vacancy defect concentrations of up to 0.5 defects/nm 3 . Increasing defect concentrations resulted in a decrease in ultimate tensile strength and fracture surface energies, driven by increased surface roughness rather than changes in the amount of surface area. Brine–salt surface energies of the fractured surfaces were 0.22 to 0.26 J/m 2 , significantly higher than values reported for atomically flat (100) surfaces at the same brine composition. This change in surface energy increased the brine–salt dihedral angle by ~27°. The dihedral angle threshold for percolation in salt is 60°, and a 27° increase due to rough fracture surfaces identifies a reduction in porosity percolation and a decrease in salt permeability. Therefore, bedded salt and salt domes may be even more stable than those previously predicted from dihedral angle calculations.

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

Fracture mechanisms of sodium silicate glasses

Reactive classical molecular dynamics simulations of sodium silicate glasses, xNa 2 O–(100 – x)SiO 2 (x = 10–30), under quasi-static loading, were performed for the analysis of molecular scale fracture mechanisms. Mechanical properties of the sodium silicate glasses were consistent with experimentally reported values, and the amount of crack propagation varied with reported fracture toughness values. The most crack propagation occurred in NS20 systems (20-mol% Na 2 O) compared with the other simulated compositions. Dissipation via two mechanisms, the first through sodium migration as a lower activation energy process and the second through structural rearrangement as a higher activation energy process, was calculated and accounted for the energy that was not stored elastically or associated with the formation of new fracture surfaces. Here, a correlation between crack propagation and energy dissipation was identified, with systems with higher crack propagation exhibiting less energy dissipation. Sodium silicate glass compositions with lower energy dissipation also exhibited the most sodium movement and structural rearrangement within 10 Å of the crack tip during loading. Therefore, high sodium mobility near the crack tip may enable energy dissipation without requiring formation of structural defects. Therefore, the varying mobilities of the network modifiers near crack tips influence the brittleness and the crack growth rate of modified amorphous oxide systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A New Simplified Discrete Fracture Model for Shearing of Intersecting Fractures and Faults

Abstract Shearing of fractures and faults is important because it can result in permeability change or even induce seismicity—both are keys for efficient and safe energy recovery and storage in Earth systems. Quantitative analysis of shearing of intersecting fractures and faults is challenging because it can involve dynamic frictional contacts that are complicated by deformation of the rock matrix. To predict the shearing of intersecting fractures/faults, we attempt to answer the question of how intersections impact the shearing of a fracture network and whether we can simplify the description as compared to classical discrete fracture network (DFN) models. To answer these questions, we conducted a series of numerical simulations on scenarios for variable numbers of intersecting fractures. All these examples yield consistent results: the results of using DFNs are consistent with those of using hypothetical major paths. This leads to a new model, which we name simplified discrete fracture network model , to analyze shearing of intersecting fractures/faults using major path(s). We found that the intersections of fractures do not fundamentally change the shearing of two intersecting fractures if the intersecting angles are small. Furthermore, increasing the number of fractures/faults may relax the stress as more fractures/faults become available for shearing and distributing the stress. The simplified DFN model, which can capture efficiently the shearing behavior of each major paths from a large number of intersecting fractures/faults, will be a promising conceptual model that is complementary to existing equivalent continuum and discrete fracture models to analyze shearing of intersecting fractures/faults.

58 GEOSCIENCES↗
DOE OSTI2021

An asymptotically compatible treatment of traction loading in linearly elastic peridynamic fracture

Meshfree discretizations of state-based peridynamic models are attractive due to their ability to naturally describe fracture of general materials. However, two factors conspire to prevent meshfree discretizations of state-based peridynamics from converging to corresponding local solutions as resolution is increased: quadrature error prevents an accurate prediction of bulk mechanics, and the lack of an explicit boundary representation presents challenges when applying traction loads. Herein, we develop a reformulation of the linear peridynamic solid (LPS) model to address these shortcomings, using improved meshfree quadrature, a reformulation of the nonlocal dilatation, and a consistent handling of the nonlocal traction condition to construct a model with rigorous accuracy guarantees. In particular, these improvements are designed to enforce discrete consistency in the presence of evolving fractures, whose a priori unknown location render consistent treatment difficult. In the absence of fracture, when a corresponding classical continuum mechanics model exists, our improvements provide asymptotically compatible convergence to corresponding local solutions, eliminating surface effects and issues with traction loading which have historically plagued peridynamic discretizations. When fracture occurs, our formulation automatically provides a sharp representation of the fracture surface by breaking bonds, avoiding the loss of mass. We provide rigorous error analysis and demonstrate convergence for a number of benchmarks, including manufactured solutions, free-surface, nonhomogeneous traction loading, and composite material problems. Finally, we validate simulations of brittle fracture against a recent experiment of dynamic crack branching in soda-lime glass, providing evidence that the scheme yields accurate predictions for practical engineering problems.

42 ENGINEERING↗
DOE OSTI2026

Spontaneous crack healing in calcite reveals the influence of dynamic strain evolution and surface chemistry

The mechanics of fracture healing in calcite remain poorly constrained yet are fundamental to managing fluid transport in geothermal reservoirs and hydrocarbon systems. Here, we apply microfocused synchrotron Laue X-ray diffraction and infrared spectroscopy to investigate subcritical crack healing in a 1 mm-thick calcite crystal subjected to controlled loading in a double-torsion device. Over a 44-hour period following load removal, we map the evolution of residual strain fields surrounding the crack tip and observe a progressive increase in compressive strain perpendicular to the crack plane accompanied by infrared spectroscopic signatures that reveal enhanced accumulation of water at the healed interface. The correlation between strain evolution and surface chemistry suggests that spontaneous crack healing in calcite is driven by dynamic anelastic relaxation coupled with irreversible fluid-mineral interactions. These findings offer insight into time-dependent crack closure processes in carbonates and highlight the role of chemically-mediated plasticity in subsurface fracture evolution.

materials science↗
DOE OSTI2022

Investigation of fracture in porous materials: a phase-field fracture study informed by ReaxFF

Microscopic features (e.g, pore shapes, sizes, and distribution) in porous material substantially affect the overall mechanical properties such as stiffness and strength. In turn, these material properties determine the macroscopic behaviors of fracture in the porous material. In certain cases, macroscopic properties can be derived from the porous skeleton and void ratio (i.e., porosity), but in many other cases, such derivation is a challenging task. This paper presents a numerical investigation of microporosity and micropore shapes effect on the macrofracture behavior in porous amorphous silica. For this study, we extend the recently-proposed combined molecular dynamic (MD) and phase-field (PF) fracture modeling approach by including different pore shapes in the atomistic domain. In the MD simulations, we adopt ReaxFF to evaluate the material properties, where four different micropore cases are considered. Based on the material properties derived from MD simulations, the macrofracture propagation of porous media is studied using hybrid PF simulation. In the characterization of the pore structure, the concept of pore ligament is proposed to relate the pore shape and the critical energy release rate. Two classical fracture problems were used to evaluate the effect of pore shape on the macrofracture behavior. The results of the case study show that although the micropore shapes change the macrofracture behaviors, these effects vary with the geometry and loading conditions of macroscopic boundary value problems. Furthermore, the case study also shows that the influence of micropore structure can be captured at the macroscopic level through the material properties derived from the MD simulations.

42 ENGINEERING↗
DOE OSTI2022

Role of the Mineral in the Self-Healing of Cracks in Human Enamel

Human enamel is an incredibly resilient biological material, withstanding repeated daily stresses for decades. The mechanisms behind this resilience remain an open question, with recent studies demonstrating a crack-deflection mechanism contributing to enamel toughness and other studies detailing the roles of the organic matrix and remineralization. Here, we focus on the mineral and hypothesize that self-healing of cracks in enamel nanocrystals may be an additional mechanism acting to prevent catastrophic failure. To test this hypothesis, we used a molecular dynamics (MD) approach to compare the fracture behavior of hydroxyapatite (HAP) and calcite, the main minerals in human enamel and sea urchin teeth, respectively. In this work, we find that cracks heal under pressures typical of mastication by fusion of crystals in HAP but not in calcite, which is consistent with the resilience of HAP enamel that calcite teeth lack. Scanning transmission electron microscopy (STEM) images of structurally intact ("sound") human enamel show dashed-line nanocracks that resemble and therefore might be the cracks healed by fusion of crystals produced in silico . The fast, self-healing mechanism shown here is common in soft materials and ceramics but has not been observed in single crystalline materials at room temperature. The crack self-healing in sound enamel nanocrystals, therefore, is unique in the human body and unique in materials science, with potential applications in designing bioinspired materials.

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

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗
DOE OSTI2026

Insights into elevated temperature tensile deformation mechanisms and kink banding in additively manufactured tungsten

Additive manufacturing (AM) potentially enables fabrication and repair of plasma facing components (PFCs) of tungsten. However, deformation mechanisms in AM-fabricated tungsten under service-relevant thermomechanical conditions remain unknown majorly due to challenges in achieving crack-free W. This study reveals the deformation mechanisms operative under tension at 800 °C and 1200 °C in the electron beam melting powder bed fusion (EBM-PBF) fabricated W. Textured columnar grains with mixed <001 >/<111 >|| build-direction were observed. Strong anisotropy in tensile properties persisted across test temperatures, and all specimens exhibited extensive deformation-induced banding phenomena. Grain boundary character analysis of the deformed specimens indicated that these deformation bands were kink bands with tilt grain boundaries prominently present at the band-matrix interface. Interestingly, the material within the kink bands rotated toward higher resolved shear stress, providing insights into the origins of kink band formation in body-centered cubic refractory metals. Intragranular misorientation axis analysis revealed that plasticity was dominated by {110} <111 > slip at 800 °C, whereas additional (213)[11⁢̄1] and (112)[11⁢̄1] slip systems activated at 1200 °C. A dense low-angle boundary networks observed near fractured surface indicated the onset of dynamic recrystallization. Results reveal plasticity governing mechanisms in AM-fabricated W under power plant-relevant conditions, and provide insights into kink band attributes in refractory metals.

Mayes, Riley [ORNL] (ORCID:000900089307010X)↗
DOE OSTI2026

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗
DOE OSTI2025

Deep Learning–Assisted Multiobjective Optimization of Geological CO 2 Storage Performance under Geomechanical Risks

In geological CO 2 storage, designing the optimal well control strategy for CO 2 injection to maximize CO 2 storage while minimizing the associated geomechanical risks is not trivial. This challenge arises due to pressure buildup, CO 2 plume migration, the highly nonlinear nature of geomechanical responses to rock-fluid interaction, and the high computational cost associated with coupled flow and geomechanics simulations. In this paper, we introduce a novel optimization framework to address these challenges. The optimization problem is formulated as follows: maximize total CO 2 storage while minimizing geomechanical risks by adjusting the injection schedules within bounded constraints. The geomechanical risks are primarily driven by injection-induced pressure build-up, which is characterized by ground displacement and the induced microseismicity. We used the Fourier neural operator (FNO)-based deep learning model to construct surrogate models, replacing the time-consuming coupled flow and geomechanics simulations for evaluating the aforementioned objective functions. The developed surrogate models have been incorporated into a multiobjective optimization framework through a genetic algorithm to reduce the computational burden. The proposed optimization framework reduces the computational cost from approximately 2,400 hours, when using objective function evaluations based on physics-based simulations, to around 20 minutes. A set of Pareto-optimal solutions of the proposed workflow yields nontrivial optimal decisions, reducing the microseismicity potential and the vertical displacement. This Pareto front highlights the optimal trade-offs between CO 2 storage amount, safety, and ground displacement, emphasizing the need for careful optimization and management of injection strategies to achieve a balanced outcome. The novelty of this work is twofold. First, we demonstrate the importance of incorporating the minimization of the geomechanical risks as objective functions into the CO 2 storage optimization workflow to mitigate the potential risk of induced microseismicity and ground displacement. Second, we leverage the FNO-based surrogate models to optimize a real-field CO 2 storage operation.

42 ENGINEERING↗
DOE OSTI2025

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗