Engineering PapersSearch

SEARCH · Engineering Papers

Results for “thermal”

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 775 records · Page 43

MOOSE Thermal-Hydraulics Module - MOOSE workshop

The MOOSE Thermal Hydraulics Module (THM) is designed to facilitate the development of thermal hydraulic system models. It provides the capability to assemble networks of coupled components such as pipes, junctions, valves, turbomachinery, and heat exchangers. Its library of components supports a single-phase, compressible flow model based on a variable-area formulation of the Euler equations of gas dynamics and discretized using a finite volume scheme. THM offers a flexible system for specifying closures such as friction factors or heat transfer coefficients, allowing the user to choose from built-in correlations or define their own in the input file. A control logic system can be used to control input parameters, necessary for implementing transient scenarios and mirroring real control systems in thermal hydraulic systems. THM can be coupled with other MOOSE-based applications for multiphysics calculations. This training will give an introduction to the capabilities of THM and provide some examples of its usage and validation.

97 - MATHEMATICS AND COMPUTING

Latent heat thermal energy storage performance maps enabling fast & accurate building energy simulations

Thermal energy storage (TES) using phase change materials (PCMs) has gained attention as an effective approach to manage energy demand fluctuations and shift peak building loads. PCM embedded heat exchangers (PCM-HXs) offer high energy storage density and low temperature variation during phase change, being suitable for load-shifting applications. However, this component is typically evaluated using computationally expensive methods, which present significant challenges when the ultimate goal is to assess the performance of PCM-HX integrated thermal energy storage systems in the full building context. In this paper, we present a methodology to generate highly accurate and computationally efficient PCM-HX performance maps which can be easily integrated into building energy simulation tools to analyze the feasibility of space conditioning systems with latent heat PCM-based TES. The performance maps are generated using a computationally efficient PCM-HX simulation tool based on a Generalized Resistance-Capacitance Model (GRCM) which can simulate arbitrary PCM-HXs with high accuracy and significantly less computational effort compared to full CFD simulations. The methodology was verified for a case study considering a 5-ton (~17.5 kW) air-to-water heat pump-thermal energy storage system (HP-TES), which was co-simulated in Modelica for a DOE prototype small-office building in Vienna, Austria, using Spawn of EnergyPlus™. The TES performance maps provided accurate predictions of PCM-HX behavior when used as Modelica component, with deviations within 2-4% while also achieving at least 103 computational time reduction. Leveraging this faster prediction capability, four PCMs with different melting temperatures for cooling (12°C, 16°C) and heating (31°C, 36°C) were assessed to investigate their impact on system performance. This work highlights the importance of robust PCM-HX models for efficient and high-fidelity building-level simulations, presenting new opportunities for advanced control strategy development and parametric analysis of TES configurations in a computationally efficient manner

Modelica Building Simulations

A Physics-Informed Reinforcement Learning Framework for Economic-Thermal Co-Optimization of Crypto Mining Data Centers: Preprint

The rapid expansion of cryptocurrency mining has created a new class of high-density data centers characterized by extreme thermal flux and high sensitivity to volatile economic markets. Traditional thermal management strategies, typically reliant on rule-based control, maintain static setpoints that fail to account for fluctuating electricity prices and cryptocurrency values - factors critical to mining profitability. To address this, we present a physics-informed reinforcement learning (PIRL) framework for economic-thermal co-optimization in crypto mining data centers. This framework consists of a proximal policy optimization (PPO) agent, a virtual testbed powered by high-fidelity physics-based models, and an interactive frontend dashboard. The PPO agent is trained using the virtual testbed and strict hardware safety limits. This physics-informed approach allows the agent to learn a stochastic policy that dynamically balances mining revenue against operational costs by co-optimizing HVAC cooling setpoints and IT computational hashrate. The simulation results demonstrate that the integrated framework achieved an 8.62% increase in net operational profit compared to traditional baseline strategies while strictly adhering to safety-critical temperature constraints (coolant supply temperature < 32 degrees C). This work provides a scalable template for the deployment of reinforcement learning in mission critical facilities where economic volatility and physical safety must be managed simultaneously.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI

Scaling behavior and giant field enhancement of the thermal conductivity in the honeycomb antiferromagnet BaCo 2⁢ (AsO 4 ) 2

The layered honeycomb material BaCo 2 (AsO 4 ) 2 is of topical interest because its magnetic state is related to that of the Kitaev magnet α-RuCl 3 . Using thermal transport to probe how magnetic excitations interact with phonons in the magnetically disordered regime, we have uncovered an unusually large enhancement of the thermal conductivity κ xx in an in-plane magnetic field H. Just above the Néel temperature T N , a field of 13 T increases κ xx by a factor of ∼ 211, which is very large compared to other magnetic insulators. Interestingly, κ xx (H, T) exhibits a scaling behavior in the entire magnetically disordered region that surrounds the ordered zigzag state. The ratio Δκ xx (H, T)/κ xx (13, T), measured throughout the disordered region, collapses to a one-parameter scaling function exp(−1/gx) (where x = μ B B/k B T and g is a constant).

36 MATERIALS SCIENCE

Early Progress Toward the Feasibility of the Centrifugal Nuclear Thermal Rocket

The Centrifugal Nuclear Thermal Rocket (CNTR) is a Nuclear Thermal Propulsion (NTP) concept designed to heat propellant directly by the reactor fuel. The primary difference between the CNTR concept and traditional NTP systems is that rather than using traditional solid fuel elements, the CNTR uses liquid fuel with the liquid contained in rotating cylinders by centrifugal force. If the concept can be successfully realized, the CNTR would have a high specific impulse (~1800 seconds) at high thrust, which may enable (i) viable near-term human Mars exploration by reducing round-trip times to 420 days and (ii) direct injection orbits for scientific missions to the Solar System outer planets and potentially Kuiper Belt objects. The CNTR could also use storable propellants such as ammonia, methane, or hydrazine at an Isp of ~900 seconds, enabling long-term in-space storage of a dormant system. Significant engineering challenges must be addressed to establish the technical viability of the CNTR. Research is presently underway to determine resolutions for these engineering challenges. In particular, research has begun on the analytical modeling and simulation of the two-phase heat transfer between the liquid metallic uranium fuel and the gaseous propellant. A paper was presented at the 2021 IAC which described these challenges and the study plan to address them. This paper describes the analytical and experimental progress to date toward resolving these challenges and establishing the engineering feasibility of the CNTR technology.

Nuclear Thermal Propulsion

The role of thermal stratification on the co‐spectral properties of momentum transport above an Amazonian forest

The influence of thermal stratification on the turbulent kinetic energy balance has been widely studied; however, its influence on the turbulent stress remains less explored in the presence of tall vegetated canopies and less ideal micrometeorological conditions. Here, the impact of thermal stratification on turbulent momentum flux is considered in the roughness sublayer (RSL) and the atmospheric surface layer (ASL) using the Amazon Tall Tower Observatory (ATTO) in Brazil. A scalewise co‐spectral budget (CSB) model is developed using standard closure schemes for the pressure–velocity decorrelation. The CSB revealed that the co‐spectrum ${F}_{wu}\left({k}_x\right)$ between longitudinal (u') and vertical (w') velocity fluctuations is impacted by the energy spectrum of the vertical velocity ${E}_{ww}\left({k}_x\right)$ and the much less studied longitudinal heat‐flux co‐spectrum ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$, where ${\theta}_{\mathrm{v}}^{\prime }$ are temperature fluctuations and ${k}_x$ is the longitudinal wavenumber. Under stable, very stable, and dynamic–convective conditions, the scaling exponent ${F}_{wu}\left({k}_x\right)$ in for the inertial subrange (ISR) scales is dominated by ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$ instead of ${E}_{ww}\left({k}_x\right)$. A near ${k}_x^{-7/3}$scaling in ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$ robust to large variations in thermal stratification is found, whereas the Kolmogorov ISR scaling for ${E}_{ww}\left({k}_x\right)\sim {k}_x^{-5/3}$ is not found. The scale‐dependent decorrelation time between u' and w' is dominated by ${\epsilon}^{-1/3}{k}_x^{-2/3}$ in the ISR, but is nearly constant for eddies larger than the vertical velocity integral scale, regardless of stability. Implications of these findings for generalized stability correction functions that are based on the turbulent stress budget instead of the turbulent kinetic energy budget are discussed.

canopy turbulence

Separation of infrared and bulk in thermal QCD

A new thermal regime of QCD, featuring decoupled scale-invariant infrared glue, has been proposed to exist both in pure-glue (N f =0) and “real-world” (N f =2+1 at physical quark masses) QCD. In this IR phase, elementary degrees of freedom flood the infrared, forming a distinct component independent from the bulk. This behavior necessitates non-analyticities in the theory. In pure-glue QCD, such non-analyticities have been shown to arise via Anderson-like mobility edges in Dirac spectra (λ IR = 0, ±λ A ≠ 0), as manifested in the dimension function d IR (λ). Here, we present the first evidence, based on lattice QCD calculation at a=0.105 fm, that this mechanism is also at work in real-world QCD, thus supporting the existence of the proposed IR regime in nature. An important aspect of our results is that, while at T = 234 MeV we find a dimensional jump between zero modes and lowest near-zero modes very close to unity (d IR = 3 to d IR ≃ 2), similar to the IR phase of pure-glue QCD, at T = 187 MeV we observe a continuous λ-dependence. This suggests that thermal states just above the chiral crossover are non-analytically (in T) connected to thermal state at T = 234 MeV, supporting the key original proposition that the transition into the IR regime occurs at a temperature strictly above the chiral crossover.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Mechanical and Thermal Characterization of Additively Manufactured Carbon/Nylon 12 and Carbon/PEEK Composites

This study explores additive manufacturing of carbon fiber-reinforced thermoplastic composites using the Composite-Based Additive Manufacturing (CBAM) process. Carbon/Nylon 12 and Carbon/PEEK composites were fabricated and evaluated through mechanical (compression, tensile, flexural, and impact) and thermal (DSC and TGA) tests. Carbon/PEEK exhibited superior mechanical performance, with 97.5% higher tensile strength, 79.8% higher elastic modulus, and 59.6% higher flexural strength compared to Carbon/Nylon 12. Thermal testing showed that Carbon/PEEK had higher thermal stability, beginning degradation at 350 °C versus 298 °C for Carbon/Nylon. These results indicate that CBAM-fabricated Carbon/PEEK composites are suitable for applications requiring high strength and temperature resistance.

Additive manufacturing

Characterization of SiO 2 Thermally Grown Oxide Stress Evolution of EBCs with Al-Containing Dopants

SiC/SiC ceramic matrix composites (CMCs) are desired for use in combustion environments to achieve higher turbine operating temperatures. However, CMCs require environmental barrier coatings (EBCs) for protection from the gas environment. EBC systems are known to primarily fail through coating delamination via growth of a thermally grown oxide (TGO) at the EBC—silicon bond coating interface when exposed to steam, which accelerates the TGO growth rate. The TGO undergoes a phase transformation during thermal cycling, which results in stresses that may encourage EBC spallation. Yb-silicate EBCs with mullite and yttrium aluminum garnet (YAG) dopant additions were deposited on SiC substrates with a Si intermediate bond coating and exposed to thermal cycling in steam at 1350 °C. The impact of Al dopant additions on the TGO growth rate and the SiO 2 phase transformation was assessed. Photo-stimulated luminescence spectroscopy (PSLS) was used to characterize the Al-containing phases and to measure stress evolution in the EBC following exposure using the stress-induced peak shift of the R-lines of mullite. Raman microscopy was used to map the stresses in the Si bond coating following exposure. It was found that the TGO phase transformation upon cooling increased compressive stress in the Si bond coating within 15 µm of the TGO.

Building Materials

Vector-level feedforward control of LPBF melt pool area using a physics-based thermal model

Laser powder bed fusion (LPBF) is an additive manufacturing technique that has gained popularity thanks to its ability to produce geometrically complex, fully dense metal parts. However, these parts are prone to internal defects and geometric inaccuracies, stemming in part from variations in the melt pool. Here, this paper proposes a novel vector-level feedforward control framework for regulating melt pool area in LPBF. By decoupling part-scale thermal behavior from small-scale melt pool physics, the controller provides a scale-agnostic prediction of melt pool area and efficient optimization over it. This is done by operating on two coupled lightweight models: a finite-difference thermal model that efficiently captures vector-level temperature fields and a reduced-order, analytical melt pool model. Each model is calibrated separately with minimal single-track and 2D experiments, and the framework is validated on a complex 3D geometry in both Inconel 718 and 316L stainless steel. Results showed that feedforward vector-level laser power scheduling reduced geometric inaccuracy in key dimensions by 62%, overall porosity by 16.5%, and photodiode root-mean-squared deviation by 38.5% on average. Overall, this modular, data-efficient approach demonstrates that proactively compensating for known thermal effects can significantly improve part quality while remaining computationally efficient and readily extensible to other materials and machines.

Additive manufacturing

MURR LEU structural and thermal hydraulics analyses: Part I – Preliminary irradiation thermo-mechanical behavior

The University of Missouri Research Reactor (MURR) is expected to be converted from highly enriched uranium (HEU, ≥ 20 wt% U-235) U-Al x dispersion fuel to low-enriched uranium (LEU, < 20 wt% U-235) with U-10Mo monolithic fuel. This work introduces high-fidelity irradiation thermo–mechanical (T-M) analysis of the MURR LEU focusing on changes in coolant channel gap thickness. Three-dimensional (3D) finite element (FE) models were developed to simulate the irradiation T-M behavior of the MURR LEU element with all 23 curved fuel plates, the two side plates, and the combs. It was shown that channel gap thickness changes were influenced not only by plate thickness variations due to fuel swelling and creep but also by the radial displacement of consecutive MURR LEU plates. Modeling the fuel element assembly captured side plate displacements, which were shown to reduce radial fuel plate displacements towards the convex side. The maximum local radial displacement in the element was predicted at the end of life (EOL) as 23.7 mil (602.0 µm) on the lateral centerline of plate 23 towards the convex side. The maximum stripe-averaged reduction in channel gap thickness, particularly relevant for thermal hydraulics (TH) safety analysis, was calculated as 15.9 mil (403.9 µm) in single-side heated channel 24 (the outermost channel). These results account for the thermal resistance from the oxide build-up on cladding surfaces which was shown to be up to 0.82 mil (20.8 µm) thick. It was demonstrated that accounting for oxide layer thermal resistance led to a 10 °C higher peak fuel temperature and a 4.4 mil (111.8 µm) greater maximum local radial displacement. The impact of the calculated channel gap thickness changes on the MURR LEU TH safety analysis is evaluated in Part II.

Irradiation thermo-mechanical analysis

Deep reinforcement learning control for co-optimizing energy consumption, thermal comfort, and indoor air quality in an office building

With the recent demand for decarbonization and energy efficiency, advanced HVAC control using Deep Reinforcement Learning (DRL) becomes a promising solution. Due to its flexible structures, DRL has been successful in energy reduction for many HVAC systems. However, only a few researches applied DRL agents to manage the entire central HVAC system and control multiple components in both the water loop and the air loop, owing to its complex system structures. Moreover, those researches have not extended their applications by incorporating the indoor air quality, especially both CO2 and PM2.5concentrations, on top of energy saving and thermal comfort, as achieving those objectives simultaneously can cause multiple control conflicts. What's more, DRL agents are usually trained on the simulation environment before deployment, so another challenge is to develop an accurate but relatively simple simulator. Therefore, we propose a DRL algorithm for a central HVAC system to co-optimize energy consumption, thermal comfort, indoor CO2 level, and indoor PM2.5 level in an office building. To train the controller, we also developed a hybrid simulator that decoupled the complex system into multiple simulation models, which are calibrated separately using laboratory test data. The hybrid simulator combined the dynamics of the HVAC system, the building envelope, as well as moisture, CO2, and particulate matter transfer. Three control algorithms (rule-based, MPC, and DRL) are developed, and their performances are evaluated on the hybrid simulator environment with a realistic scenario (i.e., with stochastic noises). The test results showed that, the DRL controller can save 21.4 % of energy compared to a rule-based controller, and has improved thermal comfort, reduced indoor CO2 concentration. The MPC controller showed an 18.6 % energy saving compared to the DRL controller, mainly due to savings from comfort and indoor air quality boundary violations caused by unmeasured disturbances, and it also highlights computational challenges in real-time control due to non-linear optimization. Finally, we provide the practical considerations for designing and implementing the DRL and MPC controllers based on their respective pros and cons.

Guo, Fangzhou

A scalable framework for efficient coupling of thermal and microstructural simulations in additive manufacturing

Predicting microstructure evolution in metal additive manufacturing (AM) is important for process optimization, but spatiotemporal scale disparities between thermal transport and microstructure evolution create significant challenges for efficient data transfer between simulation codes. To address this, we present Stork, a scalable framework for coupling thermal and microstructural simulations. Stork uses a sparse data representation to identify and store active solidification sub-volumes, enabling highly parallel quad-linear interpolation from coarse thermal grids to fine microstructure grids without large intermediate storage. We demonstrate the framework by coupling the semi-analytic heat transfer code 3DThesis with the time-parallel cellular automata code Toucan. This approach achieves over two orders of magnitude reduction in data generation time and file size compared to prior workflows. Numerical studies show that quad-linear interpolation preserves grain morphology and crystallographic texture in laser powder bed fusion (LPBF) simulations for coarsening ratios up to 16. Overall, Stork provides a scalable pathway for high-throughput, component-scale AM simulations on modern high-performance computing systems.

36 MATERIALS SCIENCE

A dynamic 2D Borehole Thermal Energy Storage (BTES) model for enhanced computational efficiency

Progressing toward a future increasingly reliant on renewable energy sources, the development of effective, durable energy storage solutions becomes essential to balance supply and demand fluctuations. Borehole Thermal Energy Storage (BTES) is a long-duration thermal energy storage technology that captures excess heat generated from renewable energy sources and stores it underground for later use, enabling the efficient utilization of sustainable energy. This approach is particularly valuable in district energy networks when integrated with Ground Source Heat Pumps (GSHP) to provide stable heating and cooling. However, traditional three-dimensional (3D) numerical models of BTES systems demand extensive computational resources, limiting their practicality for real-time and large-scale applications. This study introduces a novel two-dimensional (2D) modeling approach that reduces computational costs while maintaining high accuracy. By employing a radial ring-based discretization method, the model simulates heat injection, retention, and retrieval dynamics over seasonal cycles. A new thermal-mass weighted-average temperature parameter is introduced to evaluate the performance of BTES systems. Model validation against FEFLOW simulations demonstrates a 17-fold improvement in computational speed compared to traditional Computational Fluid Dynamics (CFD) models while achieving a mean absolute percentage error (MAPE) of 2 % during charging and 4 % during discharging. Additionally, a trade-off analysis between computational efficiency and accuracy is conducted, ensuring the model's applicability for real-world scenarios. The findings of this research contribute to the development of computationally efficient BTES models, facilitating better optimization, control, and integration into renewable energy systems. This work provides a foundation for further studies in techno-economic analysis, multi-year performance evaluation, and real-time operational strategies for BTES applications, supporting a more sustainable energy future.

2D modeling

Thermal and structural analysis of a new WEST reciprocating probe head with Titanium-Zirconium-Molybdenum armour material

The Tore Supra tokamak was transformed into an X-point divertor fusion device in the frame of the WEST (W-for tungsten-Environment in Steady-state Tokamak) project, launched in support of the ITER tungsten divertor strategy. WEST began operation in 2017. A new reciprocating probe head, equipped with a Titanium-Zirconium-Molybdenum (TZM) heat shield, was recently built and operated on WEST during the C9 Campaign in early 2024. TZM was chosen as armour material for its very good mechanical properties between 700 °C and 1400 °C. This material was already extensively used for probe heads and other plasma facing components on the Alcator C-Mod tokamak. Here, this paper presents the design of the new WEST probe head. It is equipped with tungsten filaments that can be heated to emit electrons and provide a direct measurement of the plasma potential in the scrape-off layer, and it is also equipped with a pair of Langmuir probes. The probe plunges vertically into the plasma, reaches its maximum dive depth, then returns to its protected resting position in about 0.25 s. Thermal loads are calculated using as input the thermal properties of TZM and radial profiles of the heat flux measured in the past by other reciprocating probes, and scaled upwards to the expected heat loads in the WEST scrape-off layer. An analytical estimation of Electromagnetic (EM) loads will also be presented, as well as a thermal and structural analysis of the probe for 15 cycles at maximum plunge depth. The probe is observed by an infrared camera. Comparisons between the heat flux deduced by the Langmuir probe and the infrared radiance of the probe housing during plasma exposure are finally presented.

Diagnostic

Comparative characterization of mixed spectra and thermal neutron shielded irradiated tungsten

The effect of mixed spectra and thermal neutron shielded irradiation on tungsten was evaluated with plasma exposure in the tritium plasma experiment followed by thermal desorption spectroscopy, X-ray photoelectron spectroscopy, and transmission electron microscopy. The two different irradiation campaigns were performed at the High Flux Isotope Reactor to 0.39–0.74 displacement per atom (dpa) in the 894–1379 K temperature range. A neutron spectrum influence on the void size and void number density was not observed. However, a strong correlation was found between void size and void number density with temperature, but not with dpa in the limited dpa range of this study. Thermal neutron shielding significantly reduced the transmutation to Re+Os. Higher irradiation temperature will lead to larger voids with lower number density, which reduces deuterium retention. In conclusion, grain growth was also observed for high-temperature irradiation of over ~1300 K within the limited grains visible in the transmission electron microscopy specimens.

36 MATERIALS SCIENCE

Defect formation and transmutation behaviors in irradiated tungsten under thermal-neutron shielded and mixed spectrum conditions

Here, this study investigates the microstructural evolution of pure tungsten irradiated under thermal-neutron shielded and mixed spectrum conditions in the High Flux Isotope Reactor (HFIR). Four samples were irradiated at temperatures from 570 °C to 1130 °C up to 0.73 dpa. Neutron spectrum significantly influenced the accumulation of transmutation products, with Re+Os content estimated at ∼0.3–0.6% under thermal-neutron shielded conditions and ∼5.2% under the mixed spectrum condition. Irradiation temperature strongly influences tungsten’s microstructure, with dislocation loops and fine voids forming at lower temperatures and only larger voids and Re/Os segregation observed at higher temperature. Under thermal-neutron shielded conditions, dislocation loops and voids were observed at 570 °C and 790 °C. At the highest irradiation temperature (1130 °C), dislocation loops were no longer observed, while larger but less dense voids remained. Re and Os segregation to void surfaces was evident at 790 °C and 1130 °C, though no precipitation was observed. In contrast, under the mixed-spectrum condition, both spherical and needle-like Re/Os-rich precipitates were observed, frequently accompanied by large voids. Dislocation loops were not observed, but loop-like contrast within the precipitates suggests they may have nucleated on pre-existing loops. Irradiation-induced hardening was assessed for the shielded samples at 570 °C and 790 °C. Dispersed barrier hardening (DBH) analysis, based on TEM-resolved defects, revealed that voids were the dominant contributors to hardening, consistent with literature results. A schematic model is proposed to describe defect and precipitate evolution in tungsten under fusion-relevant transmutation-to-dpa conditions.

Dislocation loops

Microstructural stability of irradiated yttrium hydride under thermal cycles

Hydrides that retain their hydrogen in service are critical for advanced microreactor and space nuclear systems. Hydrogen retention is affected by evolving microstructure under temperature and irradiation extremes. This study investigated microstructural features in neutron-irradiated yttrium hydride and the stability of these features under thermal cycles. Differential scanning calorimetry was used to determine hydrogen desorption and temperatures at which significant phase changes are occurring, and transmission electron microscopy was used for microstructural analysis after each thermal cycle. Cyclic heating led to crystallization and epitaxial growth of surface oxidation and dehydriding of the matrix. Other features such as the bulk matrix crystal structure, irradiation-induced cavities, and matrix precipitation remained constant after thermal cycling. Results allow us to design better hydride alloys through microstructure tailoring and potential irradiation conditioning treatments.

36 MATERIALS SCIENCE