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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reduced Order Models Generation for HTGRs Pebble Shuffling Procedure Optimization Studies

This report provides an initial study for producing reduced-order models (ROMs) of pebble-bed high temperature gas reactor (HTGR) models for the purposes of design optimization. As an initial study, this work is meant to be exploratory---identifying useful workflows and methods for ROM generation---and not meant to be a catch-all analysis of HTGR ROM generation and usage for optimization. This report summarizes three tasks performed in Fiscal Year 2022: 1) the creation of HTGR model, 2) the sensitivity analysis of model design parameters, and 3) an introduction to ROM generation techniques. The representative HTGR model created in this work is a multiphysics equilibrium-core using the BlueCRAB (comprehensive reactor analysis bundle) reactor analysis application, coupling four physical phenomena: neutronics, streamline depletion, porous flow thermal hydraulics, and pebble heat conduction. Part of the model creation was identifying some design parameters and quantities of interest that are relevant in an optimization analysis and adjustable in the model. The sensitivity analysis utilized a polynomial chaos expansion methodology to compute global sensitivity metrics. This analysis showed that thermal hydraulics parameters and quantities of interest had a relatively small impact on simulation results. Finally, the ROM generation work involved exploring three different ROM methodologies: polynomial regression, a Gaussian process, and artificial neural networks. Using a cross-validation technique to characterize ROM performance, the Gaussian process and single-layer artificial neural networks showed the most promising results. Overall, this study was insightful and the lessons learned will be invaluable for the eventual development of an HTGR design optimization workflow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Snow thermal conductivity and conductive flux in the Central Arctic: Estimates from observations and implications for models

During the Arctic winter, the conductive heat flux through the sea ice and snow balances the radiative and turbulent heat fluxes at the surface. Snow on sea ice is a thermal insulator that reduces the magnitude of the conductive flux. The thermal conductivity of snow, that is, how readily energy is conducted, is known to vary significantly in time and space from observations, but most forecast and climate models use a constant value. This work begins with a demonstration of the importance of snow thermal conductivity in a regional coupled forecast model. Varying snow thermal conductivity impacts the magnitudes of all surface fluxes, not just conduction, and their responses to atmospheric forcing. Given the importance of snow thermal conductivity in models, we use observations from sea ice mass balance buoys installed during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition to derive the profiles of thermal conductivity, density, and conductive flux. From 13 sites, median snow thermal conductivity ranges from 0.33 W m −1 K −1 to 0.47 W m −1 K −1 with a median from all data of 0.39 W m −1 K −1 from October to February. In terms of surface energy budget closure, estimated conductive fluxes are generally smaller than the net atmospheric flux by as much as 20 W m −2 , but the average residual during winter is −6 W m −2 , which is within the uncertainties. The spatial variability of conductive heat flux is highest during clear and cold time periods. Higher surface temperature, which often occurs during cloudy conditions, and thicker snowpacks reduce temporal and spatial variability. These relationships are compared between observations and the coupled forecast model, emphasizing both the importance and challenge of describing thermodynamic parameters of snow cover for modeling the Arctic as a coupled system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The MOOSE Thermal Hydraulics Module

The Multiphysics Object-Oriented Simulation Environment (MOOSE) is an open-source object-oriented finite element framework written in C++ (Lindsay et al., 2022). The Thermal Hydraulics Module (THM) is an optional MOOSE physics module that provides capabilities for studying thermal hydraulic systems. Its core capability lies in assembling a network of coupled components, for instance, pipes, junctions, and valves. THM provides several new systems to MOOSE to enable and facilitate thermal hydraulic simulations, most notably the Components system, which provides a higher-level syntax to MOOSE’s lower-level objects. This system is extensible by the user, but the current library primarily includes components based on a one-dimensional, single-phase, variable-area, compressible flow model, as well as heat conduction.

42 ENGINEERING↗
DOE OSTI2023

Improved Particle Heat Transfer by way of Bimodal Particle Distributions for High Temperature Solar Thermal Energy

High temperature solar thermal facilities are looking to increase operating temperatures through novel heat transfer media, one such being solid particles. These particles operating at high temperatures will require transferring their thermal energy into another working fluid like supercritical carbon dioxide which can be used in advanced power cycles. Achieving high heat transfer between the particles and supercritical carbon dioxide is essential to high efficiency and low-cost operation. Therefore, optimizing the thermal conductivity of these particles is one potential way to ensure high performance. Traditionally, unimodal particle distributions have been employed in high temperature particle solar power plants. However, ambient temperature testing of bimodal particle distributions has revealed a superior thermal conductivity when compared to its unimodal counterpart at the same temperature. This data was obtained by certified, off-the-shelf instruments that can effectively simulate the conditions a particle would be exposed to in a high temperature solar thermal system. Data obtained in this way suggests that the increased thermal conductively imputed by a bimodal particle distribution is significant at working temperatures in solar facilities. Furthermore, the thermal conductivity of these bimodal particle distributions peaks when the best combination of large and small particles is applied. At high temperatures, binary particle distributions are compared to monodispersed distributions of larger particles where heat transfer is more prolific due to the increased surface radiation. Various thermal conductivity, porosity and heat exchanger models are explored in conjunction with data acquired up to 700 C.

Stout, Dallin (ORCID:0009000294586091)↗
DOE OSTI2023

Modelling and Simulation of a Thermoelectric Heat Pump with Micro-Channel Heat Transfer

Thermoelectric heat pumps (TEHPs) have advantages of modularity and simple design for heating and cooling in some applications. Models of thermoelectric (TE) heat pumps are widely studied. Nevertheless, most existing modeling work focuses on the TE material or air-sourced TE modules. TEHP modeling, especially the relationship between modules and heat exchangers, has not been comprehensively conducted. This work presents a water-to-water TEHP modeling framework that combines Goldsmid’s approach for TE material performance, Gnielinski’s correlation for convective heat transfer, and thermal balance theory for heat exchange networks. This combined framework provides an accurate theoretical analysis of the water-to-water TEHP system. Subsequently, the framework was used to empirically determine TE material properties (electric resistivity, thermal conductivity, and Seebeck coefficient) that minimize modeling errors versus experimentally observed values from the literature. Finally, an additional set of experimental TEHP data was used to validate the model, all with relative absolute deviations of approximately 10% when predicting heating capacity and 10%–20% when forecasting cooling capacity at a 30 K temperature lift. For future work, people can further develop models of TEHPs with an air-based heat sink on one side and water channels on the other side.

Wan, Hanlong↗
DOE OSTI2020

FAST-1.0 Software Acceptance Testing Report

The purpose of this document is to detail the software testing of FAST-1.0 through unit, integration, and assessment tests. More than 400 tests were designed to provide coverage of the requirements for FAST-1.0. These requirements were developed in the NRC Statement of Work (SOW) for NRC Agreement Number NRC-HQ-25-14-D-001 and are transcribed to this document as part of the Software Quality Assurance Plan for the FAST Code System, PNNL-28767. FAST 1.0 represents the merger of the FRAPCON and FRAPTRAN codes and as such FAST-1.0 performs steady state and transient fuel performance calculations described below. FAST-1.0 calculates the steady state response of light-water reactor fuel rods during long-term burnup. The code calculates temperature, pressure, and deformation of a fuel rod as functions of time-dependent fuel rod power and coolant boundary conditions. The phenomena modeled by the code include: 1) heat conduction through fuel and cladding to the coolant; 2) cladding elastic and plastic deformation; 3) fuel-cladding mechanical interaction; 4) fission gas release from the fuel and rod internal pressure; and 5) cladding oxidation. FAST is used to perform independent calculations for regulatory evaluations of fuel performance under normal operation and anticipated operational occurrences (AOOs). FAST calculates the temperature and deformation history of a fuel rod as a function of time-dependent fuel rod power and coolant boundary conditions. The phenomena modeled by FAST include: 1) heat conduction; 2) heat transfer from cladding to coolant; 3) elastic-plastic fuel and cladding deformation; 4) cladding oxidation; 5) fission gas release; and 6) fuel rod gas pressure. FAST code assessment, development and maintenance drive a significant portion of the NRC fuel research activities and the tools are used in a substantial number of regulatory products. Given the centrality of the FAST code to the effectiveness of fuel research, it is critical to assess, develop and maintain this tool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Semi-Annual Report for Horizontal Compact High Temperature Gas Reactor (HC-HTGR) Development during Performance Period October 2022-March 2023

Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed HC-HTGR in 3 years and support its commercialization as a safe, low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. Additionally, Argonne is providing analysis of the primary heat transport system to ensure temperatures in the reactor systems, structures and components with significant safety margins during normal operation and design basis accident scenarios. This third semi-annual report summarized the progress made at Argonne on the two tasks during the first half of FY23. As a part of the RCCS design task, a design process for the water panel and the system configuration was performed to improve the thermal performance of the RCCS. The updated water panel design had achieved enhanced thermal performance with 0.96 MWth capability with major design updates made in structural interfaces with the RPV and initial configuration of the water tanks. The loop configuration of the RCCS has been proposed to have two independently working loops for system redundancy. A water panel material study was performed focusing on the use of carbon steel in water systems. Additional modeling strategies of primary system thermal fluids analyses were developed to meet modeling needs that are not well suited for the 1D-3D assembly level approach. The first of these is a reduced order assembly model, in which fuel centered unit cells are used to represent a fuel assembly. The 2D approach used in this model is much more computationally efficient, making this model useful for full core transient scenarios where fuel to coolant heat transfer is still the dominant flow path. The next model is the 3D core conduction model to be used to analyze decay heat removal in loss of primary system flow scenarios. Because these scenarios require a large domain to be modeled, a homogenized core model is being pursued to reduce the required computational costs. To accurately model decay heat scenarios it is necessary to couple a RCCS model to the 3D core conduction model. A simplified case is presented to demonstrate how the coupling methodology will be applied to the full core 3D conduction model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗
DOE OSTI2022

Performance and power modeling and prediction using MuMMI and 10 machine learning methods

Energy-efficient scientific applications require insight into how high performance computing system features impact the applications' power and performance. This insight can result from the development of performance and power models. Here, in this article, we use the modeling and prediction tool MuMMI (Multiple Metrics Modeling Infrastructure) and 10 machine learning methods to model and predict performance and power consumption and compare their prediction error rates. We use an algorithm-based fault-tolerant linear algebra code and a multilevel checkpointing fault-tolerant heat distribution code to conduct our modeling and prediction study on the Cray XC40 Theta and IBM BG/Q Mira at Argonne National Laboratory and the Intel Haswell cluster Shepard at Sandia National Laboratories. Our experimental results show that the prediction error rates in performance and power using MuMMI are less than 10% for most cases. By utilizing the models for runtime, node power, CPU power, and memory power, we identify the most significant performance counters for potential application optimizations, and we predict theoretical outcomes of the optimizations. Based on two collected datasets, we analyze and compare the prediction accuracy in performance and power consumption using MuMMI and 10 machine learning methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

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.

42 ENGINEERING↗
DOE OSTI2021

Modeling of the High Temperature Test Facility Using RELAP-7

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is an electrically heated, helium cooled, experimental facility. The HTTF was modeled using RELAP-7 in both 1-D and 3-D for validation. Due to core symmetry, 1/6th of the HTTF core was modeled in 3-D using the coupled heat conduction and forced convection capabilities of RELAP-7. This served as a validation case for RELAP-7 and its capabilities to model advanced nuclear technologies such as high temperature gas-cooled reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Modeling of the High Temperature Test Facility PG-26 Transient Using RELAP-7

The High Temperature Test Facility (HTTF) is an electrically heated, helium cooled, experimental facility located at Oregon State University (OSU). Modeling of the HTTF PG-26 transient using a multi-app approach utilizing MOOSE and RELAP-7 was conducted. A 3-D model of the HTTF core, reflector, core barrel, and RPV was generated, and 3-D heat conduction throughout the structure was coupled to 1-D fluid flow results from RELAP-7. The model was able to accurately predict peak ceramic core temperatures experienced during the transient which helps validate the capability of RELAP-7 to model advanced gas-cooled nuclear reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Joint Geophysical and Numerical Insights of the Coupled Thermal‐Hydro‐Mechanical Processes During Heating in Salt

Abstract Salt offers an optimal medium for the permanent isolation of heat‐producing radioactive waste due to its impermeability, high thermal conductivity, and ability to close fractures through creep. A thorough understanding of the thermal‐hydrological‐mechanical (THM) processes, encompassing brine migration, is fundamental for secure radioactive waste disposal within salt formations. At the Waste Isolation Pilot Plant (WIPP), we conducted joint in situ geophysical monitoring experiments during active heating to investigate brine migration near excavations. This experiment incorporated electrical resistivity tomography (ERT) alongside high‐resolution fiber‐optic‐based distributed temperature sensing within a controlled heating experiment. Additionally, discrete element model (DEM) based numerical simulations were conducted to simulate THM processes during heating, providing a more mechanistic understanding of the coupled processes leading to the observed changes in geophysical measurements. During heating, resistivity shifts near the heater were reasonably explained by temperature effects. However, in more distant, cooler regions, the resistivity decrease exceeded predictions based solely on temperature. DEM simulations highlighted brine migration, propelled by pore pressure gradients, as the likely primary factor contributing to the additional resistivity decline beyond temperature effects. The comparison between the predicted ERT responses and observations was much improved when considering the effects of brine migration based on the DEM simulations. These geophysical and simulation findings shed light on brine migration in response to salt heating, enhancing our understanding of the coupled THM processes in salt for safe radioactive waste disposal.

58 GEOSCIENCES↗
DOE OSTI2026

Heat Transfer in Nuclear Waste Glasses: Measurements and Modeling of Thermal Radiation Properties

We measured and modeled near infrared extinction of nuclear waste glasses from 300 °C to 1150 °C to enable predictive radiation heat transfer and thermal conductivity estimates. A composition and redox informed model resolved contributions from key chromophores (Fe+2-O-Fe+3, V+4, free and bonded ?Si-OH groups) and, when present, spinel particles that can cause strong scattering. The model reproduced measured absorption from room temperature up to 1150 °C, with minor discrepancies near 1 µm (likely due to possible trace impurities) and 2.5 µm (linked to uncertainty in hydroxy groups). Spectra showed silicate melts were semitransparent mainly in the 0.5–4.0 µm window, responsible for radiation thermal conductivity that generally increases with increasing temperature. We quantified the dependence of effective thermal conductivity on dissolved water and provided distributions across >100 LAW/HLW/DFHLW melts at 1150 °C, supporting improved melter heat transfer modeling.

Ferkl, Pavel↗