SEARCH · Engineering Papers

Results for “thermomechanical modeling”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Analysis of shell-type structures subjected to time-dependent mechanical and thermal loading

A general mathematical model and solution methodologies are being developed for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic, or static thermomechanical loads. Among the system responses, which were associated with these load conditions, were thermal buckling, creep buckling, and ratcheting. Thus, geometric as well as material-type nonlinearities (of high order) can be anticipated and must be considered in the development of the mathematical model. Furthermore, this must also be accommodated in the solution process.

Simitses, G. J.↗

DOE OSTI2021

Microstructure-dependent rate theory model of defect segregation and phase stability in irradiated polycrystalline LiAlO 2

We report gamma lithium aluminate (LiAlO 2 ) is a breeder material for tritium and is one of key components in a tritium-producing burnable absorber rod (TPBAR). Dissolution and precipitation of second phases such as LiAl 5 O 8 and voids are observed in irradiated LiAlO 2 . Such microstructure changes cause the degradation of thermomechanical properties of LiAlO 2 and affect tritium retention and release kinetics, and hence, the TPBAR performance. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed for investigating the accumulation of species and phase stability in polycrystalline LiAlO 2 structures under irradiation. Three sublattices (i.e. [Li, Al, V] I [O, V o ] II [Li i , Al i , O i , V i ] III ), and concentrations of six diffusive species (i.e. Li; vacancy of Li or Al at [Li, Al, V] I sublattice, O vacancy at [O, V o ] II sublattice, and Li, Al and O interstitials at [Li i , Al i , O i , V i ] III interstitial sublattices; are used to describe spatial and temporal distributions of defects and chemistry. Microstructure-dependent thermodynamic and kinetic properties including the generation, reaction, and chemical potentials of defects and defect mobility are taken into account in the model. The parametric studies demonstrated the capability of the developed RIS model to assess the effect of thermodynamic and kinetic properties of defects on the segregation and depletion of species in polycrystalline structures and to explain the phase stability observed in irradiated LiAlO 2 samples. The developed RIS model will be extended to study the precipitation of LiAl 5 O 8 and voids and tritium retention by integrating the phase-field method.

36 MATERIALS SCIENCE↗

NASA NTRS1983

Thermomechanical loading of multilayered cylindrical geometries in thermal cycling from 300 to 1300 K

The principal of multiple material layering is well known as an effective method of reducing heat transfer; however, thermal gradients can impose significant loads and lead to delamination and subsequent component failure. An analysis is developed and experimental data are discussed for the thermomechanical effects of multilayered materials on a heat sink substrate of cylindrical geometry subject to thermal cycling. The geometry is heated in cross-flow by a high-velocity flame and cooled in crossflow by ambient-temperature air from a critical flow orifice. Each layer of material possesses a threshold beyond which small changes in temperature or mechanical loading greatly influence the life in thermal cycling of the layered materials. Comparisons are made between the thermomechanical loads predicted by various numerical codes for the linear case and by a simplified analytic model.

Hendricks, R. C.↗

NASA NTRS1974

A thermal, thermoelastic, and wear analysis of high-energy disk brakes

A thermomechanical investigation of the sliding contact problem encountered in high-energy disk brakes is described. The analysis includes a modelling, using the finite element method of the thermoelastic instabilities that cause transient changes in contact area to occur on the friction surface. In order to include the effect of wear at the contact surface, a wear criterion is proposed that results in the prediction of wear rates for disk brakes that are quite close to experimentally determined wear rates. The thermal analysis shows that the transient temperature distribution in a disk brake assembly can be determined more accurately by use of this thermomechanical analysis than by a more conventional analysis that assumes constant contact conditions. It also shows that lower, more desirable, temperatures in disk brakes can be attained by increasing the volume, the thermal conductivity, and, especially, the heat capacity of the brake components.

Kennedy, F. E., Jr.↗

DOE OSTI2025

AGR TRISO Fuel Performance Modeling in FY-25

This report summarizes the activities performed in FY-25 to support fuel performance modeling for the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program. This includes implementation of a new uranium oxycarbide (UCO) kernel swelling rate model, inner pyrolytic carbon (IPyC) cracking behavior and failure predictions, development of new UCO and silicon carbide (SiC) cesium diffusion parameters, and the thermomechanical behavior of particle layers using experimental data from micro-tensile strength testing.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗

NASA NTRS1986

Thermal barrier coating life prediction model development

A methodology is established to predict thermal barrier coating life in a environment similar to that experienced by gas turbine airfoils. Experiments were conducted to determine failure modes of the thermal barrier coating. Analytical studies were employed to derive a life prediction model. A review of experimental and flight service components as well as laboratory post evaluations indicates that the predominant mode of TBC failure involves thermomechanical spallation of the ceramic coating layer. This ceramic spallation involves the formation of a dominant crack in the ceramic coating parallel to and closely adjacent to the topologically complex metal ceramic interface. This mechanical failure mode clearly is influenced by thermal exposure effects as shown in experiments conducted to study thermal pre-exposure and thermal cycle-rate effects. The preliminary life prediction model developed focuses on the two major damage modes identified in the critical experiments tasks. The first of these involves a mechanical driving force, resulting from cyclic strains and stresses caused by thermally induced and externally imposed mechanical loads. The second is an environmental driving force based on experimental results, and is believed to be related to bond coat oxidation. It is also believed that the growth of this oxide scale influences the intensity of the mechanical driving force.

Demasi, J. T.↗

DOE OSTI2022

ExaAM: Metal additive manufacturing simulation at the fidelity of the microstructure

Additive manufacturing (AM), or 3D printing, of metals is transforming the fabrication of components, in part by dramatically expanding the design space, allowing optimization of shape and topology. However, although the physical processes involved in AM are similar to those of welding, a field with decades of experimental, modeling, simulation, and characterization experience, qualification of AM parts remains a challenge. The availability of exascale computational systems, particularly when combined with data-driven approaches such as machine learning, enables topology and shape optimization as well as accelerated qualification by providing process-aware, locally accurate microstructure and mechanical property models. We describe the physics components comprising the Exascale Additive Manufacturing simulation environment and report progress using highly resolved melt pool simulations to inform part-scale finite element thermomechanics simulations, drive microstructure evolution, and determine constitutive mechanical property relationships based on those microstructures using polycrystal plasticity. We report on implementation of these components for exascale computing architectures, as well as the multi-stage simulation workflow that provides a unique high-fidelity model of process–structure–property relationships for AM parts. In addition, we discuss verification and validation through collaboration with efforts such as AM-Bench, a set of benchmark test problems under development by a team led by the National Institute of Standards and Technology.

3D printing↗

NASA NTRS1985

Constitutive modeling for isotropic materials (HOST)

This report presents the results of the second year of work on a problem which is part of the NASA HOST Program. Its goals are: (1) to develop and validate unified constitutive models for isotropic materials, and (2) to demonstrate their usefulness for structural analyses of hot section components of gas turbine engines. The unified models selected for development and evaluation are that of Bodner-Partom and Walker. For model evaluation purposes, a large constitutive data base is generated for a B1900 + Hf alloy by performing uniaxial tensile, creep, cyclic, stress relation, and thermomechanical fatigue (TMF) tests as well as biaxial (tension/torsion) tests under proportional and nonproportional loading over a wide range of strain rates and temperatures. Systematic approaches for evaluating material constants from a small subset of the data base are developed. Correlations of the uniaxial and biaxial tests data with the theories of Bodner-Partom and Walker are performed to establish the accuracy, range of applicability, and integability of the models. Both models are implemented in the MARC finite element computer code and used for TMF analyses. Benchmark notch round experiments are conducted and the results compared with finite-element analyses using the MARC code and the Walker model.

Lindholm, U. S.↗

DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗

NASA NTRS1984

Assessment of the operating characteristics of the SSME LOX turbopump pump-end bearing

A bearing/shaft model of the SSME LOX turbopump was developed using the SHABERTH bearing/shaft math modeling computer code. A previously developed bearing/shaft thermal model of the SSME LOX turbopump turbine and bearing was used in conjunction with SHABERTH to evaluate the thermomechanical operating characteristics of the LOX turbopump end bearings. Results show that for the two unmounted diametrical clearances evaluated (4.0 mils and 6.3 mils), the inboard pump end bearing supports about 81% of the isolator load for the small clearance and 77% of the isolator load for the larger clearance. Bearing clearance changes due to thermal effects were 40% for the 4.0 mil diametrical clearance case and 19% for the 6.3 mil clearance case evaluated. The thermal analysis included evaluation of bearing temperatures for a subcooled case and a saturated case. Results indicate that no drastic temperature change occurred between the two cases. Since the rolling element and race surfaces of the subcooled case were at temperatures sufficiently high enough to be vapor blanketed, exceeding saturation temperature at the bearing inlet did not increase surface temperatures greatly.

New, L. S.↗

NASA NTRS1991

Thermo-mechanical response predictions for metal matrix composite laminates

An analytical micromechanical model is employed for prediction of the stress-strain response of metal matrix composite laminates subjected to thermomechanical loading. The predicted behavior of laminates is based upon knowledge of the thermomechanical response of the transversely isotropic, elastic fibers and the elastic-viscoplastic, work-hardening matrix. The method is applied to study the behavior of silicon carbide/titanium metal matrix composite laminates. The response of laminates is compared with that of unidirectional lamina. The results demonstrate the effect of cooling from a stress-free temperature and the mismatch of thermal and mechanical properties of the constituent phases on the laminate's subsequent mechanical response. Typical results are presented for a variety of laminates subjected to monotonic tension, monotonic shear and cyclic tensile/compressive loadings.

Aboudi, J.↗

DOE OSTI2026

Advanced Power Electronics and Electric Machines

NLR's advanced power electronics and electric machines (APEEM) research group develops world-class experimental and modeling capabilities to design and evaluate efficient and reliable power electronics and electric machines, as well as innovative thermal management systems to maximize their performance. They also design, fabricate, and characterize advanced power electronics packaging and are developing state-of-health monitoring techniques.

25 ENERGY STORAGE↗

NASA NTRS2008

Measurement and Prediction of the Thermomechanical Response of Shape Memory Alloy Hybrid Composite Beams

An experimental and numerical investigation into the static and dynamic responses of shape memory alloy hybrid composite (SMAHC) beams is performed to provide quantitative validation of a recently commercialized numerical analysis/design tool for SMAHC structures. The SMAHC beam specimens consist of a composite matrix with embedded pre-strained SMA actuators, which act against the mechanical boundaries of the structure when thermally activated to adaptively stiffen the structure. Numerical results are produced from the numerical model as implemented into the commercial finite element code ABAQUS. A rigorous experimental investigation is undertaken to acquire high fidelity measurements including infrared thermography and projection moire interferometry for full-field temperature and displacement measurements, respectively. High fidelity numerical results are also obtained from the numerical model and include measured parameters, such as geometric imperfection and thermal load. Excellent agreement is achieved between the predicted and measured results of the static and dynamic thermomechanical response, thereby providing quantitative validation of the numerical tool.

Davis, Brian↗

DOE OSTI2024

Temperature-dependent mechanical properties and crystal plasticity parameters for additively manufactured Haynes-214 alloy: Experiments and numerical modeling

Our experimental mechanical testing data demonstrated that the additively manufactured (AM) laser powder bed fusion (L-PBF) Haynes-214 alloy exhibits non-linear mechanical properties as the temperature rises from ambient to 870 °C. Crystal plasticity (CP) simulations provide an effective approach to gaining deeper insights into microstructure-property linkages under thermomechanical loading. This method can reduce the need for costly high-temperature mechanical testing while accounting for the effects of crystallographic texture and grain morphology on the mechanical behavior of AM materials. However, calibrating a CP model is time-consuming because individual simulations are computationally expensive and hundreds (or more) of iterations over parameter sets may be required. To address this issue, we have designed a machine learning-differential evolution (ML-DE) CP framework that can accurately interpolate the tensile properties of AM L-PBF Haynes-214 alloy across a wide temperature range from ambient to 870 °C, with minimal reliance on experimental data. The framework uses electron backscatter diffraction (EBSD) measurements to generate statistically equivalent microstructural volume elements to serve as inputs to the CP modeling framework. Stress–strain curves were generated from 1000 CP simulations, which serve as the training data set for the three ML regression algorithms explored: linear, extra-trees, and multi-layer perceptron. These three regression models were independently evaluated to compare their efficiency and identify the most suitable algorithm for the given problem. Results revealed that the extra-trees ML regressor outperforms the other models in both qualitative and quantitative aspects with an R 2 of 0.98. Subsequently, the differential evolution optimization approach is employed to calibrate the ML-based CP material parameters with experimental results obtained at various temperatures. Finally, temperature-dependent CP material parameters are formulated. The effectiveness and efficiency of the designed framework are validated through comparison with experimental results, demonstrating a high degree of agreement. These calibrated parametric constitutive equations enable further use of the CP model to study the deformation behavior of this alloy under a wide range of thermo-mechanical loading conditions.

36 MATERIALS SCIENCE↗

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 OSTI2025

Step-loaded creep testing of Zircaloy-4 cladding at higher temperatures in the α-phase

A refined understanding of zirconium-based cladding thermomechanical performance during rapid transients is essential for enhancing the safety and operation of light-water reactors. Traditional models for zirconium alloys under accident conditions generally assume that creep dominates fuel cladding performance. Here, these historic models have largely remained unchanged and serve as the basis for safety criteria development. As the U.S. nuclear industry pursues higher burnup levels, the increased release of fission gases during transients raises the risk of cladding failure in the low-temperature hcp α-phase, making the fidelity of these models of greater importance. Creep testing was conducted from 550–700°C with 25–120 MPa applied hoop stresses to investigate Zircaloy-4 deformation at accident-relevant temperatures in the α-phase. Step-loading was employed to capture creep behavior across a wide stress range from a single sample. The stress-strain rate data at higher temperatures (650 and 700°C) were well-described by isotropic versions of the Erbacher and Kaddour models, while the lower temperature data (550 and 600°C) were underpredicted by both anisotropic and isotropic model variants. Greater strain rates during the initial loading step at 650 and 700°C were attributed to recrystallization and grain growth of sub-micron crystallites. Yet, texture analysis revealed the basal split texture remained after testing. These observations produced results suggesting Zircaloy-4 claddings experience higher creep rates across the α-phase than previously thought, possibly related to dynamic anisotropy due to temperature dependent activation of deformation mechanisms, effects of biaxial loading, and variation in material condition between the current testing used in previous model development.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Identifying multiple synergistic factors on the susceptibility to stress relaxation cracking in variously heat-treated weldments

The 347H austenitic stainless steel has been widely used for pressure vessels and pipeline (PVP) applications due to its excellent creep and corrosion resistance, which fit ideally to the harsh conditions in petrochemical industries, fossil fuel or nuclear power plants, and modern energy storages. However, a failure mode has been commonly observed with cracks emerging at the heat affected zone (HAZ) of weldments during post-weld heat treatment (PWHT) or under intermediate to high temperature service conditions. This phenomenon is termed as Stress Relaxation Cracking (SRC) since the purpose of PWHT is to relieve the welding-induced residual stress fields, or as Stress Age Cracking (SAC) if failure happens during service. A leading literature explanation of this failure suggests that the residual stress relaxation and the precipitation dissolution and/or re-precipitation occur in the same temperature range, which can lead to locally high strains and thus to crack at the grain boundaries. Since in situ spatial measurements of residual stress fields, microstructural evolution, and failure processes are nearly infeasible, this work recourses to a micromechanical finite element framework that models the high temperature failure as the nucleation and growth of grain boundary cavities, whereas various parameters such as thermomechanical loading history and its evolution, the competition of grain-interior dislocation creep and grain-boundary diffusion in failure lifetime, and microstructural heterogeneities (such as the precipitate free zone near grain boundaries) can be quantitatively incorporated. It can be concluded from these microstructure-explicit simulations that an accurate knowledge of residual stress evolution and a carefully calibrated set of material constitutive parameters are the essential prerequisites for lifetime predictions. The understanding of individual governing factors also leads to a mechanistic interpretation of the observed SRC susceptibility C-curves. In conclusion, these results suggest that the criticality of residual stress evolution, but not the precipitation-induced local strains, be the leading factor for SRC.

347H stainless steel weldments↗