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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS1987

Structural Integrity and Durability of Reusable Space Propulsion Systems

A two-day conference on the structural integrity and durability of reusable space propulsion systems was held on May 12 and 13, 1987, at the NASA Lewis research Center. Aerothermodynamic loads; instrumentation; fatigue, fracture, and constitutive modeling; and structural dynamics were discussed.

Source record↗
NASA NTRS1999

Delamination Modeling of Composites for Improved Crash Analysis

Finite element crash modeling of composite structures is limited by the inability of current commercial crash codes to accurately model delamination growth. Efforts are made to implement and assess delamination modeling techniques using a current finite element crash code, MSC/DYTRAN. Three methods are evaluated, including a straightforward method based on monitoring forces in elements or constraints representing an interface; a cohesive fracture model proposed in the literature; and the virtual crack closure technique commonly used in fracture mechanics. Results are compared with dynamic double cantilever beam test data from the literature. Examples show that it is possible to accurately model delamination propagation in this case. However, the computational demands required for accurate solution are great and reliable property data may not be available to support general crash modeling efforts. Additional examples are modeled including an impact-loaded beam, damage initiation in laminated crushing specimens, and a scaled aircraft subfloor structures in which composite sandwich structures are used as energy-absorbing elements. These examples illustrate some of the difficulties in modeling delamination as part of a finite element crash analysis.

Fleming, David C.↗
NASA NTRS1999

Delamination Modeling of Composites for Improved Crash Analysis

Finite element crash modeling of composite structures is limited by the inability of current commercial crash codes to accurately model delamination growth. Efforts are made to implement and assess delamination modeling techniques using a current finite element crash code, MSC/DYTRAN. Three methods are evaluated, including a straightforward method based on monitoring forces in elements or constraints representing an interface; a cohesive fracture model proposed in the literature; and the virtual crack closure technique commonly used in fracture mechanics. Results are compared with dynamic double cantilever beam test data from the literature. Examples show that it is possible to accurately model delamination propagation in this case. However, the computational demands required for accurate solution are great and reliable property data may not be available to support general crash modeling efforts. Additional examples are modeled including an impact-loaded beam, damage initiation in laminated crushing specimens, and a scaled aircraft subfloor structures in which composite sandwich structures are used as energy-absorbing elements. These examples illustrate some of the difficulties in modeling delamination as part of a finite element crash analysis.

Fleming, David C.↗
NASA NTRS2001

Microstructure, Mechanical Properties, Hot-Die Forming, and Joining of 47XD Gamma TiAl Rolled Sheets

The microstructure and mechanical properties, along with the hot-die forming and joining of Ti-47Al-2Nb-2Mn-0.8 vol% TiB, sheets (known as 47XD), produced by a low-cost rolling process, were evaluated. A near-gamma microstructure was obtained in the as-rolled condition. The microstructures of heat-treated sheets ranged from a recrystallized equiaxed near-gamma microstructure at 1,200 to 1,310 C, to a duplex microstructure at 1,350 C, to a fully lamellar microstructure at 1,376 C. Tensile behavior was determined for unidirectionally rolled and cross-rolled sheets for room temperature (RT) to 816 C. Yield stress decreased gradually with increasing deformation temperature up to 704 C; above 704 C, it declined rapidly. Ultimate tensile strength exhibited a gradual decrease up to 537 C before peaking at 704 C, followed by a rapid decline at 816 C. The modulus showed a gradual decrease with temperature, reaching approximately 72% of the RT value at 816 C. Strain to failure increased slowly from RT to 537 C; between 537 C and 704 C it exhibited a phenomenal increase, suggesting that the ductile-brittle transition temperature was below 704 C. Fracture mode changed from transgranular fracture at low temperature, to a mixture of transgranular and intergranular fracture at intermediate temperature, to ductile fracture at 816 C, coupled with dynamic recrystallization at large strains. Creep rupture response was evaluated between 649 and 816 C over the stress range of 69 to 276 MPa. Deformation parameters for steady-state creep rate and time-to-rupture were similar: activation energies of approximately 350 kJ/mol and stress exponents of approximately 4.5. Hot-die forming of sheets into corrugations was done at elevated temperatures in vacuum. The process parameters to join sheets by diffusion bonding and brazing with TiCuNi 70 filler alloy were optimized for test coupons and successfully used to fabricate large truss-core and honeycomb structures. Nondestructive evaluation methods, e.g., ultrasonic C-scans and thermography along with metallography, were used to characterize bond quality. Microstructural evaluation during heat treatment, identification of phases at the braze/matrix interface, determination of shear strengths of brazed joints, and deformation mechanisms during tensile and creep processes will be discussed.

Das, G.↗
DOE OSTI2021

Investigating Fracture Network Deformation Using Noble Gas Release

We investigate deformation mechanics of fracture networks in unsaturated fractured rocks from subsurface conventional detonation using dynamic noble gas measurements and changes in air permeability. We dynamically measured the noble gas isotopic composition and helium exhalation of downhole gas before and after a large subsurface conventional detonation. These noble gas measurements were combined with measurements of the subsurface permeability field from 64 discrete sampling intervals before and after the detonation and subsurface mapping of fractures in borehole walls before well completion. We saw no observable increase in radiogenic noble gas release from either an isotopic composition or a helium exhalation point of view. Large increases in permeability were observed in 13 of 64 discrete sampling intervals. Of the sampling intervals which saw large increases in flow, only two locations did not have preexisting fractures mapped at the site. Given the lack of noble gas release and a clear increase in permeability, we infer that most of the strain accommodation of the fractured media occurred along previously existing fractures, rather than the creation of new fractures, even for a high strain rate event. These results have significant implications for how we conceptualize the deformation of rocks with fracture networks above the percolation threshold, with application to a variety of geologic and geological engineering problems.

Gardner, W. Payton↗
DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗
DOE OSTI2023

Effects of particle size, shape and loading rate on the normal compaction of an advanced granular ceramic

Compaction behavior of granular materials is influenced by strain rate, particle size, and shape. In this report, boron carbide powders with different particle sizes under uni-axial strain conditions are studied using quasi-static compression, dynamic Kolsky bar experiments and normal plate impact. A rounded powder is compacted to investigate the effect of particle shape. The normal plate impact technique is an excellent tool in characterization of powder compaction behavior up to strain rates of ~10 5 s -1 . From our experiments, granular boron carbide shows a highly compressible behavior with significant volume compaction. Constitutive responses are obtained for four powders. Particle fracture is identified as key deformation mechanism. Dynamic loading introduces more particle fragmentation than quasi-static loading. Morphological characterization of particle shapes shows that the deformation from powder compaction alters the particle shape distribution, which is also rate-dependent. Particle size, shape and strain rate effects on the normal stress are discussed accordingly.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Revealing complex subsurface dynamics with continuous seismic monitoring: Observations using distributed acoustic sensing and surface orbital vibrators during hydraulic fracturing

Understanding hydraulic fracturing is crucial to improving the stimulation of unconventional reservoirs and increasing fluid production. This study develops a novel seismic monitoring technology using distributed acoustic sensing (DAS) and surface orbital vibrators (SOV) to capture fracture seismic response and mechanical properties at high temporal intervals. We analyze continuous time-lapse vertical seismic profiling (VSP) data acquired every hour during the first nine days of treatment of an unconventional reservoir in the Austin Chalk/Eagle Field Laboratory. The VSP data contain clear seismic signals scattered from the activated fractures. The spatiotemporal changes of the fracture reflectivity revealed by the SOV/DAS data correlate well with the observations of fracture locations inferred from low-frequency DAS data. These results capture the fracture opening and closure processes, as well as highlight potential prestage activations of the fractures due to hydraulic connectivity with preexisting fracture systems. Therefore, analysis of the presented data set provides a unique opportunity to understand fracture initiation and subsequent evolution, not only in the context of unconventional resources but also in enhanced geothermal systems.

Correa, Julia↗
DOE OSTI2023

Gravity-driven controls on fluid and carbonate precipitation distributions in fractures

Abstract Many challenges related to carbon-dioxide ( $$\hbox {CO}_2$$ CO 2 ) sequestration in subsurface rock are linked to the injection of fluids through induced or existing fracture networks and how these fluids are altered through geochemical interactions. Here, we demonstrate that fluid mixing and carbonate mineral distributions in fractures are controlled by gravity-driven chemical dynamics. Using optical imaging and numerical simulations, we show that a density contrast between two miscible fluids causes the formation of a low-density fluid runlet that increases in areal extent as the fracture inclination decreases from 90 $$^\circ$$ ∘ (vertical fracture plane) to 30 $$^\circ$$ ∘ . The runlet is sustained over time and the stability of the runlet is controlled by the gravity-driven formation of 3D vortices that arise in a laminar flow regime. When homogeneous precipitation was induced, calcium carbonate covered the entire surface for horizontal fractures (0 $$^\circ$$ ∘ ). However, for fracture inclinations greater than 10 $$^\circ$$ ∘ , the runlet formation limited the areal extent of the precipitation to less than 15% of the fracture surface. These insights suggest that the ability to sequester $$\hbox {CO}_2$$ CO 2 through mineralization along fractures will depend on the fracture orientation relative to gravity, with horizontal fractures more likely to seal uniformly.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Benchmark study of a new simplified DFN model for shearing of intersecting fractures and faults

It is challenging to quantitatively predict shearing of intersecting fractures/faults because of dynamic frictional contacts accompanied by possible nonlinear rock deformation. To address such challenges, a new conceptual model—the simplified DFN model—was proposed and validated by Hu et al. 46 to use major paths (MPs) to represent complicated DFNs for calculation of shearing. In this work, we conducted a benchmark study for three examples that involve different levels of complexity of intersecting fractures, and correspondingly different numbers of MPs. The codes and software that were used in the benchmark cover a range of continuum, discontinuum and hybrid numerical methods: NMM (LBNL), FLAC3D (LBNL), GBDEM (KIGAM), FRACOD (DynaFrax), and CASRock (CAS). The general consistency between DFN and MP cases as predicted by all the codes/software demonstrates that major paths can be used to simplify the geometry of DFNs in a wide range of software. Disagreement in results made by some software and potential future improvements are discussed. We show that (1) shearing of one or multiple major fractures can be reduced if there are multiple smaller intersecting fractures in that area, which is a useful basis for understanding and controlling induced seismicity and merits further analysis, and (2) the agreement achieved in the benchmark examples provide confidence that the simplified DFN model is a promising conceptual model that can be used for different types of numerical approaches and software for simplifying the analysis of the shearing of intersecting fractures and faults.

58 GEOSCIENCES↗
DOE OSTI2024

Role of water in fracture of modified silicate glasses

Decarbonizing the glass industry requires alternative melting technology, as current industrial melting practices rely heavily on fossil fuels. Hydrogen has been proposed as an alternative to carbon-based fuels, but the ensuing consequences on the mechanical behavior of the glass remain to be clarified. A critical distinction between hydrogen and carbon-based fuels is the increased generation of water during combustion, which raises the equilibrium solubility of water in the melt and alters the behavior of the resulting glass. A series of five silicate glasses with 80% silica and variable [Na 2 O]/([H 2 O] + [Na 2 O]) ratios were simulated using molecular dynamics to elucidate the effects of water on fracture. Several fracture toughness calculation methods were used in combination with atomistic fracture simulations to examine the effects of hydroxyl content on fracture behavior. Here, this study reveals that the crack propagation pathway is a key metric to understanding fracture toughness. Notably, the fracture propagation path favors hydrogen sites over sodium sites, offering a possible explanation of the experimentally observed effects of water on fracture properties.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Multiscale Modeling of Fracture Strength in Fibrous Thermal Protection System Materials

This work presents a multiscale modeling approach to predict the fracture strength of fibrous Thermal Protection System (TPS) materials. The model assumes that system failure is initiated at the joints between individual fibers. We investigated three distinct TPS compositions: amorphous silica, alumina and aluminosilicate fibers. Molecular dynamics (MD) simulations were employed to determine the fracture strength values of these fiber joints for both material systems. These fracture strength values were then integrated into simulations of 3D randomly populated fiber structures, where tensile load transfer occurs through the fiber joints. These microscale properties are upscaled through a renormalization approach [1] to predict macroscale tensile strength of 3D random fiber networks, accounting for joint-dominated failure and effective load-bearing area. The study concludes by demonstrating the resulting strength variation as a function of material composition, fiber density, and morphology. We also show validation of results by comparing them against explicit fiber finite element (FE) modeling [2] where fiber joint fracture is represented by cohesive elements.

Jaehyun Cho↗