Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal model”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Mechanical, Electrochemical & Thermal Modeling of Electric Vehicle Batteries for Crash Simulation (CRADA CRD-19-00811 Final Report)

Under the proposed effort in partnership with Hyundai Motor Company (HMC), the National Laboratory of the Rockies (NLR) will cooperate with HMC to develop mathematical models for battery cells and modules for simulating abuse response in batteries subject to the type of mechanical crushing that can occur in a full motor vehicle crash.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Mechanical, Electrochemical & Thermal Models - Training (CRADA CRD-19-00798 Final Report)

Under the proposed effort in partnership with Hyundai Motor Company (HMC), the National Renewable Energy Laboratory (NLR) will host Dr. Jaeyoung Lim from HMC for a period of one year to jointly develop mathematical models for battery cells and modules subject to mechanical crush. NLR will assist with the development of mathematical models that Dr. Lim will incorporate into his research effort on new concepts of mobility with electric vehicles.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Improved internal short circuit models for thermal runaway simulations in lithium-ion batteries

Thermal runaway (TR) modeling is one of the primary tools that can be used to overcome challenges associated with lithium-ion battery (LIB) safety. Among all LIB accidents that have occurred over the past decade, Internal Short Circuit (ISC) remains the most common trigger mechanism. Many available models in the literature either use a simplified approach to simulate ISC or completely ignore its contribution. The aim of this study is to understand the nature of the heat released for different types of ISC scenarios, including aluminum-anode, anode–cathode, and copper-cathode ISC. We study ISC behavior using a coupled electrochemical–thermal model with an integrated TR chemical kinetics solver built in the COMSOL Multiphysics framework. The time duration of heat release and the magnitude of the peak ISC current are studied as functions of parameters such as the size of the penetrating filament and the capacity of the cell. The numerical results are used to build an empirical model validated against the published experimental TR propagation data. Our model can be successfully used as a viable low-cost substitute in lower order (lumped) TR simulations to enable TR prevention and mitigation.

Singh, Bakhshish Preet (ORCID:0000000264751992)↗
DOE OSTI2026

Machine learning-enhanced hybrid modeling approach for better identification of a building thermal network model and improved prediction

The gray-box modeling approach, which uses a semi-physical thermal network model, has been widely used in building prediction applications, such as model predictive control (MPC). However, unmeasured disturbances, such as occupants, lighting, and in/exfiltration loads, make it challenging to apply this approach to practical buildings. In this word, we propose a hybrid modeling approach that integrates the gray-box model with a model for unmeasured disturbance. After reviewing several system identification approaches, we systematically designed the unmeasured disturbance model with a model selection process based on statistical tests to make it robust. We generated data based on the building model calibrated by real operational data and then trained the hybrid model for two different weather conditions. The hybrid model approach demonstrates an RMSE reduction of approximately 0.2–0.9 °C and 0.3–2 °C on 1-day ahead temperature prediction compared to the Conventional approach for mild (Berkeley, CA) and cold (Chicago, IL) climates, respectively. In addition, this approach was applied to experimental data obtained from the laboratory building to be used for the MPC application, showing superior prediction performances.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Multiscale Modeling of Silicon Carbide Cladding for Nuclear Applications: Thermal Performance Modeling

The complex multiscale and anisotropic nature of silicon carbide (SiC) ceramic matrix composite (CMC) makes it difficult to accurately model its performance in nuclear applications. The existing models for nuclear grade composite SiC do not account for the microstructural features and how these features can affect the thermal and structural behavior of the cladding and its anisotropic properties. In addition to the microstructural features, the properties of individual constituents of the composites and fiber tow architecture determine the bulk properties. Models for determining the relationship between the individual constituents’ properties and the bulk properties of SiC composites for nuclear applications are absent, although empirical relationships exist in the literature. Here, a hierarchical multiscale modeling approach was presented to address this challenge. This modular approach addressed this difficulty by dividing the various aspects of the composite material into separate models at different length scales, with the evaluated property from the lower-length-scale model serving as an input to the higher-length-scale model. The multiscale model considered the properties of various individual constituents of the composite material (fiber, matrix, and interphase), the porosity in the matrix, the fiber volume fraction, the composite architecture, the tow thickness, etc. By considering inhomogeneous and anisotropic contributions intrinsically, our bottom-up multiscale modeling strategy is naturally physics-informed, bridging constitutive law from micromechanics to meso-mechanics and structural mechanics. The effects that these various physical attributes and thermo-physical properties have on the composite’s bulk thermal properties were easily evaluated and demonstrated through the various analyses presented herein. Since silicon carbide fiber-reinforced SiC CMCs are also promising thermal–structural materials with a broad range of high-end technology applications beyond nuclear applications, we envision that the multiscale modeling method we present here may prove helpful in future efforts to develop and construct reinforced CMCs and other advanced composite nuclear materials, such as MAX phase materials, that can service under harsh environments of ultrahigh temperatures, oxidation, corrosion, and/or irradiation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Helium Incorporation into Scandium Fluoride, a Model Negative Thermal Expansion Material

Scandium trifluoride is a model negative thermal expansion (NTE) material. Its simple structure can be described as an A-site vacant perovskite, and it shows isotropic NTE over a very wide temperature range (up to ~1100 K), due to transverse vibrational motion of the fluoride. Like many framework NTE materials, it undergoes a phase transition at low pressures, adopting a rhombohedral (R3̅c) structure at >0.7 GPa and 300 K in commonly used nonpenetrating pressure media, such as silicone oil. High pressure X-ray diffraction data and gas uptake/release measurements indicate that, on compression in helium above ~200 K, helium is inserted into ScF 3 to form the defect perovskite He x ScF 3 . The incorporation of helium stiffens the structure and changes its phase behavior. At room temperature, complete filling of the structure with helium does not occur until >1.5 GPa. On compression, a cubic perovskite structure is maintained until ~5 GPa. As the pressure was increased to ~9.5 GPa, a further transition occurred at ~7 GPa. The first transition at ~5 GPa is likely to a tetragonal (P4/mbm) perovskite, but the detailed structure of the perovskite phase formed on compression above ~7 GPa is unclear. Cooling down from 300 to 100 K in helium at ~0.4 GPa leads to an approximate composition of He 0.1 ScF 3 . High pressure neutron diffraction measurements, in the temperature range 15–150 K show that the incorporation of helium increases the pressure at which the cubic (Pm3̅m) to rhombohedral (R3̅c) putative quantum structural phase transition occurs from close to 0 GPa to ~0.2 GPa at 0 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Heat Transfer in Nuclear Waste Glasses: Measurements and Modeling of Thermal Radiation Properties

We measured and modeled near infrared extinction of nuclear waste glasses from 300 °C to 1150 °C to enable predictive radiation heat transfer and thermal conductivity estimates. A composition and redox informed model resolved contributions from key chromophores (Fe+2-O-Fe+3, V+4, free and bonded ?Si-OH groups) and, when present, spinel particles that can cause strong scattering. The model reproduced measured absorption from room temperature up to 1150 °C, with minor discrepancies near 1 µm (likely due to possible trace impurities) and 2.5 µm (linked to uncertainty in hydroxy groups). Spectra showed silicate melts were semitransparent mainly in the 0.5–4.0 µm window, responsible for radiation thermal conductivity that generally increases with increasing temperature. We quantified the dependence of effective thermal conductivity on dissolved water and provided distributions across >100 LAW/HLW/DFHLW melts at 1150 °C, supporting improved melter heat transfer modeling.

Ferkl, Pavel↗
DOE OSTI2025

Thermal 𝑅⁢𝐶 circuit model for resolving the thermal paradox

Thermal measurements of heat capacity and thermal conductivity in a wide range of insulators and superconductors exhibit a “thermal paradox”: a large linear specific heat reminiscent of neutral Fermi surfaces (associated with fractionalized quasiparticles) in nonmetallic samples that exhibit no corresponding linear temperature coefficient to the thermal conductivity. At first sight, these observations appear to support the formation of a continuum of thermally localized many-body excitations, a form of many-body localization that would be fascinating in its own right. Here, in this work, by mapping thermal conductivity measurements onto thermal 𝑅⁢𝐶 circuits, we argue that the development of extremely long thermal relaxation times, a “thermal bottleneck,” is likely in systems with either many-body localization or neutral Fermi surfaces due to the large ratio between the electron and phonon specific heat capacities. We present a reevaluation of thermal conductivity measurements in materials exhibiting a thermal paradox that can be used in future experiments to deliberate between these two exciting alternatives.

heat transfer↗
DOE OSTI2026

A New 1D Model for Thermal Mixing and Stratification in Advanced Reactor Transients

Thermal mixing and stratification in large pools and enclosures play a critical role in the safety and performance of pool-type nuclear reactors, particularly during transient scenarios involving significant temperature differences between incoming and bulk coolant. Accurate modeling of these phenomena is essential for predicting system behavior and supporting passive safety features such as natural circulation. Here, this paper presents a new 1D model for thermal mixing and stratification, developed and implemented in the SAM code. The model represents a large pool as 1D coolant jet channels and zero-dimensional bulk pool volumes, enabling the simulation of a wide range of flow configurations, including hot and cold jet interactions, stratified layers, and the influence of complex geometries such as ceilings, free surfaces, and internal obstacles. Heat exchange between jet and pool regions is governed by closure relations calibrated against 3D computational fluid dynamics (CFD) simulations. The model improves upon earlier approaches by incorporating time-dependent jet characteristics and capturing the associated delay effects more accurately. Code-to-code comparisons and validation against experimental data from the Thermal Stratification Test Facility demonstrate the model’s accuracy and flexibility. This work offers two key contributions: (1) an efficient and robust method for simulating thermal mixing and stratification at the system level, eliminating the need for external coupling between system analysis codes and CFD, and (2) a significant enhancement of SAM’s capabilities to analyze thermal stratification phenomena in advanced reactor systems.

SAM↗
DOE OSTI2024

Counter Data Paucity through Adversarial Invariance Encoding: A Case Study on Modeling Battery Thermal Runaway

Lithium-ion batteries, widely used for their durability and high energy storage, face the risk of internal short circuits leading to catastrophic thermal runaway events. These events, triggered by external stimuli like mechanical loads, pose safety concerns in applications such as electric vehicles. Detecting and understanding thermal runaway events is crucial, but physics-driven models struggle to explain the non-linear evolution of battery temperature during these events, considering factors like material composition and state-of-charge. Due to the rarity of these events and the cost of data collection, we propose a deep learning (DL) model to predict battery temperature responses during thermal runaway. The challenge lies in the scarcity of data, making traditional DL models prone to overfitting and learning low-quality representations of the complex process.Our approach introduces a novel few-shot architecture that incorporates an adversarially governed invariant encoding process. This architecture aims to distill "invariant" relationships by addressing distributional shifts in data across various battery properties, facilitating the detection of thermal runaway events. Specifically, our results demonstrate that deep learning models conditioned on these "invariant" representations outperform state-of-the-art baselines, achieving a remarkable 96.8% performance improvement in terms of the popular metric MAPE. This framework presents a promising direction for enhancing battery safety modeling, particularly in the context of rare and complex events like thermal runaway. Our code and code and dataset used for the paper are public1.

Tabassum, Anika [ORNL] (ORCID:0000000254600955)↗
DOE OSTI2025

GridLAB-D Technical Support Document: Residential Equivalent Thermal Parameter Model

GridLAB-D is a power distribution systems simulation and analysis tool developed by Pacific Northwest National Laboratory. Although GridLAB-D and its underlying residential load model has been used in a variety of power systems analyses, comprehensive documentation of the derivation of the residential house model has not been formally published. The purpose of this technical support document is to serve as a comprehensive resource for the thermal dynamics modeling and house definition implemented in GridLAB-D. It covers model derivation, parameter definitions, and implementation steps in the codebase.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Hypertriton puzzle in relativistic heavy-ion collisions

The yields of hadrons and light nuclei in relativistic collisions of heavy-nuclei at a center of mass energy of $\sqrt{s_{NN}}$ = 2.6 TeV can be described remarkably well by a thermal distribution of an ideal gas of hadrons and light nuclei interacting only via the decay of resonances. Given the particularly small binding energy of hypertritons relative to the temperature describing the yields (about 156 MeV), one might naturally expect hypertrions to dissociate in medium, making the agreement of hypertriton yields with thermal predictions highly puzzling. The puzzle is compounded by the fact that small binding energy is associated with the large size of the hypertriton. This size is on a similar scale to the overall size of the fireball and much larger than the length scale over which temperatures in the fireball vary over phenomenologically relevant amounts. Here, this paper quantifies the tension this effect causes and shows that it is sufficiently large to render the thermal model inconsistent: its natural assumptions are in conflict with its outputs. The possibility that hypertritons are formed at freeze out as compact objects, quark droplets, that subsequently evolve into hypertritons is considered as a way to resolve the puzzle. It is noted that beyond making the assumption that compact quark droplets form, additional detailed dynamical assumptions which have not been justified are needed to make the thermal model work. The issue of why, despite these issues, the hypertriton is well described by a simple statistical description at freeze out is unresolved. Resolving the hypertriton puzzle is important as it may clarify whether the phenomenological success of the simple thermal model for yields accurately reflects the simple picture of the underlying physics on which it is based.

hydrodynamic models↗
DOE OSTI2024

Thermal-Mechanic Modeling of Fusion Components Using The MOOSE Framework

With fusion energy rapidly developing, there is a great need for an open source tool that can rapidly model different fusion tokamak designs. The open-source codes Multiphysics Object Oriented Simulation Environment (MOOSE), Tritium Migration Analysis Program Version 8 (TMAP8) , and Fusion ENergy Integrated Multiphysi-X (FENIX) are being used to demonstrate thermal-mechanic modeling of fusion components. In this study, an analysis is done on the fusion nuclear science facility (FNSF) blanket, and ITER toroidal magnet and diverter mono block. This study will lay the groundwork for fusion component analysis in the MOOSE framework and FENIX program.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗
DOE OSTI2026

Physics vs structure: A systematic benchmark of learning strategies for multi-zone building thermal dynamics

Recent advances in physics-informed and data-driven machine learning promise improved thermal models for advanced building control, yet there is limited quantitative evidence on when added physics structure and architectural complexity are beneficial. Here, this work presents a systematic benchmark of five representative system identification methods for modeling multi-zone building thermal dynamics: linear state-space models, multi-layer perceptrons, neural state-space models, neural ordinary differential equations, and physically-consistent neural networks. The methods are evaluated across multiple data regimes and zone coupling strategies. Using a high-fidelity multi-zone commercial building emulator, we examine short-term and long-term prediction accuracy, computational efficiency, and ease of development. Our results reveal critical trade-offs between prediction performance, model complexity, and physical consistency. We demonstrate that decoupled, nonlinear black-box models consistently outperform coupled physics-constrained architectures in both predictive accuracy and out-of-distribution robustness in majority of the test cases for the building type considered in the study. Our findings quantify the cost of complexity in building thermal modeling and provide concrete, actionable, scenario-based guidelines for selecting model classes for control-oriented applications.

Building thermal modeling↗
DOE OSTI2025

High-temperature seals for supercritical carbon-dioxide (sCO 2 ) turbines (Final Report)

This is the final report for project DE-FE0031924 titled “High-temperature seals for supercritical carbon-dioxide (sCO 2 ) turbines.” The report provides a summary of the entire project efforts from October 2020 through December 2024 including the high-temperature commercial dry gas seal (DGS) tests and thermal modeling of Task 2, as well as the high-temperature, large-diameter seal design and high-temperature tests of large-diameter seals in Task 3. A key outcome of Task 2 was the testing completion of specially instrumented commercial DGS in the GE-SwRI Apollo sCO 2 compressor (27,000 rpm). Test data from the DGS showed elevated temperatures upwards of 350 o F, which are close to the higher operating temperature limit of the DGS. The temperature measurements provide insight into the expected thermal loads on DGS operating in high-speed sCO 2 compressor and provided test data for validation of an in-house thermal model of the compressor/seal. Under Task 2.0, this report also presents the development of a steady-state conjugate heat-transfer model of the DGS operating in the sCO 2 compressor – a first of its kind model for modeling heat transfer of sCO 2 in an actual operating compressor. The findings of the thermal model show a reasonable match between temperature predictions of the model and the measured temperature data, also pointing out the validity of the approach and assumptions made in modeling the flows, heat transfer coefficients and windage modeling in the rig. Under Task 3.0, this report presents the preliminary design of a large-diameter hybrid face seal (14 inch and 26-inch diameter) for field testing in a land-based GE turbine. The preliminary seal design effort presented in this report under project DE-FE0031924 builds on the development and successful laboratory testing for such large diameter hybrid face seal under the prior DE-FE0024007 project. Key aspects of seal fluid analyses with CFD, mechanical design considerations and assembly considerations in a land-based turbine are presented. Finally, under Task 3.0, this report also presents the continued high-temperature testing of the 14-inch diameter hybrid face seal developed previously under the DE-FE0024007 program. Specifically, test data demonstrating successful non-contact seal operation and seal effective leakage of 0.001-inch with seal inlet temperatures above 700 o F are presented in this report. Successful hybrid seal operation in a laboratory environment for a large diameter (14-inch) seal at temperatures above 700 o F is a major technological milestone for this technology.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Bayesian Optimized Deep Ensemble for Uncertainty Quantification of Deep Neural Networks: a System Safety Case Study on Sodium Fast Reactor Thermal Stratification Modeling

Deep neural networks (DNNs) are increasingly important to scientific computing and engineering system simulations. Accurate uncertainty quantification (UQ) for DNNs is critical in safety-sensitive engineering domains. Traditional Deep Ensemble (DE) methods, while easy to implement, frequently suffer from poorly calibrated uncertainty estimates and limited predictive accuracy due to reliance on fixed architectures with varied weight initializations. To address these issues, we introduce a workflow that combines Bayesian Optimization (BO) and DE. The workflow is modular, scalable, and integrates parallel BO initialized with Sobol sequences to individually optimize the hyperparameters of each ensemble member. This method enhances ensemble diversity, improves predictive accuracy, and provides reliable uncertainty estimates. We evaluate the proposed BODE approach in a sodium fast reactor thermal stratification modeling case study, where we used a densely connected convolutional neural network to predict turbulent viscosity during the reactor transient with consideration of data noise. We benchmark its performance against several optimization approaches, including baseline deep ensemble, evolutionary algorithm-optimized ensemble, ensemble formed via random search combined with greedy selection, and a BO ensemble using random initialization. Here, our results demonstrate superior performance of the developed BODE approach. In noise-free scenarios, BODE notably reduces incorrect aleatoric uncertainty and significantly enhances predictive accuracy. Under conditions of 5% and 10% Gaussian noise, BODE adaptively quantifies uncertainty proportional to data noise, achieving up to an 80% reduction in root mean square error compared to baseline methods and producing well-calibrated prediction intervals.

Bayesian optimization↗
DOE OSTI2024

Latent Pitfalls in Microstructure-Based Modeling for Thermally Aged 9Cr-1Mo-V Steel (Grade 91)

A case study was conducted on a mechanistic model development that predicted tensile strength deterioration with thermal aging of 9Cr-1Mo-V steel in supporting the 60-year design life expected for advanced nuclear reactors. For property prediction beyond practical testing times, mechanistic modeling is highly desired, as it taps into the physics of structure–property relationships and therefore can generate reliable results for extrapolation. Meanwhile, as mechanistic models are often complicated, reflecting the intricacy of microstructure and strengthening mechanisms, pitfalls that are difficult to detect often exist. Here, this paper discusses latent pitfalls that are common in mechanistic modeling or specific in this 9Cr-1Mo-V case development through using the American Society of Mechanical Engineers verification and validation in computational solid mechanics (ASME V&V 10) standard for evaluating credibility of modeling in materials engineering. Suggestions are also made for enhancing reliability of microstructure-based modeling.

36 MATERIALS SCIENCE↗