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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

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↗
DOE OSTI2025

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↗
NASA NTRS1974

NASTRAN variance analysis and plotting of HBDY elements

The NASTRAN Thermal Analyzer has been intended to do variance analysis and plot the thermal boundary elements. The objective of the variance analysis addition is to assess the sensitivity of temperature variances resulting from uncertainties inherent in input parameters for heat conduction analysis. The plotting capability provides the ability to check the geometry (location, size and orientation) of the boundary elements of a model in relation to the conduction elements. Variance analysis is the study of uncertainties of the computed results as a function of uncertainties of the input data. To study this problem using NASTRAN, a solution is made for both the expected values of all inputs, plus another solution for each uncertain variable. A variance analysis module subtracts the results to form derivatives, and then can determine the expected deviations of output quantities.

Harder, R. L.↗
NASA NTRS1975

Design of experiments for measuring heat-transfer coefficients with a lumped-parameter calorimeter

A theoretical investigation was conducted to determine optimum experimental conditions for using a lumped-parameter calorimeter to measure heat-transfer coefficients and heating rates. A mathematical model of the transient temperature response of the calorimeter was used with the measured temperature response to predict the heat-transfer coefficient and the rate of heating. A sensitivity analysis was used to determine the optimum transient experiment for simultaneously measuring the heat addition during heating and the convective heat-transfer coefficient during heating and cooling of a lumped-parameter calorimeter. Optimum experiments were also designed for measuring the convective heat-transfer coefficient during both heating and cooling and cooling only.

Vanfossen, G. J., Jr.↗
NASA NTRS1975

Vibration and flutter analysis of reusable surface insulation panels

A study has been conducted of the aeroelastic behavior of heat shield panels for the thermal protection system of the space shuttle. Questions of structural and aerodynamic modeling are considered along with the analysis of a single isotropic tile and an isotropic main substructure. Representative flutter results are presented, taking into account a rigid and a flexible metallic substructure.

Dowell, E. H.↗
NASA NTRS1979

Review of the metallographic cooling rates of meteorites and a new model for the planetesimals in which they formed

The cooling rates of meteorites through approximately 900 -650 K, as read from their metal alloy compositions, are reviewed. Metallographic cooling rates are compared with the cooling rates that appear to be required by the K/Ar and Ar-40/Ar-39 ages of five meteorite classes, and discrepancies are found in all cases. Either (1) the metallographic cooling rates (and also Pu-244 fission cooling rates) are systematically in error, being too slow by a factor of approximately 6; or (2) the traditional thermal model for parent meteorite planets (having constant dimension and uniform physical properties) is oversimplified and the Ar closure temperatures for chondrites derived by Turner et al. (1978) are too low. An alternative parent planet model is proposed and numerically modeled, in which accretion of thermally insulating particulate matter, heat generation by Al-26 decay, melting or sintering of the particulate matter into conductive rock, and establishment of the properties of the meteorites occurred concurrently. Meteorite chronologies are somewhat easier to understand in this context, since the initially small, hot (thus sintered and conductive) bodies would have cooled rapidly to isotopic closure, but later cooling might have been much slower as a result of the continued accretion of insulating particulate matter.

Wood, J. A.↗
NASA NTRS1982

Heat transfer analysis of the Bridgman-Stockbarger configuration for crystal growth. Part 1: Analytical treatment of the axial temperature distribution

All first order effects on the axial temperature distribution in a solidifying charge in a Bridgman-Stockbarger configuration for crystal growth are analyzed on the basis of a one dimensional model whose validity can be verified through comparison with published finite difference ana;uses of two dimensional models. The model presented includes an insulated region between axially aligned heat pipes and considers the effects of charge diameter, charge motion, thickness, and thermal conductivity of a confining crucible, thermal conductivity change at the crystal-melt interface, generation of latent heat at the interface, and finite charge length. Results are primarily given in analytical form and can be used without recourse to computer work for both improve furnace design and optimization of growth conditions in a given thermal configuration.

Jasinski, T. J.↗
NASA NTRS1983

Numerical study of the effect of an embedded surface-heat source on the separation bubble of supersonic flow

A numerical study of the conjugated problem of a separated supersonic flow field and a conductive solid wall with an embedded heat source is presented. Implicit finite-difference schemes were used to solve the two-dimensional time-dependent compressible Navier-Stokes equations and the time-dependent heat-conduction equation for the solid in both general coordinate systems. A detailed comparison between the thin-layer and Navier-Stokes models was made for steady and unsteady supersonic flow and showed insignificant differences. Steady-state and transient cases were computed and the results show that a temperature pulse at the solid-fluid interface can be used to detect the flow direction near the wall in the vicinity of separation without significant distortion of the flow field.

Degani, D.↗
NASA NTRS1985

Crustal radiogenic heat production and the selective survival of ancient continental crust

It is pointed out that the oldest terrestrial rocks have so far revealed no evidence of the impact phase of Earth evolution. This observation suggests that processes other than impact were dominant at the time of stabilization of these units. However, a use of the oldest terrestrial rocks as a sample of the early terrestrial crust makes it necessary to consider the possibility that these rocks may represent a biased sample. In the present study, the global continental heat flow data set is used to provide further evidence that potassium, uranium, and thorium abundances are, on the average, low in surviving Archean crust relative to younger continental crust. An investigation is conducted of the implications of relatively low crustal radiogenic heat production to the stabilization of early continental crust, and possible Archean crustal stabilization models are discussed.

Morgan, P.↗
NASA NTRS1985

Crustal radiogenic heat production and the selective survival of ancient continental crust

It is pointed out that the oldest terrestrial rocks have so far revealed no evidence of the impact phase of earth evolution. This observation suggests that processes other than impact were dominant at the time of stabilization of these units. However, a use of the oldest terrestrial rocks as a sample of the early terrestrial crust makes it necessary to consider the possibility that these rocks may represent a biased sample. In the present study, the global continental heat flow data set is used to provide further evidence that potassium, uranium, and thorium abundances are, on the average, low in surviving Archean crust relative to younger continental crust. An investigation is conducted of the implications of relatively low crustal radiogenic heat production to the stabilization of early continental crust, and possible Archean crustal stabilization models are discussed.

Morgan, P.↗
NASA NTRS1993

Flow instability in particle-bed nuclear reactors

The particle-bed core offers mitigation of some of the problems of solid-core nuclear rocket reactors. Dividing the fuel elements into small spherical particles contained in a cylindrical bed through which the propellant flows radially, may reduce the thermal stress in the fuel elements, allowing higher propellant temperatures to be reached. The high temperature regions of the reactor are confined to the interior of cylindrical fuel assemblies, so most of the reactor can be relatively cool. This enables the use of structural and moderating materials which reduce the minimum critical size and mass of the reactor. One of the unresolved questions about this concept is whether the flow through the particle-bed will be well behaved, or will be subject to destructive flow instabilities. Most of the recent analyses of the stability of the particle-bed reactor have been extensions of the approach of Bussard and Delauer, where the bed is essentially treated as an array of parallel passages, so that the mass flow is continuous from inlet to outlet through any one passage. A more general three dimensional model of the bed is adopted, in which the fluid has mobility in three dimensions. Comparison of results of the earlier approach to the present one shows that the former does not accurately represent the stability at low Re. The more complete model presented should be capable of meeting this deficiency while accurately representing the effects of the cold and hot frits, and of heat conduction and radiation in the particle-bed. It can be extended to apply to the cylindrical geometry of particle-bed reactors without difficulty. From the exemplary calculations which were carried out, it can be concluded that a particle-bed without a cold frit would be subject to instability if operated at the high temperatures desired for nuclear rockets, and at power densities below about 4 megawatts per liter. Since the desired power density is about 40 megawatts per liter, it can be concluded that operation at design exit temperature but at reduced power could be hazardous for such a reactor. But the calculations also show that an appropriate cold frit could very likely cure the instability. More definite conclusions must await calculations for specific designs.

Kerrebrock, Jack L.↗
NASA NTRS2016

Electrically Driven Liquid Film Boiling Experiment

This presentation presents the science background and ground based results that form the basis of the Electrically Driven Liquid Film Boiling Experiment. This is an ISS experiment that is manifested for 2021. Objective: Characterize the effects of gravity on the interaction of electric and flow fields in the presence of phase change specifically pertaining to: a) The effects of microgravity on the electrically generated two-phase flow. b) The effects of microgravity on electrically driven liquid film boiling (includes extreme heat fluxes). Electro-wetting of the boiling section will repel the bubbles away from the heated surface in microgravity environment. Relevance/Impact: Provides phenomenological foundation for the development of electric field based two-phase thermal management systems leveraging EHD, permitting optimization of heat transfer surface area to volume ratios as well as achievement of high heat transfer coefficients thus resulting in system mass and volume savings. EHD replaces buoyancy or flow driven bubble removal from heated surface. Development Approach: Conduct preliminary experiments in low gravity and ground-based facilities to refine technique and obtain preliminary data for model development. ISS environment required to characterize electro-wetting effect on nucleate boiling and CHF in the absence of gravity. Will operate in the FIR - designed for autonomous operation.

Electrohydrodynamics↗
NASA NTRS2017

Mars 2020 Mobility Actuator Thermal Testing and Model Correlation

This paper describes the thermal testing and model correlation of a mobility actuator planned for use on the Mars 2020 Rover. The mobility actuator is identical to those which have been successfully flown and operated on the Mars Science Laboratory (MSL) Curiosity rover since its successful landing on Mars in August of 2012. The actuator consists of a motor, brake, and encoder paired with a four stage planetary gear box. In this thermal test, the actuator was instrumented with a number of thermocouples on both the interior and exterior of the gearbox. Heaters and a cold plate were used to generate thermal gradients across the actuator in vacuum, low pressure GN2, and low pressure CO2 environments in an effort to correlate a thermal model and develop a better understanding of how heat flows through the mechanism. This testing resulted in the successful model correlation of a simplified thermal model, and yielded important insights regarding the conductance of ball bearings and gear-to-gear contact. Ball bearing thermal conductance in a low pressure environment can be estimated by using correlations for vacuum thermal conductance along with a multiplier to account for increased gas conduction, and gear-to-gear conductance can be estimated by accounting for gas and grease conduction between two gears.

Novak, Keith↗
DOE OSTI2023

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↗
DOE OSTI2023

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↗
NASA NTRS2002

Characteristics of Pool Boiling on Graphite-Copper Composite Surfaces

Nucleate pool boiling performance of different liquids on graphite-copper composite (Gr-Cu) surfaces has been experimentally studied and modeled. Both highly wetting fluids, such as freon-113 and pentane, and a moderately wetting fluid (water) were tested on the Gr-Cu surfaces with different graphite-fiber volume fractions to reveal the enhancement effects of the composite surfaces on the nucleate pool boiling. Results of the experiments show that the graphite-fiber volume fraction has an optimum value. The Gr-Cu composite surface with 25 percent graphite-fiber volume (f=0.25) has a maximum enhancement effect on the nucleate boiling heat transfer comparing to the pure copper surface. For the highly wetting fluid, the nucleate boiling heat transfer is generally enhanced on the Gr- Cu composite surfaces by 3 to 6 times shown. In the low heat flux region, the enhancement is over 6 times, but in the high heat flux region, the enhancement is reduced to about 40%. For the moderately wetting fluid (water), stronger enhancement of nucleate boiling heat transfer is achieved on the composite surface. It shown the experimental results in which one observes the nucleate boiling heat transfer enhancement of 5 to 10 times in the low heat flux region and an enhancement of 3 to 5 times in the high heat flux region. Photographs of bubble departure during the initial stage of nucleate boiling indicate that the bubbles detached from the composite surface are much smaller in diameter than those detached from the pure copper surface. Typical photographs are presented.It shows that the bubbles departed from the composite surface have diameters of only O(0.1) mm, while those departed from the pure copper surface have diameters of O(1) mm. It is also found that the bubbles depart from the composite surface at a much higher frequency, thus forming vapor columns. These two phenomena combined with high thermal conductivity of the graphite fiber are considered the mechanisms for such a significant augmentation in nucleate boiling heat transfer on the composite surfaces. A physical model is developed to describe the phenomenon of bubble departure from the composite surface: The preferred site of bubble nucleation is the fiber tip because of higher tip temperature than the surrounding copper base and poor wettability of the graphite tip compared with that of the base material (copper). The high evaporation rate near the contact line produces the vapor cutback due to the vapor recoil pushing the three-phase line outwards from the fiber tip, and so a neck of the bubble is formed near the bubble bottom. Evaporation and surface tension accelerate the necking process and finally result in the bubble departure while a new small bubble is formed at the tip when the surface tension pushes the three-phase line back to the tip. The process is schematically shown. The proposed model is based on and confirmed by experimental results.

Zhang, Nengli↗
DOE OSTI2024

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↗
DOE OSTI2023

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↗