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.
Development of the PROTEUS and ANSYS Coupled System for Simulating Heat Pipe Cooled Micro Reactors
The objective of this project is to develop the coupled system of PROTEUS/ANSYS for simulating a heat pipe cooled micro reactor. As an initial effort, while figuring out controlling ANSYS Mechanical with Workbench, the coupled system of PROTEUS/FLUENT in the CORBA environment were tried to demonstrate the coupled simulation for a transient problem of a heat pipe cooled micro reactor. Finally, the coupled system of PROTEUS/ANSYS was constructed using APDL commands in the CORBA environment, instead of ANSYS Workbench, since we found that ANSYS with the current version of Workbench Python scripting did not provide capabilities necessary for simulating a transient problem.
Thermo-hydraulic steam pipe models for district heating simulations: Simplifications to balance accuracy and simulation speed
Steam piping networks are essential for optimizing performance in industrial processes and district heating systems. However, dynamic models that balance thermo-hydraulic accuracy with computational efficiency remain limited. In response, this paper presents a new discretized steam pipe model based on the plug flow approach, capturing key thermo-hydraulic behaviors while simplifying steam phase change processes. Implemented in Modelica, the model accurately calculates temperature and pressure distributions along steam pipelines. To improve computational efficiency for district-scale simulations, five model simplifications are introduced: lumped thermo-hydraulic functions, empirical correlations, fluid state approximations, steady-state dynamics and inclusion of flow derivatives. These simplified models achieve 85%-98% accuracy in predicting pressure drop and condensation losses, including dynamic condensate behavior during pipe warm-up—a factor often overlooked in existing models. The models support diverse network configurations, scaling effectively to systems with multiple distribution pipes and connected building loads. Discrete models provide detailed insights but exhibit a cubic increase in simulation time as the network scales by N connected building O(N 2.42 ). In contrast, lumped models simulate 10–28 times faster than discrete, offering quadratic scaling of simulation time O(N 1.73 ). However, they still require 6 times more computation time than a lossless network, highlighting the inherent computational challenges of modeling compressible fluid flow. In conclusion, the steady-state lumped variant, with its near-linear scalability in computational time O(N 1.01 ), emerges as an efficient solution for preliminary design evaluations and extensive parametric studies.
Strain age cracking in a simulated heat affected zone of Inconel 740H laser-powder bed fusion components
Inconel® 740H components produced via laser-powder-bed fusion (L-PBF) additive manufacturing are ideal for supercritical CO2 primary heat exchangers. Arc welding is often needed, and subsequent post-weld heat treatment aging at 790–840 °C is required to improve strength via gamma prime (γ’) precipitation; however, strain age cracking (SAC) can occur in the heat affected zone (HAZ) during this process. This study uses stress relaxation testing at 800 °C on simulated HAZ specimens from vertically and horizontally built L-PBF IN740H to assess SAC susceptibility across heating rates of 40–3480 °C/h and weld-induced strains of 3–10 %. Vertical builds require greater mechanical energy input and exhibit longer times to fracture than horizontal builds, and time to fracture occurs sooner at slower heating rates. Creep voids were observed in γ’-denuded regions along grain boundaries, and cracking propagated along migrated grain boundaries and at interfaces between elongated secondary phases (γ’ or carbides) and the γ-matrix.
HEAT simulation and IR data comparison for ST40 plasma-facing components
Not Available
Experimental and simulated heating in nanocrystalline diamond charge exchange injection films using an electron beam to mimic SNS conditions
Nanocrystalline diamond (NCD) foils have been in use at the Spallation Neutron Source as primary stripper foils for charge exchange injection since the machine’s inception. NCD’s low atomic number leads to reduced beam scattering while still being easy to handle, as opposed to even lighter elements. NCD foils also have the benefit of being free-standing, and rigid when grown under the right conditions to minimize residual stresses. This study overviews a method to mimic the Spallation Neutron Source (SNS) beam and characterize thermally driven failure of NCD foils. To do this, a foil test stand was developed with in situ diagnostics that tracks signs of foil sublimation and thinning. The foil test stand’s electron beam is equated to the SNS beamline by comparing peak deposited energy densities for both beams. Post-mortem testing is also conducted to help elucidate changes the NCD foils undergo during exposure to the beam. Further, a COMSOL simulation was also developed as a method to assess potential future changes to the NCD film and predict how changes to beam conditions will effect foil temperatures. This method allows for examination of the NCD foil’s thermal limits and can be used to assess future changes to SNS stripper foils.
Demonstrating UPC++/Kokkos Interoperability in a Heat Conduction Simulation (Extended Abstract)
We describe the replacement of MPI with UPC++ in an existing Kokkos code that simulates heat conduction within a rectangular 3D object, as well as an analysis of the new code’s performance on CUDA accelerators. The key challenges were packing the halos in Kokkos data structures in a way that allowed for UPC++ remote memory access, and streamlining synchronization costs. Additional UPC++ abstractions used included global pointers, distributed objects, remote procedure calls, and futures. We also make use of the device allocator concept to facilitate data management in memory with unique properties, such as GPUs. Our results demonstrate that despite the algorithm’s good semantic match to message passing abstractions, straightforward modifications to use UPC++ communication deliver vastly improved performance and scalability in the common case. We find the one-sided UPC++ version written in a natural way exhibits good performance, whereas the message-passing version written in a straightforward way exhibits performance anomalies. We argue this represents a productivity benefit for one-sided communication models.
Development and validation of non-axisymmetric heat flux simulations with 3D fields using the HEAT code
A new comprehensive module to simulate heat fluxes from three-dimensional (3D) magnetic fields has been implemented in the HEAT code. Especially compact tokamaks like SPARC require tools to predict and manage large heat fluxes. Existing release versions of HEAT can only simulate axisymmetric heat flux on 3D plasma facing components. The new module uses an M3D-C1 perturbed equilibrium and the MAFOT code to trace field lines of the perturbed 3D magnetic field. Heat flux is then assigned to the resulting footprints via a 3D layer model. The model distinguishes between the scrape-off layer, the magnetic lobes and the private flux region, and employs only 0D parameters like the layer width, diffusive spread and the last closed flux surface position in the perturbed edge to generate a heat flux profile. The magnitude is normalized to the total input power. Resulting heat flux simulations are compared and validated against infrared measurements in the DIII-D tokamak with applied 3D fields; good agreement is found for several cases. The new module can now be applied to the SPARC tokamak; a preliminary result for applied rotating 3D fields is shown.
Effects of the heat source on the steady-state transport in gradient-driven global gyrokinetic simulations
Heat sources are routinely used in gradient-driven global gyrokinetic simulations to prevent profile relaxation and ensure steady-state turbulence and transport. To have any predictive capability, gradient-driven simulations need to be insensitive to the particular choice of a source. Several forms of the heat source are investigated using gyrokinetic δf particle-in-cell simulations. The ion heat flux is found to be sensitive to the source strength when the source rate is below a fraction (∼1/10) of the linear growth rate.
A fast and accurate modeling approach for water and steam thermodynamics with practical applications in district heating system simulation
In U.S. district heating (DH) systems, steam is the most common heat transport medium. Industry demand for new advanced modeling capabilities of complete steam DH systems is increasing; however, the existing models for water/steam thermodynamics are too slow for large system simulations because of computationally expensive algebraic loops that require the solution to nonlinear systems of equations. For practical applications, this work presents a novel split-medium approach that implements numerically efficient liquid water models alongside various water/steam models, breaking costly algebraic loops by decoupling mass and energy balance equations. New component models for steam DH systems are also presented. We implemented the models in the equation based Modelica language and evaluated accuracy and computing speed across multiple scales: from fundamental thermodynamic properties to complete districts featuring 10 to 200 buildings. Compared to district models with the IF97 water/steam model and equipment models from the Modelica Standard Library, the new implementation improves the scaling rate for large districts from cubic to quadratic with negligible compromise to accuracy. Additionally, for an annual simulation with 180 buildings, this translates to a computing time reduction from 33 to 1-1.5 h. These results are critically important for industry practitioners to simulate steam DH systems at large scales.
Coupled Multiphysics Simulations of Heat Pipe Microreactors Using DireWolf
DireWolf is a multiphysics software driver application designed to simulate heat pipe–cooled nuclear microreactors. Developed under the U.S. Department of Energy, Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the DireWolf software application’s objective is to provide the nuclear community with a design and safety analysis simulation capability. Based upon the NEAMS program Multiphysics Object-Oriented Simulation Environment (MOOSE) computational framework, DireWolf tightly couples nuclear microreactor physics, reactor physics, radiation transport, nuclear fuel performance, heat pipe thermal hydraulics, power generation, and structural mechanics to resolve the interdependent nonlinearities. DireWolf is capable of simulating both steady and transient normal reactor operation and several postulated failure scenarios. We will present the fundamental physics of heat pipe–cooled nuclear microreactors and the MOOSE-based software employed in DireWolf. Both steady and transient results for coupled reactor physics, radiation transport, and nuclear fuel performance are demonstrated.
Model Form Error Correction for a Black-Box Thermal Battery Heat Transfer Simulation
Thermal batteries are crucial for supplying power to high-consequence engineering applications such as rockets. Computational simulations have been developed to predict thermal battery behavior, but these simulations often suffer from modeling errors, including model form uncertainty. Addressing this uncertainty can be achieved by quantifying either the model discrepancy in the output or the model form error (MFE) in the governing equation. MFE is particularly valuable as it can be better extrapolated beyond observed outputs, which is essential for predictions involving changes in external system loading, system configuration and geometry, or output quantities. This paper employs a state estimation approach to estimate MFE using experimental data and then utilizes machine learning (ML) to model its relationship with state variables. A nonintrusive technique is used to estimate MFE in a black-box thermal battery heat transfer simulation. The trained machine learning model for MFE is then applied to correct simulation predictions under extrapolated initial conditions and battery configurations. In conclusion, the methodology's performance is evaluated using additional experimental data, demonstrating its effectiveness in improving prediction accuracy.
Simulations of Heat Transfer Using Tight-Fitting Twisted Tape Inserts for First Wall Cooling in Molten Salt Breeder Blankets
One of the major components in fusion energy systems is the fusion blanket, which has a vacuum vessel to contain the plasma. As part of the fusion blanket/vacuum vessel, the first wall and plasma-facing components require sufficient cooling to prevent material degradation during operation from the superheated plasma. Most fusion blanket concepts involve first wall and divertor coolant channels with heat transfer enhancements (HTEs) that are intended to withstand the incident high heat fluxes of 1 to 5 MW/m 2 . Twisted tape inserts are a proposed HTE that have been investigated previously for first wall cooling and monoblock divertor cooling channels and in other nonfusion heat transfer components. By inserting twisted tapes into straight pipes, the amount of turbulence in the system can be increased at lower Reynolds numbers by swirling the flow. This results in better heat transfer characteristics with marginal increases in frictional pressure losses. In particular, simulations of high-Prandtl-number fluids such as the proposed molten salt FLiBe in twisted tapes, which is prototypic to liquid immersion blankets, have not been previously explored. Here, in this study, we simulate various Prandtl numbers in pipes with twisted tape inserts using large eddy simulations to determine the effects of increasing Prandtl numbers on heat transfer performance. The quantities of particular interest are the Nusselt number and the friction factor, which were recovered using data reduction techniques to determine impacts on heat transfer and pressure losses. This work serves as a starting point for determining the feasibility of twisted tape inserts for liquid immersion blanket concepts.
Numerical simulations of membrane distillation systems with actively heated membranes
Membrane distillation (MD) is a thermal desalination process that is gaining attention for treating hypersaline brines. One recent approach uses composite membranes with a thermally conductive layer that delivers heat to the membrane–feed interface, where it drives evaporation. Though this increases single-pass recovery, it comes with the challenge of promoting lateral heat conduction through the thin membrane. Herein, we develop a 3D, computational fluid dynamics (CFD) model that simulates conjugate heat, mass, and momentum transport in the feed channel and composite membrane. We then use the CFD to verify a simpler numerical model that approximates the feed velocity field analytically. We validate the numerical model experimentally, and use it to investigate the impacts of the conductive layer thickness, feed channel geometry, and operating conditions on temperature and concentration polarization, permeate production, and specific energy consumption. Overall, we find that lateral heating increases permeate production at the expense of increased concentration polarization. In extreme cases, the concentration increases more than four-fold along the membrane surface. Furthermore, we show however, that the feed flow rate and conductive layer thickness can be tailored to mitigate concentration polarization, for only a small reduction in permeate production.
Prototype HB650 Cryomodule Heat Loads Simulations
During the design stages of the PIP-II cryomodules, many analytical calculations and FEA have been performed on simpler geometry in order to estimate the heat loads and also to optimize the design. To better analyze the cryomodule cold tests, simulations have been performed with MATLAB to determine the temperature of the main components during cool down and to determine the heat loads of the cryomodule. These simulations have been applied to the High Beta 650 MHz prototype cryomodule design and compared to the cold tests performed on it.
A time-parallel method for scalable heat transfer simulations of additive manufacturing
Here, a major challenge in simulating the thermal behavior in additive manufacturing processes is the disparate length and time scales between transport phenomena occurring in the melt pool and the component. A common simulation approach relies on spatial decomposition for parallel computing, but due to the nature of heat transfer in AM, where most of the computational expenditure is localized near the melt pool, the computational speedup from spatial parallelization saturates quickly. Therefore, additional parallelism by means of time-domain decomposition is needed to fully take advantage of high-performance computing (HPC) resources. This work introduces a time-parallel method to improve the computational scalability of additive manufacturing simulations on HPC systems, while maintaining high temporal resolution of heat transfer near the melt pool. The method, inspired by the nonlinear paraexp formalism, performs an iterative superposition of nonlinear solutions to the initial value problem, integrating the heat equation across overlapping time-parallel intervals. For a single layer of the NIST AMB2018–01 L7 benchmark problem, the method achieves a 38.51x speedup in wall-clock time with a maximum error in the global temperature solution of 0.99%. This reduces the total solution time from 196.72 min to 5.11 min on 128 nodes of the ORNL Frontier supercomputer. The tradeoff between accuracy and total wall-clock time is investigated and recommendations for time-parallel deployment for AM problems are made.
Near-field acoustical holography and acoustic power analysis of a simulated, highly heated supersonic jet
Although near-field acoustical holography (NAH) and acoustic intensity analysis have previously been used to investigate the apparent jet noise sources produced by military aircraft, explicit connections to supersonic jet characteristics cannot be made due to a lack of information about the exhaust plume. Here, to begin to bridge this gap and better understand the source information yielded by NAH, the current study instead applies NAH to a virtual measurement of the near-field pressures of a highly heated laboratory-scale supersonic jet generated by large-eddy simulation (LES). The holographic reconstructions of the pressure, particle velocity, and acoustic intensity are found to match the LES-generated acoustic field well and are used to calculate the acoustic power of the jet. The jet's calculated overall acoustic power is compared to the free-stream mechanical power, resulting in an acoustic efficiency of 1.5%. Ray-tracing of the acoustic intensity to the jet centerline generates an axial distribution of the acoustic power origin, showing that almost all the power originates from the supersonic portion of the flow and with the distribution peak upstream of the potential core tip. Holographic reconstruction of the pressures along the nozzle lipline captures the general spectral shape of the LES-generated pressures, though it underestimates the amplitude.