Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture modelling”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Scale-bridging with the extended/generalized finite element method for linear elastodynamics

This paper presents an extended/generalized finite element method for bridging scales in linear elastodynamics in the absence of scale separation. More precisely, the GFEMgl framework is expanded to enable the numerical solution of multiscale problems through the automated construction of specially-tailored shape functions, thereby enabling high-fidelity finite element modeling on simple, fixed finite element meshes. Furthermore, this introduces time-dependencies in the shape functions in that they are subject to continuous adaptation with time. The temporal aspects of the formulation are investigated by considering the Newmark-β time integration scheme, and the efficacy of mass lumping strategies is explored in an explicit time-stepping scheme. This method is demonstrated on representative wave propagation examples as well as a dynamic fracture problem to assess its accuracy and flexibility.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Increased Rates of Hip and Spine Fractures in Astronauts are Associated with Longer Spaceflight Durations

Given the accelerated declines in bone mineral density reported in US astronauts serving on the International Space Station, there is a possibility that a prolonged exposure to spaceflight could put astronauts at increased fracture risk during subsequent planetary explorations or after return to Earth. To date, there have been no fractures in astronauts during spaceflight. Recently-reported data from a flight study suggested that intake of the bisphosphonate alendronate combined with resistive exercise (on the Advanced Resistive Exercise Device – ARED) preserved more astronauts at their preflight skeletal status after a ~6-mos spaceflight than with ARED exercise alone. However, there is low acceptance of alendronate, or the recommended zoledronic acid, to suppress in-flight bone loss largely due to concern for unintended medication side effects and a perceived absence of astronaut fractures that can be directly attributable to spaceflight exposure. Hence, a recent survey of all astronaut fractures was analyzed to assess whether the rate of postflight fractures could be associated with extended spaceflight exposures. The fracture survey was of 262 astronauts across 8433.6 person-years (PY) and grouped according to previous exposure to a “Long Duration (LD)” spaceflight (3 mos or greater, but typically 6-mos). The number of fractures was modeled using count regression adjusting for i) the repeated measures within individuals across the groups, ii) the differing length of surveillance for each individual and iii) age at start of surveillance for each of the groups. Astronaut Group by Spaceflight Exposure #Individuals PY Covered No Exposure to Spaceflight 262 1882.7 No Exposure to LD spaceflight “Non-LD” 232 6121.3 Exposure to LD spaceflight 42 429.6 Analysis revealed an increased rate of hip and spine fractures for surveillance periods compared between LD Exposure and Non-LD Exposure: Hip – Relative Rate: 3.4 (95%CI: 1.0, 12.0); Spine– Relative Rate: 3.4 (95%CI: 1.0, 11.6). For all other skeletal sites, rates of fracture were similar for LD and Non-LD exposures. Data comparisons were conducted between groups of astronauts to preclude the interpretation of greater physical activity in astronauts as a contributing factor to fracture rates. This analysis of fracture data is the first report associating postflight hip and spine fractures in astronauts with exposure to extended duration spaceflights. Further analyses integrating co-morbidities for bone loss are planned.

skeleton↗
NASA NTRS2013

Micromechanics Fatigue Damage Analysis Modeling for Fabric Reinforced Ceramic Matrix Composites

A micromechanics analysis modeling method was developed to analyze the damage progression and fatigue failure of fabric reinforced composite structures, especially for the brittle ceramic matrix material composites. A repeating unit cell concept of fabric reinforced composites was used to represent the global composite structure. The thermal and mechanical properties of the repeating unit cell were considered as the same as those of the global composite structure. The three-phase micromechanics, the shear-lag, and the continuum fracture mechanics models were integrated with a statistical model in the repeating unit cell to predict the progressive damages and fatigue life of the composite structures. The global structure failure was defined as the loss of loading capability of the repeating unit cell, which depends on the stiffness reduction due to material slice failures and nonlinear material properties in the repeating unit cell. The present methodology is demonstrated with the analysis results evaluated through the experimental test performed with carbon fiber reinforced silicon carbide matrix plain weave composite specimens.

Min, J. B.↗
DOE OSTI2023

Phase Field Modeling of Corrosion for Design of Next-Generation Magnesium-Aluminum Vehicle Joints

In this study, WPI and its partners sought to understand galvanic corrosion in Mg-Al friction stir weld (FSW) joints. Unlike fusion welds, FSW is done in the solid state. This limits (but does not eliminate) the formation of brittle Al-Mg intermetallics, resulting in better joint strength and corrosion resistance. This study aimed to understand corrosion behavior of Mg-Al FSW joints by using a Cahn-Hilliard phase field electrochemistry model based on the work of Pongsaksawad et al. Its context challenge problem was an ultra-light door designed by Magna with 6061-Al sheet outer and ZEK100-Mg sheet inner panels, joined by an adhesive and FSW. Unlike most FSW joints, which weld through the softer material into the harder one, the Al sheet’s hem joint around Mg required the FSW tool to through the harder Al into the softer Mg, which was itself a challenge. Despite this challenge, PNNL succeeded in consistently making very strong FSW lap joints, by using a triflute tool and power control. Characterization showed mostly separate Al and Mg regions in the weld, usually with a hook protruding from Al into Mg, which likely contributed to high joint strength. Nanohardness mapping showed higher hardness in thin Mg swirls into the Al nugget, likely due to the intermetallics. The study used the industry standard SAE J2334 Cyclic Corrosion Test (CCT) to simulate corrosion conditions over the life of a vehicle, as well as linear polarization testing and ASTM G71 pitting corrosion testing. Corrosion reactions were: H 2 O + ½O 2 + 2e - → 2OH - at the Al cathode, and either Mg → Mg 2+ + 2e - or Mg + 2OH - → Mg(OH) 2 + 2e - at the Mg anode, with electron transfer through the joint. Mass loss increased with Al section length, supporting the hypothesis that galvanic corrosion in Al-Mg joints is limited by the cathodic reaction. Both pitting corrosion (first anode reaction) and hydroxide film formation (second reaction) were observed. The J2334 CCT test method was slightly modified, adding a chromium solution cleaning step after each week of corrosion testing, in order to remove Mg oxides and hydroxides and accurately measure mass loss over time. That said, FSW joints treated in this way did not did not exhibit significant reduction in lap shear strength of the joints vs. newly welded samples which had not undergone corrosion. Corrosion appeared to consist mainly of pitting in Mg in a way which did not directly affect the bond between the materials where the joints failed. A Cahn-Hilliard phase field model coupled with electrical potential described the electronically mediated galvanic corrosion electrochemical reactions described above in the four-component Al-Mg-H 2 O-(H 2 O+½O 2 ) system. This 2-D model successfully predicted Mg(OH)₂ formation at the anode in some circumstances, and Mg pitting corrosion in others, and showed the correct electric field directions in both cases. The model predicted Mg pit corrosion depth within a factor of two of measured pits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Localization from Hilbert space shattering: From theory to physical realizations

We show how a finite number of conservation laws can globally “shatter” Hilbert space into exponentially many dynamically disconnected subsectors, leading to an unexpected dynamics with features reminiscent of both many-body localization and quantum scars. A crisp example of this phenomenon is provided by a “fractonic” model of quantum dynamics constrained to conserve both charge and dipole moment. We show how the Hilbert space of the fractonic model dynamically fractures into disconnected emergent subsectors within a particular charge and dipole symmetry sector. This shattering can occur in arbitrary spatial dimensions. A large number of the emergent subsectors, exponentially many in system volume, have dimension one and exhibit strictly localized quantum dynamics—even in the absence of spatial disorder and in the presence of temporal noise. Other emergent subsectors display nontrivial dynamics and may be constructed by embedding finite-sized nontrivial blocks into the localized subspace. While “fractonic” models provide a particularly clean realization, the shattering phenomenon is more general, as we discuss. We also discuss how the key phenomena may be readily observed in near term ultracold atom experiments. In experimental realizations, the conservation laws are approximate rather than exact, so the localization only survives up to a prethermal timescale that we estimate. We comment on the implications of these results for recent predictions of Bloch/Stark many-body localization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2013

A Comparison of Materials Issues for Cermet and Graphite-Based NTP Fuels

This paper compares material issues for cermet and graphite fuel elements. In particular, two issues in NTP fuel element performance are considered here: ductile to brittle transition in relation to crack propagation, and orificing individual coolant channels in fuel elements. Their relevance to fuel element performance is supported by considering material properties, experimental data, and results from multidisciplinary fluid/thermal/structural simulations. Ductile to brittle transition results in a fuel element region prone to brittle fracture under stress, while outside this region, stresses lead to deformation and resilience under stress. Poor coolant distribution between fuel element channels can increase stresses in certain channels. NERVA fuel element experimental results are consistent with this interpretation. An understanding of these mechanisms will help interpret fuel element testing results.

Nuclear Propulsion↗
NASA NTRS2007

Risk Assessment of Bone Fracture During Space Exploration Missions to the Moon and Mars

The possibility of a traumatic bone fracture in space is a concern due to the observed decrease in astronaut bone mineral density (BMD) during spaceflight and because of the physical demands of the mission. The Bone Fracture Risk Module (BFxRM) was developed to quantify the probability of fracture at the femoral neck and lumbar spine during space exploration missions. The BFxRM is scenario-based, providing predictions for specific activities or events during a particular space mission. The key elements of the BFxRM are the mission parameters, the biomechanical loading models, the bone loss and fracture models and the incidence rate of the activity or event. Uncertainties in the model parameters arise due to variations within the population and unknowns associated with the effects of the space environment. Consequently, parameter distributions were used in Monte Carlo simulations to obtain an estimate of fracture probability under real mission scenarios. The model predicts an increase in the probability of fracture as the mission length increases and fracture is more likely in the higher gravitational field of Mars than on the moon. The resulting probability predictions and sensitivity analyses of the BFxRM can be used as an engineering tool for mission operation and resource planning in order to mitigate the risk of bone fracture in space.

Lewandowski, Beth E.↗
NASA NTRS2008

Risk Assessment of Bone Fracture During Space Exploration Missions to the Moon and Mars

The possibility of a traumatic bone fracture in space is a concern due to the observed decrease in astronaut bone mineral density (BMD) during spaceflight and because of the physical demands of the mission. The Bone Fracture Risk Module (BFxRM) was developed to quantify the probability of fracture at the femoral neck and lumbar spine during space exploration missions. The BFxRM is scenario-based, providing predictions for specific activities or events during a particular space mission. The key elements of the BFxRM are the mission parameters, the biomechanical loading models, the bone loss and fracture models and the incidence rate of the activity or event. Uncertainties in the model parameters arise due to variations within the population and unknowns associated with the effects of the space environment. Consequently, parameter distributions were used in Monte Carlo simulations to obtain an estimate of fracture probability under real mission scenarios. The model predicts an increase in the probability of fracture as the mission length increases and fracture is more likely in the higher gravitational field of Mars than on the moon. The resulting probability predictions and sensitivity analyses of the BFxRM can be used as an engineering tool for mission operation and resource planning in order to mitigate the risk of bone fracture in space.

Lewandowski, Beth E.↗
DOE OSTI2023

Identifying crack tip position and stress intensity factors from displacement data

Fracture prognosis and characterization efforts require knowledge of crack tip position and the Stress Intensity Factors (SIFs) acting in the vicinity of the crack. Here, in this work, we present an efficient numerical approach to infer both of these characteristics under a consistent theoretical framework from noisy, unstructured displacement data. The novel approach utilizes the separability of the asymptotic linear elastic fracture mechanics fields to expedite the search for crack tip position and is particularly useful for noisy displacement data. The manuscript begins with an assessment of the importance of accurately locating crack tip position when quantifying the SIFs from displacement data. Next, the proposed separability approach for quickly inferring crack tip position is introduced. Comparing to the widely used displacement correlation approach, the performance of the separability approach is assessed. Cases involving both noisy data and systematic deviation from the asymptotic linear elastic fracture mechanics model are considered, e.g. inelastic material behavior and finite geometries. An open source python implementation of the proposed approach is available for use by those doing field and laboratory work involving digital image correlation and simulations, e.g. finite element, discrete element, molecular dynamics and peridynamics, where the crack tip position is not explicitly defined.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Elastic-Plastic J-Integral Solutions or Surface Cracks in Tension Using an Interpolation Methodology

No closed form solutions exist for the elastic-plastic J-integral for surface cracks due to the nonlinear, three-dimensional nature of the problem. Traditionally, each surface crack must be analyzed with a unique and time-consuming nonlinear finite element analysis. To overcome this shortcoming, the authors have developed and analyzed an array of 600 3D nonlinear finite element models for surface cracks in flat plates under tension loading. The solution space covers a wide range of crack shapes and depths (shape: 0.2 less than or equal to a/c less than or equal to 1, depth: 0.2 less than or equal to a/B less than or equal to 0.8) and material flow properties (elastic modulus-to-yield ratio: 100 less than or equal to E/ys less than or equal to 1,000, and hardening: 3 less than or equal to n less than or equal to 20). The authors have developed a methodology for interpolating between the goemetric and material property variables that allows the user to reliably evaluate the full elastic-plastic J-integral and force versus crack mouth opening displacement solution; thus, a solution can be obtained very rapidly by users without elastic-plastic fracture mechanics modeling experience. Complete solutions for the 600 models and 25 additional benchmark models are provided in tabular format.

Allen, P. A.↗
NASA NTRS2013

Elastic-Plastic J-Integral Solutions or Surface Cracks in Tension Using an Interpolation Methodology. Appendix C -- Finite Element Models Solution Database File, Appendix D -- Benchmark Finite Element Models Solution Database File

No closed form solutions exist for the elastic-plastic J-integral for surface cracks due to the nonlinear, three-dimensional nature of the problem. Traditionally, each surface crack must be analyzed with a unique and time-consuming nonlinear finite element analysis. To overcome this shortcoming, the authors have developed and analyzed an array of 600 3D nonlinear finite element models for surface cracks in flat plates under tension loading. The solution space covers a wide range of crack shapes and depths (shape: 0.2 less than or equal to a/c less than or equal to 1, depth: 0.2 less than or equal to a/B less than or equal to 0.8) and material flow properties (elastic modulus-to-yield ratio: 100 less than or equal to E/ys less than or equal to 1,000, and hardening: 3 less than or equal to n less than or equal to 20). The authors have developed a methodology for interpolating between the goemetric and material property variables that allows the user to reliably evaluate the full elastic-plastic J-integral and force versus crack mouth opening displacement solution; thus, a solution can be obtained very rapidly by users without elastic-plastic fracture mechanics modeling experience. Complete solutions for the 600 models and 25 additional benchmark models are provided in tabular format.

Allen, Phillip A.↗
DOE OSTI2020

Stress effects on flow and transport in three-dimensional fracture networks

We investigate the effects of various external stress regimes on fracture apertures, fluid flow, and solute transport in three-dimensional fracture networks. We use well-established geomechanics equations coupled with discrete fracture network modeling to characterize changes in primary flow paths within a complex network as a function of stress magnitude and orientation. These changes manifest in the alterations of the fluid flow field and are measured in terms of Eulerian and Lagrangian flow observables including solute transport, which is a key problem in many hydrologic applications. Changes in primary flow paths affect the solute transport in the network by promoting anomalously early arrival or long tailing behavior. However, early time arrival is not ubiquitous in anisotropically stressed networks and in most cases there is a delayed arrival of solute, which is attributed to (i) the presence of low-velocity zones normal to the flow direction, (ii) the angle between flow direction and major compressive principal stress directions, and (iii) changes in the primary flow paths (i.e., increases in tortuosity and active network structure). Overall, flows become more channelized in anisotropically stressed fracture networks than in unstressed and isotropically stressed networks.

58 GEOSCIENCES↗
NASA NTRS1987

A nondestructive technique for predicting the strength remaining in filament wound composites subjected to low-level impact

A model for predicting the fracture strength of homogeneous materials is proposed. Impacted FWC samples were evaluated using ultrasonic testing and an X-ray dye penetration method. The ability of the model to measure fracture strength was also examined. The relation between attenuation and velocity measurements is studied. It is observed that the X-ray method is not useful for predicting fracture strength because the dye could not penetrate the matrix. It is noted that fracture strength predictions derived from the fracture mechanical model and the ultrasonic measurements correlate well with actual measured fracture strengths.

Madaras, E. I.↗
NASA NTRS1991

On representation of fracture profiles by fractional integrals of a Wiener process

An enhancement of the crack diffusion model of brittle fracture is proposed. In the model, certain probabilities are evaluated as averages over the set of all virtual crack paths, and the latter are modeled as realization of a Wiener process. However, commonly observed fracture profiles are smoother. In fact, one can show that their fractal dimension is smaller than the value 3/2 for a Wiener process. Accordingly, it is proposed to represent fracture profiles by graphs of fractional integrals of a Wiener process. This allows one to employ additional fractographic information obtained from fracture profiles, e.g. their fractal dimension, in a predictive formalism, evaluating probabilities of failure.

Kunin, B.↗
DOE OSTI2025

Phase-field modeling of orientation-dependent crack growth in ductile single crystals with anisotropic elasticity

Crack growth in ductile single crystals (DuSCs) is orientation dependent due to the anisotropies of crystal plasticity and elastic tensor. This study develops a phase-field model incorporating both crystal plasticity and crack growth and proposes a general method to decompose the elastic energy into compressive and tensile parts to prevent crack growth under compression in the phase-field description. The phase-field model, in combination with three Euler angles, is employed to simulate orientation-dependent crack growth in DuSCs. The contributions from crystal plasticity and anisotropic elasticity are compared, and the former is found to dominate in the anisotropy of crack growth in copper single crystals. Furthermore, the simulation results demonstrate that crystal orientation strongly affects the heterogeneous distribution of plastic strain and the interaction between plastic strain and crack growth. High-throughput phase-field simulations are performed with exhaustive crystal orientations, and the results are explained based on the anisotropy of the Taylor factor.

Computational Solid Mechanics↗
DOE OSTI2023

Understanding Subsurface Fracture Evolution Dynamics Using Time‐Lapse Full Waveform Inversion of Continuous Active‐Source Seismic Monitoring Data

Abstract Predicting the behavior, geometry, and flow properties of subsurface fractures remains a challenging problem. Seismic models that can characterize fractures usually suffer from low spatiotemporal resolution. Here, we develop a correlative double‐difference time‐lapse full waveform inversion of continuous active source seismic monitoring data for determining high‐spatiotemporal‐resolution time‐lapse V p models of in‐situ fracture evolution at a shallow contamination site in Wyoming, USA. Assisted by rock physics modeling, we find that (a) rapidly increasing pore pressure initializes and grows the fracture, increasing the porosity slightly (from ∼13.7% to ∼14.6%) in the tight clay formation, thus decreasing V p (∼50 m/s); (b) the fluid injection continues decreasing V p , likely through the introduction of gas bubbles in the injectate; and (c) final V p reductions reach over ∼150 m/s due to a posited ∼4.5% gas saturation. Our results demonstrate that high‐resolution V p changes are indicative of mechanical and fluid changes within the fracture zone during hydrofracturing.

58 GEOSCIENCES↗
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↗