Search NASA⌕ Search

SEARCH · Search NASA

Results for “fractured domains”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A probabilistic-phase field model for the fracture of brittle materials

Abstract We develop a computational method to determine the failure probability of brittle materials under general mechanical loading conditions. The method is a combination of two parts: (1) numerical simulations of materials with multiple cracks using phase field theory, where the complete fracture process is viewed as ‘damage percolation’ along critical paths or clusters of cracks, rather than the traditional weak-link failure mechanism of Weibull, and (2) an extension of the Batdorf statistical theory of fracture to finite domains, where it is implemented within the finite element framework. The results of phase-field simulations at the ‘percolation threshold’ are used as failure data in the Batdorf theory to determine the overall probability of failure. The input to this approach is the size distribution of cracks in a pristine material. An example is shown, where alumina samples that were previously tested by Abe and coworkers (Abeet al2003J. Am. Ceram. Soc.861019–21) in four-point loading are compared to the results of our numerical simulations. The approach developed here has the advantage of being extendable to more complex thermomechanical loading.

Materials Science↗
DOE OSTI2022

Heterogeneous Ice Growth in Micron-Sized Water Droplets Due to Spontaneous Freezing

Understanding how ice nucleates and grows into larger crystals is of crucial importance for many research fields. The purpose of this study was to shed light on the phase and structure of ice once a nucleus is formed inside a metastable water droplet. Wide-angle X-ray scattering (WAXS) was performed on micron-sized droplets evaporatively cooled to temperatures where homogeneous nucleation occurs. We found that for our weak hits ice grows more cubic compared to the strong hits that are completely hexagonal. Due to efficient heat removal caused by evaporation, we propose that the cubicity of ice at the vicinity of the droplet’s surface is higher than for ice formed within the bulk of the droplet. Moreover, the Bragg peaks were classified based on their geometrical shapes and positions in reciprocal space, which showed that ice grows heterogeneously with a significant population of peaks indicative of truncation rods and crystal defects. Frequent occurrences of the (100) reflection with extended in-planar structure suggested that large planar ice crystals form at the droplet surface, then fracture into smaller domains to accommodate to the curvature of the droplets. Planar faulting due to misaligned domains would explain the increased cubicity close to the droplet surface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

Evidence for divergent plate-boundary characteristics and crustal spreading on Venus

Detailed examination of the topography and morphology of western Aphrodite Terra reveals numerous features that are similar to terrestrial divergent plate-boundary characteristics. Individual domains between fracture-zone-like discontinuities contain a variety of bilaterally symmetrical topographic elements that suggest that topographic features have been created at rise crests, rifted and separated, and moved laterally to their present symmetrical positions. The topographic and morphologic similarities, together with strikingly similar mirror-image map patterns on both sides of the rise axis, suggest that western Aphrodite Terra shares the characteristics of oceanic divergent plate boundaries, and is the site of crustal spreading on Venus. Topographic profiles are consistent with spreading rates of the order of several centimeters per year.

Head, James W., III↗
NASA NTRS1992

Phase transformations in 40-60-GPa shocked gneisses from the Haughton Crater (Canada): An Analytical Transmission Electron Microscopy (ATEM) study

In order to better understand phase transformations, chemical migration, and isotopic disequilibrium in highly shocked rocks, we have performed a microprobe and an ATEM study on gneisses shocked up to 60 GPa from the Haughton Crater. This study reveals the following chemical and structural characteristics: (1) SiO2 dominant areas are formed by a mixture of pure SiO2 polycrystalline quartz identified by electron diffraction pattern and chemical analysis and a silica-rich amorphous phase containing minor amounts of aluminium, potassium, and iron; (2) Areas with biotitelike composition are formed by less than 200-nm grains of iron-rich spinels embedded in a silica-rich amorphous phase that is very similar to the one described above; (3) Layers with feldsparlike composition are constituted by 100-200-nm-sized alumina-rich grains (the indexation of the crystalline structure is under progress) and the silica-rich amorphous phase; (4) Zones characterized by the unusual Al/Si ratio close to 1 are formed by spinel grains (200-nm-sized) embedded in the same silica-rich amorphous phase; and (5) The fracturated sillimanites contain domains with a lamellar structure, defined by the intercalation of 100-nm-wide lamellae of mullite crystals and of a silica-rich amorphous phase. These mullite crystals preserved the crystallographical orientation of the preshock sillimanite. All compositional domains, identified at the microprobe scale, can thus be explained by a mixture in different proportion between the following phases: (1) a silica-rich amorphous phase, with minor Al and K; (2) quartz crystals; (3) spinel crystals and alumina-rich crystals; (4) sillimanite; and (5) mullite. Such mixtures of amorphous phases and crystals in different proportions explain disturbed isotope systems in these rocks and chemical heterogeneities observed on the microprobe.

Martinez, I.↗
DOE OSTI2020

Discrete Fracture Network Modeling to Estimate Upscaled Parameters for the Topopah Spring, Lava Flow, and Tiva Canyon Aquifers at Pahute Mesa, Nevada National Security Site

This report describes the results of Discrete Fracture Network (DFN) simulations for the Topopah Spring Aquifer (TSA), Lava Flow Aquifer, and Tiva Canyon Aquifer (TCA), at Pahute Mesa on the Nevada National Security Site (NNSS), formerly the Nevada Test Site. The research focuses on calculating upscaled groundwater flow and contaminant transport parameters using DFNs generated according to fracture characteristics observed in the TSA, LFA and TCA at Pahute Mesa. The highly fractured and heterogeneous nature of these aquifers makes them candidates for stochastic DFN modeling of radionuclide transport on a small scale with subsequent upscaling. One hundred independent DFN realizations are generated for each aquifer, and the upscaled parameters for continuum simulations of subsurface flow and transport in fractured media at Pahute Mesa are calculated. Our goal is to implement a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, such as channelization of flow and transport along a few well connected, large fractures. Additionally, to simulate advective and advective-diffusive transport through the fracture networks, the Time Domain Random Walk (TDRW) approach is applied to account for matrix diffusion into a finite half-space. Moreover, a novel approach to calculate dynamic (active) fracture surface area to reflect flow channeling is implemented. This work will improve the representation of radionuclide transport processes in largescale, regulatory-focused models by providing estimates of hard-to-measure flow and contaminant transport parameters at large scales. In this report, we (1) show recent results of flow and transport simulations on multiple DFN realizations of the TSA, LFA, TCA; (2) discuss the resulting distributions of estimated upscaled parameters; (3) describe the estimation of upscaled parameters for an equivalent parallel-plate continuum model and (4) present a comparison between simulated transport from the equivalent continuum model and an actual DFN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Particle-Based Conditional Sampling Scheme for the Simulation of Transport in Fractured Rock With Diffusion Into Stagnant Water and Rock Matrix

In situ experiments and field-scale characterization studies have pointed out that, in fractured crystalline media, groundwater flow is highly channelized. This implies that, at the scale of a single fracture, only part of the fracture surface area is in contact with flowing water, while the rest of in-plane water is essentially stagnant and can be accessed by solutes via molecular diffusion. Despite their importance for contaminant retention, to date, there are no numerical or analytical approaches that could be used to assess the implication of stagnant water zones on solute transport in realistic large-scale Discrete Fracture Network-based models. In this study, we present an efficient and flexible algorithm for the simulation of transport in fractured rock with diffusion into stagnant water and rock matrix. The algorithm is a generalization of a previously developed numerical framework for time domain particle tracking in sparsely fractured rock. The key of the generalization is that total time in fracture (τ f ) is first evaluated using a Monte Carlo sampling and then a second sampling is performed conditioned on τ f . The algorithm has been successfully validated against existing independent solutions and the implication of diffusion into stagnant water and secondary diffusion into the matrix has been assessed for a realistic modeling scenario. The results show that, due to diffusion into stagnant water, contaminants are more strongly retarded. This increased retention is more significant for sorbing species, as a larger number of sorption sites is accessible. A high sensitivity to the flowing channel/stagnant water zone geometry has also been observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mapdfn.py

SAND2021-4824 O mapdfn.py takes the output of dfnWorks describing a discrete fracture network (DFN) and converts the DFN to an equivalent porous medium representation. This is done by looping through cells in a 3-D model domain, tagging those intersected by fractures, and calculating effective properties. Coordinate transformations are imported from an open source script (transformations.py) developed at the University of California. An associated driver script, mapdfn2pflotran.py, calls the functions in mapdfn.py and writes the mapdfn.py output to hdf5 files for use with PFLOTRAN. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-4824 O

Stein, Emily↗
DOE OSTI2020

A Method to Represent a Well in a Three‐Dimensional Discrete Fracture Network Model

Abstract In discrete fracture network (DFN) modeling, fractures are randomly generated and placed in the model domain. The rock matrix is considered impermeable. Small fractures and isolated fractures are often ignored to reduce computational expense. As a result, the rock matrix between fractures could be large and intersections may not be found between a well introduced in the model and the hydraulically connected fracture networks (fracture backbones). To overcome this issue, this study developed a method to conceptualize a well in a three‐dimensional (3D) DFN using two orthogonal rectangular fractures oriented along the well's axis. Six parameters were introduced to parameterize the well screen and skin zone, and to control the connectivity between the well and the fracture backbones. The two orthogonal fractures were discretized using a high‐resolution mesh to improve the quality of flow and transport simulations around and along the well. The method was successfully implemented within dfnWorks 2.0 (Hyman et al. 2015) to incorporate a well in a 3D DFN and to track particles leaving an injection well and migrating to a pumping well. Verification of the method against MODFLOW/MODPATH found a perfect match in simulated hydraulic head and particle tracking. Using three examples, the study showed that the method ensured the connectivity between wells and fracture backbones, and honored the physical processes of flow and transport along and around wells in DFNs. Recommendations are given for estimating the values of the six introduced well parameters in a real‐world case study.

Pham, Hai (ORCID:0000000310164394)↗
DOE OSTI2022

Peridynamics and surrogate modeling of pressure-driven well stimulation

In this work we use the peridynamics theory of solid mechanics to simulate fracture in an annular rock domain subject to an in-situ stress and create surrogate models that predict the area of the resulting cracks. Peridynamics is a non-local formulation of continuum mechanics that naturally accommodates material discontinuities. Furthermore, unlike other fracture modeling techniques there is no need to provide information about the crack path. We utilize the peridynamics code Peridigm and take a two-stage approach to fracture modeling. First an implicit solve is performed to compute the in-situ stress state. We then execute an explicit solve where a pressure loading designed to emulate fluid-driven hydraulic fracture is applied at the borehole and transmitted to the pre-stressed rock. We present results from polynomial and single and multi-level Gaussian process surrogate models constructed from a sampling study of the peridynamics model. The surrogates predict crack area given a measure of the in-situ stress anisotropy and rise time and amplitude of the pressure loading. These surrogates take a minuscule fraction of peridynamics model's running time to evaluate and are a step towards enabling advanced optimization and uncertainty quantification workflows that require many model evaluations.

42 ENGINEERING↗
DOE OSTI2024

A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction

This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

fracture↗
DOE OSTI2025

Review of High‐Speed Digital Image Correlation: Advancements and Good Practices

This paper reviews the current state of the art in high‐speed (HS) and ultrahigh‐speed (UHS) digital image correlation (DIC) techniques, emphasizing their critical role in experimental research across various scientific domains. HS and UHS DIC have evolved significantly, driven by advancements in camera systems, image processing algorithms and experimental methodologies. These developments have opened new avenues for capturing and analysing dynamic events with unprecedented temporal and spatial detail, but not without introducing challenges such as optical distortions, motion blur and lighting issues that can affect measurement quality. This review advocates for standardized reporting practices in HS/UHS DIC methodologies to improve reproducibility and reliability across studies, drawing on guidelines from the International Digital Image Correlation Society (iDICs). Through a comprehensive analysis of over 150 articles, this review identifies key advancements in imaging technology and their application in six research domains: material characterization, test development, fracture mechanics, model validation, ballistic and explosive phenomena assessment and measurement uncertainties. Distinctions between two‐dimensional (2D) and stereo‐DIC applications are explored, offering insights into their practical implementation and trade‐offs. Good practices for HS/UHS DIC applications are proposed along with suggestions for future directions for this evolving field, highlighting the indispensable role of technological innovation in expanding the capabilities of optical metrology.

experimental mechanics↗
DOE OSTI2023

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase-field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase-field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy-based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase-field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase-field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. In conclusion, the comparison of numerical results for two-dimensional test cases with existing analytical solutions demonstrates that the proposed phase-field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three-dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Towards machine-learning a fully-coupled constitutive model for thermal-hydraulic fracture in geothermal systems: phase I (Final Report)

This project, entitled “Towards machine-learning a fully-coupled constitutive model for thermal-hydraulic fracture in geothermal systems: phase I,” addresses challenges in understanding and controlling subsurface fracture networks, which are crucial for applications like deep geothermal heat mining and deep-crustal minerals/metals/hydrogen extraction. The research focuses on advancing the understanding of coupled thermal-hydro-mechanical-chemical (THMC) processes in geologic materials, particularly under the high temperature and pressure conditions found in the deep crust. This seed grant focused specifically on thermal cracking and the development of new constitutive models. Significant progress was made in both experimental and theoretical domains. To study micro-scale fracture formation, the project demonstrated the ability to create thermal cracking under stress in granite samples using a Paterson Gas-medium Deformation Apparatus.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Phase instability-coupled fracture behavior in garnet LLZO solid electrolytes: a machine learning-enabled atomistic study

Fracture in the garnet-type solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) poses a critical threat to both the performance and safety of solid-state batteries. To unravel the coupled chemomechanical processes that govern fracture in LLZO under external loading, we carry out large-scale molecular-dynamics simulations with a validated machine-learning force field. Our results demonstrate that triaxial stresses at crack flanks trigger a localized cubic-to-tetragonal phase transformation, which is accompanied by Li-ion rearrangement. The emergent tetragonal domains feature lattice contraction normal to the fracture plane, imposing coherent misfit strains that provide an additional driving force for further crack propagation. Crucially, introducing Li deficiencies stabilizes the cubic phase, postponing the phase transition and thereby delaying fracture initiation. These findings highlight the role of intrinsic phase instability in dictating LLZO's fracture resistance and its critical connection to local Li concentration. This chemomechanical coupling points toward targeted strategies to enhance the mechanical robustness of garnet electrolytes, including tuning Li content, ensuring dopant homogeneity, and refining processing protocols.

Materials science↗
DOE OSTI2020

A numerical-homogenization based phase-field fracture modeling of linear elastic heterogeneous porous media

Most porous media, such as geomaterials and biomaterials are highly heterogeneous in nature, and they contain large variations of microscopic pore structures, such as pore sizes, pore distribution, and pore shapes. The oscillation of microscopic structures is a substantial challenge in theoretical characterization and is usually ignored in continuous modeling. However, mechanical behavior of porous media such as deformation and failure, are essentially impacted by the microscopic heterogeneity which needs to be considered in modeling a porous media. Here, this research proposes a numerical modeling framework with a capability to investigate the effect of microscopic heterogeneity on the macroscopic fracture behavior in porous media by using a numerical homogenization technique, combined with the phase-field fracture modeling method. This numerical modeling strategy computes a homogenized elasticity tensor based on microscopic heterogeneous pore structures heterogenous porous domain by solving boundary value problems at microscopic domain. The strain energy and subsequent propagation of macroscopic fractures will be updated using homogenized stiffness information. Using this numerical scheme, the microscopic pore structure’s impact on the fracture behavior through the homogenized elastic tensor will be taken into account. This multiscale technique is benchmarked against classical problems. Finally, the results highlight the importance of the underlying pore structure and reveal that both fracture strength and propagation path can be influenced by the microscopic heterogeneity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗