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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The formation of an radial edge electric field due to finite ion orbit width effects is the possible root cause of the H-mode edge

Abstract Full orbit-following simulations of thermal ions show that finite ion-orbit width effects create charge separation near the last closed flux surface (LCFS) which generates a localized radial electric field. Experimentally, edge electric fields are observed in H-mode plasmas and they are necessary for the edge turbulence suppression via the E × B flow shear mechanism. Confined trapped (and to a lesser extent co-passing) ions near the plasma edge form a positive charge distribution outside the LCFS, while thermal electrons are tied more tightly to field lines owing to their small mass and are poorly confined outside the LCFS, hence charge neutrality is violated outside the LCFS. A large number of reported observations from spherical and conventional tokamaks support the results from the simulations although the simulations were not performed fully self consistently. The results suggest ways to lower the H-mode power threshold and optimize the H-mode plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Which way does the dendrite grow? Competition among epitaxy, preferred growth direction, and thermal gradients in powder bed fusion additive manufacturing

The as-processed microstructure of metal alloy parts manufactured through laser powder bed fusion (LPBF) is heavily derived from the cellular dendritic solidification. The growth direction of dendrites within the melt pool is determined through competition among epitaxial growth, preferred growth directions, and maximum thermal gradients. However, the dominant factor and the specific role of each in developing melt pool microstructures remain unknown. Here, in this study, we performed single laser track scans on an SS316L single crystal substrate and combined experimental characterization of microstructure and crystal orientations with Computational Fluid Dynamics simulations of thermal gradients to evaluate the role of each factor in determining dendritic growth direction and evolution. Our results reveal that epitaxial growth dominates microstructure development by preferentially growing along a single 〈100〉 variant of the single crystal substrate adjacent to the melt pool boundary. Under LPBF’s highly curved and rapidly evolving thermal field, this preferential dendrite variant selection and its continued growth from the melt pool boundary to the centerline are governed by the local temperature gradient magnitude at the solid-liquid interface, rather than by the instantaneous maximum temperature gradient direction alone. Using these findings, we successfully predict changes in the dendrite growth direction with changing laser scan direction on a single crystal substrate, and show that the geometric melt pool centerline can deviate from the microstructural centerline because asymmetric local temperature gradient magnitudes transiently limit growth, resulting in different dendrite travel distances on each side of the melt pool.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantifying the thermal effect and methyl radical production in nanosecond repetitively pulsed glow discharges applied to a methane-air flame

In this work, we investigated non-equilibrium plasma produced by nanosecond repetitively pulsed glow discharges applied across a lean premixed methane-air flame. The flame is stationary, axisymmetric, and laminar. The discharges are applied on the symmetry axis crossing the reactant gases, flame front, and product gases, allowing phase-locked averaged measurements and comparisons with axisymmetric numerical simulations. The thermal effect and methyl radical production are quantified in the discharge in the reactant gas region. One-dimensional, two-beam, hybrid, femtosecond-picosecond, coherent anti-Stokes Raman scattering is used to acquire spatial and temporal profiles of temperature and oxygen-to-nitrogen concentration ratio. Photo-fragmentation laser-induced fluorescence is used to acquire quantitative two-dimensional profiles of methyl radicals in the discharge providing the first quantitative imaging of methyl produced ahead of a flame by plasma-induced methane dissociation. The spatial profiles of temperature and oxygen-to-nitrogen concentration ratio are in steady state, indicating that individual discharges have an insignificant heating effect. Upper and lower bounds of the produced mole fraction of methyl radicals in the plasma are obtained due to uncertainties in the collisional quenching rates of excited state methylidyne radicals in the plasma. The discharges produce a maximum of 600–1100 ppm of methyl radicals upstream of the flame front within 25 ns. This amount is similar to the predicted methyl mole fraction for the flame without plasma and thus represents a significant chemical perturbation to the reactants upstream of the flame front. The produced methyl follows an exponential decay in the first microsecond after the discharge with a decay constant of 8 µs close to the flame, and 0.8 µs further from the flame. The decay then deviates from the exponential curve and the methyl persists for tens of microseconds. The results suggest that for the tested configuration, the thermal effect of individual discharges through fast gas heating is negligible, while active chemical species are produced in large quantities in the reactant gases, upstream of the flame front.

hybrid fs-ps CARS↗
DOE OSTI2021

Design and performance evaluation of a dual-circuit thermal energy storage module for air conditioners

We present experimental results and a validated numerical model of a dual-circuit phase-change thermal energy storage module for air conditioners. The module incorporates a phase-change material encapsulated in compressed expanded natural graphite foam. We used n-tetradecane as the PCM with a transition temperature (~4.5 °C) suitable for air-conditioning applications. Heat exchange to and from the module is accomplished through two fluid loops operating as a heat source and sink embedded inside multiple slabs of the composite material. This dual-circuit design enables easier integration with air-conditioning equipment and provides enhanced flexibility in system operation as compared to the state-of-the-art thermal storage systems. When integrated with an air-conditioner, this design will enable peak-load shaving and enhances operational efficiency. The thermal storage device was designed for a nominal storage capacity of ~ 3.5 kWh. We evaluated the heat transfer and energy storage performance of this device using standalone heat transfer experiments to estimate key thermal resistances and identify design improvements before integration with an air conditioner. The numerical model of the heat exchanger uses a combination of discretized and lumped parameter approaches to maintain a balance between accuracy and computational expense. Our analyses show that the geometric features and integration of fluid tubes are key contributors to the thermal contact resistance between the fluid and the thermal storage material, and consequently, to the overall performance of the thermal storage module. Our standalone experiments also identified important operating scenarios in which this thermal storage module can be used for air-conditioning in buildings.

25 ENERGY STORAGE↗
DOE OSTI2026

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

WEC Irradiation Can Thermal Testing Report

With renewed interest for the next generation of small modular reactors and space-based nuclear reactors, metal hydrides have been a focus of research for moderating materials to reduce the need for highly enriched fuel to maintain small core sizes. Yttrium hydride (YH) is an attractive candidate for these metal hydrides due to its ability to retain hydrogen at elevated temperatures as well as the low neutron absorption cross section of yttrium metal. Given the necessity for reactors to operate at elevated temperatures for extended periods of time, understanding hydrogen retention in moderator channels in simulated operating thermal conditions is paramount in predicting reactor performance and the long-term stability of YH. Additionally, coating strategies on the moderator encapsulating will be explored to ascertain if it aids in retention of hydrogen within the capsule and moderator.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Random circuit block-encoded matrix and a proposal of quantum LINPACK benchmark

The LINPACK benchmark reports the performance of a computer for solving a system of linear equations with dense random matrices. Although this task was not designed with a real application directly in mind, the LINPACK benchmark has been used to define the list of TOP500 supercomputers since the debut of the list in 1993. We propose that a similar benchmark, called the quantum LINPACK benchmark, could be used to measure the whole machine performance of quantum computers. The success of the quantum LINPACK benchmark should be viewed as the minimal requirement for a quantum computer to perform a useful task of solving linear algebra problems, such as linear systems of equations. We propose an input model called the Random Circuit Block-Encoded Matrix (RACBEM), which is a proper generalization of a dense random matrix in the quantum setting. The RACBEM model is efficient to be implemented on a quantum computer and can be designed to optimally adapt to any given quantum architecture, with relying on a black-box quantum compiler. Besides solving linear systems, the RACBEM model can be used to perform a variety of linear algebra tasks relevant to many physical applications, such as computing spectral measures, time series generated by a Hamiltonian simulation, and thermal averages of the energy. We implement these linear algebra operations on IBM Q quantum devices as well as quantum virtual machines, and demonstrate their performance in solving scientific computing problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Fabrication and Performance Evaluation of Double-Sided Copper Nanowire Arrays as Thermal and Electrical Interfacial Layers

Reducing contact interface thermal and electrical resistances is in great demand across various industries, particularly in the semiconductor industry. This study introduces an approach using double-sided copper nanowire (Cu NW) arrays on copper sheets as both thermal and electrical interfacial layers, designed to effectively accommodate the topographical inconsistencies between contact surfaces. Further, experimental outcomes reveal a significant reduction in thermal contact resistance (TCR), with a value of 2.5 mm 2 K W –1 , thereby exceeding the efficiency of reported nanostructural thermal interface materials (TIMs). Additionally, when utilized as an electrical interfacial layer, these double-sided Cu NWs arrays dramatically reduced electrical contact resistance (ECR), outperforming traditional conductive grease in applications necessitating separable bonding, though showing comparable performance to costly silver-based conductive epoxies required for permanent, inseparable bonds. The promising results of the double-sided Cu NWs arrays in reducing both TCR and ECR, confirmed by finite element simulation, highlight their substantial potential in advancing TIMs and electrical interconnection applications across various sectors.

42 ENGINEERING↗
DOE OSTI2022

Design of a neutral thermal scattering (NeTS) module for hydrogen in light water

The accurate representation of thermal scattering law (TSL) data is integral to the design and characterization of many modern nuclear systems, particularly those using light water as a moderator/coolant. As a material-dependent distribution over energy-momentum phase space, the TSL may exhibit a variety of unique and relevant conditional (e.g., temperature, pressure) and compositional (e.g., porosity, stoichiometry, radiation damage) dependencies. Currently, there are various approaches to incorporate temperature dependence, which is especially important in coupled neutronic-thermal hydraulic simulations. In each approach, there is an inherent tradeoff between memory consumption and accuracy. Some techniques require tens to hundreds of MBs or more, while others fail to reproduce the underlying data to within 10% error over the considered input domain, despite having a reduced storage burden. This work aims to address both sides of the tradeoff simultaneously by implementing a novel deep learning (DL) approach to TSL representation. The neural thermal scattering (NeTS) concept, which is amenable to an arbitrary number of dependencies, is demonstrated via the inclusion of temperature dependence into a highly compact, highly accurate functional form of the multi- variate TSL (i.e., S(α, β, T)) for hydrogen in light water. Resulting storage requirements are on the order of 100 kB, and median and maximum percent deviations are on the order of 0.1% and 1%, respectively. These measures represent a step improvement over previous techniques. Notably, the developed neural network and feature methodology build on those employed in prior work on beryllium oxide. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Analysis of thermal response of new diagnostic probe in TCV

Thermal response of the new multi-purpose diagnostic probe inside the TCV (Tokamak à Configuration Variable) tokamak is studied. The New Probe Head (NPH) was constructed at the University of Innsbruck and developed in collaboration with various European laboratories within EUROfusion program. The NPH can be used in different MST devices by applying corresponding adapters. Its design, followed by the prototype, has to go through a number of tests (numerical and experimental), before it can be used in the tokamak environment, where it has to withstand repeated exposure to the Scrape-Off-Layer (SOL) plasma. Presented numerical analyses of the NPH thermal response inside the TCV tokamak are performed by the Jožef Stefan Institute.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Role of scan strategies and heat treatment on grain structure evolution in Fe-Si soft magnetic alloys made by laser-powder bed fusion

A major goal in printing soft magnetic Fe-Si steels using additive manufacturing is to take advantage of the potential for complex geometric designs and site-specific grain control. One major step in the processing of these alloys is understanding how processing parameters might impact how the as-built microstructure responds to annealing (i.e. the annealing response). The impact of scan strategy on the annealing response for thin wall geometries is specifically explored. Two scan strategies were explored for a thin wall geometry that produced a strongly columnar grain structure and equiaxed grain structure. Additionally, samples from both scan strategies annealed at 1200 °C showed a marked difference in annealing response with the more equiaxed sample seeing full recrystallization and grain growth, while the more columnar grain structure saw little change in microstructure. After analysis through characterization techniques and thermal-mechanical simulations Differences in internal energy within the grains were ruled out because calculated GND density values were similar for both samples. The formation of secondary particles was ruled out as a contributing factor due to the type of oxide formations and their size. It was concluded that the contributing factor to the difference in the annealing response were a difference in the resulting grain size and the density of high angle grain boundaries. These two differences were largely attributed to differences in the thermal gradient conditions due to grains preferentially growing in the direction of the steepest thermal gradient.

36 MATERIALS SCIENCE↗
DOE OSTI2025

stor4build

The EnergyPlus simulation engine supports modeling and simulation of thermal energy storage (TES) systems in several ways, including using the Python-EMS feature, which extends the operation of the engine with custom code written in Python. Creation of models using this feature can be tedious and error prone, with the connection of the model components to the Python code a particularly troublesome area. The stor4build Python package simplifies this process by modifying an input model to add a selected TES technology (implemented with the Python-EMS feature) and runs the simulation. The package leverages the OpenStudio middleware software development kit to automate this process as much as possible, eliminating potential errors and simplifying usage of EnergyPlus. The package provides objects, functions, and OpenStudio measures that implement the necessary operations to automate the creation of EnergyPlus models that integrate TES technologies with building systems. In addition, two user interfaces are provided: a command line interface and a web application programming interface. The automated process implemented by the package greatly simplifies the modeling and simulation process, allowing for parametric studies to be executed much more efficiently and effectively. The OpenStudio-based workflow is also very flexible and will allow for future additions of new technologies.

DeGraw, JasonWilliam [Oak Ridge National Laborator↗
DOE OSTI2025

Multifrequency eddy-current detection of fast transient thermal signatures for in situ monitoring applications

Eddy-current (EC) nondestructive evaluation has a long history of use in a variety of ex situ defect monitoring applications because of its exquisite sensitivity to local material variations. Due to the relationship between a material's conductivity and its temperature, EC methods have also been used to investigate quasistatic, long-range temperature variations in casting applications. However, these techniques remain underutilized for the measurement of rapidly varying, spatially nonuniform temperature distributions. In this work, we construct a model system capable of generating repeatable temperature transients in steel plates and measure real-time eddy-current signals with millisecond time resolution and spatial resolution of the order of 1 mm. Using a combination of Multiphysics simulations, fast thermal imaging, and time-resolved holographic interferometry, we tease apart contributions to the eddy-current signals arising from temperature variations and transient plate deformation. Finally, we perform a systematic study in which we vary the plate thickness and the eddy-current excitation frequency to demonstrate that eddy-current techniques can provide information about a three-dimensional, time-varying, subsurface thermal distribution, which is inaccessible to the traditional thermal imaging techniques.

Rosenberg, Ethan R. [Lawrence Livermore National L↗
DOE OSTI2025

Thermal and electric multidomain dynamic model for integration of power grid distribution with behind-the-meter devices

As renewable energy sources like solar and wind power become more integrated into the grid, coordinated control of behind-the-meter devices is crucial for enhancing grid flexibility and reliability and for meeting cost targets, with standardized models being developed to support this transition. The increasing flexibility and uncertainty of integrated renewable energy grids, along with interactions between various subsystems, make traditional steady-state modeling insufficient to capture transient and dynamic behaviors. Current models (e.g., composite load and battery equivalent models) focus on thermodynamic or electrical characteristics but overlook critical electromechanical interactions. This limits the ability to share performance information for grid services and hampers fast dynamic simulations. In addition, motor stalling is usually triggered by a fault event and attributed to the characteristics of the mechanical torque of the motor, resulting in absorption of a large amount of reactive power during the stalling period. Further, this significant withdrawal of reactive power will deteriorate the dynamic voltage stability of power grids and cause delayed voltage recovery. Therefore, an in-depth modeling of the thermodynamics or mechanical torque is essential to study the impacts of the realistic torque characteristics of those behind-the-meter devices on power system voltage stability. This study developed a dynamic multidomain model for building HVAC systems, such as air-source heat pumps, to simulate their thermal and electrical responses to grid transients. The model can accurately predict power metrics with a mean absolute percentage error of 10%, by validating against with power system computer-aided design performance data. Case studies demonstrate the model capability of capturing the transient response to sudden voltage changes, rapid load fluctuations, and system shutdowns respectively. During a sudden voltage drop (30% for 0.1s), a fully loaded heat pump’s motor speed dropped, continued declining, and shut down after 3.6s, with severe power oscillations and a torque spike. A partially loaded unit experienced temporary oscillations but stabilized. Under higher building loads, compressor speed increased from 64% to 100%, with power and torque rising before stabilizing. In safety-triggered shutdowns, power decreased after minor fluctuations, and torque briefly spiked before dropping to zero.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Nuclear fission reaction simulations in compact stars

Type-Ia supernovas are powerful stellar explosions that provide important distance indicators in cosmology. Recently, we proposed a new Type-Ia supernova mechanism that involves a nuclear fission chain-reaction in an isolated white dwarf [Phys. Rev. Lett. 126, 131101 (2021)]. Here, in this work, we perform novel reaction network simulations of the actinide-rich first solids in a cooling white dwarf. The network includes neutron-capture and fission reactions on a range of U and Th isotopes with various possible values for 235 U enrichment. We find, for modest 235 U enrichments, neutron capture on 238 U and 232 Th can breed additional fissile nuclei so that a significant fraction of all U and Th nuclei may fission during the chain reaction. Finally, we compute the energy release from the fission chain reaction for various uranium enrichments; a novel result that is a necessary input for thermal diffusion simulations of carbon ignition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Machine learning from RANS and LES to inform coarse grid simulations

Nuclear system thermal hydraulic analysis has historically relied on computationally inexpensive 1D codes. However, such tools are unable to capture multiscale multidimensional effects in large nuclear reactor enclosures. On the other hand, simulations with higher fidelity can be too expensive for such purposes. One of the ways to reduce computational cost is to perform simulations on a coarse grid, which, unfortunately, introduces large discretization errors. In this paper, two high-to-low data-driven approaches are investigated: (1) a coarse grid turbulence model to predict eddy viscosity and (2) correction of errors in coarse grid velocity fields. The approaches aim to reduce grid- and turbulence model-induced errors in coarse grid Reynolds-averaged Navier–Stokes (RANS) simulations. Two sources of high-fidelity data, RANS and large eddy simulations (LES), are explored. To extract the eddy viscosity from the LES data, an inverse optimization problem is solved. However, the LES eddy viscosity is shown to be comparable to the RANS eddy viscosity in terms of error reduction. Therefore, the directly available RANS eddy viscosity was used to develop a coarse grid data-driven turbulence model. Additionally, error correction in velocity is used to reduce the remaining uncertainties and bring the results closer to reality. In conclusion, the performance of the frameworks is demonstrated for a scaled upper plenum of a gas-cooled reactor facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗