Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inelastic deformations”

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 55 records · Page 3

A Multiphysics Multiscale Simulation Platform for Damage, Environmental Degradation, and Life Prediction of CMCs in Extreme Environments

This project successfully developed a multiphysics, multiscale computational framework to enhance the design and development of CMCs, with a focus on modeling highly nonlinear, time-dependent damage mechanisms and material degradation under extreme conditions, such as those experienced in turbine service environments. The project made significant advances in improving our understanding of progressive damage, oxidative degradation, and time-dependent inelastic deformation in CMCs, with particular attention to the role of uncertainties in predictions. Key outcomes include the integration of advanced material characterization, uncertainty quantification, and multiphysics constitutive models to predict the behavior of CMCs over their service life. A novel multiscale methodology was employed, which integrated microscale constituent behaviors with structural-scale responses, enabling the manufacturing defects in the microstructure that are prone to damage nucleation. Through the development of DL algorithms, the project advanced the prediction of damage initiation and crack propagation, taking into account the defect morphology and statistical variations across multiple scales. The framework was rigorously validated using thermomechanical experiments, which tested CMCs under various mechanical loadings at elevated temperatures, further enhancing the model's predictive capability. Overall, the research outcomes have provided a more accurate, reliable method for predicting CMC component life, significantly advancing material design, and improving component reliability in extreme environments. This work has strong implications for the optimization of turbine components and other high-performance applications where CMCs are used.

03 NATURAL GAS↗

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗

High-burnup Experiments in Reactivity Initiated Accidents (HERA)

HERA is a joint experimental program (JEEP) operating within the Nuclear Energy Agency’s (NEA’s) framework for irradiation experiments (FIDES). HERA is dedicated to the understanding of light water reactor (LWR) fuel performance at high burnup under reactivity-initiated accidents (RIA). During RIAs the fuel cladding can be breached by one of two mechanisms; pellet-cladding mechanical interaction or high temperature balloon and rupture. Pellet-cladding mechanical interaction (PCMI) failure occurs in the early phase of an RIA transient where the fuel pellet expands rapidly into the cladding prior to any meaningful heat transfer from the pellet to the cladding. This results in high stresses in the cladding resulting in crack propagation through the base Zircaloy metal and axially along the cladding tube. Failure by PCMI is a brittle failure mode and occurs at low cladding strains, approximately 2% hoop strain. If the cladding survives the PCMI phase of the transient, it may still fail due to ballooning and rupture induced by extended exposure to high temperature. Thermal transport from the fuel pin to the coolant is significantly reduced if a transition to film boiling on the cladding surface occurs. This result in both a rapid increase in fuel temperature, which drives fission gas release, and a spike in cladding temperature, which results in loss of mechanical strength. If the internal pressure of the fuel rod is greater than the system pressure the cladding can inelastically deform and eventually rupture. Failure by swelling and rupture occurs at larger strains, generally greater than 5% hoop strain. While PCMI failures are generally seen as more limiting, the dominance of one failure mode over the other can be affected by both the cladding material conditions, and the transient evolution.

(HERA)↗

High-burnup Experiments in Reactivity Initiated Accidents (HERA)

HERA is a joint experimental program (JEEP) operating within the Nuclear Energy Agency’s (NEA’s) framework for irradiation experiments (FIDES). HERA is dedicated to the understanding of light water reactor (LWR) fuel performance at high burnup under reactivity-initiated accidents (RIA). During RIAs the fuel cladding can be breached by one of two mechanisms; pellet-cladding mechanical interaction or high temperature balloon and rupture. Pellet-cladding mechanical interaction (PCMI) failure occurs in the early phase of an RIA transient where the fuel pellet expands rapidly into the cladding prior to any meaningful heat transfer from the pellet to the cladding. This results in high stresses in the cladding resulting in crack propagation through the base Zircaloy metal and axially along the cladding tube. Failure by PCMI is a brittle failure mode and occurs at low cladding strains, approximately 2% hoop strain. If the cladding survives the PCMI phase of the transient, it may still fail due to ballooning and rupture induced by extended exposure to high temperature. Thermal transport from the fuel pin to the coolant is significantly reduced if a transition to film boiling on the cladding surface occurs. This result in both a rapid increase in fuel temperature, which drives fission gas release, and a spike in cladding temperature, which results in loss of mechanical strength. If the internal pressure of the fuel rod is greater than the system pressure the cladding can inelastically deform and eventually rupture. Failure by swelling and rupture occurs at larger strains, generally greater than 5% hoop strain. While PCMI failures are generally seen as more limiting, the dominance of one failure mode over the other can be affected by both the cladding material conditions, and the transient evolution.

(HERA)↗

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets—Part I: Mesoscale simulations

Recently a mesoscale model was developed that includes thermomechanically consistent constitutive equations with dissipation due to porous compaction, inelastic distortional deformation rate, and slipping on persistent dry joint surfaces. In this Part I of a two-part paper, we describe enhancements of the mesoscale model which include: normal compliance with hysteresis, non-persistent joints with finite area, joint sets with angular deviation, and the weakening effect of fluid pressure in the joints. Examples study the influences of these effects on average strength. These computationally expensive simulations are used to inspire functional forms in a continuum model, described in Part II, which is needed for large scale simulations of real applications.

42 ENGINEERING↗

Review: Inelastic Constitutive Modeling: Polycrystalline Materials

This article provides a literature review that details the development of inelastic constitutive modeling as it relates to polycrystalline materials. This review distinguishes between inelastic constitutive models that account for nonlinear behavior at the microstructural level, time-independent classic plasticity models, and time-dependent unified models. Particular emphasis is placed on understanding the underlying theoretical framework for unified viscoplasticity models where creep and classical plasticity behavior are considered the result of applied boundary conditions instead of separable rates representing distinct physical mechanisms. This article establishes a clear understanding of the advantages of the unified approach to improve material modeling. This review also discusses recent topics in constitutive modeling that offer new techniques that bridge the gap between the microstructure and the continuum.

36 MATERIALS SCIENCE↗

Phase-Field Model Incorporating Large Inelastic Strain with Application to the Oxidation of High-Temperature Coating Systems

Structural metallic and ceramic materials at high temperatures often suffer from creep, fatigue and oxidation, which typically involves local large inelastic strain. The large strain is usually approached by Lagrangian-description-based finite deformation theory. Here we develop an efficient, Eulerian-description-based, Incremental Realization of Inelastic Deformation (IRID) algorithm to deal with large inelastic strain with small elastic strain. We apply this algorithm to simulate high temperature oxidation with severe volume expansion using a phase-field approach. The oxidation and deformation of bond coat can lead to spallation and failure of the thermal barrier coatings (TBC) and environmental barrier coatings (EBC), which is a critical concern in design of these coatings. We simulate the growth of thermally grown oxide (TGO) in these coating systems, with phase transformations and microstructure evolution incorporated. The simulation results are compared with relevant experimental reports in the literature.

Cheng, Tianle↗

Lifetimes of low-lying levels in 158 Gd

The low-lying structure of the well-deformed nucleus 158 Gd has been revisited to elucidate the nature of the low-lying states in 158 Gd. Earlier (p, t) studies identified numerous 0 + states below 4.3 MeV, prompting questions about whether these states correspond to collective vibrations or shape coexistence. New and previously reported (n, n' γ) measurements are combined, including γ-γ coincidences, excitation functions, and angular distributions, to extract lifetimes and transition probabilities for 44 excited states up to 2.7 MeV, including 32 previously unmeasured levels. Our results confirm or revise γ-ray placements and provide detailed transition strengths, revealing both weakly collective and strongly enhanced B(E2) and B(E1) transition probabilities. In particular, a tentative 0 + state at 2437.8 keV exhibits a strong interband B(E2) transition, which may be a candidate for a possible two-phonon (ββ) excitation. Systematic comparisons with neighboring Gd isotopes, Hartree–Fock–Bogoliubov, and interacting-boson model predictions suggest that the first excited 0 + state in 158 Gd is predicted to be a β-vibration, although it is weakly collective.

158Gd↗

Observing the effects of numbers of valence nucleons on 0$^{+}_{𝑔⁡𝑠}$ → 2$^{+}_{1}$ transitions in deformed nuclei by comparing proton and neutron transition matrix elements

We examined the ratios of neutron and proton transition matrix elements, 𝑀 𝑛 /𝑀 𝑝 , for the 0$^{+}_{𝑔⁡𝑠}$ → 2$^{+}_{1}$ transitions in 48 even-even stable nuclei with 𝑁 > 20 for which electromagnetic matrix elements were compiled by Pritychenko et al. and for which high-quality inelastic proton-scattering data were available. Several deformed rare-earth nuclei have (𝑀 𝑛 /𝑀 𝑝 )/(𝑁/𝑍) values significantly below 1.0, which is not consistent with a simple liquid-drop picture. However, this phenomenon can be explained using a schematic picture in which 𝑀 𝑝 reaches a maximum at proton midshell (𝑍 = 66) and 𝑀𝑛 reaches its maximum at neutron midshell (𝑁 = 104). Several midmass vibrational nuclei have 𝑀 𝑛 /𝑀 𝑝 values significantly below 𝑁/𝑍, which is not consistent with the expectation that 𝑀 𝑛 /𝑀 𝑝 = 𝑁/𝑍 in such nuclei. As a result, a shell-model investigation of these observations might yield insights about this behavior.

Collective levels↗

Athermal glass work at the nanoscale: Engineered electron-beam-induced viscoplasticity for mechanical shaping of brittle amorphous silica

Amorphous silica deforms viscoplastically at elevated temperatures, which is common for brittle glasses. The key mechanism of viscoplastic deformation involves interatomic bond switching, which is thermally activated. Here, for this study, we precisely control the mechanical shaping of brittle amorphous silica at the nanoscale via engineered electron–matter interactions without heating. We observe a ductile plastic deformation of amorphous silica under a focused scanning electron beam with low acceleration voltages (few to tens of kilovolts) during in-situ compression studies, with unique dependence on the acceleration voltage and beam current. By simulating the electron–matter interaction, we show that the deformation of amorphous silica depends strongly on the volume where inelastic scattering occurs. The electron–matter interaction via e-beam irradiation alters the Si–O interatomic bonds, enabling the high-temperature deformation behavior of amorphous silica to occur athermally. Finally, by systematically controlling the electron–matter interaction volume, it is possible to mechanically shape the brittle amorphous silica on a small scale at room temperature to a level comparable to glass shaping at high temperatures. The findings can be extended to develop new fabrication processes for nano- and microscale brittle glasses.

36 MATERIALS SCIENCE↗

Simultaneous inference of the compressibility and inelastic response of tantalum under extreme loading

We study the deformation of tantalum under extreme loading conditions. Experimental velocity data are drawn from both ramp loading experiments on Sandia’s Z-machine and gas gun compression experiments. The drive conditions enable the study of materials under pressures greater than 100 GPa. We provide a detailed forward model of the experiments including a model of the magnetic drive for the Z-machine. Utilizing these experiments, we simultaneously infer several different types of physically motivated parameters describing equation of state, plasticity, and anelasticity via the computational device of Bayesian model calibration. Characteristics of the resulting calculated posterior distributions illustrate relationships among the parameters of interest via the degree of cross correlation. The calibrated velocity traces display good agreement with the experiments up to experimental uncertainty as well as improvement over previous calibrations. Examining the Z-shots and gun-shots together and separately reveals a trade-off between accuracy and transferability across different experimental conditions. Implications for model calibration, limitations from model form, and suggestions for improvements are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Near-Threshold Dipole Strength in 10 Be with Isoscalar Character

Isoscalar dipole transitions are a distinctive fingerprint of cluster structures. A 1 − resonance at 7.27(10) MeV, located just below the 𝛼-emission threshold, has been observed in the deuteron inelastic scattering reactions off 10 Be. The deformation lengths of the excited states in 10 Be below 9 MeV have been inferred from the differential cross sections using coupled channel calculations. This observed 1 − resonance has isoscalar characteristics and exhausts approximately 5%–15% of the isoscalar dipole energy-weighted sum rule, providing evidence for pronounced 𝛼 cluster structure in 10 Be. The Gamow coupled channel approach supports this interpretation and suggests the near-threshold effect might be playing an important role in this excitation energy domain. Here, the 𝛼 +𝛼 +𝑛 +𝑛 four-body calculation reproduces the observed enhanced dipole strength, implying that the four-body cluster structure is essential to describe the 1 − states in 10 Be.

Cluster models↗

Sample problems for Section III, Division 5 design by inelastic analysis of Grade 91 components

This report works through two sample design analyses of representative Grade 91 components using the ASME Boiler & Pressure Vessel Code Section III, Division 5, Appendix HBB-T rules for the strain and deformation and creep-fatigue design criteria using a design by inelastic analysis approach. The component geometries and loadings were selected to span a wide variety of potential applications. The sample components are an endplate geometry, which is fairly representative of typical pressure vessel construction, and a core block, representing a critical component in a heat pipe microreactor, which has a complicated geometry and complex thermal stress history. These sample problems serve several functions. First, they are complete worked examples of how to apply the current design by inelastic analysis rules. Second, they serve as test cases to evaluate the current rules and provide feedback on how they could be optimized to deal with complicated component geometries or to be better automated in modern finite element analysis software. Finally, the analysis here uses an inelastic constitutive model for Grade 91 developed by Argonne National Laboratory and proposed for inclusion in the ASME Boiler & Pressure Vessel Code Section III, Division 5 rules as a reference constitutive model. Using this constitutive model in the sample problem analyses tests the model in the context of the ASME Code rules. The final conclusions of this report are that the material model and the current Code rules are satisfactory, but that the Code rules could be optimized to simplify automation and reduce over conservatism. Future work should address the specific optimization topics discussed in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental study of compaction localization in carbonate rock and constitutive modeling of mechanical anisotropy

Abstract Sedimentary rocks are inherently anisotropic and prone to strain localization. While the influence of rock anisotropy on the brittle/dilative regime has been studied extensively, its influence on the ductile/compactive regime is much less explored. This paper discusses the anisotropic behavior of a high‐porosity carbonate rock from central Europe (the Maastricht Tuffeau). A set of triaxial tests with concurrent x‐ray tomography has been performed at different confining pressures. The anisotropic characteristics of this rock have been investigated by testing samples cored at different inclinations of the bedding, thus revealing non‐negligible effects of the coring direction on yielding and compaction behavior. Specifically, samples cored perpendicular to bedding display higher strength and longer stages of post‐yielding deformation before manifesting re‐hardening. Despite such alterations of the inelastic response, Digital Image Correlation has revealed that the strain localization mode is independent of the coring direction, thus being primarily affected by the confinement level. To capture the observed interaction between material anisotropy and compaction behavior at the continuum‐scale, an elastoplastic constitutive law has been proposed. For this purpose, a set of tensorial bases has been introduced to replicate how the oriented rock fabric modulates the yielding and plastic flow characteristics of the material. The analyses show that the impact of the coring direction on yield function and plastic flow rule is fundamentally different, thus requiring the use of distinct projection strategies (a strategy here defined heterotopic mapping). The performance of the model, studied through parametric analyses and by calibrating the experimental results, illustrates the improved capability of the proposed constitutive approach when applied to strongly anisotropic porous rocks.

Shahin, Ghassan↗

A Physics-Based Data-Driven Approach for Modeling of Environmental Degradation in Elastomers

Abstract Elastomers are now commonly used in a number of industries, including aerospace, structure, transportation, shipbuilding, and automotive, due to their excellent workability, formability, and flexibility. During their activity, elastomers are subjected to harsh environmental conditions, which decreases their resilience. False predictions made early in their lives can have major financial and environmental implications. Elastomers’ performance and properties, such as strength, durability, and density, are influenced by chemical changes in these materials, known as degradation, which occurs over time. This process can alter the morphology of a polymer matrix as well as cause chain scission and cross-linking, resulting in different behaviors than that of the unaged material. To demonstrate the effect of thermaloxidative aging on the mechanical behavior of elastomers, several experimental and theoretical models have been proposed. In view of the large volume of experimental data available on micro-structural evolution in the course of aging, we propose a physics-based data-driven approach to overcome the shortcomings of both phenomenological and micro-mechanical models. This work presents a novel thermodynamically consistent, multiagent machine-learned model for predicting the constitutive behavior of cross-linked elastomers during environmental aging, such as thermo-oxidative and hydrolytic aging for various states of deformation. Single mechanism degradation changes the polymer matrix over time where it is causing chain scission, reduction of cross-links, and morphology change. To capture the idealized Mullins effect and permanent set due to the effect of single aging mechanisms on nonlinear mechanical responses of elastomers, we propose a data-driven model for simulating inelastic elements in a polymer matrix. By using a sequential order reduction, we were able to reduce the 3D stress-strain tensor mapping problem to a small number of super-constrained 1D mapping problems. To systematically classify such mapping problems into a few categories, an assembly of multiple replicated conditional neural network learning agents (L-agents) is used based on our recent work. Each category is represented by a different type of agent. The effect of deformation history, aging time, and aging temperature is captured by this model. The model is validated using a broad collection of data, ranging from our experimental results to data from the literature. In addition, thermodynamic consistency and frame independence are investigated. The most significant achievements of this model are its precision, simplicity, and prediction of inelasticity under various states of deformation. The model’s accuracy and simplicity make it a good option for commercial and industrial applications. Conveniently, due to the model modular nature, it can be expanded in the future to include viscoelasticity and non-isotropic formation for better precision.

Ghaderi, Aref↗

Hexadecapole strength in the rare isotopes 74,76 Kr

In the Ge-Sr mass region, isotopes with neutron number N ≤ 40 are known to feature rapid shape changes with both nucleon number and angular momentum. To gain new insights into their structure, inelastic proton scattering experiments in inverse kinematics were performed on the rare isotopes 74,76 Kr. This work focuses on observables related to the J π = $4$$^{+}_{1}$ states of the Kr isotopes and, in particular, on the hexadecapole degree of freedom. By performing coupled-channels calculations, hexadecapole deformation parameters β 4 were determined for the J π = $4$$^{+}_{1}$ states of 74,76 Kr from inelastic proton scattering cross sections. Two possible coupled-channels solutions were found. A comparison to predictions from nuclear energy density functional theory, employing both non-relativistic and relativistic functionals, clearly favors the large, positive β 4 solutions. These β 4 values are unambiguously linked to the well deformed prolate configuration. Given the β 2 - β 4 trend, established in this work, it appears that β 4 values could provide a sensitive measure of the nuclear shell structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Mechanism of the wurtzite to rocksalt phase transformation in cadmium sulfide single crystals shock compressed along the 𝑐-axis to elastic impact stresses of ∼ 5 GPa

Cadmium sulfide (CdS), which exhibits a wurzite (WZ) to rocksalt (RS) phase transformation at elevated stresses, is an ideal system to address the role of deformation on phase transformation mechanisms and kinetics. CdS has been shown to have very different elastic-inelastic behavior for shock propagation along the WZ 𝑐- and 𝑎-axes, and also exhibits significant differences in the time scale for the phase transformation when shocked along these different axes. As an important first step in examining the role of deformation in shock-induced phase transformations, here we present in situ, time-resolved x-ray diffraction (XRD) measurements in single-event, shock wave experiments on CdS single crystals shocked along the WZ 𝑐-axis to elastic impact stresses near or above ∼ 5 GPa, where a marked increase in the kinetics of this transformation has been reported. The XRD measurements are compared with forward diffraction simulations to evaluate different proposed phase transformation mechanisms and are shown to be consistent with a shearing mechanism that results in the RS [001] direction being aligned with the original WZ 𝑐-axis and the RS [110] direction being aligned with the original WZ 𝑎-axis. In order to reproduce all observed RS diffraction spots the forward diffraction simulations required three different RS crystalline domains, each with a [110] direction aligned with one of the three original WZ 𝑎-axes, and a significant mosaicity, particularly about the RS [001] direction (the direction of shock propagation). In conclusion, this large mosaicity is likely due to variations in the possible shearing mechanisms that connect the WZ and RS structures, which can create crystallites with up to ∼ ±10° of relative rotations about the RS [001] direction, as well as further rotation of these crystallites as they grow, merge, and coalesce.

Crystal structure↗

Inelastic peridynamic model for molecular crystal particles

The peridynamic theory of solid mechanics is applied to modeling the deformation and fracture of micrometer-sized particles made of organic crystalline material. A new peridynamic material model is proposed to reproduce the elastic–plastic response, creep, and fracture that are observed in experiments. The model is implemented in a three-dimensional, meshless Lagrangian simulation code. In the small deformation, elastic regime, the model agrees well with classical Hertzian contact analysis for a sphere compressed between rigid plates. Under higher load, material and geometrical nonlinearity is predicted, leading to fracture. Finally, the material parameters for the energetic material CL-20 are evaluated from nanoindentation test data on the cyclic compression and failure of micrometer-sized grains.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗