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
NASA NTRS1992

Tensile behavior of the L(1)2 compound Al67Ti25Cr8

Temperature-related variations in tensile yield strength and ductility were studied on cast, homogenized and isothermally forged Al67Ti25Cr8. Yield strength dropped discontinuously between 623 K and 773 K and then decreased gradually with increasing temperature. Below 623 K, fracture occurred prior to macroscopic yielding. Ductility decreased from 0.2 percent at 623 K to zero at 773 K, but increased again at higher temperatures. At 1073 K, an elongation of 19 percent was obtainable, and ultimate tensile strength and localized necking were observed. Fracture surfaces and deformed microstructures were examined. The 1073 K tensile specimen that exhibited 19 percent elongation showed grain boundary serrations and some evidence of recrystallization (likely dynamic) although fracture occurred predominantly via an intergranular mode.

Kumar, K. S.↗
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 OSTI2023

MEMPHIS (Mesoscle Multiphysics Phase Field Simulator v.0

SAND2024-01427O Mesoscale Multiphysics Phase Field Simulator (MEMPHIS) is used for solving interfacial problems and can be applied to many evolutionary problems. Some examples are solidification dynamics, viscous fingering, fracture mechanics, hydrogen embrittlement, vesicle dynamics, electro-mechanical degradation of batteries, aging of microstructures, and segregation. MEMPHIS is written in Fortran 90 in a modular fashion, where model and numerical solvers are decoupled from one another. 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.

Dingreville, Remi↗
DOE OSTI2023

Using Low-Field Nuclear Magnetic Resonance and X-Ray Computed Microtomography Imaging to Explore Potential of Microbially-Induced Calcium Carbonate Precipitation Treatment to Seal Shale Fractures

Microbially-induced calcium carbonate precipitation (MICP) is a biological process in which microbially produced urease enzymes convert urea and calcium into solid calcium carbonate (CaCO 3 ) deposits. Studies have shown that MICP can be used to seal fractures in shale, raising the possibility of applying this technology to restimulate fracking wells by plugging underperforming fractures. For this and other applications to become a reality, non-invasive tools are needed to determine how effectively MICP seals shale fractures under subsurface conditions. In this study, a 2.54 cm wide and 5.08 cm long Marcellus shale core with a single, ~1 mm wide fracture held open by sand "proppant" underwent MICP-treatment at 60°C until reaching three orders of magnitude permeability reduction. Low-field nuclear magnetic resonance (LF-NMR) and X-Ray computed microtomography (μ-CT) techniques were used to assess the extent of biomineralization within the fracture. These tools revealed that while CaCO 3 precipitation occurred throughout the fracture, there was preferential precipitation around proppant, and the core sealed at the effluent end before filling most of the fracture. Both tools were able to independently calculate of the amount of solid biomineral formed inside the fracture. Furthermore, this study found that the distribution of proppant within the shale fracture was an important parameter controlling the degree of biomineralization.

clastic rock↗
DOE OSTI2023

Comparative Analysis of Imaging and Measurements of Micrometer-Scale Fracture Aperture Fields Within a Heterogeneous Rock Using PET and X-ray CT

Knowledge of the spatial distribution of fracture apertures is essential for reliable characterization of flow and transport processes in fractured systems and for better understanding of physicochemical matrix–fracture interactions. Here, we propose and test two image-based methods, thereby extending the current experimental capabilities to characterize aperture size distribution in structurally heterogeneous geologic porous media noninvasively. The first approach utilizes an inversion method based on the dataset acquired from positron emission tomography (PET) and the second approach considers an extension of the classic missing attenuation technique that relies on clinical X-ray computed tomography (X-ray CT). Independent sets of imaging experiments are conducted on a fractured basalt core with heterogeneous matrix properties and aperture distributions to compare the two methodologies. A repeat of each experiment is conducted to verify the proposed workflows. The performance of these two imaging techniques is systematically evaluated through the analysis of signal-to-noise ratio, minimum fracture size detectability, and measurement errors. While both approaches provide a reliable estimation of fracture aperture distributions, PET yields a signal-to-noise ratio that is substantially higher than the corresponding X-ray CT measurements. Furthermore, uncertainties of the aperture values for PET are considerably lower ($\bar {\sigma}_{\text {d}} = 15\%$) compared to those obtained from X-ray CT ($\bar {\sigma}_{\text {d}} = 29\%$), allowing for the detection of minimum aperture sizes of 20 $\mu$m with 70% confidence level. Finally, these approaches provide key experimental tools for better understanding dynamic hydromechanical fracture properties in geologic systems.

58 GEOSCIENCES↗
DOE OSTI2021

DEM analysis of the effect of interface transition zone on dynamic splitting tensile behavior of high-strength concrete based on multi-phase model

Highlights: • Crushable aggregates with the real geometry shape are generated by using clump-cluster method. • Flat-Joint Model is adopted as the constitutive law. • The tensile strength is almost not affected by ITZ strength when its strength exceeds 0.8 times that of mortar. • A modified formula that can reflect the functional relationship between the ITZ strength and the TDIF is proposed. The tensile strength is one of the main parameters in the design of concrete structures and has an important influence on the cracking performance and durability of concrete materials. As the weak region in concrete, the strength of the interface transition zone (ITZ) directly affects the dynamic mechanical response of concrete. Meanwhile, considering that the tensile performance of concrete is far weaker than its compressive performance, the study of the strength of ITZ is of great significance to the dynamic tensile fracture characteristics and mechanical properties of concrete. In this paper, a method combining 3D scanning technology and Clump-Cluster was proposed to generate a real geometry shape model of crushable aggregate. On this basis, the 3D discrete element model (DEM) of high-strength concrete was constructed using the Flat-Joint Model (FJM). The effects of the ITZ strength on the failure mode, number of microcracks, tensile strength and strain rate effect of high-strength concrete under strain rates from 0.01/s to 100/s were studied by splitting tensile test. The simulation results indicate that the number of microcracks generated decreases with the increase of ITZ strength, and the proportion of ITZ cracks decreases obviously. When the ITZ strength exceeds the critical value of 0.8 times that of mortar, the tensile strength is almost not affected by ITZ strength. In addition, based on the numerical simulation results, a modified formula that can reflect the functional relationship between the ITZ strength and the dynamic increasing factor of the concrete tensile strength (TDIF) is proposed, and the rationality of the TDIF modified formula is verified.

36 MATERIALS SCIENCE↗
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 OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗
DOE OSTI2023

Peridynamic elastic waves in two-dimensional unbounded domains: Construction of nonlocal Dirichlet-type absorbing boundary conditions

The focus of this paper is on application of peridynamics (PD) to propagation of elastic waves in unbounded domains. We construct absorbing boundary conditions (ABCs) derived from a semi-analytical solution of the PD governing equation at the exterior region. This solution is made up of a finite series of plane waves, as fundamental solutions (modes), which satisfy the PD dispersion relations. The modes are adjusted to transmit the energy from the interior region (near field) to the exterior region (far field). The corresponding unknown coefficients of the series are found in terms of the displacement field at a layer of points adjacent to the absorbing boundary. This is accomplished through a collocation procedure at subregions (clouds) around each absorbing point. The proposed ABCs offer appealing advantages, which facilitate their application to PD. They are of Dirichlet-type, hence their implementation is relatively simple as no derivatives of the field variables are required. They are constructed in the time and space domains and thus application of Fourier and Laplace transforms, cumbersome for nonlocal models, is not required. At the discrete level, the modes satisfy the same numerical dispersion relations of the near field, which makes the far-field solution compatible with that of the near field. We scrutinize the performance of the proposed ABCs through several examples. So our investigation shows that the proposed ABCs perform stably in time with an appropriate level of accuracy even in problems characterized by highly-dispersive propagating waves, including crack propagation in semi-unbounded brittle solids.

42 ENGINEERING↗
DOE OSTI2022

Fragmentation analysis of a bar with the Lip-field approach

The Lip-field approach was introduced in Moës and Chevaugeon (2021) as a new way to regularize softening material models. It was tested in 1D quasistatic in Moës and Chevaugeon (2021) and 2D quasistatic in Chevaugeon and Moës (2021): this paper extends it to 1D dynamics, on the challenging problem of dynamic fragmentation. The Lip-field approach formulates the mechanical problem to be solved as an optimization problem, where the incremental potential to be minimized is the non-regularized one. Spurious localization is prevented by imposing a Lipschitz constraint on the damage field. Here, the displacement and damage field at each time step are obtained by a staggered algorithm, that is the displacement field is computed for a fixed damage field, then the damage field is computed for a fixed displacement field. Indeed, these two problems are convex, which is not the case of the global problem where the displacement and damage fields are sought at the same time. The incremental potential is obtained by equivalence with a cohesive zone model, which makes material parameters calibration simple. A non-regularized local damage equivalent to a cohesive zone model is also proposed. It is used as a reference for the Lip-field approach, without the need to implement displacement jumps. These approaches are applied to the brittle fragmentation of a 1D bar with randomly perturbed material properties to accelerate spatial convergence. Both explicit and implicit dynamic implementations are compared. Favorable comparison to several analytical, numerical and experimental references serves to validate the modeling approach.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Numerical modeling and experimental validation of low velocity impact of woven GFRP/CFRP composites

Low-velocity impact of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) composite laminates was studied experimentally and numerically. Hybrid laminates containing blocked layers of GFRP/CFRP/GFRP with all plies oriented at 0° were investigated. Relatively high impact energies were used to obtain full perforation of the laminate in a low-velocity impact setup. Numerical simulations were carried out using the in-house transient dynamics finite element code, Sierra/SM, developed at Sandia National Laboratories. A three-dimensional continuum damage model was used to describe the response of a woven composite ply. Two methods for handling delamination were considered and compared: (1) cohesive zone modeling and (2) continuum damage mechanics. The reduced model size achieved by omission of the cohesive zone elements produced acceptable results at reduced computational cost. Further, the comparison between different modeling techniques can be used to inform modeling decisions relevant to low velocity impact scenarios. The modeling was validated by comparing with the experimental results and showed good agreement in terms of predicted damage mechanisms and impactor velocity and force histories.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Advanced Simulation and Experiments of Strongly Coupled Geomechanics and Flow for Gas Hydrate Deposits: Validation and Field Application. Final Report

The objectives of the proposed research are (1) to investigate geomechanical responses induced by depressurization experimentally and numerically; (2) to enhance the current numerical simulation technology in order to simulate complex physically coupled processes by depressurization and (3) to perform in-depth numerical analyses of two selected potential production test sites: one based on the deposits observed at the Ulleung basin UBGH2-6 site; and the other based on well-characterized accumulations from the westend Prudhoe Bay. To these ends, the recipient will have the following specific objectives: 1). Information obtained from multi-scale experiments previously conducted at the recipient’s research partner (the Korean Institute of Geoscience and Mineral Resources (KIGAM)) that were designed to represent the most promising known Ulleung Basin gas hydrate deposit as drilled at site UBGH2-6 will be evaluated (Task 2). These findings will be further tested by new experimental studies at Lawrence Berkeley National Laboratory (LBNL) and Texas A&M (TAMU) (Task 3) that are designed capture complex coupled physical processes between flow and geomechanics, such as sand production, capillarity, and formation of secondary hydrates. The findings of Tasks 2 and 3 will be used to further improve numerical codes. 2) Develop (in Tasks 4 through 6) an advanced coupled geomechanics and non-isothermal flow simulator (T+MAM) to account for large deformation and strong capillarity. This new code will be validated using data from the literature, from previous work by the project team, and with the results of the proposed experimental studies. The developed simulator will be applied to both Ulleung Basin and Prudhoe Bay sites, effectively addressing complex geomechanical and petrophysical changes induced by depressurization (e.g., frost-heave, strong capillarity, cryo-suction, induced fracturing, and dynamic permeability).

02 PETROLEUM↗
NASA NTRS1981

Fundamental aspects in the quantitative ultrasonic determination of fracture toughness: General equations

The problem of establishing a theoretical groundwork for experimentally-found correlations between ultrasonic and fracture toughness factors in polycrystalline metals is discussed. It is noted that the link between these material properties and ultrasonic factors are the microstructural parameters that interact with stress wave propagation during deformation and fracture. The dynamic response of material inhomogeneities and the strains and displacements they undergo under incident stress waves are considered. Dynamic strains and displacements inside and outside scatterers are treated. The underlying approach, the formulation and governing equations for the eigenstrains, and the determination of the energy due to the presence of inhomogeneities are presented. The stress wave interaction problem is presented in terms of the dynamic eigenstrain concept.

Fu, L. S.↗
NASA NTRS1985

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 June 4 and 5, 1985, at the NASA Lewis Research Center. Presentations were made by industry, university, and government researchers organized into four sessions: aerothermodynamic loads; structural dynamics, fatigue, fracture, and constitutive modeling; and instrumentation. The principal objectives of the conference were to disseminate research results to date and future plans in each of the four areas. This publication contains the extended abstracts and the visual material presented during the conference.

Stanley J Marsik↗
NASA NTRS1990

Activities of the Structures Division, Lewis Research Center

The purpose of the NASA Lewis Research Center, Structures Division's 1990 Annual Report is to give a brief, but comprehensive, review of the technical accomplishments of the Division during the past calendar year. The report is organized topically to match the Center's Strategic Plan. Over the years, the Structures Division has developed the technology base necessary for improving the future of aeronautical and space propulsion systems. In the future, propulsion systems will need to be lighter, to operate at higher temperatures and to be more reliable in order to achieve higher performance. Achieving these goals is complex and challenging. Our approach has been to work cooperatively with both industry and universities to develop the technology necessary for state-of-the-art advancement in aeronautical and space propulsion systems. The Structures Division consists of four branches: Structural Mechanics, Fatigue and Fracture, Structural Dynamics, and Structural Integrity. This publication describes the work of the four branches by three topic areas of Research: (1) Basic Discipline; (2) Aeropropulsion; and (3) Space Propulsion. Each topic area is further divided into the following: (1) Materials; (2) Structural Mechanics; (3) Life Prediction; (4) Instruments, Controls, and Testing Techniques; and (5) Mechanisms. The publication covers 78 separate topics with a bibliography containing 159 citations. We hope you will find the publication interesting as well as useful.

Source record↗
DOE OSTI2022

Dynamic aspects of graphene deformation and fracture from approximate density functional theory

Graphene is one of the most intriguing two-dimensional carbon materials. Its mechanical strength and failure are key concerns for materials engineering and applications. Despite the success of fracture mechanics, the mechanism of how pristine materials fail remains an elusive problem. While many theoretical studies based on molecular dynamics using empirical forcefields have tried to address this question, atomic-scale mechanics are not clearly understood. Especially, a widely employed bond-breaking approach based on the critical bond length has not been rigorously tested. Here, utilizing molecular dynamics simulations with density functional based tight binding, we investigate how the failure of the pristine material initiates. The Wiberg bond order (W BO ) to estimate the change of chemical bonds shows a transition from sp 2 (W BO ~ 1.33) to sp 3 (WBO < 1.0) during the deformation. However, it reveals that a single threshold value for either the WBO or bond length is insufficient to decide failure of pristine material without free edges or defects. Instead, collective behaviors of the local atomic group govern the fracture initiation of pristine graphene. Overall, our study provides dynamic mechanical responses based on quantum mechanics, which have not been captured by empirical forcefields, opening opportunities to design properties by precisely coupling the mechanics and quantum mechanics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2006

Molecular-dynamics Simulation-based Cohesive Zone Representation of Intergranular Fracture Processes in Aluminum

A traction-displacement relationship that may be embedded into a cohesive zone model for microscale problems of intergranular fracture is extracted from atomistic molecular-dynamics simulations. A molecular-dynamics model for crack propagation under steady-state conditions is developed to analyze intergranular fracture along a flat 99 [1 1 0] symmetric tilt grain boundary in aluminum. Under hydrostatic tensile load, the simulation reveals asymmetric crack propagation in the two opposite directions along the grain boundary. In one direction, the crack propagates in a brittle manner by cleavage with very little or no dislocation emission, and in the other direction, the propagation is ductile through the mechanism of deformation twinning. This behavior is consistent with the Rice criterion for cleavage vs. dislocation blunting transition at the crack tip. The preference for twinning to dislocation slip is in agreement with the predictions of the Tadmor and Hai criterion. A comparison with finite element calculations shows that while the stress field around the brittle crack tip follows the expected elastic solution for the given boundary conditions of the model, the stress field around the twinning crack tip has a strong plastic contribution. Through the definition of a Cohesive-Zone-Volume-Element an atomistic analog to a continuum cohesive zone model element - the results from the molecular-dynamics simulation are recast to obtain an average continuum traction-displacement relationship to represent cohesive zone interaction along a characteristic length of the grain boundary interface for the cases of ductile and brittle decohesion. Keywords: Crack-tip plasticity; Cohesive zone model; Grain boundary decohesion; Intergranular fracture; Molecular-dynamics simulation

Yamakov, Vesselin I.↗