Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Damage evolution”

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 217 records · Page 12

Incorporation of Failure Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in various coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

Impact↗

Implementation of a Tabulated Failure Model Into a Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

polymer matrix composites↗

ZrN coating as diffusion barrier in U(Mo) dispersion fuel systems

The control of the interaction between the U(Mo) fuel phase and the Al matrix is one of the challenges of dispersion fuel plate development for research reactors. Given the specific properties of this interaction layer, larger amounts of it in the meat could lead to a reduction of the plate mechanical integrity and thermal conductivity, eventually leading to pillowing. The SELENIUM project showed that by depositing a ZrN coating on the surface of the U(Mo) fuel particles, the amount of formed U(Mo)-Al interaction layer is limited but still present. Microstructural analysis performed on the as fabricated coating and fresh fuel plates containing ZrN coated U(Mo) dispersed in an Al matrix, revealed that the coating gets damaged during plate production. The post irradiation examinations (PIE) of the ZrN coated U(Mo) fuel plates, from the SELENIUM and SEMPER FIDELIS experiments, show how the U(Mo)-Al interaction layer is formed -only at those locations where the coating is missing or damaged -and the evolution of coating microstructure during irradiation. Finally, as a remedy, to further reduce the amount of interaction layer formed, the use of an Al-Si matrix was proposed based on the higher affinity of Si for U compared to the affinity of Al for U. PIE of a fuel plate consisting of ZrN coated U(Mo) dispersed in an Al-Si matrix irradiated in the SEMPER FIDELIS experiment, clearly demonstrates the benefit of adding Si to the matrix.

36 MATERIALS SCIENCE↗

Pit growth kinetics in aluminum: effects of salt loading and relative humidity

Abstract The growth kinetics of localized corrosion, e.g. pits, in corrosive environments often controls the service life of metallic components. Yet, our understanding of these kinetics is largely based on coupon-level, e.g. mass-loss, studies which provide limited insights into the evolution of individual damage events. It is critical to relate observed cumulative loss trends, such as links between changing humidity and mass loss rates, to the growth kinetics of individual pits. Towards this goal, we leverage in-situ X-ray computed tomography to measure the growth rates of over sixty pits in aluminum in four different humid, chloride environments over ≈3 days of exposure. Pit growth rates and final volumes increased with increasing droplet volume, which was observed to increase with increasing humidity and salt loading. Two factors, droplet spreading and oxide jacking, dramatically increased pit growth rates and final volumes.

Noell, Philip J. (ORCID:0000000346998999)↗

Mitigating X-ray-induced damage for accurate analysis of battery cathodes using transmission X-ray microscopy

Here, we investigate the impact of X-ray-induced damage on the evolution of battery cathodes using synchrotron transmission X-ray microscopy. In an operando coin cell, exposure to the X-ray beam can induce bubble formation, potentially rendering cathode particles electrochemically inactive and displacing particles from their initial position. However, utilizing external pressure on a pouch cell can mitigate bubble formation and its associated effects.

25 ENERGY STORAGE↗

Effect of local structure and stoichiometry on the dynamic behavior of bi-metal interfaces

Molecular dynamics (MD) simulations are performed to investigate the failure behavior of Al/Al, Al/Fe, and Al/AlFe bi-layers under high strain-rate shock loading. Additional simulations are performed to invesigate solely the tensile response of these systems and to isolate the effects of loading history on damage and failure. Specifically, the compression stage of the shock is excluded from the tensile simulations. This study shows that local structure variation, including the introduction of serrations, plays different roles in controlling damage nucleation depending on the specific loading condition. Under shock loading, the stress for void nucleation and the resulting void distribution is insensitive to an interface structure, whereas under pure tensile loading, the opposite is true, even though the strain rate is comparable. The above difference can be explained based on the corresponding change in the total deformation prior to damage nucleation under shock loading, which is missing from the pure tensile loading. Therefore, it is concluded that whether local interface structure matters under high strain-rate loading depends on the specific loading history and the stress state: the local structure variation matters only when it alters the deformation behavior, and subsequently, damage nucleation and evolution.

36 MATERIALS SCIENCE↗

Impact of trapping on tritium self-sufficiency and tritium inventories in fusion power plant fuel cycles

The dynamic analysis of fusion power plant (FPP) fuel cycles highlights the challenge of achieving tritium self-sufficiency in future FPPs. While state-of-the-art fuel cycle models offer valuable insights into the necessary design parameters for attaining tritium self-sufficiency, none of these models currently consider the impact of tritium trapping within fuel cycle components. However, detailed analysis of individual components reveals that substantial amounts of tritium can be trapped within the first wall, divertors, and breeding blanket systems, suggesting that tritium trapping may significantly influence the FPP ability to achieve self-sufficiency. The compounded effects of additional tritium traps generated by irradiation effects and component replacements further exacerbate this challenge. The novelty of this work is the integration of an explicit, physics-based model for tritium trapping, evolution of damage-induced traps, and component replacements into a dynamic, system-level model of a fuel cycle. The results show an increase of a factor 10 3 – 10 4 of tritium inventory in the first wall and vacuum vessel of an ARC-class FPP when accounting for the aforementioned phenomena. This, coupled with the replacement of components subject to significant tritium trapping, slows down fuel cycle dynamics, resulting in an extended tritium doubling time (50% increase), higher start-up inventory (30% increase), and higher required tritium breeding ratio (2%–5%) compared to a scenario without tritium trapping.

fuel cycle↗

Analyzing the effect of pressure on the properties of point defects in $\gamma\mathrm{U–Mo}$ through atomistic simulations

Uranium–molybdenum (U–Mo) alloys in monolithic fuel foil are the primary candidate for the conversion of high-performance research reactors in the USA. Monolithic fuel is utilized in a plate-type design with a zirconium diffusion barrier and aluminum cladding. These fuel types are unique in that they contain no plenum for the release of fission gases, which, in conjunction with the aluminum cladding, can lead to large stress states within the fuel. The nature of how fundamental processes of radiation damage, including the evolution of point defects, under such stresses occur is unknown. In this work, we present molecular dynamics simulations of the formation energy of point defects under applied stress. We report this work will allow for the implementation of stress-dependent microstructural evolution models of nuclear fuels, including those for both fission gas bubble growth and creep, which are critical to ensure the stable and predictable behavior of research reactor fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Atomic Oxygen (ATOX) simulation of Teflon FEP and Kapton H surfaces using a high intensity, low energy, mass selected, ion beam facility

A high intensity (greater than 10(exp 15) ions/sq cm) low energy (down to 5 eV) mass selected ion beam (MSIB) facility was used to study the effects of ATOX on two polymers commonly used for space applications (Kapton H and Teflon FEP). The polymers were exposed to O(+) and Ne(+) fluences on 10(exp 15) - 10(exp 19) ions/sq cm, using 30eV ions. A variety of analytical methods were used to analyze the eroded surfaces including: (1) atomic force microscopy (AFM) for morphology measurements; (2) total mass loss measurements using a microbalance; (3) surface chemical composition using x-ray photoelectron spectroscopy (XPS), and (4) residual gas analysis (RGA) of the released gases during bombardment. The relative significance of the collisional and chemical degradation processes was evaluated by comparing the effects of Ne(+) and O(+) bombardment. For 30 eV ions it was found that the Kapton is eroded via chemical mechanisms while Teflon FEP is eroded via collisional mechanisms. AFM analysis was found very powerful in revealing the evolution of the damage from its initial atomic scale (roughness of approx. 1 nm) to its final microscopic scale (roughness greater than 1 micron). Both the surface morphology and the average roughness of the bombarded surfaces (averaged over 1 micron x 1 micron images by the system's computer) were determined for each sample. For 30 eV a non linear increase of the Kapton roughness with the O(+) fluence was discovered (a slow increase rate for fluences phi less than 5 x 10(exp 17) O(+)/sq cm, and a rapid increase rate for phi greater than 5 x 10(exp 17) O(+)/sq cm). Comparative studies on the same materials exposed to RF and DC oxygen plasmas indicate that the specific details of the erosion depend on the simulation facility emphasizing the advantages of the ion beam facility.

Vered, R.↗

Chapter 21. Atomic-Level Dislocation Dynamics in Irradiated Metals

Primary damage and microstructure evolution in structural nuclear materials operating under conditions of a high flux of energetic atomic particles and high temperature and stress lead to the formation of a high concentration, non-homogeneous distribution of defect clusters in the form of dislocation loops, voids, gas-filled bubbles and radiation-induced precipitates of nanometer scale. They cause changes in many material properties. Being obstacles to dislocation glide, they strongly affect the mechanical properties, with an increase in yield and flow stresses and a reduction in ductility. Atomic-scale computer simulations can provide details of how these effects are influenced by the obstacle structure, applied stress, strain rate and temperature. Processes such as obstacle cutting, transformation, absorption and drag are observed. Some recent results for body-centered and face-centered cubic metals are described in this review and, where appropriate, comparisons are drawn with predictions based on the elasticity theory of crystal defects. Perspectives on how to use this information at higher scales and in particular in mesoscale, dislocation dynamics simulations are also discussed.

Osetskiy, Yury N.↗

Incorporation of Damage and Failure into an Orthotropic Elasto-Plastic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in the composite impact models currently available in LS-DYNA(Registered Trademark) is under development. In particular, the material model, which is being implemented as MAT 213 into a tailored version of LS-DYNA being jointly developed by the FAA and NASA, incorporates both plasticity and damage within the material model, utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength), and is able to analyze the response of composites composed with a variety of fiber architectures. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. The capability to account for the rate and temperature dependent deformation response of composites has also been incorporated into the material model. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, in composites it has been found that loading in one direction can lead to damage in multiple coordinate directions. To account for this phenomena, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. The onset of material failure, and thus element deletion, is being developed to be a function of the stresses and plastic strains in the various coordinate directions. Systematic procedures are being developed to generate the required input parameters based on the results of experimental tests.

Finite Element Method↗

Origin and limit of the recovery of damaged information by time reversal

Recently it was found that scrambled information can be partially recovered by a time-reversed evolution, even after being damaged by an intruder. Here we reconsider the origin of the information recovery, and argue that the presence of classical chaos does not preclude it and only leads to a quantitative reduction of the recovery ratio. We also show how decoherence (i.e., entanglement with the intruder) limits the recovery, by proving an upper bound on the recovery ratio in terms of the entangling power of the intruder's action.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗

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↗

The constitutive representation of high-temperature creep damage

The elastic-viscoplastic constitutive equations of Bodner-Partom were applied to modeling creep damage in a high temperature Ni-alloy, B1900 + Hf. Both tertiary creep in bulk materials and creep crack growth in flawed materials were considered. In the latter case, the energy rate line integral was used for characterizing the crack driving force, and the rate of crack extension was computed using a local damage formulation that assumed fracture was controlled by cavitation occurring within the crack-tip process zone. The results of this investigation were used to assess the evolution equation for isotropic damage utilized in the Bodner-Partom constitutive equations.

Chan, K. S.↗

Risk Analysis of Radiological Release from Pu-238 Targets During Manual Handling

Pu-238 Isotope Production Targets are routinely installed in the Advanced Test Reactor (ATR) core, transferred, and stored in the spent fuel canal. These evolutions involve manual handling and manipulation of the targets underwater using long handled tools. The ATR Safety Analysis Report (SAR) postulates a design basis accident which results in damage from manual manipulation of targets, and radiological consequences must be determined for receptors inside the reactor facility, as well as public receptors. This presentation presents the analysis used to determine the radiological consequences due to potential target damage in the ATR canal. The analysis considered radionuclide release fractions, damage ratios for handling evolutions, and entrainment of radionuclides in the canal water.

Advanced Test Reactor↗

Primary Radiation Damage Formation in Solids

The processes that give rise to changes in the microstructure and the physical and mechanical properties of solid metals and ceramics exposed to energetic particles are initiated by essentially elastic collisions between atoms in what has been called an atomic displacement cascade. These processes are similar in metals and ceramics, with some differences related to the ionic and covalent bonding in ceramics compared to metallic bonding. The formation and evolution of primary radiation damage mechanisms are described to provide an overview of how stable defects are formed by displacement cascades, as well as the nature and morphology of the defects themselves. The impact of the primary variables cascade energy and irradiation temperature are discussed, along with a range of secondary factors that can influence damage formation.

Zarkadoula, Eva↗

Development and Characterization of a Rate-Dependent Three-Dimensional Macroscopic Plasticity Model Suitable for Use in Composite Impact Problems

Several key capabilities have been identified by the aerospace community as lacking in the material/models for composite materials currently available within commercial transient dynamic finite element codes such as LS-DYNA. Some of the specific desired features that have been identified include the incorporation of both plasticity and damage within the material model, the capability of using the material model to analyze the response of both three-dimensional solid elements and two dimensional shell elements, and the ability to simulate the response of composites composed with a variety of composite architectures, including laminates, weaves and braids. In addition, a need has been expressed to have a material model that utilizes tabulated experimentally based input to define the evolution of plasticity and damage as opposed to utilizing discrete input parameters (such as modulus and strength) and analytical functions based on curve fitting. To begin to address these needs, an orthotropic macroscopic plasticity based model suitable for implementation within LS-DYNA has been developed. Specifically, the Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The coefficients in the yield function are determined based on tabulated stress-strain curves in the various normal and shear directions, along with selected off-axis curves. Incorporating rate dependence into the yield function is achieved by using a series of tabluated input curves, each at a different constant strain rate. The non-associative flow-rule is used to compute the evolution of the effective plastic strain. Systematic procedures have been developed to determine the values of the various coefficients in the yield function and the flow rule based on the tabulated input data. An algorithm based on the radial return method has been developed to facilitate the numerical implementation of the material model. The presented paper will present in detail the development of the orthotropic plasticity model and the procedures used to obtain the required material parameters. Methods in which a combination of actual testing and selective numerical testing can be combined to yield the appropriate input data for the model will be described. A specific laminated polymer matrix composite will be examined to demonstrate the application of the model.

Impact↗