SEARCH · Search NASA

Results for “Thermomechanical behavior”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Characterization of Fuel-to-Coolant Heat Transfer During Reactivity-Initiated Accidents Using Tightly Coupled Thermal Hydraulics and Fuel Thermomechanics

The reactivity-initiated accident (RIA) is a complex scenario with several tightly interacting physical phenomena. Accurately predicting fuel behavior during these transients is difficult due to limitations in the modeling of fuel-to-coolant heat transfer. Common approaches to simulate RIAs involve standalone calculations using either a fuel performance code or a thermal-hydraulic code. The complex interdependencies of thermal-hydraulic and fuel mechanical behavior suggest that a tight coupling between these codes may provide more accurate predictions of fuel-to-coolant heat transfer and cladding mechanical response. Here, RELAP5-3D and BISON are coupled in this paper to simulate RIAs, and a sensitivity analysis is performed to rank key thermal properties and two-phase heat transfer parameters relevant for fuel-to-coolant heat transfer and cladding failure mechanisms in UO 2 –Zircaloy-4 systems. Gas gap conductance, film boiling heat transfer uncertainty, pulse width, fuel-specific heat capacity, and cladding-specific heat capacity were identified as important parameters. Variations in figures of merit resulting from changes to pulse width and the material thermal properties indicate that time-dependent heat transfer rates are significant for safety-relevant mechanical parameters due to the time dependence of cladding ductility and pellet-cladding mechanical interaction loading. The results suggest that the thermal-hydraulic factors have a nonnegligible influence on the thermomechanical solution and vice versa. Tight coupling of both sets of physics is recommended to improve prediction of fuel behavior during RIAs. Highlights include the following: 1. The RELAP5-3D thermal-hydraulic code and the BISON fuel performance code are tightly coupled for simulation of RIA transients with energy depositions at the Zircaloy-4 cladding failure threshold. 2. Departure from nucleate boiling occurred for all simulated cases. Due to the ductility of fresh fuel, substantial ballooning occurred in most cases. 3. Gas gap conductance, fuel-specific heat capacity, cladding-specific heat capacity, transient pulse width, and film boiling heat transfer were the dominant thermal factors impacting the safety figures of merit at energy depositions.

Critical Heat Flux (CHF)↗

DOE OSTI2025

Recrystallization driven softening and heating rate dependencies of FeCrAl nuclear fuel cladding during accident transients

A refined understanding of FeCrAl cladding behavior during rapid transients is critical for its potential deployment in light-water reactors. Current assessments focus on transient burst testing metrics such as balloon geometry, burst temperature, and hoop stress, often used as proxies for simpler conventional tensile properties. However, directly correlating isothermal tensile and creep data with accident transient scenarios remains a challenge, although it is essential for high fidelity model development. Recent modeling based on tensile tests up to 800 °C, conducted with both immediate loading and a 10-minute soak, showed that immediate loading better predicts experimental burst temperatures, indicating a thermal softening effect. Building upon this observation, the current study connects transient performance, microstructural evolution, and high-temperature tensile properties by leveraging results from C26M claddings burst tests performed at heating rates of 1–50 °C/s and hoop stresses from 25 to 100 MPa. At 25 MPa, rupture temperatures varied by only 6 °C, but at 100 MPa, the difference reached 116 °C, with faster heating yielding higher burst temperatures. Microstructural analysis identified recrystallization as the primary cause of heating rate-dependent softening, eliminating prior cold-working. In-situ thermomechanical data linked ballooning onset to localized instabilities, similar to ultimate tensile strength behavior in conventional tensile tests. High heating rates correlated with immediate loading tensile data, while lower rates matched soaked data. Furthermore, by linking burst performance to microstructural evolution and tensile properties, this work provides a foundation for more accurate modeling of FeCrAl claddings and potentially other Fe-based materials under accident conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Thermomechanical conversion in metals: dislocation plasticity model evaluation of the Taylor-Quinney coefficient

Using a partitioned-energy thermodynamic framework which assigns energy to that of atomic configurational stored energy of cold work and kinetic-vibrational, in this study we derive an important constraint on the Taylor-Quinney coefficient, which quantifies the fraction of plastic work that is converted into heat during plastic deformation. Associated with the two energy contributions are two separate temperatures – the ordinary temperature for the thermal energy and the effective temperature for the configurational energy. We show that the Taylor-Quinney coefficient is a function of the thermodynamically defined effective temperature that measures the atomic configurational disorder in the material. Finite-element analysis of recently published experiments on the aluminum alloy 6016-T4 [1], using the thermodynamic dislocation theory (TDT), shows good agreement between theory and experiment for both stress-strain behavior and temporal evolution of the temperature. The simulations include both conductive and convective thermal energy loss during the experiments, and significant thermal gradients exist within the simulation results. Computed values of the differential Taylor-Quinney coefficient are also presented and suggest a value which differs between materials and increases with increasing strain.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Texture and microstructure evolution in thermomechanically processed Mg-Ca and Mg-Zn-Ca alloys

Mg-Zn-Ca alloys have the potential for producing Mg alloy sheets with weaker basal textures and therefore improved formability. Thermomechanical processing (TMP) also plays a substantial role in determining the final sheet texture. The interplay of these variables, TMP and alloying, complicates comparisons of alloys across the literature. This work systematically explores the texture evolution and recrystallization behavior in Mg-Ca and Mg-Zn-Ca alloys during plane strain compression (PSC) using a Gleeble thermomechanical simulator. It is demonstrated that the basal texture intensity and texture characteristics change significantly during post-deformation annealing, particularly in the ternary alloys. It is also shown that careful selection of the TMP processing variables used during PSC is essential to producing weak textures. In particular, it is important to limit the recrystallization that occurs between compressive hits. This is achieved by adjusting solute content, strain rate, and the duration of the soak between passes.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Incorporating dynamic recrystallization into a crystal plasticity model for high-temperature deformation of Ti-6Al-4V

During hot deformation of (α + β) titanium alloys, the simultaneous action of strain and temperature in the (α + β) regime facilitates dynamic recovery, dynamic recrystallization (DRX), and phase transformations via non-equilibrium paths. DRX is manifested in the form of fine recrystallized α or β grains. Here, in the present study, a two-phase crystal plasticity finite element framework (CP-DRX) has been developed which incorporates DRX kinetics into the crystal plasticity (CP) model to predict the flow characteristics of Ti-6Al-4V alloys during thermo-mechanical processing. The CP slip system parameters, as well as elastic properties from both α and β phases of Ti, are calibrated for different strain rate conditions. An EBSD-informed two-phase microstructure representation has been utilized in the CP-DRX framework to explore the time-dependent evolution of DRX microstructure and crystal orientation for different strain rate conditions. The proposed CP-DRX can capture the evolution of crystal orientation and plastic flow stress–strain response of polycrystalline Ti-6Al-4V during the deformation process. Furthermore, the proposed model is able to capture the softening behavior, observed in average stress–strain response from experiments performed using a Gleeble thermomechanical simulator and predict the recrystallization texture.

(α+ β)-Titanium (Ti) Alloys↗

DOE OSTI2020

Modeling The Effects of Loading Scenario and Thermal Expansion Coefficient on Potential Failure of Cryo-compressed Hydrogen Vessels

A multiscale thermomechanical model for a simplified Type-3 cryogenic compressed hydrogen (H2) storage vessel is developed in this paper. The model accounts for the temperature-dependent elastic-plastic behavior of the vessel carbon/epoxy composite overwrap and aluminum alloy liner. The homogenized thermo-elastic-plastic behavior for the individual laminas of the vessel layup is obtained by an incremental Eshelby-Mori-Tanka approach associated with a micromechanical failure criterion to predict lamina failure while a standard elastic-plastic constitutive model is used to describe the behavior of a typical aluminum alloy assumed for the liner. The vessel response to external loadings is achieved by a finite element method. Four loading scenarios representing four thermomechanical cycles applied to the vessel are analyzed to evaluate constituent and lamina stresses as well as the associate failure criterion during the cycle according to these scenarios. The model can provide helpful guidance to mitigate thermal stresses by an adequate selection of loading scenario, optimizing the layup and by tailoring thermomechanical properties of the resin matrix.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Arctic Critical Infrastructure: Assessing and Predicting the Risk to Critical Permafrost Infrastructure from Climate Change: A New Thermomechanical Approach

This study presents the development of a computational framework designed to predict the interaction between permafrost and infrastructure, addressing potential failure modes and mitigation strategies in the context of climate change. The framework, rooted in advanced modeling and simulation (mod/sim) techniques, integrates thermomechanical coupling to account for the complex interplay between heat flow, ice content, and mechanical behavior in permafrost. Existing models fail to fully capture these dynamics, particularly as they relate to the effects of ice saturation on structural integrity. Our innovative Arctic Coastal Erosion (ACE) framework fills this gap by coupling thermal and mechanical models to accurately simulate subsidence and deformation in permafrost environments. We applied the ACE framework to a representative runway, demonstrating its capability to predict settlement due to rising temperatures and subsequent permafrost thaw. This proof-of-concept showcases the potential of the framework to evaluate risks to Arctic infrastructure, which supports over four million people and 70% of existing permafrost-based structures. By simulating various infrastructure types and environmental conditions, our research offers insights into failure mechanisms and evaluates structural solutions to mitigate risk. The anticipated deliverables, including a prototype runway exemplar, position this project as a critical advancement in permafrost infrastructure modeling, with applications in national security and resilience planning.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A CALPHAD-informed approach to modeling constituent redistribution in Zr-based metallic fuels using BISON

Here, a CALPHAD-informed (Computer Coupling of Phase Diagrams and Thermochemistry) constituent redistribution model was developed for Zr-based metallic fuels and incorporated into the BISON fuel performance code. Three uncertain model parameters associated with β and γ phase kinetics were calibrated using integral test data from U-Zr fuel elements irradiated in Experimental Breeder Reactor II. The calibrated constituent redistribution model was shown to predict the behavior of U-Zr fuels with excellent accuracy. Model predictions for U-Pu-Zr fuels were physically reasonable but less accurate. Reduction of uncertainties in the ternary phase transition temperatures and collection of kinetic data for the ζ phase are expected to improve the model’s ternary predictions. Finally, the new model was coupled to existing thermomechanics models in BISON to simulate irradiation of an entire U-Zr fuel element, demonstrating its ability to accurately predict the behavior of U-Zr fuels with realistic geometries and mesh resolutions at the engineering scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Modeling discontinuous dynamic recrystallization containing second phase particles in magnesium alloys utilizing phase field method

Pre-existing second phase particles can induce particle stimulated nucleation (PSN) and pinning effect on grain boundaries during DRX process of magnesium alloys. The interaction among pinning effect, PSN mechanism and strain-induced grain boundary migration (SIBM) nucleation mechanism imposes challenge in quantitative study of microstructure evolution. A phase field model that describes discontinuous dynamic recrystallization process considering second phase particles (PF-DDRX-SP) which mainly distribute along grain boundaries is developed. The second phase order parameter with diffusion interface is introduced to the free energy function. With considering the pinning effect of second phase particles on grain boundaries, the grain boundary energy couples the interaction between second phase and grains. The effect of second phase on dislocation density evolution is considered in the model, and the PSN and SIBM mechanism are both executed by limiting the nucleation position to the grain boundary and the second phase interface. For the validation, the initial topology consistent with initial microstructure is applied to the simulation of DRX kinetics and flow behavior of AZ80 magnesium alloy, and great reliability and accuracy of developed model is indicated. Furthermore, through the model, the microstructure evolution with different particle size and with different volume fraction of round particle during DRX process were simulated. Finally, the results show that particles with higher dispersion reduce the number of DRX nuclei, which indicate that second phase particles suppress the nucleation at grain boundaries because of the pinning effect. More importantly, the proposed model could also quantitatively predict the relationships between the parameters of second phase particles and DRX behaviors, and enable to optimize the initial second phase structure in a uniform grain structure during thermomechanical process.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Intense cyclic heating effects on thermo-fracture and thermal shock of solid tungsten and open-cell tungsten foam

Here we investigate here the effects of transient (cyclic) arc-jet plasma and laser heating on fracture behavior of W-foam and solid tungsten. The two key parameters that control the foam thermomechanical response are its density and mean cell size expressed in Pores Per Inch (PPI). Tungsten foam samples were fabricated with Chemical Vapor Deposition (CVD) with variety of PPI and relative density. These were tested under two types of qualitatively different conditions: (1) high-enthalpy arc-jet, and (2) high-power cyclic laser heating. None of the foam samples showed macroscopic through-thickness cracks. However, distributed micro-cracks were observed on ligaments and their triple junctions. Under the same loading conditions, W-foam and solid tungsten showed similar crack network pattern and characteristic length-scale. However, Crack Opening Displacement (COD) was twice as large in solid W as compared to W-foam. Foam samples that have been previously exposed to a low-pressure helium plasma showed significant changes in their surface forming nano-texture fuzz which was removed by subsequent testing in the arc-jet. Extensive fracture and re-crystallization were observed in the thin solid W disk that was fully-constrained from expansion. Thicker and fully-constrained solid W disks did not display recrystallization, grain growth, and extensive cracking. However, thicker disks that were free to expand showed some recrystallization and extensive through-thickness cracks due to less effective cooling and thus higher temperatures. Laser beam testing showed no visible damage formation at 0.19 GW/m 2 and 0.38 GW/m 2 for both low-density (23%) and high-density (43%) foams at low pulses (100-1000). Micro-cracks were observed after 10,000 pulses at 0.19 GW/m 2 in both foams, and in low-density foam after 100,000 at 0.38 GW/m 2 . The nature of thermomechanical damage in W-foam exposed to extreme power (GW/m 2 ) short-duration laser pulses was found to be qualitatively similar to that of high power (MW/m 2 ) long-duration arc-jet.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Recycling of natural fiber composites: Challenges and opportunities

Natural fibers have been widely used for reinforcing polymers attributed to their sustainable nature, excellent stiffness to weight ratio, biodegradability, and low cost compared with synthetic fibers like carbon or glass fibers. Thermoplastic composites offer an advantage of recyclability after their service life, but challenges and opportunities remain in the recycling of natural fiber reinforced polymer composites (NFRPCs). This article summarized the effects of reprocessing/recycling on the material properties of NFRPCs. The material properties considered include mechanical performance, thermal properties, hygroscopic behavior, viscoelasticity, degradation, and durability. NFRPCs can generally be recycled approximately 4–6 times until their thermomechanical properties change. After recycling 7 times, the tensile strength of NFRPCs can decrease by 17%, and the tensile modulus can decrease by 28%. The mitigation approaches to overcome degradation of material properties of NFRPCs such as adding functional additives and virgin plastics are also discussed. The main challenges in these approaches such as degradation and incompatibility are discussed, and an effort is made to provide a rationale for reprocessing/recyclability assessment. Finally, future applications of NFRPCs such as additive manufacturing and automotive part use are discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2020

TMS XTRRA experiments on NIF

The study of x-ray induced rapid heating and the resulting material responses - (thermomechanical shock (TMS) and thermo-structural response (TSR)) - is important to our understanding of material behavior in extreme environments. A series of high conversion efficiency (XRCE) x-ray sources have been developed for the NIF laser facility ranging in fluence and spectral content. In conjunction with source development, we have developed a new XTRRA (X-ray Transport and Radiation Response Analysis) test cassette to hold a set of samples, each at equal distance from target chamber center (TCC) where the x-ray source is placed. Design and performance details of the x-ray sources is discussed elsewhere. Here, we describe the development and use of this TMS XTRRA test cassette. The XTRRA test cassette is a snout mounted on a NIF Diagnostic Instrument Manipulator (DIM) designed to position six samples equidistant from the x-ray source to assure uniform irradiation of all samples. Photon Displacement interferometry (PDI) is used to measure the displacement of the rear surface of each sample resulting from x-ray exposure on the front of the sample. This will help determine the x-ray generated impulse on the sample. The samples are retrieved from the snout after the shot for inspection and post-shot analysis including the measurement of sample ablated mass and engineering tests examining changes in the material structure. The data from experiments are compared to predictions from LLNL’s ASC codes and used to inform models implemented in these codes. Thus, XRCE x-ray sources with the XTRRA snout at the National Ignition Facility (NIF) laser at the Lawrence Livermore National Laboratory [1,2] provide the ability to expose a broad range of material samples directly to high x-ray fluence, measure shock transits, and recover the sample for further analysis post-shot.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Thermo-Mechanical Shock, X-Ray Transport and Radiation Response Analysis Experiments on NIF

The study of x-ray induced rapid heating and the resulting material responses - (thermomechanical shock (TMS) and thermo-structural response (TSR)) - is important to our understanding of material behavior in extreme environments. A series of high conversion efficiency (XRCE) x-ray sources have been developed for the NIF laser facility ranging in fluence and spectral content. In conjunction with source development, we have developed a new XTRRA (X-ray Transport and Radiation Response Analysis) test cassette to hold a set of samples, each at equal distance from target chamber center (TCC) where the x-ray source is placed. Design and performance details of the x-ray sources is discussed elsewhere. Here, we describe the development and use of this TMS XTRRA test cassette. The XTRRA test cassette is a snout mounted on a NIF Diagnostic Instrument Manipulator (DIM) designed to position six samples equidistant from the x-ray source to assure uniform irradiation of all samples. Photon Displacement interferometry (PDI) is used to measure the displacement of the rear surface of each sample resulting from x-ray exposure on the front of the sample. This will help determine the x-ray generated impulse on the sample. The samples are retrieved from the snout after the shot for inspection and post-shot analysis including the measurement of sample ablated mass and engineering tests examining changes in the material structure. The data from experiments are compared to predictions from LLNL’s ASC codes and used to inform models implemented in these codes. Thus, XRCE x-ray sources with the XTRRA snout at the National Ignition Facility (NIF) laser at the Lawrence Livermore National Laboratory provide the ability to expose a broad range of material samples directly to high x-ray fluence, measure shock transits, and recover the sample for further analysis post-shot.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

An Improved Fracture Mechanics-Informed Multiscale Thermomechanical Damage Model for Ceramic Matrix Composites

This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.

Skinner, Travis↗

DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗

DOE OSTI2025

Influence of helium bubbles on the bulk equation of state of gold under static compression

Nuclear materials often evolve into two-phase systems comprising a bulk matrix with dispersed inert-gas bubbles. The presence of these bubbles can have consequences to the thermomechanical response of materials and is a key life-limiting factor in some nuclear fuel forms. Understanding the behavior of these two-phase, bubble-matrix systems is, thus, important to improved predictive models and frameworks for many nuclear materials applications. While temperature excursions of these two-phase systems have been characterized, fewer studies have focused on the evolution of inert-gas bubbles under pressure. Here, in this paper, we use x-ray tools to interrogate a He-implanted gold foil to determine the pressure-dependent evolution of the individual components (Au matrix + bubbles), and we compare that total pressure dependence to theoretical equation-of-state descriptions based on mixing rules.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

36 MATERIALS SCIENCE↗