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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A detailed study of pre-heating effects in electron beam melting powder bed fusion process

Metal-based additive manufacturing processes, such as powder bed fusion with electron beam (PBF-EB) process, also referred to as electron beam melting (EBM), can produce high-density parts with minimal residual stresses due to the uniform and coherent preheating of the powder bed. However, understanding and controlling the multiple stages of preheating is required to enable the production of high-quality, consistent parts of various materials. This work presents a large-scale, multi-layer, three-dimensional numerical analysis focused on studying the preheating stages for predicting thermal history during the PBF-EB process. The model follows a continuous multi-stage cyclic process, that incorporates all the main stages of the PBF-EB process for 316 L stainless steel. This includes the gradual deposition of a new powder layer, the first and second preheating levels of the powder bed, and the energy deposition during melting (excluding the actual melt-pool behavior simulation). The model employs an adaptive time-scaling approach that automatically adjusts the energy deposition for each solution time-increment. This allows for localized changes in time-resolution over an otherwise computationally expensive multi-layer procedure. The material property variations are also taken into account, with an emphasis on the subtle irreversible changes in powder effective thermal conductivity after the two requisite preheating stages of the powder bed. This effect is studied using simplified conductivity models from the literature for partially sintered powder, validated by a dedicated experiment and numerical simulation. The large-scale model is then used to estimate the actual temperatures during first and second preheating levels for 316 L steel, which is not yet fully supported commercially for PBF-EB. Model predictions are corroborated by experiments, using and analyzing IR images, taken at the completion of each layer by the machine’s built-in infrared camera. The current model also incorporates a qualitative assessment for the effects of conductivity change during pre-heating, as well as evaluates the applicability of the time-scaling approach.

36 MATERIALS SCIENCE↗

NASA NTRS1976

Oceanic lithosphere and asthenosphere - Thermal and mechanical structure

A coupled thermomechanical subsolidus model of the oceanic lithosphere and asthenosphere is developed which includes vertical heat conduction, a temperature-dependent thermal conductivity, heat advection by a horizontal and vertical mass flow that depends on depth and age, contributions of viscous dissipation or shear heating, a linear or nonlinear deformation law relating shear stress and strain rate, as well as a temperature- and pressure-dependent viscosity. The model requires a constant horizontal velocity and temperature at the surface, but zero horizontal velocity and constant temperature at great depths. The depth- and age-dependent temperature, horizontal and vertical velocities, and viscosity structure of the lithosphere and asthenosphere are determined along with the age-dependent shear stress in those two zones. The ocean-floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of ocean-floor age; seismic velocity profiles which exhibit a marked low-velocity zone are constructed from the age-dependent geotherms and assumed values of the elastic parameters. It is found that simple boundary-layer cooling determines the thermal structure at young ages, while effects of viscous dissipation become more important at older ages.

Schubert, G.↗

NASA NTRS1972

Effects of overlaps, stitches, and patches on multilayer insulation.

The effects of spacecraft installation-type discontinuities on aluminized Mylar multilayer insulation have been measured using an electrically heated cylindrical calorimeter. The experimental results have been reproducible within plus or minus 10% or better for identical blankets and less than plus or minus 5% for the same blanket. Good blanket performance can be obtained using short interleaving for joints, eliminating stitched seams, and using small patches, particularly over the hot side of stitches. The heat loss through stitching is excessive but highly localized. Analytical models matched to thermocouple test results are presented. The manner in which radiation and conduction participate in heat transfer adjacent to discontinuities is indicated. High lateral heat transfer is also established as an important phenomenon in the vicinity of a discontinuity.

Stimpson, L. D.↗

DOE OSTI2021

Modeling Molten Salt Spreading and Heat Transfer using MELTSPREAD – An Uncertainty Analysis

This report summarizes spreading and heat transfer results for molten salt pouring onto a flat stainless steel substrate that were calculated using MELTSPREAD to gain insight into bulk salt behavior during a molten salt spill accident. MELTSPREAD was developed at Argonne to model the one-dimensional flow and freezing behavior of molten corium and is being applied to model molten salts for the first time as part of this work. An uncertainty analysis of thermophysical properties and spill conditions was performed using MELTSPREAD to determine which properties and conditions have the greatest impact on model outcome. The model was run with and without the inclusion of decay heat. Eutectic FLiNaK was used as the salt composition because it is well-characterized and appropriate for model development. Small salt volumes were used for the initial model runs described in this report so that the model results can be compared to the results from ongoing experiments on molten salt spreading and heat transfer currently being conducted at Argonne at a benchtop scale. The spreading of a small volume of FLiNaK (50 mL) was found to be limited by the balance between the molten salt surface tension and gravity and not by freezing for the pour conditions and assumed mechanisms of heat transfer in the model. Heat transfer from the salt was highly inefficient to due to its low thermal conductivity, high heat capacity, and high heat of fusion and varying these properties had no effect on the spreading behavior. Changing the viscosity and density of the salt slightly affected the spreading behavior. Setting the salt surface tension to approximately zero to simulate the salt wetting the substrate had the greatest effect on the spreading behavior and this increased both the spreading velocity and overall spreading area. Including decay heat in the model had a significant effect on the long-term cooling behavior of the salt but had no effect on the overall spreading behavior. Future effort should involve modeling a larger volume of molten salt to represent a reactor case spill scenario because the spreading behavior of large volumes of salt may differ from the small volumes that were the focus of this document.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1975

Thermal analytic model of 30 cm engineering model mercury ion thruster

A lumped parameter thermal nodal network was developed for a 30 cm Engineering Model Mercury Ion Thruster. The network consists of approximately 100 nodes coded in SINDA format for use on the Univac 1106/1108 computer. This model takes into account internal dissipation, radiation, and conduction as well as environmental heating. A series of tests were performed to simulate a wide range of thermal environments on an operating 30 cm thruster, instrumented to measure the temperature distribution within the thruster. The results of these tests were used to calibrate the analytical model. The analytical model along with comparisons between analytical and experimental results for the various operating conditions are presented.

Oglebay, J. C.↗

DOE OSTI2024

Two-Fluid and Discrete Element Modeling of a Parallel Plate Fluidized Bed Heat Exchanger for Concentrating Solar Power

A novel high-temperature particle solar receiver is developed using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a fluidized bed at the lower section, leaving the system from the bottom at a constant flowrate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built by collaborators at the Colorado School of Mines and their data are used for model validation. In this experimental setup, near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using the discrete element method and a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A new continuum conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2 K.

concentrating solar power↗

NASA NTRS1985

Preliminary Results from Numerical Experiments on the Summer 1980 Heat Wave and Drought

During the summer of 1980, a prolonged heat wave and drought affected the United States. A preliminary set of experiments has been conducted to study the effect of varying boundary conditions on the GLA model simulation of the heat wave. Five 10-day numerical integrations with three different specifications of boundary conditions were carried out: a control experiment which utilized climatological boundary conditions, an SST experiment which utilized summer 1980 sea-surface temperatures in the North Pacific, but climatological values elsewhere, and a Soil Moisture experiment which utilized the values of Mintz-Serafini for the summer, 1980. The starting dates for the five forecasts were 11 June, 7 July, 21 July, 22 August, and 6 September of 1980. These dates were specifically chosen as days when a heat wave was already established in order to investigate the effect of soil moistures or North Pacific sea-surface temperatures on the model's ability to maintain the heat wave pattern. The experiments were evaluated in terms of the heat wave index for the South Plains, North Plains, Great Plains and the entire U.S. In addition a subjective comparison of map patterns has been performed.

Wolfson, N.↗

NASA NTRS1985

Cellular interface morphologies in directional solidification. III - The effects of heat transfer and solid diffusivity

The shape and stability of two-dimensional finite-amplitude cellular interfaces arising during directional solidification are compared for several solidification models that account differently for latent heat released at the interface, unequal thermal conductivities of melt and solid, and solute diffusivity in the solid. Finite-element analysis and computer-implemented perturbation methods are used to analyze the families of steadily growing cellular forms that evolve from the planar state. In all models a secondary bifurcation between different families of finite-amplitude cells exists that halves the spatial wavelength of the stable interface. The quantitative location of this transition is very dependent on the details of the model. Large amounts of solute diffusion in the solid retard the growth of large-amplitude cells.

Ungar, Lyle H.↗

NASA NTRS2016

Optical Probe of the Superconducting Normal Mixed State in a Magnetic Penetration Thermometer

Using ultraviolet photon pulses, we have probed the internal behavior of a molybdenum-gold Magnetic Penetration Thermometer (MPT) that we designed for x-ray microcalorimetry. In this low-temperature detector, the diamagnetic response of a superconducting MoAu bilayer is used to sense temperature changes in response to absorbed photons. We have previously described an approximate model that explains the high responsivity of the detector to temperature changes as a consequence of a Meissner transition of the molybdenum-gold film in the magnetic field applied by the superconducting circuit used to bias the detector. We compare measurements of MPT heat capacity and thermal conductance, derived from UV photon pulse data, to our model predictions for the thermodynamic properties of the sensor and for the electron cooling obtained by quasiparticle recombination. Our data on electron cooling power is also relevant to the operation of other superconducting detectors, such as Microwave Kinetic Inductance Detectors.

superconductivity↗

DOE OSTI2021

Ultra-High Temperature Thermal Conductivity Measurements of a Reactive Magnesium Manganese Oxide Porous Bed Using a Transient Hot Wire Method

Pelletized magnesium manganese oxide shows promise for high temperature thermochemical energy storage. It can be thermally reduced in the temperature range between 1250 °C and 1500 °C and re-oxidized with air at typical gas-turbine inlet pressures (1–25 bar) in the temperature range between 600 °C and 1500 °C. The combined thermal and chemical volumetric energy density is approximately 2300 MJ/m3. The rate at which a thermochemical storage module can be charged is limited by heat transfer inside the solid packed bed. Hence, the effective thermal conductivity of packed beds of magnesium-manganese oxide pellets is a crucial parameter for engineering Mg-Mn-O redox storage devices. We have measured the effective thermal conductivity of a packed bed of 3.66 ± 0.516 mm sized magnesium manganese oxide (Mn to Mg molar ratio of 1:1) pellets in the temperature range of 300–1400 °C. Since the material is electrically conductive at temperatures above 600 °C, the sheathed transient hot wire method is used for measurements. Raw data is analyzed using the Blackwell solution to extract the bed thermal conductivity. The effective thermal conductivity standard deviation is less than 10% for a minimum of three repeat measurements at each temperature. Experimental results show an increase in the effective thermal conductivity with temperature from 0.50 W/m °C around 300 °C to 1.81 W/m °C close to 1400 °C. We propose a dual porosity model to express the effective thermal conductivity as a function of temperature. This model also considers the effect of radiation within the bed, as this is the dominant heat transfer mode at high temperatures. The proposed model accounts for microscale pellet porosity, macroscale bed porosity, pellet size, solid thermal conductivity (phonon transport), and radiation (photon transport). The coefficient of determination between the proposed model and the experimental results is greater than 0.90.

Engineering↗

NASA NTRS1993

A finite element analysis of the freeze/thaw behavior of external artery heat pipes

A two-dimensional finite element model was used to determine the freeze/thaw characteristics of an external artery heat pipe. During startup, the working fluid, which was located in the liquid channel and the circumferential wall grooves, experienced a phase transformation from a solid to a liquid state. The transient heat conduction equations with moving interfacial conditions were solved using the appropriate initial boundary conditions. The modelling results include the cross-sectional temperature distribution and the interfacial or melt front position as a function of time. A fixed grid approach was adopted in the model for the phase-change process during thawing of frozen working fluid. The interfacial position between the liquid and solid regions was found by balancing the latent heat caused by interfacial movement with the heat addition or extraction at the related grid points.

Lu, X. J.↗

DOE OSTI2025

Thermal 𝑅⁢𝐶 circuit model for resolving the thermal paradox

Thermal measurements of heat capacity and thermal conductivity in a wide range of insulators and superconductors exhibit a “thermal paradox”: a large linear specific heat reminiscent of neutral Fermi surfaces (associated with fractionalized quasiparticles) in nonmetallic samples that exhibit no corresponding linear temperature coefficient to the thermal conductivity. At first sight, these observations appear to support the formation of a continuum of thermally localized many-body excitations, a form of many-body localization that would be fascinating in its own right. Here, in this work, by mapping thermal conductivity measurements onto thermal 𝑅⁢𝐶 circuits, we argue that the development of extremely long thermal relaxation times, a “thermal bottleneck,” is likely in systems with either many-body localization or neutral Fermi surfaces due to the large ratio between the electron and phonon specific heat capacities. We present a reevaluation of thermal conductivity measurements in materials exhibiting a thermal paradox that can be used in future experiments to deliberate between these two exciting alternatives.

heat transfer↗

DOE OSTI2026

Concept of Dynamic Heat Insulation for Rotating Detonation Engines

This work introduces a new class of materials concept to dynamically reduce instantaneous heat fluxes in Rotating Detonation Engine (RDE) combustor chamber walls. The high-frequency and high-amplitude surface heat fluxes observed in RDEs arise from large instantaneous temperature differences between the detonation shockwave and chamber wall surface. These temperature gradients drive substantial energy loss and reduce the chamber gas pressure, ultimately limiting the cycle’s thermodynamic efficiency. This work introduces a concept for dynamically insulating the combustion chamber surfaces using surface layers or coatings with low thermal time scale. With such coatings, the surface temperature may follow the fluctuations of the cyclic detonation wave temperature, thus reducing the instantaneous heat flux therefore cycle-mean heat flux. To analyze these cyclic thermal phenomena, a one-dimensional analytical conduction solver was utilized with the capability to handle multilayered structures. Parametric modeling was performed using transient heat flux boundary conditions representative of a hydrogen–air RDE across a broad range of coating thermal properties and engine conditions. The coating effectiveness scaled with the product of thermal time constant and detonation wave frequency and the results were non-dimensionalized to guide future materials development. This strategy may offer benefits in increasing material survivability, reducing cooling requirements, and enhancing pressure gain.

heat transfer↗

DOE OSTI2022

Dominant Energy Carrier Transitions and Thermal Anisotropy in Epitaxial Iridium Thin Films

High aspect ratio metal nanostructures are commonly found in a broad range of applications such as electronic compute structures and sensing. The self-heating and elevated temperatures in these structures, however, pose a significant bottleneck to both the reliability and clock frequencies of modern electronic devices. Any notable progress in energy efficiency and speed requires fundamental and tunable thermal transport mechanisms in nanostructured metals. Here, in this work, time-domain thermoreflectance is used to expose cross-plane quasi-ballistic transport in epitaxially grown metallic Ir(001) interposed between Al and MgO(001). Thermal conductivities ranges from roughly 65 (96 in-plane) to 119 (122 in-plane) W m -1 K -1 for 25.5–133.0 nm films, respectively. Further, low defects afforded by epitaxial growth are suspected to allow the observation of electron–phonon coupling effects in sub-20 nm metals with traditionally electron-mediated thermal transport. Via combined electro-thermal measurements and phenomenological modeling, the transition is revealed between three modes of cross-plane heat conduction across different thicknesses and an interplay among them: electron dominant, phonon dominant, and electron–phonon energy conversion dominant. The results substantiate unexplored modes of heat transport in nanostructured metals, the insights of which can be used to develop electro-thermal solutions for a host of modern microelectronic devices and sensing structures.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Modeling of Additively Manufactured Ceramic Heat Exchangers with Semi-Elliptical Cross-Section Flow Channels

An approximate and easily applied analytical model was developed for heat transfer calculations of heat exchangers consisting of multiple rows and columns of heat transfer fluid flow channels with semi-elliptical cross sections. Heat exchangers of this type are being developed by using ceramic material and additive manufacturing for high temperature and pressure-concentrating solar electric power plants. Calculations using the model require only the geometrical dimensions and flow conditions of the heat exchanger. Comparisons of modeling predictions, with both simulation results and experimental data, were conducted to verify the viability of the model. The results showed good agreement where almost all the modeling predictions were within 20% of the simulation results or the experimental data. Finally, the proposed modeling approach is more generally applicable to heat transfer analysis of heat exchangers with similar flow channel configurations to those considered in this study.

Ceramic↗

NASA NTRS1992

TOPEX Microwave Radiometer - Thermal design verification test and analytical model validation

The testing of the TOPEX Microwave Radiometer (TMR) is described in terms of hardware development based on the modeling and thermal vacuum testing conducted. The TMR and the vacuum-test facility are described, and the thermal verification test includes a hot steady-state segment, a cold steady-state segment, and a cold survival mode segment totalling 65 hours. A graphic description is given of the test history which is related temperature tracking, and two multinode TMR test-chamber models are compared to the test results. Large discrepancies between the test data and the model predictions are attributed to contact conductance, effective emittance from the multilayer insulation, and heat leaks related to deviations from the flight configuration. The TMR thermal testing/modeling effort is shown to provide technical corrections for the procedure outlined, and the need for validating predictive models is underscored.

Lin, Edward I.↗

NASA NTRS1995

A Network Model for the Effective Thermal Conductivity of Rigid Fibrous Refractory Insulations

A procedure is described for computing the effective thermal conductivity of a rigid fibrous refractory insulation. The insulation is modeled as a 3-dimensional Cartesian network of thermal conductance. The values and volume distributions of the conductance are assigned to reflect the physical properties of the insulation, its constituent fibers, and any permeating gas. The effective thermal conductivity is computed by considering the simultaneous energy transport by solid conduction, gas conduction and radiation through a cubic volume of model insulation; thus the coupling between heat transfer modes is retained (within the simplifications inherent to the model), rather than suppressed by treating these heat transfer modes as independent. The model takes into account insulation composition, density and fiber anisotropy, as well as the geometric and material properties of the constituent fibers. A relatively good agreement, between calculated and experimentally derived thermal conductivity values, is obtained for a variety of rigid fibrous insulations.

Marschall, Jochen↗

DOE OSTI2026

Design and performance assessment of a dual-mode latent thermal storage integrated heat pump across multiple climate zones

The electrification of space cooling and heating systems risks overloading the existing electrical grid during peak hours. Heat pumps integrated with thermal energy storage (HP-TES) offer a promising solution by shifting peak loads to off-peak hours, reducing grid strains. This work presents the design and assessment of a dual-mode HP-TES that uses a single 22°C phase-change material (PCM) as a heat source or sink in heating and cooling modes, respectively. Performance assessment was conducted using Modelica-based transient models, and full-year simulations were conducted across eight US climate zones using typical weather data. Two HP-TES performance metrics: peak energy reduction and recharge energy increases, were defined by comparing HP-TES energy consumption with base HP peak consumptions. Annual heating demand reductions (40–65%) exceeded cooling demand reductions (15–20%) due to the elimination of peak-hour backup heating. Cooling mode recharge energy requirements (0–10%) were reduced in locations with lower summer nighttime temperatures, while heating recharge requirements were lower than the high base-peak energy demand from backup heating. Simplified rapid methods that were 10 7 times faster than annual simulations were developed to predict seasonal cooling and heating energy reductions, and recharge energy requirements, with maximum deviations of ±2.5% and ±3.5% points. These methods enabled rapid parametric studies that identified optimal location-specific PCM temperatures between 22°C and 27°C, highlighting the need to consider both discharge and recharge energy requirements to achieve sustainable, energy-efficient peak load shifting across various climate zones.

25 ENERGY STORAGE↗