Search NASA⌕ Search

SEARCH · Search NASA

Results for “heat simulation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

High Performance Heat Pipe Power Transient Testing at SPHERE Facility

Microreactors are being researched, designed, and built at Idaho National Laboratory (INL). Microreactors are small reactors defined at less than 20MW of power. These reactor concepts are also being looked at throughout industry for various applications. An important aspect of these reactor designs is economic feasibility i.e. lower overnight capital cost. The driving factors for implementing microreactors are quick setup and takedown, minimal operators, and the ability to manufacture them readily and to fit in mid-sized containers for transport. A specific area of research to aid in successful integration of these factors within the designs is passive heat removal of the core’s thermal power. Interest in heat pipes to achieve this passive heat removal has been shown across multiple industry partners. Because of this interest, INL has developed a test facility to facilitate experimental tests for sodium filled heat pipes. INL has developed the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility to run experiments on high performance, sodium filled heat pipes. As mentioned above, heat pipes are passive heat transfer devices. Radially, heat pipes are broken up into an outer wall, a small annular gap, a wick structure, and a centerline gap. They function by utilizing latent heat transfer. Heat pipes are traditionally separated into three regions, an evaporator (heat input), an adiabatic region, and finally a condenser region (heat removal). As heat is being applied to the evaporator, the working fluid undergoes a phase change to a vapor. This phase change causes a differential pressure across the axial length of the pipe driving flow down the center gap of the heat pipe. The vapor flows down past the adiabatic region to the condenser where the heat is removed. This heat removal forces the working fluid to phase change back to a liquid. The wick structure is then utilized to drive the flow back towards the evaporator by capillary forces. This backflow is aided by the annular gap. Because this heat transfer mechanism functions with latent heat transfer, the heat pipe is close to isothermal down the axial length. Heat pipes can operate under a wide range of working fluids. Considerations for these working fluids are primarily driven by operating temperatures amongst other important factors based around overall performance. Sodium filled heat pipes operate from 450°C up to 900°C. This temperature range works well for the current microreactor designs. In conjunction with this experimental capability, INL has developed a modeling software to simulate heat pipe physics within reactor cores. This modeling software is called Sockeye and functions under the established INL Multiphysics Object Oriented Simulation Environment (MOOSE). SPHERE also supports Sockeye development by providing the modeling team with experimental data on an array of setups and operating parameters to support validation efforts. A power transient experiment was performed utilizing the SPHERE facility to continue to aid with Sockeye development. The testing followed a proposed test plan to ramp up and down the temperature of the heat pipe. Sockeye models steady state heat pipe operation with high accuracy, the data provided by the power transient testing aims to assist with the validation efforts and further enhance transient modeling capability of the tool [2].

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Heat Equation Neural Simulation

This code solves a steady state heat equation problem on a wire with random walks implemented using a spiking neural algorithm. SAND2020-12194 M

Reeder, Leah↗
DOE OSTI2020

Preliminary Reporting of Thermophysical Property Measurements for the Ghareb Formation

Accurate knowledge of thermophysical properties of rock is vital to develop meaningful models of high level nuclear waste emplacement scenarios. The Israel Atomic Energy Commission is considering storing high level nuclear waste in the Ghareb formation, a porous kerogen bearing chalk. Sandia is supporting this effort with an evolving lab- based geomechanics testing program. We have completed measurements of thermal properties up to 275C and room temperature hydrostatic compaction measurements. We report thermal conductivity, thermal diffusivity, specific heat, and mass loss from our thermal measurements, and we report bulk moduli and porosity loss from our compaction measurements. These values are crucial for the numerical models to simulate heat transfer and formation compressibility around a heat generating repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Machine learning surrogate of physics-based building-stock simulator for end-use load forecasting

Building energy models are used to simulate heat and mass transfer and estimate end-use load in buildings. With the proliferation of solar photovoltaics on residential and commercial buildings, increasingly, buildings are expected to provide grid services, for which accurate and computationally efficient building energy simulations and end-use load prediction are imperative. Existing building energy simulation tools, however, have significant computational overhead that make them less practical in real-time deployment for optimization, design, uncertainty quantification and control in building energy management systems. Here this article presents a data-driven machine learning model based on light gradient boosting method (LightGBM) as a surrogate for a physics-based simulator for residential buildings to predict end-use load. The machine learning based surrogate model accounts for time-series related variables, seasonality and trend component of end-use load, and history of end-use load. The accuracy of the surrogate model is assessed on the prediction of the load profiles of 100 different houses in Cook County, Illinois, USA. The LightGBM surrogate model is shown to reduce the root-mean-squared error by 53% relative to a reference decision tree (DT) based model reported previously in the literature. Moreover, the model predicts the load spikes and high-ramp rate events throughout the year which are often the Achilles heel of other models in the literature. The machine learning based surrogate model is demonstrated to be computationally efficient, with a ten-fold reduction in the computational time compared to a physics-based building energy simulation, and suitable for uncertainty analysis and real-time control of building characteristics in response to uncertainty.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

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↗
DOE OSTI2022

Gyrokinetic simulations compared with magnetic fluctuations diagnosed with a Faraday-effect radial interferometer-polarimeter in the DIII-D pedestal

Abstract Experimental data on electromagnetic fluctuations in DIII-D, made available by the Faraday-effect radial interferometer-polarimeter (RIP) diagnostic Chen (2016 Rev. Sci. Instrum. 87 11E108), is examined in comparison with detailed gyrokinetic simulations using gyrokinetic electromagnetic numerical experiment (GENE). The diagnostic has the unique capability of making internal measurements of fluctuating magnetic fields ∫ n e δ B r d R ∫ n e d R . Local linear simulations identify microtearing modes (MTMs) over a substantial range of toroidal mode numbers (peaking at n = 15) with frequencies in good agreement with the experimental data. Local nonlinear simulations reinforce this result by producing a magnetic frequency spectrum in good agreement with that diagnosed by RIP. Simulated heat fluxes are in the range of experimental expectations. However, magnetic fluctuation amplitudes are substantially lower than the experimental expectations. Possible sources of this discrepancy are discussed, notably the fact that the diagnostics are localized at the mid-plane—the poloidal location where the simulations predict the fluctuation amplitudes to be smallest. Despite some discrepancies, several connections between simulations and experiments, combined with general criteria discriminating between potential pedestal instabilities, strongly point to MTMs as the source of the observed magnetic fluctuations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

EXPERIMENTAL AND MODELLING ANALYSIS OF A LARGE-SCALE TWO-PHASE LOOP THERMOSYPHON

Liquid pumping requires large quantities of electrical energy, including about 7% of the energy of building heating, ventilation, and air conditioning (HVAC) systems. To replace pumped condenser-cooling tower water loops with a passive alternative system, we implemented a commercial-scale two-phase loop thermosyphon (TPLT). The unit consists of a 13 m riser integrated with a commercially available cooling tower and circulation water heater that simulates heat loads up to 25 kW. In addition to providing passive cooling capabilities, the cooling tower unit is also maintenance free, reliable, and can operate in both dry and wet modes. This study characterizes the performance (temperature difference between the evaporator and condenser) and the two-phase flow behavior of the loop under various refrigerant charges. Sight glasses installed throughout the loop are used to identify the operating flow regimes in the riser and downcomer. Over the range of operating conditions, we identified that there is an optimal refrigerant charge range for a specific heat load at which near-isothermal operation can be achieved. We further developed a model to predict the operating flow rate and gravitational height in the downcomer and compared it with the experimental data. The results show that the model agrees with the experimental data, in particular the threshold gravitational head height at which there will be subcooled liquid at the condenser exit, indicating that it can be used as a predictive tool for identifying the optimal loop charge for a given heat load.

Loop thermosyphon, passive two phase↗
DOE OSTI2026

A single-zone zero-dimensional study of HCCI combustion of methanol dehydration products to enable ignition of direct-injected methanol

Methanol is an alternative fuel gaining traction in the maritime sector. Its direct adoption, however, is accompanied by a unique set of technical challenges, such as low cetane number and high latent heat of vaporization. An approach to overcome these challenges is being developed at the US Department of Energy’s Oak Ridge National Laboratory, where onboard generation of dimethyl ether (DME) via catalytic dehydration of methanol can be used to assist in the mixing controlled combustion of direct-injected (DI) methanol. The generated mixture from this dehydration process can be premixed with intake air to condition the cylinder via. homogeneous charge compression ignition (HCCI) for subsequent DI methanol. In this preliminary work, various catalyst or reactor conversion efficiencies were simulated (using bottles) at constant DME and water flow at low load on a single-cylinder marine-variant of a CAT® C18 18 L engine with a 145 mm bore. To substantiate the experimental findings, a zero-dimensional engine model was developed in Cantera using a DME mechanism with 79 species and 658 reactions. Results presented include experimental and simulation heat release rate comparisons, species evolution information, and constant volume ignition delay (ID) for DI methanol with and without background species from HCCI of the premixed products from different reactor efficiencies. The results suggest that thermal effects dominate the DI methanol ignition process, and this work provides a chemical kinetic foundation or guideline for developing future control schemes.

Tyrewala, Daanish [ORNL] (ORCID:0000000208599324)↗
DOE OSTI2026

Overview of the physics basis for the ARC fusion power plant

Commonwealth Fusion Systems plans to build ARC as the first fusion power plant at a site in Chesterfield County, Virginia, USA by the early 2030s. We present an overview of analysis comprising the physics basis of the ARC V3A design, a high-magnetic-field tokamak with 𝐵 0 =11.4 T, 𝐼 𝑝 =12.0 MA, 𝑅 0 =4.62 m, 𝑎 =1.18 m. ARC V3A is designed to produce 𝑃 𝑓⁡𝑢⁢𝑠 ≈1.13 GW DT fusion power and deliver ⩾400 MW net electric power to the grid. This overview includes quantitative analysis of fundamental issues for design of and operational plasma scenarios for a tokamak power plant, and lays out the design targets and strategic choices for ARC, including empirical fusion performance projections, assessment of H-mode access, ion cyclotron resonance heating simulations, alpha particle physics and time-dependent full-pulse simulations. This is complemented by topical papers on fusion performance and transport, disruption physics, boundary physics and magnetohydrodynamic stability. Critically, these studies identify key model uncertainties and physics risks to be retired through SPARC operation. Due to the modular nature of ARC, early results from SPARC can be incorporated into the design of the first ARC as well as subsequent replacements of the ARC vacuum vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Interfacial water asymmetry at ideal electrochemical interfaces

Controlling electrochemical reactivity requires a detailed understanding of the charging behavior and thermodynamics of the electrochemical interface. Experiments can independently probe the overall charge response of the electrochemical double layer by capacitance measurements and the thermodynamics of the inner layer with potential of maximum entropy measurements. Relating these properties by computational modeling of the electrochemical interface has so far been challenging due to the low accuracy of classical molecular dynamics (MD) for capacitance and the limited time and length scales of ab initio MD. Here, we combine large ensembles of long-time-scale classical MD simulations with charge response from electronic density functional theory to predict the potential-dependent capacitance of a family of ideal aqueous electrochemical interfaces with different peak capacitances. We show that while the potential of maximum capacitance varies, this entire family exhibits an electrode charge of maximum capacitance (CMC) between -2.9 and -2.2 μC/cm 2 , regardless of the details in the electronic response. Simulated heating of the same interfaces reveals that the entropy peaks at a charge of maximum entropy (CME) of -5.1 ± 0.6 μC/cm 2 , in agreement with experimental findings for metallic electrodes. The CME and CMC both indicate asymmetric response of interfacial water that is stronger for negatively charged electrodes, while the difference between CME and CMC illustrates the richness in behavior of even the ideal electrochemical interface.

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

Investigation of core impurity transport in DIII-D diverted negative triangularity plasmas

Abstract Tokamak operation at negative triangularity has been shown to offer high energy confinement without the typical disadvantages of edge pedestals (Marinoni et al 2021 Nucl. Fusion 61 116010). In this paper, we examine impurity transport in DIII-D diverted negative triangularity experiments. Analysis of charge exchange recombination spectroscopy reveals flat or hollow carbon density profiles in the core, and impurity confinement times consistently shorter than energy confinement times. Bayesian inferences of impurity transport coefficients based on laser blow-off injections and forward modeling via the Aurora package (Sciortino et al 2021 Plasma Phys. Control. Fusion 63 112001) show core cross-field diffusion to be higher in L-mode than in H-mode. Impurity profile shapes remain flat or hollow in all cases. Inferred radial profiles of diffusion and convection are compared to neoclassical, quasilinear gyrofluid, and nonlinear gyrokinetic simulations. Heat transport is observed to be better captured by reduced turbulence models with respect to particle transport. State-of-the-art gyrokinetic modeling compares favorably with measurements across multiple transport channels. Overall, these results suggest that diverted negative triangularity discharges may offer a path to a highly-radiative L-mode scenario with high core performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Infrared thermography NDT for in-situ defect detection in sandwich composite panel manufacturing

Composite manufacturing presents numerous challenges, as defects can arise from various sources throughout the process. In sandwich composite structures, the integration of a foam core introduces additional complexity and increases the likelihood of defect formation like delamination. To mitigate these issues and reduce the risk of future structural failures, in-situ monitoring during manufacturing is essential. This study investigates infrared (IR) thermography as a non-destructive technique for detecting manufacturing defects in foam-core sandwich composite panels under thermally excited conditions representative of in-situ processing. A stationary FLIR A8590 IR camera (640 × 512 pixels, 30Hz, 17mm lens, 9 ft stand-off distance) was used to monitor prefabricated panels subjected to controlled external heating simulating compression molding and resin cure exotherm. Interlaminar delamination defects with characteristic sizes ranging from 0.25 × 0.25in² to 5 × 5in² produced measurable surface temperature depressions of approximately 4–10°C during transient cooling, exceeding the effective noise floor of the camera by more than two standard deviations. Thicker laminates exhibited prolonged defect detectability windows due to increased thermal diffusion time. In contrast, embedded Teflon inclusions generated weak thermal contrasts of ≤ 3°C, approaching the measurement noise floor, due to limited thermal property contrast with the surrounding glass fiber composite. These results establish quantitative detectability limits for stationary thermographic inspection of sandwich composite panels under manufacturing-representative thermal cycles.

Barakat, Abdallah [ORNL] (ORCID:0000000296141398)↗
DOE OSTI2022

PDC Modifications for Analysis of Gas-Cooled Reactors with Direct Helium Brayton Cycle

Capabilities of the Plant Dynamics Code (PDC) have been extended to allow steady-state and transient simulation of graphite-cooled reactors with direct helium Brayton cycle. On the cycle side, the most significant code modification is the addition of helium properties, in the format required by the code’s equations. Since the code was already formulated to handle more complex real gas properties, adding helium as working fluid that behaves like ideal gas was fairly straightforward. A reactor module was added to PDC to simulate a reactor cooled by the working fluid of the Brayton cycle. Two options are supported: channel type, typical for graphite gas-cooled reactors, and pin type, typical for light-water and liquid metal-cooled reactors. The reactor module is an extension of the electrical heater model and simulates heat deposition in the fuel and transfer of this heat from the fuel to the coolant through the matrix and tube materials. The new reactor module becomes the third option in PDC for modeling heat addition to the cycle, besides previously modeled heat addition heat exchanger and electrical heater. In addition to those changes, other minor code modifications and improvements were introduced during the work of expanding PDC to modeling of gas-cooled reactors. These modifications are summarized in the last chapter of this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Design Heating and Cooling Load Calculation Versus Building Load Simulation for Cold Climate Heat Pumps: Understanding the "Gap"

This report explores the differences between Manual J-equivalent block load calculations and building HVAC energy simulation results using EnergyPlus™ calculations when designing cold climate heat pump systems for residential use. This study will help HVAC researchers and advanced designers understand the impacts of oversizing heat pumps on home energy use.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Simulations of ion heating due to ion-acoustic instabilities in presheaths

Particle-in-cell, direct simulation Monte Carlo simulations reveal that ion-acoustic instabilities excited in presheaths can cause significant ion heating. Ion-acoustic instabilities are excited by the ion flow toward a sheath when the neutral gas pressure is small enough and the electron temperature is large enough. A series of 1D simulations were conducted in which neutral plasma (electrons and ions) was uniformly sourced with an ion temperature of 0.026 eV and different electron temperatures (0.1 eV–50 eV). Ion heating was observed when the electron-to-ion temperature ratio exceeded the minimum value predicted by linear response theory to excite ion-acoustic instabilities at the sheath edge (T e /T i ≈ 28). When this threshold was exceeded, the temperature equilibration rate between ions and electrons rapidly increased near the sheath so that the local temperature ratio did not significantly exceed the threshold for instability. This resulted in significant ion heating near the sheath edge, which also extended back into the bulk plasma; presumably due to wave reflection from the sheath. Furthermore, this ion-acoustic wave heating mechanism was found to decrease for higher neutral pressures, where ion-neutral collisions damp the ion-acoustic waves and ion heating is instead dominated by inelastic collisions in the presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Experiments and simulations of isochorically heated warm dense carbon foam at the Texas Petawatt Laser

An experimental and simulation study of warm dense carbon foams at ambient density (n e ~10 21 cm –3 ) is presented. This study of isochorically heated foams is motivated by their potential application in carbon-atmosphere white-dwarf envelopes, where there are modeling uncertainties due to the equation of state. The foams are heated on an approximately picosecond time scale with a laser-accelerated proton beam. The cooling and expansion of the heated foams can be modeled with appropriately initialized radiation-hydrodynamics codes; xRAGE code is used in this work. The primary experimental diagnostic is the streaked optical pyrometer, which images a narrow band of radiation from the rear surface of the heated material. Presented are xRAGE modeling results for both solid aluminum targets and carbonized resorcinol-formaldehyde foam targets, showing that the foam appears to cool slowly on the pyrometer because of partial transparency. So that simulations of cooling foam are processed properly, it is necessary to account for finite optical depth in the photosphere calculation, and the methods for performing that calculation are presented in depth.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Studies of the influence of heat pipe model and coupling strategy in heat pipe microreactor simulations using Sockeye

This work demonstrates various modeling approaches for thermally coupling heat pipes to a graphite monolithic block in a microreactor context, accounting for the presence of a gap. We show that heat pipe model selection has a significant impact on the transient behavior of a three-dimensional coupled microreactor assembly problem by comparing results obtained using different heat pipe models from the heat pipe code Sockeye. Additionally, we show that the choice of coupling strategy has significant consequences for both accuracy and performance. Furthermore, we found that our heat flux transfer strategy exhibited greater accuracy than our temperature transfer strategy, even when comparing a loose coupling of the heat flux transfer strategy to a tight coupling of the temperature transfer strategy. Although this research uses only theoretical test problems, it provides important insights in modeling thermal fluids phenomena in heat pipe microreactors.

Direwolf↗
DOE OSTI2023

Multi-scenario Extreme Weather Simulator application to heat waves: Ko’olauloa community resilience hub

Heat waves are increasing in severity, duration, and frequency. The Multi-Scenario Extreme Weather Simulator (MEWS) models this using historical data, climate model outputs, and heat wave multipliers. In this study, MEWS is applied for planning of a community resilience hub in Hau'ula, Hawaii. The hub will have normal operations and resilience operations modes. Both these modes were modeled using EnergyPlus. The resilience operations mode includes cutting off air conditioning for many spaces to decrease power requirements during emergencies. Results were simulated for 300 future weather files generated by MEWS for 2020, 2040, 2060, and 2080. Shared socioeconomic pathways 2-4.5, 3-7.0 and 5-8.5 were used. The resilience operations mode results show two to six times increase of hours of exceedance beyond 32.2 °C from present conditions, depending on climate scenario and future year. The resulting decrease in thermal resilience enables an average decrease of energy use intensity of 26% with little sensitivity to climate change. The decreased thermal resilience predicted in the future is undesirable, but was not severe enough to require a more energy-intensive resilience mode. Instead, planning is needed to assure vulnerable individuals are given prioritized access to air-conditioned parts of the hub if worst-case heat waves occur.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗