Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal 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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Hydraulic and Thermal Stimulation Program at Raft River Idaho, A DOE EGS

A Department of Energy Enhanced Geothermal System (EGS) stimulation program has injected over 254 million gallons of water into the well RRG-9 ST1 since the summer of 2013. Three major stimulations have been conducted during the program increasing injection flow rates from less than 20 gpm to 550 gpm. Geologic, water chemistry, microseismic activity, and borehole imaging data have been used to develop a conceptual model describing possible flow paths of this injected water. This model contains two major fracture zones one of which intersects the RRG-9 ST1 wellbore. Modified Hall and injectivity index plots constructed using injection flow rates, surface temperatures, and wellhead pressures show steady improvement in the injectivity of the well. Here, the injectivity index has risen from 0.15 gpm/psi to 2.0 gpm/ psi. A pressure falloff test conducted on April 28, 2015 indicates a reservoir permeability of 1,220 md and -5.38 skin factor. The well stimulation program was simulated numerically using an Idaho National Laboratory reservoir simulation code, FALCON. These simulations show a significant increase in the permeability of connecting fracture pathways after each stimulation event.

Enhanced geothermal system↗
DOE OSTI2025

Improving Thermal Management Strategies for Data Centers: A Physical Testbed Incorporating Small Modular Reactor and Microreactor Technology

This study aims to accelerate the demonstration of various thermal management systems for data centers using nuclear-generated heat to enhance energy and grid reliability. Utilizing mobile containerized and stationary test beds at INL's High Performance Computing (HPC) facility, this project integrates with various nuclear-related energy systems testing facilities. Key components include immersion cooling apparatus, absorption chillers, and adjustable thermal management simulators. Tasks involve acquiring necessary hardware, sensors, and cooling apparatus, engaging with data center industry stakeholders, and providing a testing platform for algorithms, models, tools, and software. The objective is to expedite the deployment of nuclear-powered data centers, thereby improving energy reliability and affordability.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2022

Dispatch optimization of a concentrating solar power system under uncertain solar irradiance and energy prices

The integration of thermal energy storage into a concentrating solar power system allows for mitigating some of the risk associated with uncertain solar irradiance and uncertain energy prices. We solve a 48 h dispatch optimization model with continually updated conditional point forecasts of both direct normal irradiance (DNI) and electricity prices with a rolling-horizon scheme at hourly resolution over the course of a year. Joint, conditional forecasts for DNI and prices are formed using an autoregressive moving-average time series model with exogenous weather predictors. We guide dispatch using a mixed-integer programming model, but in order to evaluate performance we use the System Advisor Model (SAM) of the National Renewable Energy Laboratory. SAM is a techno-economic simulation model that accounts for plant thermodynamics with higher fidelity. Our conditional DNI forecasts improve annual revenue by 4%–12% over using historical forecasts based on data from previous years. Conditional price forecasts improve annual revenue by 6%–19% in the real-time market over analogous historical forecasts. Updating these forecasts every six hours, rather than every 24 h, further improves annual revenue by 5%–6%. Here, we also investigate a method that values terminal inventory in our dispatch optimization model, again when used in a rolling-horizon scheme.

14 SOLAR ENERGY↗
DOE OSTI2025

Mapping Thermal Conductivity at the Atomic Scale: A Step toward the Thermal Design of Materials

We describe a spatial decomposition of the thermal conductivity, termed site-projected thermal conductivity, which quantifies the thermal conduction activity at each atomic site—a critical parameter for the thermal design of materials. The method is based on the Green–Kubo formula and the harmonic approximation and requires the force-constant and dynamical matrices, as well as a relaxed structural model. Throughout the study, it uses high quality models previously tested and compared to many experiments. It discusses the method and underlying approximations for amorphous silicon, carries the detailed analysis for amorphous silicon, and then examines an amorphous-crystal silicon interface and representative carbon materials. This study identifies the sites and local structures that reduce heat transport, and quantify these (estimate the spatial range) over which these “thermal defects” are effective. It identifies filamentary structures in the amorphous silicon network which impact heat transport and electronic structure (the Urbach edge) and electronic transport.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Numerical Simulation of Electron Magnetohydrodynamics with Landau-quantized Electrons in Magnetar Crusts

Abstract In magnetar crusts, magnetic fields are sufficiently strong to confine electrons into a small to moderate number of quantized Landau levels. This can have a dramatic effect on the crust's thermodynamic properties, generating field-dependent de Haas–van Alphen oscillations. We previously argued that the large-amplitude oscillations of the magnetic susceptibility could enhance the ohmic dissipation of the magnetic field by continuously generating small-scale, rapidly dissipating field features. This could be important to magnetar field evolution and contribute to their observed higher temperatures. To study this, we performed quasi-3D numerical simulations of electron MHD in a representative volume of neutron star crust matter, for the first time including the magnetization and magnetic susceptibility resulting from Landau quantization. We find that the potential enhancement in the ohmic dissipation rate due to this effect can be a factor ∼3 for temperatures of the order of 10 8 K, and ∼4.5 for temperatures of the order of 5 × 10 7 K, depending on the magnetic field configuration. The nonlinear Hall term is crucial to this amplification: without it, the magnetic field decay is only enhanced by a factor ≲2 even at 5 × 10 7 K. These effects generate a high wavenumber plateau in the magnetic energy spectrum associated with the small-scale de Haas–van Alphen oscillations. Our results suggest that this mechanism could help explain the magnetar heating problem, though due to the effect's temperature-dependence, full magneto-thermal evolution simulations in a realistic stellar model are needed to judge whether it is viable explanation.

Rau, Peter B. (ORCID:0000000152209277)↗
DOE OSTI2020

As a Matter of State: The Role of Thermodynamics in Magnetohydrodynamic Turbulence

Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Thermomechanical Modeling and Analysis of a High-Temperature Light Trapping Planar Cavity Receiver

Solar energy harnessed through concentrating solar power (CSP) systems offers a promising path to sustainable energy production, with the efficiency and longevity of these systems relying on key components like solar receivers. This study analyzes the thermomechanical behavior of an innovative enclosed light-trapping solar receiver optimized for particle heating applications. The receiver utilizes sheet metal alloys to form enclosed cavities that reflect and trap incoming solar flux, as well as enclosed channels that contain fluidized particle beds absorbing solar heat. Finite element analysis (FEA) is applied to predict the receiver's thermomechanical performance under extreme solar flux conditions. Temperature distributions from a thermal model simulating a multi-panel assembly at steady state are input into the FEA thermomechanical model for stress analysis. A key aspect of the analysis focuses on evaluating creep-fatigue damage, with a design target of achieving a 30- year service life. Various stress relief techniques are also proposed to extend the receiver's service life. The results highlight the significant impact of the particle-to-wall heat transfer coefficients (HTCs), ranging from 800 W/m2*K to 1400 W/m2*K. The 800 W/m2*K case shows a maximum von Mises stress of 164 MPa, while the 1400 W/m2*K case reduces it to 150 MPa. The creep life increases from 4,000 hrs in the 800 W/m2*K case to over 100,000 hrs in the 1400 W/m2*K case with Inconel 740H used, indicating that higher HTCs reduce stress and extend lifespan. This research advances the design of high-efficiency, low-stress solar receivers for particle-based thermal energy storage in CSP and industrial heating applications.

concentrating solar power↗
DOE OSTI2021

Renewable Thermal Hybridization Framework for Industrial Process Heat Applications

Renewable thermal energy systems (RTES) could be hybridized with different renewable options (e.g., flat plate collectors (FPCs) with parabolic collectors), or combined with existing heat supplies (e.g., fossil fuels), to give options for targeted solar IPH applications, industrial decarbonization and the reduction of fuel consumption. Buildings and industrial thermal energy applications require different temperature ranges, quantities, and rates of thermal energy, and as such require flexible, cost competitive solutions, able to provide heat over various temperature ranges. Hybrid solutions and thermal energy storage (TES) will be important for the dispatch of heat at optimal times needed by the demand side of the buildings and industrial applications. A variety of tools and platforms, such as the System Advisor Model (SAM), can provide hourly thermal yield simulations from single renewable energy (RE) technology options, including FPCs or CSP for solar IPH. We have investigated a variety of approaches to hybrid system modeling for RTES at different temperatures or combinations of technologies and developed an initial framework. The hybridization framework starts by creating a heat stream and raising the temperature of that heat stream by various combinations of RE technologies and other sources such as fossil fuels, new fuels or electric heating in multiple stages, with options for TES and/or waste heat recovery (WHR).

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

G-Function Library for Modeling Vertical Bore Ground Heat Exchanger

This library contains g-functions (thermal response functions) for standard, regularly spaced vertical borehole ground heat exchangers. In total, it contains 34, 321 configurations. To permit interpolation, each configuration has g-functions for heights of 24, 48, 96, 192, and 384 m. All the g-functions were calculated with burial depths of 2m, and borehole diameters of 15 to 17.5 cm, depending on height. In configurations with uniform spacing, the spacing between the boreholes is set to 5m, though it can be scaled to other horizontal spacings.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Analysis of Design Alternatives of Actively Cooled RF Window for MPEX

The Materials Plasma Exposure eXperiment (MPEX) has been designed as a linear plasma divertor simulator in order to address plasma material interaction (PMI) science for next-generation fusion devices. It will have the capability to test neutron irradiated samples with plasma fluxes of greater than 10 24 m -2 s -1 . It is expected to operate steady state for up to 10 6 s to consider PMI affects through reactor end of life. The conceptual design of MPEX was completed in 2019, with preliminary design having begun in 2020. The plasma source for MPEX is a helicon antenna, where the energized helical antenna sits outside of the vacuum in order to minimize impurities in the plasma. It is expected to receive up to 200 kW of continuous power, and so the antenna and the window must be actively cooled. The water-cooled copper antenna has been operated at full power on the Proto-MPEX device (which is a test facility to demonstrate the plasma source and heating systems). The water-cooled window, however, is a novel component that must meet numerous competing requirements. It requires a low dielectric loss to allow the Radio Frequency (RF) power to create the plasma within the vacuum boundary. It must be structurally robust to handle the significant heat flux from the plasma and any heat from dielectric coupling. It must be compatible with the coolant (preferably water). It requires a vacuum seal that minimizes impurities into the plasma and does not compromise the structural integrity of the window. Two window designs have been tested. Results from these tests, where temperatures are measured and heat fluxes inferred from infrared camera data, have been correlated with thermal-structural simulations. When these simulations are extrapolated to the full power steady-state heat fluxes that are expected in MPEX, the designs do not appear to have the necessary structural robustness. This study explores design alternatives for the MPEX helicon antenna window, presents analysis results for several of the alternatives, and shows a viable solution that satisfies the requirements for MPEX operation.

Tipton Jr, Joseph↗
DOE OSTI2026

Consumer safety-oriented scheduling of rotating power outages during heat waves

Extreme heat events have widespread effects on power systems, reducing available generation capacity, limiting transmission capabilities, and causing unusual demand patterns on the consumer side. As these combined effects expose bulk transmission systems to potential large-scale blackouts, utilities may be required to schedule and apply rotating outages, by temporarily and alternately disconnecting distribution substations to reduce overload. However, utilities lack mechanisms to inform these events, exacerbating the negative effects of heat waves on affected communities. This paper introduces a novel framework for scheduling rotating outages during heat waves while considering impacts on consumers’ safety. Instead of random sequential load shedding, we propose a methodology to rotate power outages considering a metric that quantifies the indoor overheating risk of groups of consumers during a power outage. The overheating risk is derived from a detailed building simulation using CityBES, where the buildings are modeled based on available data—use type, year built, floor area, number of stories, location—while presence of air conditioning and occupancy are calibrated from smart meter data. Based on the metric, an algorithm to schedule the rotating outages is applied to prioritize feeders for disconnection at each hour according to their overheating risk to meet a utility load reduction target. Applied to two substations and seven feeders in the Portland General Electric territory, the results show that this approach effectively leads to the lowest overheating risk during the resulting outage schedules, with an average 10.1% lower overheating compared to uninformed schedules.

Building thermal simulation↗
DOE OSTI2024

Creation of Multi-Principal Element Alloy NiCoCr Nanostructures via Nanosecond Laser-Induced Dewetting

The multi-principal element alloy nanoparticles (MPEA NPs), a new class of nanomaterials, present a highly rewarding opportunity to explore new or vastly different functional properties than the traditional mono/bi/multimetallic nanostructures due to their unique characteristics of atomic-level homogeneous mixing of constituent elements in the nanoconfinements. Here, in this paper, the successful creation of NiCoCr nanoparticles, a well-known MPEA system is reported, using ultrafast nanosecond laser-induced dewetting of alloy thin films. Nanoparticle formation occurs by spontaneously breaking the energetically unstable thin films in a melt state under laser-induced hydrodynamic instability and subsequently accumulating in a droplet shape via surface energy minimization. While NiCoCr alloy shows a stark contrast in physical properties compared to individual metallic constituents, i.e., Ni, Co, and Cr, yet the transient nature of the laser-driven process facilitates a homogeneous distribution of the constituents (Ni, Co, and Cr) in the nanoparticles. Using high-resolution chemical analysis and scanning nanodiffraction, the environmental stability and grain arrangement in the nanoparticles are further investigated. Thermal transport simulations reveal that the ultrashort (≈100 ns) melt-state lifetime of NiCoCr during the dewetting event helps retain the constituent elements in a single-phase solid solution with homogenous distribution and opens the pathway to create the unique MPEA nanoparticles with laser-induced dewetting process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simultaneously Improving Process Efficiency and Mechanical Properties in Aluminum Alloys with Applied Magnetic Fields

In the aluminum heat treatment sector, large resistive batch furnaces are used for many heat treatment processes and are generally always running, operating at approximately 20–30% efficiency. The primary avenue for energy savings lies in shortening heat treatment times, which can be 24+ hours for solutionizing (450–625 °C) and aging treatments (150–250 ℃), depending on the alloy and specification requirements. The application of high magnetic fields (<9 T) has shown to significantly reduce aluminum alloy heat treatment times, while achieving enhanced mechanical properties. Here, we explore the effects of thermomagnetic processing on the heat treatment response of a commercial aluminum alloy A390.1. Heat treatment times were reduced by >80%, while achieving ~10% improvement in microhardness, and these improvements persisted after longer thermal exposures simulating use. The microstructural evolution and resulting properties will be discussed.

Kesler, Michael [ORNL] (ORCID:0000000161423128)↗
DOE OSTI2024

Open-cycle thermochemical energy storage for building space heating: Practical system configurations and effective energy density

Salt-hydrate thermochemical materials (TCM) are promising candidates for energy storage systems for building space heating due to their high theoretical energy density and the need for low regeneration temperature. However, water vapor is required to drive the hydration process of the TCM reactor, which poses a challenge during winter when water vapor is typically scarce. Using indoor air directly lowers the building's humidity to an unconformable level in practice, while the cold outdoor air contains limited moisture. Here we consider different integration strategies for open-cycle TCM reactors in buildings and develop a model to simulate their thermal performance across diverse buildings and climates, specifically for building space heating. The potential energy densities and the levelized cost of storage of the TCM reactor are evaluated in practical scenarios to demonstrate the load-shifting potential of TCM systems for heating applications. We use a strontium chloride (SrCl 2 )-based composite as the baseline and explore the impact of various reactor and material changes to the energy density and levelized cost of storage.

25 ENERGY STORAGE↗
DOE OSTI2022

NekRS, a GPU-accelerated spectral element Navier–Stokes solver

The development of NekRS, a GPU-oriented thermal-fluids simulation code based on the spectral element method (SEM) is described. For performance portability, the code is based on the open concurrent compute abstraction and leverages scalable developments in the SEM code Nek5000 and in libParanumal, which is a library of high-performance kernels for high-order discretizations and PDE-based miniapps. Critical performance sections of the Navier–Stokes time advancement are addressed. Performance results on several platforms are presented here, including scaling to 27,648 V100s on OLCF Summit, for calculations of up to 60B gridpoints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

NOCI-F Electronic Couplings in Assemblies of Indolonaphthyridine Molecules: From Dimers to the Full Stack

Key electronic processes related to molecular excitonic states of finite stacks of indolonaphthyridine molecules are analyzed via the non-orthogonal configuration interaction with fragments (NOCI-F) method. Indolonaphthyridine is an organic chromophore that can undergo several electronic photoexcitation-related intermolecular processes, such as exciton and electron transfer. The structures studied here are noncrystalline arrangements built as either ordered stacks of indolonaphthyridine or stacks extracted from molecular dynamics simulations including thermal disorder. Taking dimers or trimers from either model, we performed CASSCF and NOCI-F calculations to quantify the intermolecular electronic couplings governing singlet fission, excited singlet and triplet diffusion, and hole and electron diffusion processes. Also, comparing the results for the different models, we studied the effect of structural disorder and distortion on these couplings. Finally, we present a newly developed, advanced postanalysis tool. It takes the NOCI-F data as input to carry out a multifragment full Hamiltonian procedure that involves the complete stack, providing physical information not available from the dimer/trimer models, hence giving access to additional insight into the material’s properties.

coupling reactions↗
DOE OSTI2021

Effect of Glacial/Interglacial Recharge Conditions on Flow of Meteoric Water Through Deep Orogenic Faults: Insights Into the Geothermal System at Grimsel Pass, Switzerland

Many meteoric-recharged, fault-hosted geothermal systems in amagmatic orogenic belts have been active through the Pleistocene glacial/interglacial climate fluctuations. The effects of such climate-induced recharge variations on fluid flow patterns and residence times of the thermal waters are complex and may influence how the geothermal and mineralization potential of the systems are evaluated. We report systematic thermal–hydraulic simulations designed to reveal the effects of recharge variations, using a model patterned on the orogenic geothermal system at Grimsel Pass in the Swiss Alps. Fault-bound circulation of meteoric water through the Grimsel fault is driven to depths of ~10 km by the high alpine topography. Simulations suggest that the current single pass flow is typical of interglacial periods, during which 1) meteoric recharge into the fault is high (above tens of centimeters per year), 2) conditions are at or somewhat below the critical Rayleigh number, and 3) the hydraulic connectivity along the fault plane is extensive (an extent of at least 10 km into increasingly higher terrain is required to explain the 10 km penetration depth). The subcritical condition constrains the bulk fault permeability to < 1e-14 m2. In contrast, the limited recharge during the numerous Pleistocene glaciation events likely induced a layered flow system with single pass flow confined to shallow depths while non-Rayleigh convection occurred in the deep fault. The same layering can be observed at low aspect ratios (length/depth) of the fault plane, when the available recharge area limits the flux through the fault.

58 GEOSCIENCES↗