Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal 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.

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

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Advancing Geothermal Research: Fiscal Year 2025 Accomplishments Report

This is a summary of geothermal work done at the National Renewable Energy Laboratory (NREL) in Fiscal Year 2025. This year brought increased attention to the geothermal industry and NREL's geothermal research portfolio. With more than 70 active projects, NREL research spanned the areas of resource exploration and characterization; conventional and next-generation geothermal technologies; subsurface thermal energy storage; heating and cooling; co-production of geothermal with critical minerals and oil and gas; modeling and analysis leveraging expertise in data science and machine learning; and more.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Demand response of loads having thermal reserves

Systems and methods are described herein that improve grid performance by smoothing demand using thermal reserves. The smoothed demand can reduce peak loads as well as the ramp rate of demand that will otherwise require the use of inefficient, expensive generation sources. These improvements are tied to the selective switching on or off electrical loads that are coupled to thermal reserves, effectively using the thermal reserves as an energy storage mechanism. Historical data of past usage can be used to create load model and ensure that effects on customer comfort are minimized while still accomplishing the beneficial effects for the overall grid, which enables grid owners to both reduce their operational cost by avoiding expensive generation and improve system reliability by achieving more predictable power demand.

Ren, Wei↗
DOE OSTI2026

Advancing Geothermal Research: Fiscal Year 2025 Accomplishments Report

This is a summary of geothermal work done at the National Laboratory of the Rockies in Fiscal Year 2025. This year brought increased attention to the geothermal industry and NLR's geothermal research portfolio. With more than 70 active projects, NLR research spanned the areas of resource exploration and characterization; conventional and next-generation geothermal technologies; subsurface thermal energy storage; heating and cooling; co-production of geothermal with critical minerals and oil and gas; modeling and analysis leveraging expertise in data science and machine learning; and more.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Unraveling electronic correlations in warm dense quantum plasmas

The study of matter at extreme densities and temperatures has emerged as a highly active frontier at the interface of plasma physics, material science and quantum chemistry with relevance for planetary modeling and inertial confinement fusion. A particular feature of such warm dense matter is the complex interplay of Coulomb interactions, quantum effects, and thermal excitations, making its rigorous theoretical description challenging. Here, we demonstrate how ab initio path integral Monte Carlo simulations allow us to unravel this intricate interplay for the example of strongly compressed beryllium, focusing on two X-ray Thomson scattering data sets obtained at the National Ignition Facility. We find excellent agreement between simulation and experiment with a very high level of consistency between independent observations without the need for any empirical input parameters. Our results call into question previously used chemical models, with important implications for the interpretation of scattering experiments and radiation hydrodynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

HFIR LEU High Density Silicide Dispersion Optimized Design Neutronics Analyses with PHAME

A high-fidelity neutronics model of the Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) with the low-enriched uranium (LEU) high-density silicide dispersion Optimized fuel design was updated and analyzed to generate reactor physics-based metrics to support follow-on thermal hydraulic and transient analyses of this design. The Python HFIR Analysis and Measurement Engine (PHAME) was also updated to enhance the automation capabilities of the framework developed and maintained to perform these reactor physics modeling and simulation efforts. The automated framework significantly increases the efficiency and reproducibility to design and thoroughly analyzes HFIR LEU core designs, changes, and uncertainties. Reactor physics metrics evaluated include but are not limited to fuel depletion, cycle length, fission rate density distributions, axial power peaking factors, kinetics data, reactivity coefficients, control element worths, heat deposition rates, and decay heat. These neutronics results provide essential input to follow-on steady state thermal, thermal hydraulic and reactor transient analyses, which are subject of other reports. The Optimized design operates at 95 MW to maintain HFIR’s current highly enriched uranium core performance level at 85 MW.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRSDate not supplied

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗
NASA NTRSDate not supplied

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗
DOE OSTI2026

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗
DOE OSTI2025

Clustering of DESI galaxies split by thermal Sunyaev-Zeldovich effect

The thermal Sunyaev-Zeldovich (tSZ) effect is associated with galaxy clusters - extremely large and dense structures tracing the dark matter with a higher bias than isolated galaxies. We propose to use the tSZ data to separate galaxies from redshift surveys into distinct subpopulations corresponding to different densities and biases independently of the redshift survey systematics. Leveraging the information from different environments, as in density-split and density-marked clustering, is known to tighten the constraints on cosmological parameters, like $\Omega_m$, $\sigma_8$ and neutrino mass. We use data from the Dark Energy Spectroscopic Instrument (DESI) and the Atacama Cosmology Telescope (ACT) in their region of overlap to demonstrate informative tSZ splitting of Luminous Red Galaxies (LRGs). We discover a significant increase in the large-scale clustering of DESI LRGs corresponding to detections starting from 1-2 sigma in the ACT DR6 + Planck tSZ Compton-$y$ map, below the cluster candidate threshold (4 sigma). We also find that such galaxies have higher line-of-sight coordinate (and velocity) dispersions and a higher number of close neighbors than both the full sample and near-zero tSZ regions. We produce simple simulations of tSZ maps that are intrinsically consistent with galaxy catalogs and do not include systematic effects, and find a similar pattern of large-scale clustering enhancement with tSZ effect significance. Moreover, we observe that this relative bias pattern remains largely unchanged with variations in the galaxy-halo connection model in our simulations. This is promising for future cosmological inference from tSZ-split clustering with semi-analytical models. Thus, we demonstrate that valuable cosmological information is present in the lower signal-to-noise regions of the thermal Sunyaev-Zeldovich map, extending far beyond the individual cluster candidates.

Astronomy data analysis↗
DOE OSTI2025

Distinctive features of fluorescence and waveguides in magnesium aluminate spinel crystals driven by structural discrepancy

Transparent polycrystalline ceramics are of significant importance for a wide range of scientific and industrial applications. Developing a deeper understanding of their thermodynamic behavior is essential for achieving their maximum output performance in technological applications. This study provides a systematic investigation into the thermal excitation-induced fluorescence kinetics and waveguide characteristics of Magnesium Aluminate Spinel (MgAl 2 O 4 ) single crystals, with a focus on the intricate relationship between structure and properties. High-temperature extreme environments through irradiation with swift heavy ions 645.0 MeV Xe and 352.8 MeV Fe ions were created; the atomic deposition energy threshold (E th ) associated with disorder morphologies was assessed between 0.91 and 0.99 eV atom –1 . The track prediction model was developed to support the theoretical prediction of track formation. Electronic energy loss (E ele ) disrupts the balance of the initial structure through the thermal spike effect, leading to the formation of absorption-related F and F + color centers. These defects enhance photoluminescence in the visible spectrum and result in an effective modulation of the intrinsic bandgap. Moreover, as ion beams penetrate into the material, the uneven damage distribution induces the formation of waveguide structures. In conclusion, these findings provide valuable insights into the fabrication of functional devices through irradiation technologies and the structural changes of MgAl 2 O 4 at high temperatures within extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Large area transparent refractory aerogels with high solar thermal performance

Application of transparent silica aerogels in low-temperature solar thermal systems has led to major improvements in performance. In high temperature concentrating solar thermal (CST) systems, aerogels have yet to demonstrate the necessary scalability, durability, and performance to support their widespread deployment. Here, large-area transparent refractory aerogel tiles are synthesized and shown to achieve a record-high receiver figure-of-merit (FOM) at high temperatures. The work leverages a scaled-up process for sol–gel synthesis to control the density of the aerogels for improved solar transmittance and adapts a previous atomic layer deposition (ALD) technique with the aid of predictive reaction-transport modeling. After aging for 10 days at 700 °C, the large-area tiles exhibit a solar-weighted transmittance of 95.6 % and a thermal emittance of 0.31, corresponding to a FOM of 80 % at 100 suns and 700 °C. The observed sintering rates at 700 °C are comparably low to earlier one-inch aerogels, suggesting long-term stability under relevant operating conditions. Furthermore, the study indicates that refractory aerogels are scalable materials for efficient photothermal conversion at high temperatures.

Aerogels↗
DOE OSTI2025

Thermal disorder and phonon softening in the ferroelectric phase transition of lead titanate

We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTi⁢O3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Additive Manufacturing with Cellulose-Based Composites: Materials, Modeling, and Applications

Recent advances in large-scale additive manufacturing (AM) with polymer-based composites have enabled efficient production of high-performance materials. Cellulose nanomaterials (CNMs) have emerged as bio-based feedstocks due to their exceptional strength and sustainability. However, challenges such as hornification and poor dispersion in polymer matrices still limit large-scale CNM–polymer composite manufacturing, requiring novel strategies. Here, this review outlines an approach starting with atomic-level simulations to link molecular composition to key parameters like bulk density, viscosity, and modulus. These simulations provide data for finite element analysis (FEA), which informs large-scale experiments and reduces the need for extensive trials. The strategy explores how atomic interactions impact the morphology, adhesion, and mechanical properties of CNM-based composites in AM processes. The review also discusses current developments in AM, along with predictions of mechanical and thermal properties for structural applications, packaging, flexible electronics, and hydrogel scaffolds. By integrating experimental findings with molecular dynamics (MD) simulations and finite element modeling (FEM), valuable insights for material design, process optimization, and performance enhancement in CNM-based AM are provided to address ongoing challenges.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhancing heat pump water heater performance with embedded phase change materials thermal energy storage: First hour rating improvement and demand response operation

The increasing global emphasis on energy efficiency and sustainability has put heat pump water heaters (HPWHs) in the spotlight as an energy-efficient alternative to traditional water heating systems. However, their widespread adoption is limited by challenges such as insufficient First Hour Rating (FHR), suboptimal control mechanisms, and limited flexibility for demand response operations. Here, to address these limitations, this study proposes an innovative HPWH system integrated with embedded phase change material (PCM)-based thermal energy storage (TES). The research introduces a novel design and control strategy that leverages optimized PCM integration to enhance thermal storage capacity, improve hot water delivery during peak demand, and increase load-shifting potential. A combination of system modeling, performance simulation, and demand response control strategy evaluation was employed to quantify the benefits of PCM integration. Results demonstrate that the proposed PCM-TES HPWH system significantly enhances FHR, with an optimal 7.0 lb. of PCM increasing FHR by over 26 %—from 62 to 78 gal—for a standard 50-gal HPWH. Additionally, under advanced demand response operation using a preheat strategy, the system reduces the percentage of control temperature out-of-band time from 65 % (conventional HPWH) to just 11.6 %, enabling a more stable and efficient hot water supply. This research contributes a novel PCM-embedded HPWH design and control framework that addresses both performance and grid-interactivity challenges. The findings offer a viable pathway to enhancing the operational efficiency, flexibility, and grid responsiveness of residential water heating systems.

Demand response control↗
DOE OSTI2025

Subvoxel Control of Fiber Orientation via Multidirectional Shearing in 3D Printing

Anisotropy, the characteristic of materials exhibiting different properties based on their direction, is widespread in nature. Conventional manufacturing techniques often fall short in recreating such complex anisotropy. While 3D printing allows for precise fiber deposition in anisotropic composites, previous studies have only achieved bulk reorientation of fibers at the voxel level in two dimensions. This limits the replication of localized 3D anisotropies found in natural materials. To address this, a novel 3D printing technique is presented that enables subvoxel control of fiber orientation in all three directions using multidirectional shearing via nozzle rotation and inclination. The fiber orientation control in these experimental tests is driven by a numerical model, enabling a fully digital approach to program microstructures within a strand. This programmability is demonstrated through mechanical and thermal tests, illustrating a localized and controllable response to external stimuli. Achieving such complex anisotropy holds potential across several fields, including wearables and biomedical implants, lightweight composite structures, and energy storage technologies such as batteries and supercapacitors.

3D printing↗
DOE OSTI2026

Machine learning surrogates for ion energy–angle distributions in thermal and RF plasma sheaths

Ion energy–angle distributions (IEADs) at material surfaces are a critical input for plasma–material interaction (PMI) studies in fusion devices, yet they are computationally expensive to obtain using particle-in-cell (PIC) simulations. In this work, we develop a machine learning surrogate based on a deep deconvolutional neural network (DDeCNN) trained on large databases generated with the hPIC2 code. The surrogate is capable of reconstructing IEADs from sheath parameters for both thermal and radio-frequency (RF) plasmas, including cases with multiple ion species. Across thousands of test cases, the model achieves high accuracy, with over 97 % of predictions classified as good or average based on standard error metrics (MAE, MSE, L2). Even in the more challenging RF and multi-species regimes, the surrogate reliably captures the multi-peak structure of PIC results. Once trained, the surrogate produces IEADs in milliseconds on a common workstation, yielding speedups of six to seven orders of magnitude compared with running a full PIC simulation. This computational gain enables dense parameter scans and direct coupling of IEAD predictions with PMI and erosion models on whole-device scales in fusion-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fabrication of Piezoelectric Polymer and Metal–Organic Framework Composite Thin Films Using Solution Shearing

Polymer-metal–organic framework (polymer-MOF) composites have garnered significant interest as polymers can enhance the processability and industrial applicability of MOFs. Thin films of these composites are particularly attractive for applications in sensing, separations, and flexible electronics. Solution shearing, a meniscus-guided coating technique, has emerged as a scalable process for fabricating thin films of MOFs, and can produce large-area films within minutes. In this study, we utilized solution shearing to fabricate composite thin films of a MOF UiO-66 and a piezoelectric polymer poly(vinylidene fluoride-trifluoroethylene) (P(VDF-TrFE)), investigating how polymer concentration during MOF synthesis and composite formation influences thin film properties, including crystallinity, surface coverage, and piezoelectric performance. Additionally, solid-state NMR spectroscopy was utilized to probe the interactions between P(VDF-TrFE) and UiO-66 in the composite. Evidence from solid-state NMR indicated polymer-MOF interactions, suggesting that the polymer strands are in close proximity to the UiO-66 pores, supporting a mixed surface coating and pore infiltration model. Furthermore, incorporating P(VDF-TrFE) enhanced the film’s areal coverage from 70% to 100%. While the thermal conductivity remained essentially unchanged, the composite film showed an improved piezoelectric effect. The composite with 91 wt % P(VDF-TrFE) exhibited the highest output voltage of 9.1 V and a sensitivity of 0.26 V/N under applied pressure. This work demonstrates the potential of solution shearing as a scalable technique for fabricating polymer-MOF composite thin films.

P(VDF-TrFE)↗
DOE OSTI2025

Electron cyclotron emission measurements and modeling in optically semi-transparent plasma on the HSX stellarator

On the Helically Symmetric Experiment (HSX), electron cyclotron emission (ECE) at the second harmonic of the electron gyrofrequency is measured with a 16-channel heterodyne receiver. An ECE radiometer diagnostic has been calibrated with a broadband frequency source, allowing independent measurements of the electron temperature. The antenna pattern has a significant impact on the measurements, and modeling its effect is necessary to produce an ECE radiation temperature spectrum. Measurements, calibration results, and modeling of the ECE data are presented in this paper. At a relatively high plasma density, HSX plasmas are thermal, and the radiation temperature measured by the radiometer is in good agreement with the electron temperature measured by a Thomson scattering diagnostic. At a relatively low plasma density, the deviation from a Maxwellian distribution is large, and asymmetry in the radiation temperature profile is measured. These measurements are supported by single-pass and multi-pass ray tracing calculations of the electron cyclotron emission using the TRAVIS ray tracing code. Thermal and nonthermal electron populations exist in HSX, and the effect of these populations on the asymmetry in the measured emission is presented.

Cyclotron radiation↗