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 325 records · Page 18

Tensorized Interior Radiative Heat Transfer for a Scalable and Calibrated Building Energy Simulator

Building energy simulation is a critical tool for developing and testing advanced control strategies, such as Reinforcement Learning (RL), to provide demand flexibility and affordable energy costs. The recently introduced Smart Buildings Control Suite (sbsim) provides a lightweight, scalable, and data-calibrated simulation environment based on a 2D finite-difference model. However, the initial model primarily focused on conductive and convective heat transfer, neglecting the significant impact of long-wave radiative heat exchange between interior surfaces. This paper presents a significant extension to the sbsim framework by incorporating a physically-grounded model for interior radiative heat transfer. Our primary contribution is the development and integration of a fully tensorized radiative heat transfer module, which preserves the computational efficiency and scalability of the original simulator. This was achieved by developing a pipeline for view factor calculation, including an algorithm to identify directly seeing surfaces within complex floor plans, and formulating the net radiation equations for efficient execution on modern hardware accelerators. We validate the numerical accuracy of our tensorized implementation by comparing its results against a traditional iterative approach, demonstrating identical outcomes. This enhancement increases the physical fidelity of sbsim, enabling more accurate training of RL agents for building energy optimization.

Ham, Sang woo↗

Evaluation of pressure reconstruction techniques for Model Order Reduction in incompressible convective heat transfer

This paper compares pressure reconstruction strategies in Model Order Reduction for incompressible flows with convective heat transfer. The Navier-Stokes equation are reduced along with the passive scalar transport equation for the temperature using the POD-Galerkin technique. Six different pressure reconstruction methods are evaluated, two of which are novel to the best of the authors’ knowledge. Accurate pressure reconstruction is key to avoid error buildup in when solving for the conservation of linear momentum at the reduced level. The six approaches are compared using Direct Numerical Simulations of convective heat exchange processes in a 3D Backward Facing Step with a heated cylinder. Additionally, when comparing time-averaged metrics, we observe that the reconstruction methods that approximate the reduced pressure field using techniques borrowed from full order models (mechanical analogy, pressure Poisson, and velocity supremizers) yield higher errors than the methods that seek to stabilize the reduced systems (reduced residual stabilization, artificial divergence, and Uzawa operator).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

Investigating Performance of Geothermal Reservoirs With Bi-Directional Heat Transfer for Seasonal Thermal Storage Using Thermal-Hydraulic Reservoir Simulations

This paper presents initial simulation results of deep subsurface seasonal thermal storage of waste heat from a natural gas-fired power plant during the summer months for providing heating during winter months at the Cornell University main campus in Ithaca, New York. An injection temperature profile is calculated with the MEnU tool and numerical reservoir simulations are conducted with the simulator FALCON. Both a porous medium and large single fracture geothermal reservoir are considered at about 3 km depth with initial reservoir temperature of 75 degrees C. Simulations indicate that a much smaller reservoir - which is technically easier and less risky to develop - is sufficient when operating the reservoir for seasonal thermal storage instead of for continuous heat production. For example, storing 13 GWh of heat seasonally in a 100 m thick porous medium based on the Galway formation requires about 65 m of well spacing to limit the yearly variation in production temperature to 10 degrees C. Operating this reservoir for continuous heat production would require a well spacing of 365 m to limit the thermal decline to 10 degrees C after 10 years. Neglecting wellbore heat losses and assuming uniform and homogeneous properties, a round-trip thermal efficiency of 100% can be obtained if the yearly averaged injection temperature matches the initial reservoir temperature.

district heating↗

Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials (Second Annual Progress Report)

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of nuclear reactor parts. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which can introduce internal material flaws, such as pores and anisotropy. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures due to exposure to high temperature, radiation and corrosive environment in a nuclear reactor. Thermal tomography (TT) provides a capability for non-destructive evaluation of sub-surface defects in arbitrary size structures. Thermal tomography is a computational method for heat diffusion-based imaging of solids, which provides 3D visualization of data from flash thermography measurements. We investigate thermal tomography imaging and nondestructive evaluation of stainless steel and nickel super alloy metallic structures produced with laser powder bed fusion (LPBF) additive manufacturing (AM) process. Metallic structures produced with LPBF contain defects, and there are limited capabilities to evaluate these structures non-destructively. Thermal tomography reconstruction of 3D apparent spatial effusivity provides information about AM structure geometry and internal material flaws. We study performance of thermal tomography in imaging of metallic structures through COMSOL computer simulations of transient heat transfer, and through reconstruction of data obtained from experimental measurements. Reconstruction of internal defects is investigated using a stainless steel 316L specimen with flat bottom hole (FBH) indentations, and Inconel 718 plate produced with laser powder bed fusion (LPBF) method, which contains imprinted hemispherical shape low density regions containing non-sintered metallic powder. The FBH’s have the same sizes as the imprinted defects in the LPBF specimens, but offer better imaging contrast. Thermal tomography reconstructions provide visualizations of internal defects, and allow for estimation of their sizes and locations. Detection sensitivity of TT is limited by noises. We investigate separation of signal from noise in thermography images using several machine learning (ML) methods, including new spatio-temporal blind source separation (STBSS) and spatio-temporal sparse dictionary learning (STSDL) methods. Performance of the ML methods is benchmarked using thermography data obtained from imaging stainless steel 316L and Inconel 718 specimens produced LPBF method with imprinted calibrated porosity defects. The ML methods are ranked by F-score and execution runtime. Finally, we investigate TT of AM stainless steel 316L specimen with imprinted internal porosity defects using relatively low-cost, small form factor infrared (IR) camera based on uncooled micro bolometer detector. Sparse coding related K-means singular value decomposition (SVD) machine learning, image processing algorithms are developed to improve quality of TT images through removal of Additive white Gaussian noise without blurring the images. Following initial qualification of an AM component for deployment in a nuclear reactor, a compact TT can also be used for in-service nondestructive evaluation (NDE). With capability to perform in-service NDE of the AM component lifecycle, TT data can be used for development of a component digital twin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal tomography 3D imaging of additively manufactured metallic structures

Thermal tomography is a computational method for heat diffusion-based imaging of solids, which provides 3D visualization of data from flash thermography measurements. We investigate thermal tomography imaging and nondestructive evaluation of stainless steel and nickel super alloy metallic structures produced with the laser powder bed fusion (LPBF) additive manufacturing (AM) process. Metallic structures produced with LPBF contain defects, and there are limited capabilities to evaluate these structures non-destructively. Thermal tomography reconstruction of 3D apparent spatial effusivity provides information about AM structure geometry and internal material flaws. We study performance of thermal tomography in imaging of metallic structures through COMSOL computer simulations of transient heat transfer and through reconstruction of data obtained from experimental measurements. Thermal tomography reconstructions of structure shape and dimensions are shown for the Inconel 718 AM structure which has variations in the horizontal plane but is uniform along the depth dimension. Reconstruction of internal defects is investigated using a stainless steel 316L specimen with flat bottom hole (FBH) indentations, and the Inconel 718 plate is produced with the LPBF method, which contains imprinted hemispherical shape low density regions containing non-sintered metallic powder. The FBHs have the same sizes as the imprinted defects in the LPBF specimens but offer better imaging contrast. Thermal tomography reconstructions provide visualizations of internal defects and allow for estimation of their sizes and locations. Results of this study demonstrate that thermal tomography can be used for visualization and quality control in AM.

36 MATERIALS SCIENCE↗

Expansion-Driven Self-Magnetization of High-Energy-Density Plasmas

Understanding plasma self-magnetization is one of the fundamental challenges in both laboratory and astrophysical plasmas. Self-magnetization can modify plasma transport properties, altering the dynamical evolution of plasmas. Multiple high-energy-density (HED) experiments have observed the formation of ion-scale magnetic filaments of megagauss strength, though their origin remains debated. Here, in this study, we conduct 2D collisional particle-in-cell (PIC) simulations with a laser ray-tracing module for a fully self-consistent simulation of the plasma ablation, expansion, and magnetization. The simulations use a planar geometry, effectively suppressing the Biermann magnetic fields, to focus on anisotropy-driven instabilities. The laser intensity is varied between 10 13 and 10 14 W/cm 2 , which is relevant to HED and inertial fusion experiments where collisions must be considered. We find that, above a critical intensity, the plasma rapidly self-magnetizes via an expansion-driven Weibel process, producing a plasma beta of 100 (𝛽 = 8⁢𝜋⁢𝑘 𝐵 ⁢𝑛 𝑒 ⁢𝑇 𝑒 /𝐵 2 ) and Hall parameter 𝜔 ce ⁢𝜏 𝑒 >1 within the first few hundred picoseconds. The magnetic field is sufficiently strong to modify plasma heat transport, and simulations with an artificially suppressed magnetic field show noticeably different temperature profiles.

Lezhnin, K. V. [Princeton Plasma Physics Laborator↗

Enhanced Models for the Simulation of Two-Phase Boiling Flow and Critical Heat Flux with the STAR-CCM+ and NEK-2P CFD Codes: Executive Summary (Final Report)

The goal of the project was to develop, verify and validate an advanced two-phase boiling flow modeling capability that enhances the existing two-phase flow modeling of the STAR-CCM+ Computational Fluid Dynamics (CFD) code. The project focused on the simulation of two phase flow and heat transfer phenomena that occur in Boiling Water Reactor (BWR) and Pressurized Water Reactor (PWR) fuel bundles under various operating conditions. The extended STAR-CCM+ two-phase flow models implemented during the project follow the approach similar to the Extended Boiling Framework [1-6] previously developed at Argonne but include more fundamental physical models of boiling phenomena. The development and validation of the extended two-phase models is supported by the U.S. DOE GAIN program and performed in collaboration with Framatome U.S. Inc.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Active Divertor Heat Flux Control using Impurity Powder Dropper

Divertor plasma-facing components (PFCs) in a tokamak are typically designed to withstand average steady-state heat loads of about 5–10 MW/m 2 , a limit that applies to both solid and liquid lithium (LL) PFCs. Exceeding these design values can result in surface damage to tungsten PFCs or excessive lithium (Li) evaporation in liquid lithium divertor (LLD) PFCs. Since exceeding the divertor heat load limits has serious consequences, it is therefore prudent to develop a tool to reduce the divertor heat load and bring the heat load to within the design limit without affecting the plasma performance. Active low Z impurity injection such as Li has been suggested as a potential solution to mitigate excess heat flux as suggested previously, given that non-coronal radiation can be quite large ~ 20–30 MJ per mole of injected Li. Li is considered desirable for reducing the edge neutral recycling helping to improve plasma energy confinement. In this paper, we model the Impurity Power Dropper (IPD) to investigate its potential of divertor heat flux control. The IPD is typically located at the top of the tokamak device and uses a vertical drift tube of a few meters. In the 2 m drift tube case, the IPD powder is accelerated to ~ 6 m/sec before reaching the plasma with the upper divertor configuration, matching the condition for the in-board side pellet injection case. By modeling the IPD geometry we determined the IPD powder deposition profile, and thus the non-coronal radiation and ionization profiles in time as well. From the enhanced radiation power loss, it is therefore possible to reduce the divertor heat load using the divertor simulation code. In conclusion, the IPD divertor heat flux control can be tested in the facilities with IPD including ST-40, DIII-D, EAST, WEST and NSTX-U.

Active lithium injection↗

Daytime urban heat stress in North America reduced by irrigation

There is considerable uncertainty regarding the impact of irrigation on heat stress, partly stemming from the choice of heat stress index. Moreover, existing simulations are at scales that cannot appropriately resolve population centres or clouds and thus the potential for human impacts. Using multi-year convection-permitting and urban-resolving regional climate simulations, we demonstrate that irrigation alleviates summertime heat stress across more than 1,600 urban clusters in North America. This holds true for most physiologically relevant heat stress indices. The impact of irrigation varies by climate zone, with more notable irrigation signals seen for arid urban clusters that are situated near heavily irrigated fields. Through a component attribution framework, we show that irrigation-induced changes in wet-bulb temperature, often used as a moist heat stress proxy in the geosciences, exhibit an opposite sign to the corresponding changes in wet bulb globe temperature—a more complete index for assessing both indoor and outdoor heat risk—across climate zones. In contrast, the local changes in both wet-bulb and wet bulb globe temperature due to urbanization have the same sign. Here, our results demonstrate a complex relationship between irrigation and heat stress, highlighting the importance of using appropriate heat stress indices when assessing the potential for population-scale human impacts.

54 ENVIRONMENTAL SCIENCES↗

The effects of multispecies Hohlraum walls on stimulated Brillouin scattering, Hohlraum dynamics, and beam propagation

We report experiments and simulations have been conducted to investigate the efficacy of Ta 2 O 5 -lined Hohlraum walls at reducing stimulated Brillouin backscattering (SBS) as well as any subsequent effects on the Hohlraum dynamics and capsule implosions in indirect drive experiments at the National Ignition Facility. Using a 1.1 MJ 400 TW, 351 nm, shaped laser pulse, we measure a 5× reduction in SBS power in the peak of the pulse from the wall on the outer 50° cone beams. The SBS spectrum indicates a reduction in the high-Z spectral signature when using multispecies wall materials. Detailed hydrodynamic simulations were performed using different heat conduction models with flux limiters. Additional simulations were performed on the plasma maps using the 3D parallel paraxial code pF3D to compare backscatter powers between the pure Au and Ta 2 O 5 -lined Hohlraums. Further analysis, using hydrodynamically equivalent plasmas, shows that the SBS reduction is clearly a result of the added ion Landau damping caused by the oxygen ions and not from differences in plasma conditions. The experimental and simulation results also show an increase in the wall plasma expansion when using the Ta 2 O 5 liner leading to a 70% more oblate implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of pre-cooling as a recommended measure to improve residential buildings’ thermal resilience during heat waves

More intense heat waves are expected to occur more frequently in the twenty-first century. During severe heat waves, cooling capacity shortfall and overheating are likely to occur in residential buildings, and this will adversely affect occupant's thermal comfort and productivity. We propose a strategy of pre-cooling the house during off-peak hours to mitigate overheating during heat waves. Additionally, simulation results of a prototype single-family house show that adopting the rule-based control (RBC) of pre-cooling thermostat setpoint schedule is effective in reducing thermal discomfort, and that the efficacy of pre-cooling depends upon several building characteristics. An optimized control (OC) of the thermostat setpoint schedule was developed based on the simulation of a prototype building. A simplified yet improved RBC (IRBC) pre-cooling schedule was then extracted from the OC schedule for practical implementation at a larger scale. The effects of the RBC schedule and IRBC schedule were evaluated in the King District of Fresno, which contains 814 residential buildings. Results show that both thermostat setpoint schedules can reduce overheating effectively and that IRBC is slightly better than RBC for most buildings. The findings support the California government's recommendation on pre-cooling to mitigate overheating, which can be further improved with an optimized thermostat setpoint schedule broadcast to residents through early alert messages before a heat wave.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Daytime convective development over land: the role of surface forcing

Water availability at the Earth surface determines the portioning of the surface heat flux into its sensible and latent components, that is, the surface heat flux Bowen ratio. The two components affect differently the surface buoyancy flux and thus the development and growth of the convective boundary layer. As a result, the Bowen ratio has a strong impact on the daytime dry and moist convection development over land. Here, we use two canonical modeling test cases, one for the shallow convection and one for the shallow-to-deep convection transition, to document the impact of the surface Bowen ratio on daytime convection development. A simple approach is used where results from simulations featuring the original setup are contrasted with simulations where surface sensible heat flux takes on values of the latent heat flux and vice versa. Such a change illustrates the key impact of the surface water availability without changing the total surface heat flux. Because of the larger surface buoyancy flux, simulations with the reversed surface heat fluxes feature faster deepening of the convective boundary layer and wider clouds once moist convection develops. Mean cloud base width of cumulus clouds increases with the boundary layer depth. For shallow-to-deep convection transition, comparison of simulations with the observed large-scale flow that feature substantial lower-tropospheric shear with simulations with no mean horizontal flow suggests that the organization of convective circulations within the boundary layer strongly affects initial size and organization of shallow convection, as well as the subsequent width increase once convection deepens and eventually transitions to deep precipitating convection. A simple explanation is provided of why deeper well-mixed convective subcloud layer results in wider clouds. The key is the larger width of boundary layer coherent updraft structures when the mixed layer is deeper.

54 ENVIRONMENTAL SCIENCES↗

Fierro

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗

Fierro Version 2.x

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three-dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗

Ground heat exchanger design tool with rowwise placement of boreholes

Simulation-based design tools have been used since the late 1980s for designing ground heat exchangers (GHE) used with ground source heat pump (GSHP) systems. The ground heat exchanger simulations used in these tools rely on thermal response functions known as g-functions. Because of the significant computational burden in computing g-functions for even a single configuration, the design tools have relied on libraries of pre-computed g-functions. These g-functions were available for standard configuration shapes, such as lines, rectangles, open rectangles, L-shapes, and U-shapes. Standard shapes are often sub-optimal. For any building on a site, the available land may preclude use of a standard shape. For large GSHP systems with significantly imbalanced annual heat rejection and extraction loads, large rectangular fields may experience significant heat build-up (or heat draw-down) in the interior of the field. This paper describes a new ground heat exchanger design tool capable of automatically selecting and sizing both standard and irregular configurations. The focus of this paper is a method for creating, selecting and sizing irregular configurations where the available land area and "no-go" zones are described as irregular polygons.

Spitler, Jeffery↗

Effect of self-generated magnetic fields on x-ray emission in Kr-filled targets at the National Ignition Facility

We examine the effects of self-generated magnetic fields in a Kr gas pipe x-ray source platform. X-ray emission from Kr plasma is dependent on the plasma conditions, as the ionization state is largely a function of temperature. Magnetic fields are known to limit heat conduction, which increases temperature. We show that the emission in simulations of the gas pipe x-ray source is dependent on how self-generated magnetic fields are modeled. The inclusion of self-generated magnetic fields in simulations more accurately captures the emission of lower energy x-ray emission (L-shell), bringing results closer to experiments. The modeled x-ray emission and self-generated magnetic fields are shown to be particularly sensitive to the inclusion of the Nernst effect in simulations. Severely limiting the Nernst effect leads to a hotter Kr plasma, which can account for the discrepancy seen in earlier studies. By modifying the Nernst effect multiplier, we can achieve better experimental agreement in x-ray emission from gas pipes; the value of the multiplier that leads to the best agreement is dependent on the laser power of the drive. Currently, the suppression factor of the Nernst effect needed for high power drives (PL>200 TW) is more restrictive than what is currently put forward by non-local models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Changes in precipitation and air temperature contribute comparably to permafrost degradation in a warmer climate

Abstract Surface energy budgets of high-latitude permafrost systems are poorly represented in Earth system models (ESMs), yet permafrost is rapidly degrading and these dynamics are critical to future carbon-climate feedback predictions. A potentially important factor in permafrost degradation neglected so far by ESMs is heat transfer from precipitation, although increases in soil temperature and thaw depth have been observed following increases in precipitation. Using observations and a mechanistic ecosystem model, we show here that increases in precipitation hasten active layer development beyond that caused by surface air warming across the North Slope of Alaska (NSA) under recent and 21st century climate (RCP8.5). Modeled active layer depth (ALD) in simulations that allow precipitation heat transfer agreed very well with observations from 28 Circumpolar Active Layer Monitoring sites (R 2 = 0.63; RMSE = 10 cm). Simulations that ignored precipitation heat transfer resulted in lower spatially-averaged soil temperatures and a 39 cm shallower ALD by 2100 across the NSA. The results from our sensitivity analysis show that projected increases in 21st century precipitation deepen the active layer by enhancing precipitation heat transfer and ground thermal conductivity, suggesting that precipitation is as important an environmental control on permafrost degradation as surface air temperature. We conclude that ESMs that do not account for precipitation heat transfer likely underestimate ALD rates of change, and thus likely predict biased ecosystem responses.

54 ENVIRONMENTAL SCIENCES↗