Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

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 379 records · Page 21

Anisotropic thermal behavior of extrusion-based large scale additively manufactured carbon-fiber reinforced thermoplastic structures

Large format additive manufacturing (AM) enables rapid manufacturing of large parts and structures with minimum waste in material and energy. Extrusion-based AM deposition processes provide parts with highly anisotropic thermal properties, which are not typically reflected in textbook values for these materials. In order to develop accurate models that describe the directionally dependent thermal behavior of these materials in processing and service, accurate measurements of specific heat capacity and thermal conductivity are required. Here, this work characterizes, documents, and analyzes the effect of the anisotropic nature of the extrusion-based deposition process on the specific heat capacity and thermal conductivity of the resulting AM products. All measurements were made over a temperature range of 20–180°C using the transient plane source technique, also referred to as the hot disk technique. Three of the most commonly used large format AM feedstock materials that utilize carbon fiber reinforcement were examined in this work: acrylonitrile butadiene styrene, polyphenylene sulfide and polyphenylsulfone. Finally, these findings can serve as a thermal design/process guideline for future large format AM applications.

36 MATERIALS SCIENCE↗

Numerical evaluation of AGR-2 fission product release

The AGR-2 experiment produced normal operation and accident condition fuel performance data for tri-structural isotropic (TRISO) particles with UCO and UO 2 kernels. Data for compacts with no failed particles and for compacts with one or more failed particles are available. To model the fission product diffusion and release from these compacts, it is important to consider the computational mesh as well as the heat conduction and mass diffusion properties of the materials in a TRISO particle. Code and solution verification studies with the Bison fuel performance code were performed to show that the code is computing correct solutions. Bison comparisons to post-irradiation examination data and PARFUME predictions were made for several sets of compacts. Furthermore, these comparisons considered silver, cesium, strontium, and krypton with and without failed particles. Bison predictions closely match those of PARFUME, while both codes generally overpredict post-irradiation examination data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Direct numerical simulations for hybrid rocket boundary layers: Performance modeling and scaling

This paper presents a comprehensive performance and scaling analysis of direct numerical simulations for reacting boundary layers, focusing on slab burner configurations. Using a PETSc-based finite volume CFD framework, the study evaluates the scalability and computational cost of flow, chemistry, and radiation evaluations across 2D and 3D simulations. Polymethyl methacrylate (PMMA) is the fuel with pure O 2 as the oxidizer, modeled using a detailed chemical kinetics mechanism with 113 species and 660 reactions. A ray-tracing-based radiation solver, designed for distributed memory applications, is implemented to model radiation heat transfer. Parallel scalability is analyzed for the coupled flow, chemistry, and radiation heat transfer processes. Weak and strong scaling studies are conducted on up to 15,000 computational ranks, revealing robust performance when flow cells exceed 200 per rank. Chemistry evaluations dominate the computational cost in large 3D simulations, accounting for approximately 40% of the total runtime, while flow processes contribute around 35%, and radiation solver contributions remain below 10% due to reduced evaluation frequencies. GPU accelerated chemistry evaluation, implemented with Zero-RK, demonstrates significant promise, achieving up to a 4x speedup for workloads exceeding 30,000 cells per GPU. However, diminishing returns are observed for smaller workloads due to CPU-GPU communication overhead. This study identifies key challenges, including memory bottlenecks and the effects of domain partitioning on flow scalability, while highlighting the potential of GPU-accelerated chemistry to reduce computational costs. In conclusion, these findings provide realizable run configurations for 2D, 3D, and GPU-accelerated cases, offering insights for optimizing reactive flow solvers.

CFD Scalability↗

Benefits assessment of cool skin and ventilated cavity skin: Saving energy and mitigating heat and grid stress

This study assessed the energy-saving and climate-adaptive potential of cool skin and ventilated cavity skin facade technologies in Seoul's high-rise apartment buildings. We created weather scenarios for historical, mid-term future, and long-term future conditions using Coordinated Regional Downscaling EXperiment (CORDEX) method. Building energy simulations were conducted on a South Korean high-rise apartment model to evaluate their performance under different weather conditions. The results indicate that cool skin and ventilated cavity skin technologies can save cooling energy during summers but lead to heating energy penalties in winters. Ventilated cavity skin outperforms cool skin, offering better cooling energy savings and reduced heating penalties. Combining both technologies yields the highest overall energy savings, with 7 %, 9 %, and 10 % cooling energy savings for cool skin, ventilated cavity skin, and the combined package, respectively. However, cool skin increases heating energy consumption by 5 %, while ventilated cavity skin has minimal impact on heating energy. These envelope technologies also reduce peak electricity demand by at least 5 %, 8 %, and 9 %, respectively. They contribute to heat stress reduction, enhance resilience, and decrease extreme heat risks for occupants during power outages by at least 18 % under various weather conditions. Considering the prevalence of aging high-rise apartments in South Korea, adopting these envelope renovation strategies can effectively reduce cooling loads, enhance thermal comfort, and boost resilience under future climates, while avoiding costly reconstruction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Visualizing Uranium Crystallization from Melt: Experiment-Informed Phase Field Modeling and Machine Learning

The focus of this project was to observe and simulate the solidification of uranium metal at the crystallographic level from its molten state. Melting experiments were conducted at two different scales to observe microstructural evolution using either a laboratory-scale induction furnace (hundreds of grams of metal) or a microscope heating stage (hundreds of milligrams of metal), respectively. Experimental parameters and characterization data were then used to inform a phase field model of gamma-U crystal growth as dendrites with or without secondary phase impurities in the form of uranium carbide particles. Finally, training datasets were generated by the phase field model as inputs to a neural network, developed with the aim of providing a faster, cheaper surrogate model for microstructural simulations within a given parameter space. Progress is reported herein for each of these task areas. Ultimately, 1) an optical microscope heating stage capability has been stood-up for uranium metal solidification studies, 2) a phase field model was advanced to simulate multiple uranium grains growing in the presence of carbide impurity particles and 3) a neural network was constructed and optimized to predict the microstructure features of individually growing uranium crystals.

36 MATERIALS SCIENCE↗

Evaluation of Mechanical and Thermomechanical Water Vapor Compression Techniques for Enabling High Temperature Lift Hydration-Based Chemical Heat Pumps

Achieving high temperature lifts (>200 K) via a chemical heat pump based on salt hydration/dehydration reactions requires the transport of water vapor from low to high pressure. Alternative compression approaches require condensing of low-pressure water vapor, pumping of liquid water, and subsequent evaporation when the low-side pressure corresponds to sub-ambient water saturation temperatures. Thus, this study compares four steam compression methods for use within a chemical heat pump system based on a reversible calcium oxide hydration/dehydration reaction with a temperature lift from 350 °C heat to >600 °C. Purely mechanical and thermochemical/mechanical compression technologies are considered. A parametric study of maximum allowable temperature, the isentropic efficiency of mechanical compressors, the effectiveness of heat exchangers, and the assumed allowable heat exchanger pressure drop is conducted to determine the mechanical and thermal energy consumed per kilogram of compressed steam. The system complexity in terms of the number of main system components, maximum pressure ratio, and maximum allowable temperature is estimated. Model results show an absorption-based steam compressor has the highest exergetic efficiency for the required chemical heat pump required conditions. As a result, this system configuration was then experimentally demonstrated to illustrate the impact of system performance on component effectiveness.

Advanced reactors↗

Analytical solution and parameter estimation for heat of wetting and vapor adsorption during spontaneous imbibition in tuff

An analytical expression is derived for the thermal response observed during spontaneous imbibition of water into a dry core of zeolitic tuff. Sample tortuosity, thermal conductivity, and thermal source strength are estimated from fitting an analytical solution to temperature observations during a single laboratory test. Here, the closed-form analytical solution is derived using Green's functions for heat conduction in the limit of “slow” water movement; that is, when advection of thermal energy with the wetting front is negligible. The solution has four free fitting parameters and is efficient for parameter estimation. Laboratory imbibition data used to constrain the model include a time series of the mass of water imbibed, visual location of the wetting front through time, and temperature time series at six locations. The thermal front reached the end of the core hours before the visible wetting front. Thus, the predominant form of heating during imbibition in this zeolitic tuff is due to vapor adsorption in dry zeolitic rock ahead of the wetting front. The separation of the wetting front and thermal front in this zeolitic tuff is significant, compared to wetting front behavior of most materials reported in the literature. This work is the first interpretation of a thermal imbibition response to estimate transport (tortuosity) and thermal properties (including thermal conductivity) from a single laboratory test.

36 MATERIALS SCIENCE↗

Improved multiphysics model of the High Temperature Engineering Test Reactor for the simulation of loss-of-forced-cooling experiments

We present a multiphysics model of the High Temperature Engineering Test Reactor for comparison with past and predict future loss-of-forced-cooling (LOFC) experiments. The approach selected combines (1) 3-D full-core superhomogenization-corrected neutronics, (2) 3-D full-core homogenized or semi-heterogeneous heat transfer (macroscale), (3) 2-D axisymmetric fuel rod heat transfer (mesoscale), and (4) 1-D thermal-hydraulics channels. Although large sensitivities remain and the available experimental data are insufficient for use in a rigorous validation exercise, the overall transient behavior is successfully reproduced. The novelty of our work includes (1) a new macroscale-mesoscale heat transfer coupling approach relying on gap conductance to drastically speed up numerical convergence by two orders of magnitude, (2) determination of a radial effective thermal conductivity that yields reasonably close results between the homogenized and semi-heterogeneous macroscale thermal models, and (3) a preliminary study of the reactor's early behavior following a LOFC event, thus enabling further assessment of numerical models against fission power measurements.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Modeling laser interactions with aluminum and tantalum targets using a hybrid atomistic-continuum model

A hybrid atomistic-continuum method can model the microstructure evolution of metals subjected to laser irradiation. This method combines classical molecular dynamics (MD) simulations with the two-temperature model (TTM) to account for the laser energy absorption and heat diffusion behavior. Accurate prediction of the temperature evolution in the combined MD-TTM method requires reliable accuracy in electron heat capacity, electron thermal conductivity, and electron–phonon coupling factor across the temperatures generated. This study uses the electronic density of states (DOS) obtained from first-principle calculations. The calculated electron temperature-dependent parameters are used in MD-TTM simulations to study the laser metal interactions in FCC and BCC metals and the phenomenon of laser shock loading and melting. This study uses FCC Al and BCC Ta as model systems to demonstrate this capability. When subjected to short pulsed laser shocks, the dynamic failure behavior predicted using temperature-dependent parameters is compared with the experimentally reported single-crystal and nanocrystalline Al and Ta systems. The MD-TTM simulations also investigate laser ablation and melting behavior of Ta to compare with the ablation threshold reported experimentally. In conclusion, this manuscript demonstrates that integrating the temperature-dependent parameters into MD-TTM simulations leads to the accurate modeling of the laser–metal interaction and allows the prediction of the kinetics of the solid–liquid interface.

36 MATERIALS SCIENCE↗

Sensitivity and vulnerability to summer heat extremes in major cities of the United States

Many cities are experiencing increases in extreme heat because of global temperature rise combined with the urban heat island effect. The heterogeneity of urban morphology also leads to fine-scale variability in potential for heat exposure. Yet, how this rise in temperature and local variability together impacts urban residents differently at exposure-relevant scales is still not clear. Here we map the Universal Thermal Climate Index, a more complete indicator of human heat stress at an unprecedentedly fine spatial resolution (1 m), for 14 major cities in the United States using urban microclimate modeling. We examined the different heat exposure levels across different socioeconomic and racial/ethnic groups in these cities, finding that income level is most consistently associated with heat stress. We further conducted scenario simulations for a hypothetical 1 °C increase of air temperature in all cities. Results show that a 1 °C increase would have a substantial impact on human heat stress, with impacts that differ across cities. The results of this study can help us better evaluate the impact of extreme heat on urban residents at decision-relevant scales.

54 ENVIRONMENTAL SCIENCES↗

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗

Operating Hardware Impact on the Heat Transfer Properties of Windows

Despite significant advancements in fenestration technology in the last two decades, the thermal transmittance of fenestration products is still significantly higher than that of walls. This corresponds to 60% of the total energy loss of a modern building envelope through the windows. Hence, further development and improvements of fenestration products are necessary. Increasingly stringent codes and standards for fenestration stimulate industry to work on improved solutions. Thus, it is crucial that assessment techniques are able to account for innovations accurately. The thermal effects of non-continuous hardware in window frames are currently ignored by international rating procedures. A preliminary investigation conducted by our team showed significant performance degradation in two of the three out-opening casement profiles caused by the presence of operating hardware. Frames with the structure made of vinyl and fiberglass consist of many air cavities that are penetrated by operating hardware made of highly conductive materials. In these frames, in order to have an accurate assessment, it may be required to employ three-dimensional modeling due to the convective nature of heat transfer within the cavities. However, in this study, we demonstrate that the three-dimensional (3D) effects of non-continuous hardware can be approximated accurately with simpler two-dimensional (2D) simulations. We then develop a simplified model based on weighted average capable of replacing the time- and computation-intensive 3D simulations with 2D simulations and validate it against market available frames and their corresponding hardware. Validation results show that our approximation technique results in discrepancies lower than 0.05 W/(m2K), or 3% of the total thermal transmittance. Thus, we conclude that simplified 2D simulation models may be used for predicting hardware impact in window frames with reasonable accuracy. As windows and glazing structures are becoming ever better thermally insulated, it is becoming even more important to be able to model the impact of the operating hardware on the total thermal performance in order to design the best windows possible and not let the operating hardware ruin an otherwise well-proven design, which is hence addressed in this study.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Gravothermal evolution of dark matter halos with differential elastic scattering

Here, we study gravothermal evolution of dark matter halos in the presence of differential self-scattering that has strong velocity and angular dependencies. We design controlled N-body simulations to model Rutherford and Møller scatterings in the halo, and follow its evolution in both core-expansion and -collapse phases. The simulations show the commonly-used transfer cross section underestimates the effects of dark matter self-interactions, but the viscosity cross section provides an accurate approximation for modeling angular-dependent dark matter scattering. We investigate thermodynamic properties of the halo, and find that the three moments of the Boltzmann equation under the fluid approximation are satisfied. We further propose a constant effective cross section, which integrates over the halo's characteristic velocity dispersion with weighting kernels motivated by kinetic theory of heat conduction. The effective cross section provides a good approximation to differential self-scattering for most of the halo evolution. It indicates that we can map astrophysical constraints on a constant self-interacting cross section to an SIDM model with velocity- and angular-dependent scatterings.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear dynamic analysis of solution multiplicity of buoyancy ventilation in a typical underground structure

Buoyancy ventilation is widely used in underground buildings, such as underground hydropower stations. Multiple solutions of buoyancy ventilation may exist in those underground structures. In this study, we developed a transient model comprising an ordinary differential equation system to describe buoyancy ventilation patterns in typical two-zone underground structures. Additionally, the accuracy of the model was validated. Nonlinear dynamical analysis was conducted to study multiple steady-state airflow. According to mathematical derivation, the configuration of one local heat source at the bottom corner introduces two stable solutions. The criterion to determine the stability and existence of solutions for more general scenarios was developed. Using this criterion, we obtained the multiple steady states of any two-zone underground buildings for different stack height ratios and the strength ratios of the heat sources. Furthermore, this criterion can be adopted for the design of buoyancy ventilation or natural smoke ventilation systems. Designers can change the height ratio of the stack or the heat ratio of two zones to induce the desired ventilation patterns. Finally, a case study was conducted with field measurements to demonstrate the use of the nonlinear dynamical analysis method to investigate the multiple steady states of buoyancy ventilation. Through the case study, we validated that the proposed criterion could produce the same result as the nonlinear dynamical analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Hydrogen Risk Assessment Models (HyRAM) (Version 3.0 Technical Reference Manual)

The HyRAM software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen infrastructure and transportation systems. HyRAM is designed to facilitate the use of state-of-the-art science and engineering models to conduct robust, repeatable assessments of hydrogen safety, hazards, and risk. HyRAM includes generic probabilities for hydrogen equipment failures, probabilistic models for the impact of heat flux on humans and structures, and computationally and experimentally validated first-order models of hydrogen release and flame physics. HyRAM integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hydrogen hazards (thermal effects from jet fires, overpressure effects from deflagrations), and assessing impact on people and structures. HyRAM is developed at Sandia National Laboratories for the U.S. Department of Energy to increase access to technical data about hydrogen safety and to enable the use of that data to support development and revision of national and international codes and standards. HyRAM is a research software in active development and thus the models and data may change. This report will be updated at appropriate developmental intervals. This document provides a description of the methodology and models contained in the HyRAM version 3.0. HyRAM 3.0 includes the new ability to model cryogenic hydrogen releases from liquid hydrogen systems, using a different property calculation method and different equations of state. Other changes include modifications to the ignition probability calculations, component leak frequency calculations, and addition of default impulse data.

08 HYDROGEN↗

Hydrogen Risk Assessment Models (HyRAM) (V.3.1) (Technical Reference Manual)

The HyRAM software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen infrastructure and transportation systems. HyRAM is designed to facilitate the use of state-of-the-art science and engineering models to conduct robust, repeatable assessments of hydrogen safety, hazards, and risk. HyRAM includes generic probabilities for hydrogen equipment failures, probabilistic models for the impact of heat flux on humans and structures, and experimentally validated first-order models of hydrogen release and flame physics. HyRAM integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hydrogen hazards (thermal effects from jet res, overpressure effects from deflagrations), and assessing impact on people and structures. HyRAM is developed at Sandia National Laboratories for the U.S. Department of Energy to increase access to technical data about hydrogen safety and to enable the use of that data to support development and revision of national and international codes and standards. HyRAM is a research software in active development and thus the models and data may change. This report will be updated at appropriate developmental intervals. This document provides a description of the methodology and models contained in HyRAM version 3.1. There have been several impactful updates since version 3.0. HyRAM 3.1 contains a correction to use the volume fraction for two-phase speed of sound calculations; this only affects cryogenic releases in which two-phase ow (vapor and liquid) is predicted in the orifice. Other changes include clarifications that inputs for tank pressure should be given in absolute pressure, not gauge pressure. Additionally, the interface now rejects invalid inputs to probability distributions, and the less accurate single-point radiative source model selection was removed from the interface.

08 HYDROGEN↗

Development of an Improved RELAP5-3D Model for the High Temperature Test Facility

High-temperature gas-cooled reactors (HTGRs) are rapidly approaching deployment. Confidence in transient analysis of these systems requires modeling and simulation tools that have been validated against data relevant to HTGR conditions. The High Temperature Test Facility (HTTF) is an integral effects thermal hydraulics test facility for prismatic HTGRs. In spring and summer of 2019, HTTF was used for a series of experiments that now serve as the basis for the Organization of Economic Cooperation and Development / Nuclear Energy Agency Thermal Hydraulic Code Validation Benchmark for High Temperature Gas-Cooled Reactors using HTTF Data (HTGR T/H Benchmark). Previous analyses as part of the HTGR T/H benchmark used a RELAP5-3D model developed at Idaho National Laboratory (INL) and demonstrated an ability to reproduce trends in the measured data but difficulties reproducing experimental values within their uncertainty. These difficulties were largely attributed to assumptions made during the development of the initial RELAP5-3D model, which predated the HTTF experiments. A significant cause of difficulty reproducing the measured temperatures may be the radial nodalization of the previous RELAP5-3D model. In this paper, we present a new RELAP5-3D model of HTTF with finer radial nodalization built to assess the impact of radial heat transfer. We describe the new model and compare it against the old one at full-power steady state and for the pressurized conduction cooldown (PCC) transient. These analyses are based on the code-to-code comparison exercise for the PCC problem of the HTGR T/H benchmark. We compare maximum block temperature as the primary figure of merit and include discussion on intracore natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling of vapor-liquid interactions in condensing ejectors

Ejectors are compact mechanical devices that utilize the expansion of a high-pressure primary fluid to entrain and compress a low-pressure secondary fluid by means of momentum transfer between the two streams of fluid. Condensing ejectors feature both momentum and heat transfer through the interaction between the vapor and liquid streams. The goal of this paper is to develop and validate a one-dimensional slug flow model to simulate the vapor–liquid interactions in two types of condensing ejectors – one with primary liquid and secondary vapor flows (Type I), and the other with primary vapor and secondary liquid flows (Type II). Control volume analysis of the mass, momentum, and energy balance in each phase and across the interface was conducted for the ejectors. The slug flow models for both types of ejectors are validated with published results in the literature. The friction coefficient on the inner wall of the ejectors and the interfacial heat transfer coefficient are identified as controlling parameters for the ejector performance in terms of pressure and temperature distributions along the axis. The detailed parametric study shows that the liquid inlet velocity and mixing tube diameter have a significant impact on the performance of the Type I ejectors, and the performance of the Type II ejectors is mainly controlled by the secondary liquid flow inlet temperature.

42 ENGINEERING↗