Search NASA⌕ Search

SEARCH · Search NASA

Results for “Phase field 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

CO 2 -Responsive Fracturing Fluids for Enhanced Geothermal Systems (Final Report)

Our group has recently developed StimuFrac, a non-toxic stimuli-responsive fracturing fluid consisting of a CO 2 -reactive polymer which has shown at the lab-scale to consistently fracture rock cores at significantly lower net pressures in a range of representative geothermal pressure/temperature conditions. However, until now the mechanism/s responsible for more effective fracturing, of critical importance to optimize fracturing performance as well as strategize injection methodologies for field deployment, was not understood. In this document, we report (1) on the two main mechanisms responsible for fracturing rock at lower net pressures with StimuFrac; (2) the phase behavior of StimuFrac/CO 2 under geothermal wellbore conditions; and (3) based on high-temperature true triaxial stimulations, detailed evidence that StimuFrac/CO 2 is the best performing stimulation fluid under EGS T/P conditions as compared to water, waterless CO 2 , and CO 2 /water fracturing fluids. This is because i) it requires significantly lower volumes of CO 2 due to its reduced leak off into the formation by the crosslinked polymer solution; ii) large fractures can be generated reproducibly at both low and high CO 2 injection flow rates, and iii) the reversible (previously reported) viscosity increase of StimuFrac could be beneficial to transport proppants when they become available for EGS. These results were particularly evident for hot nearly dry rock as well as partially and nearly fully water-saturated granitic rock. GTO requires StimuFrac to be evaluated in fully water-saturated rock to determine whether the above-described performance applies under these conditions. Since (1) GTO considers there is no enough evidence for Sa different StimuFrac formulation to work below 90C (where full water-saturation in an open system is possible) and (2) no polyaxial loading frame larger than a few centimeters that can do hydraulic fracturing tests while maintaining the rock sample fully saturated with water at 200 C exists; PNNL concludes that the only way to determine if StimuFrac represents an advanced fracturing fluid alternative for EGS, is to perform a stimulation in an actual EGS reservoir.

15 GEOTHERMAL ENERGY↗
NASA NTRS2008

Effect of Microstructure on Time Dependent Fatigue Crack Growth Behavior In a P/M Turbine Disk Alloy

A study was conducted to determine the processes which govern hold time crack growth behavior in the LSHR disk P/M superalloy. Nineteen different heat treatments of this alloy were evaluated by systematically controlling the cooling rate from the supersolvus solutioning step and applying various single and double step aging treatments. The resulting hold time crack growth rates varied by more than two orders of magnitude. It was shown that the associated stress relaxation behavior for these heat treatments was closely correlated with the crack growth behavior. As stress relaxation increased, the hold time crack growth resistance was also increased. The size of the tertiary gamma' in the general microstructure was found to be the key microstructural variable controlling both the hold time crack growth behavior and stress relaxation. No relationship between the presence of grain boundary M23C6 carbides and hold time crack growth was identified which further brings into question the importance of the grain boundary phases in determining hold time crack growth behavior. The linear elastic fracture mechanics parameter, Kmax, is unable to account for visco-plastic redistribution of the crack tip stress field during hold times and thus is inadequate for correlating time dependent crack growth data. A novel methodology was developed which captures the intrinsic crack driving force and was able to collapse hold time crack growth data onto a single curve.

Telesman, Ignacy J.↗
DOE OSTI2023

Engineering the fracture resistance of 2H-transition metal dichalcogenides using vacancies: An in-silico investigation based on HRTEM images

Vacancy engineering of 2H-transition metal dichalcogenides (2H-TMDs) has recently attracted great attention due to its potential to fine-tune the phonon and opto-electric properties of these materials. From a mechanical perspective, this symmetry-breaking process typically reduces the overall crack resistance of the material and adversely affects its reliability. However, vacancies can trigger the formation of heterogeneous phases that synergistically improve fracture properties. In this study, using MoSe2 as an example, we characterize the types and density of vacancies that can emerge under electron irradiation and quantify their effect on fracture. Molecular dynamic (MD) simulations, employing a re-parameterized Tersoff potential capable of accurately capturing bond dissociation and structural phase changes, reveal that isolated transition metal monovacancies or chalcogenide divacancies tend to arrest the crack tip and hence enhance the monolayer toughness. In contrast, isolated chalcogenide monovacancies do not significantly affect toughness. The investigation further reveals that selenium vacancy lines, formed by high electron dose rates, alter the crack propagating direction and lead to multiple crack kinking. Using atomic displacements and virial stresses together with a continuum mapping, displacement, strain, and stress fields are computed to extract mechanistic information, e.g., conditions for crack kinking and size effects in fracture events. The study also reveals the potential of specific defect patterns, “vacancy engineering,” to improve the toughness of 2H-TMDs materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coupled momentum balance and phase-field solver with fenicsx module

Code solves momentum balance and phase-field equations simultaneously. The differential equations are solved on a discretized domain with appropriate boundary and initial conditions using finite element method. Primary purpose of the code is to simulate brittle fracture under dynamic loading. Constitutive equations are that of linear elasticity with degradation of stress due to fracture. Small strain formulation is used.

Zecevic, Milovan↗
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↗
NASA NTRS2004

Impact Processes in the Solar System

Our laboratory has previously conducted impact fracture and dynamic failure tests. Polanskey and Ahrens [1990] mapped the fractures from a series of laboratory craters (Fig. 1) and Ahrens and Rubin [ 1993] inferred that the usually further extending radial cracks resulted from tensional failure during the compression of the shock propagation. The radial spreading induced by the particle velocity field caused the stresses perpendicular to the shock front to become sufficiently large and tensile. This induces "radial fractures." The concentric fractures are attributed to the tensional failure occurring after the initial compressive phase. Upon radial propagation of the stress wave the negative tension behind the stress-wave front caused failure along the quasi-spherical concentric fractures. The near-surface and spall fractures are attributed to the fractures described by Melosh [1984]. These are activated by impact and can launch relatively unshocked samples of planetary surfaces to speeds exceeding escape velocity. In the case of Mars, some of these surface samples presumably become the SNC (Mars) meteorites.

Ahrens, Thomas J.↗
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 OSTI2020

Deformation Precursors to Catastrophic Failure in Rocks

Forecasting the timing of catastrophic failure, such as crustal earthquakes, has been a central concern for centuries. Such forecasting requires identifying signals that evolve or accelerate in the precursory phase leading to failure, and the subset of signals that may be detected in the crust. We develop machine learning models to predict the proximity of catastrophic failure in synchrotron X-ray tomography triaxial compression experiments on rocks using characteristics of evolving fracture networks. We then examine the characteristics that most strongly influence the model results, and thus may be considered the best predictors of the proximity of macroscopic failure. Finally, the resulting suite of predictive parameters underscores the importance of dilation in the precursory phase leading to catastrophic failure. The results indicate that the evolution of the strain energy density field may provide more robust predictions of the proximity of failure than other existing metrics of rock deformation.

58 GEOSCIENCES↗
DOE OSTI2023

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Utah FORGE: Phase 3 Magnetotelluric (MT) Data

New high-quality tensor MT data at 122 sites, including the vertical magnetic field and utilizing ultra-remote referencing, have been acquired over the Utah FORGE project area. The results will be used to delineate the densities of faults and fractures in crystalline basement rocks, to define the heat sources, and to derive baseline 3D resistivity structure for later MT monitoring of temporal changes in resistivity structure following well stimulation in the EGS reservoir. There are three files here related to Utah FORGE magnetotelluric (MT) data acquisition and processing. The FORGE MT EDIs zip file contains the observed MT responses in industry-standard EDI format. For each site, there is an EDI response file that utilized a local independent reference for noise cancellation, and a file that utilized a distant reference cancelling noise associated with the DC transmission line of the Delta IPP passing down the west side of Milford Valley. These two site files could be merged as appropriate. The FORGE Model Cell Center file contains the model volume of the 3D Forge MT inversion for characterizing the resistivity structure in the project area. It was derived using finite element inversion methodology described in Wannamaker et al, in the attached FORGE Phase3 Geophysics paper, from the MT observation EDI files. It is ASCII format (.dat) and entries are defined at top of the file in a simple x-y-z-Rho listing in UTM coordinates. The element layers drape the topography so the Rho value layers are not purely horizontal slices. This greatly simplifies the listing.

15 GEOTHERMAL ENERGY↗
NASA NTRS1973

Tensile properties to 650 C and deformation structures in a precipitation strengthened titanium-aluminum alloy

Appreciable strength levels were retained to 650 C in a Ti-10Al-1Si alloy aged in the (alpha + alpha sub 2) phase field to yield optimum room temperature strength and ductility. The aging treatment precipitated a uniform distribution of alpha sub 2-particles such that, at room temperature, dislocations bypassed instead of shearing the particles at low strains. Specimens fractured at room temperature exhibited fine uniform dimples even for those aging conditions that imparted no macroscopic ductility. The main crack appeared to propagate through the planar slip bands that had cut through the alpha sub 2-particles. A two-step aging process produced a higher volume fraction of bimodally distributed alpha sub 2-particles that led to higher strength levels at elevated temperatures. Both for the single size and the bimodal alpha sub 2-particle distributions, elevated-temperature deformation structures consisted mainly of planar slip bands that sheared through the alpha sub 2-particles.

Mendiratta, M. G.↗
DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

Oxidation and Associated Pore Structure Modification During Experimental Alteration of Granite

Weathering plays a crucial role in a number of environmental processes, and the microstructure and evolution of multi-scale pore space is a critically important factor in the weathering of geological formations. In igneous rocks the infiltration of meteoric water into initially relatively dry material can cause the rock to begin to disaggregate, increasing porosity and surface area, and allowing further disaggregation and weathering to occur. These processes, in turn, allow biota to colonize the rock, further enhancing the weathering rate. In some igneous rocks this may be driven by oxidation of primary minerals. One such mineral, biotite, has been repeatedly mentioned as a cause of cracking during oxidation. However, the scale-dependence of the processes by which this occurs are poorly understood. We cannot, therefore, accurately extrapolate laboratory reaction rates to the field in predictive numerical models. In order to better understand the effects of oxidation and test the hypothesis that fracture and disaggregation are initiated by swelling of oxidizing biotites, we reacted granite cores in a selenic acid-rich aqueous solution at 200°C for up to 438 days. Elevated temperatures and selenic acid were used to provide relatively fast reaction rates and highly oxidizing conditions in sealed reaction vessels. These experiments were analyzed using a combination of imaging, X-ray diffraction, Mössbauer spectroscopy, and small- and ultra-small angle neutron scattering to interrogate porosity and microfracture formation. The experimental results show little observable biotite swelling, but significant transport and growth of iron oxides and/or clays along grain boundaries throughout the sample. Significant increases in porosity were also observed at the sample rim, likely associated with feldspar alteration. Fractures and transport were observed throughout the core, suggesting that stresses due to crystallization pressures caused by the growing iron phases may be the initiating factors in granite weathering, possibly followed by biotite swelling after sufficient permeability is achieved.

Anovitz, Lawrence M.↗
DOE OSTI2026

Capturing the fractocohesive length scale in elastomers through a statistical mechanics-based gradient enhanced damage model

Here, this study aims to examine modeling flaw sensitivity in elastomers. The direct incorporation of polymer chain statistical mechanics considerations into a continuum stretch-based gradient-enhanced damage formulation, in turn, allows a representation of diffuse chain damage and fracture events that align with known micromechanical mechanisms. Through a series of numerical experiments, we simulate crack propagation and extract the fracture energy as an output of the model, while keeping track of the micromechanical signatures of diffuse chain damage that accommodate fracture propagation and eventually influence flaw sensitivity. Finally, by combining the fracture toughness and the work to rupture, we identify a fractocohesive length of the material, corresponding to the full width of the damage process zone. As the damage-to-fracture cascade in the proposed GED model is influenced by the introduction of a length scale associated with network imperfection and long-range load transfer, the emerging relationship of the two length scales is discussed, providing a potential link between microscopic damage mechanisms and the observed macroscopic fracture response.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Investigating Weathering of Basaltic Materials in Gale Crater, Mars: A Combined Laboratory, Modeling and Terrestrial Field Approach

Recent observations from Gale Crater, Mars document past aqueous alteration both in the formation of the Stimson sandstone unit, as well as in the formation of altered fractures within that unit. Geochemical and mineralogical data from Curiosity also suggest Fe-rich amorphous weathering products are present in most samples measured to date. Here we interpret conditions of possible past weathering in Gale Crater using a combination of field, laboratory, and modeling work. In order to better understand secondary Fe-rich phases on Mars, we are examining formation of weathering products in high Fe and Mg and low Al serpentine soils in the Klamath Mountains, CA. We have isolated potential weathering products from these soils, and are analyzing them using synchrotron μXRF and μXRD as well as FullPat for a direct comparison to analyses from Gale Crater. In order to interpret the implications of the persistence of potential secondary Fe-containing phases on Mars, we are also measuring the dissolution rates of the secondary weathering products allophane, Fe-rich allophane, and hisingerite. Ongoing dissolution experiments of these materials suggest that they dissolve significantly more rapidly than more crystalline secondary minerals with similar chemical compositions. Finally, to quantify the specific conditions of past aqueous alteration in Gale Crater we are performing reactive transport modeling of a range of possible past environmental conditions. Specifically, we are testing the conditions under which a Stimson unit-like material forms from a parent material similar to Rocknest or Bagnold eolian deposits, and the conditions under which observed altered fracture zones form from a Stimson unit-like parent material. Our modeling results indicate that the formation of the Stimson unit is consistent with leaching of an eolian deposit with a solution of pH = 6-8, and that formation of the altered fracture zones is consistent with leaching with a very acidic (pH = 2-3) high sulfate solution containing Ca. These results suggest circumneutral pH conditions during authigenesis or early diagenesis in the Stimson formation sediments followed by diagenetic alteration by very acidic solutions along fracture zones.

Hausrath, Elisabeth↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2022

High Performance Circuit Pastes for Solid Oxide Fuel Cell Applications: Final Technical Report

Using a combination of in-plane electrical conductivity measurements, electrical contact resistance measurements, tensile fracture tests, double shear lap fracture tests, rapid thermal cycling adhesion tests, in situ wetting angle measurements, 3D X-ray tomography, reduction-oxidation (RedOx) cycling, scanning electron microscopy (SEM), energy dispersive x-ray spectroscopy (EDAX), Density Functional Theory (DFT) modeling, Molecular Dynamics (MD) modeling, and Phase Field modeling, this project determined that a new brazing technique developed by the PIs (specifically the use of nickel particles to direct the wetting and spreading of molten silver) could be used to produce well-adhered, >97% dense, electrically-conductive circuits, current collectors, and/or electrical contacts on a variety of ceramic and ceramic-coated substrates (specifically lanthanum strontium manganite, yttria stabilized zirconia, sapphire, Cr2O3-passivated stainless steel, and Al2O3-passivated stainless steel). Compared to other silver-based alternatives, Ag:Ni led to less manufacturing defects and was more tolerant of the extreme environments (isothermal high-temperature aging) and occasional abuse (rapid thermal cycling and RedOx cycling) encountered during Solid Oxide Fuel Cell and/or Solid Oxide Electrolysis Cell operation. Further, in some situations (such as when to produce Ag:Ni electrical contacts to Al2O3-protected stainless steel) the Ni particles also improved operation by acting as chemical getters that removed electrically-insulating, surface-segregating substrate impurities (i.e. by lowering the contact resistance with the underlying substrate).

36 MATERIALS SCIENCE↗