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,063 records · Page 59

A PC-based collapse factor analysis for high pressure cryogenic systems

A collapse factor simulation computer program, suitable for use on personal computers, for supercritical cryogen system pressurization is examined. The available literature is previewed, individual component models are developed, and the overall simulation program logic is examined. Extensive use is made of existing thermodynamic and transport property evaluation packages. Heat transfer is estimated using natural convection correlations, and transient wall heat conduction is considered. Examples of the simulation are demonstrated for two facilities. Partial validation is achieved by comparing the simulation output with data from an existing test facility. Preliminary results indicate that the collapse factors will seldom exceed 1.15 during supercritical operation of high-pressure cryogenic systems. Recommendations for improvements in the component models and the simulation are made.

Hodge, B. K.↗

Emission Spectroscopy and Radiometric Measurements in the NASA Ames IHF Arc Jet Facility

Plasma diagnostic measurement campaigns in the NASA Ames Interaction Heating Facility (IHF) have been conducted over the last several years with a view towards characterizing the flow in the arc jet facility by providing data necessary for modeling and simulation. Optical emission spectroscopy has been used in the plenum and in the free jet of the nozzle. Radiation incident over a probe surface has also been measured using radiometry. Plenum measurements have shown distinct radial profiles of temperature over a range of operating conditions. For cases where large amounts of cold air are added radially to the main arc-heated stream, the temperature profiles are higher by as much as 1500 K than the profiles assumed in flow simulations. Optical measurements perpendicular to the flow direction in the free jet showed significant contributions to the molecule emission through inverse pre-dissociation, thus allowing determination of atom number densities from molecular emission. This has been preliminarily demonstrated with the N2 1st Positive System. Despite the use of older rate coefficients, the resulting atom densities are reasonable and surprisingly close to flow predictions.

Winter, Michael W.↗

An Experimental Kinetics Study of Isopropanol Pyrolysis and Oxidation behind Reflected Shock Waves

Isopropanol has potential as a future bio-derived fuel and is a promising substitute for ethanol in gasoline blends. Even so, little has been done in terms of high-temperature chemical kinetic speciation studies of this molecule. To this end, experiments were conducted in a shock tube using simultaneous CO and H2O laser absorption measurements. Water and CO formation during isopropanol pyrolysis was also examined at temperatures between 1127 and 2162 K at an average pressure of 1.42 atm. Species profiles were collected at temperatures between 1332 and 1728 K and at an average pressure of 1.26 atm for equivalence ratios of 0.5, 1.0, and 2.0 in highly diluted mixtures of 20% helium and 79.5% argon. Species profiles were also compared to four modern C3 alcohol mechanisms, including the impact of recent rate constant measurements. The Li et al. (2019) and Saggese et al. (2021) models both best predict CO and water production under pyrolysis conditions, while the AramcoMech 3.0 and Capriolo and Konnov models better predict the oxidation experimental profiles. Additionally, previous studies have collected ignition delay time (τign) data for isopropanol but are limited to low pressures in highly dilute mixtures. Therefore, real fuel–air experiments were conducted in a heated shock tube with isopropanol for stoichiometric and lean conditions at 10 and 25 atm between 942 and 1428 K. Comparisons to previous experimental results highlight the need for real fuel–air experiments and proper interpretation of shock-tube data. The AramcoMech 3.0 model over predicts τign values, while the Li et al. model severely under predicts τign. The models by Capriolo and Konnov and Saggese et al. show good agreement with experimental τign values. A sensitivity analysis using these two models highlights the underlying chemistry for isopropanol combustion at 25 atm. Additionally, modifying the Li et al. model with a recently measured reaction rate shows improvement in the model’s ability to predict CO and water profiles during dilute oxidation. Finally, a regression analysis was performed to quantify τign results from this study.

09 BIOMASS FUELS↗

The development of a high-resolution Eulerian radiation-hydrodynamics simulation capability for laser-driven Hohlraums

Hohlraums are hollow cylindrical cavities with high-Z material walls used to convert laser energy into uniform x-ray radiation drives for inertial confinement fusion capsule implosions and high energy density physics experiments. Credible computational modeling of hohlraums requires detailed modeling and coupling of laser physics, hydrodynamics, radiation transport, heat transport, and atomic physics. We report on improvements to Los Alamos National Laboratory's xRAGE radiation-hydrodynamics code in order to enable hohlraum modeling. xRAGE's Eulerian hydrodynamics and adaptive mesh refinement make it uniquely well suited to study the impacts of multiscale features in hohlraums. In order to provide confidence in this new modeling capability, we demonstrate xRAGE's ability to produce reasonable agreement with data from several benchmark hohlraum experiments. We also use xRAGE to perform integrated simulations of a recent layered high density carbon capsule implosion on the National Ignition Facility in order to evaluate the potential impacts of the capsule support tent, mixed cell conductivity methodologies, plasma transport, and cross-beam energy transfer (XBT). We find that XBT, seeded by plasma flows in the laser entrance hole (LEH), causes a slight decrease in energy coupling to the capsule and that all of these impact the symmetry of the x-ray drive such that they have an appreciable impact on the capsule implosion shape.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled effects of conduction in the crystal and thermo-solutal convection in a rectangular inclined enclosure

To date modeling of crystal growth of optoelectronic materials using Physical Vapor Transport has been limited to the study of the fluid phase. To achieve it, the equations of coupled heat, mass and momentum transfer in the gas have to be solved. The first objective of this study is to examine the effect of heat conduction in the crystal on the fluid flow in the neighborhood of the interface. Heat transfer boundary conditions on both interfaces were modified to take into account the additional heat flux between gas and solid. It is proved that heat conduction does not affect the fluid flow. In the presence of gravity, density gradients in the fluid phase generate convection responsible for the problem of a nonplanar growth of the interface. The second objective is to study systematically under one-g the different possible flows in order to solve this problem. Depending on the parameters, a diffusive mode and three convective modes (thermal, solutal and thermo-solutal) are observed. The competition between thermal and solutal convections leads to a mathematical condition which can be used to achieve a planar growth. It is proven that, under the physical conditions chosen, this mathematical condition cannot be thermodynamically satisfied.

Mennetrier, Christophe↗

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Modeling solar flare conduction fronts. I - Homogeneous plasmas and ion-acoustic turbulence. II - Inhomogeneous plasmas and ambipolar electric fields

A one-dimensional, electrostatic, particle-in-cell simulation is used here to model the expansion of a heated electron population in a coronal loop during a solar flare and the characteristics of the associated X-ray emissions. The hot electrons expand outward from the localized region, creating an ambipolar electric field which accelerates a return current of cooler, ambient electrons. Ion-acoustic waves are generated by the return currents as proposed by Brown et al. (1979), but they play little or no role in containing energetic electrons and the conduction front proposed by Brown et al. does not form. The X-ray emission efficiency of the electrons is too low in the corona for them to be the source of hard X-ray bursts. The particle dynamics changes dramatically if the heated plasma is at low altitudes and expands upward into the more tenuous plasma at higher altitudes. Two important applications of this finding are the radio-frequency heating of the corona and the collisional heating of the chromosphere by precipitating energetic electrons. In both cases, the overlying plasma has a density that is too low to supply a balancing return current to the expanding hot electrons. As a result, an ambipolar electric field develops that tends to confine the energetic electrons behind a front that propagate outward at about the speed of sound.

Mckean, M. E.↗

Evaluation of Finite-Rate Gas/Surface Interaction Models for a Carbon Based Ablator

Two sets of finite-rate gas-surface interaction model between air and the carbon surface are studied. The first set is an engineering model with one-way chemical reactions, and the second set is a more detailed model with two-way chemical reactions. These two proposed models intend to cover the carbon surface ablation conditions including the low temperature rate-controlled oxidation, the mid-temperature diffusion-controlled oxidation, and the high temperature sublimation. The prediction of carbon surface recession is achieved by coupling a material thermal response code and a Navier-Stokes flow code. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and Ablation Program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting full Navier-Stokes equations using Data Parallel Line Relaxation method. Recession analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities with heat fluxes ranging from 45 to 1100 wcm2 are performed and compared with data for model validation. The ablating material used in these arc-jet tests is Phenolic Impregnated Carbon Ablator. Additionally, computational predictions of surface recession and shape change are in good agreement with measurement for arc-jet conditions of Small Probe Reentry Investigation for Thermal Protection System Engineering.

Ablation↗

Analysis of contact conditions and microstructure evolution in shear assisted processing and extrusion using smoothed particle hydrodynamics method

Shear assisted processing and extrusion (ShAPE) is a solid-phase processing technique that adds an additional shear force as compared with a conventional extrusion approach. Recently, ShAPE has demonstrated the capability of extruding high-performance aluminum alloy 7075 (AA7075) tubes at speeds up to 12.2 m/min without surface tearing. However, the relationship among the ShAPE processing parameters, thermomechanical conditions, contact conditions, heat generation, and microstructure evolution remains primarily empirical because an insightful understanding of the associated physics is still lacking. To help elucidate these relationships, this work proposes a thermomechanical meshfree model for the first time for ShAPE processing of AA7075 using the smoothed particle hydrodynamics (SPH) method. The meshfree model is first validated thoroughly by experimental data in terms of material flow, die face temperature, and extrusion force with various processing parameters. The validated model is then used to analyze the steady-state contact conditions and heat generation rates during ShAPE processing. Distributions of the average grain size of AA7075 being extruded are calculated using the SPH model output. The meshfree model results reveal that extrusions conducted at lower temperatures and higher strain rates yield more refined grains and possibly higher material strength, which is also consistent with the experimental observations.

36 MATERIALS SCIENCE↗

PINCH POINT ANALYSIS OF GAS COOLERS IN TRANSCRITICAL R1336mzz(Z) HIGH-TEMPERATURE HEAT PUMPS

The majority of high-temperature heat pumps (HTHPs) operate subcritical vapor compression cycles, where the maximum sink temperature is limited by the critical temperature of refrigerants. In HTHPs with a sink temperature of ≥ 150°C, transcritical cycles overcome this limitation by rejecting heat in the supercritical regions of refrigerant. Supercritical refrigerants in gas coolers have a large temperature glide and substantial nonlinear variation in specific heat capacity vs. temperature. Closely matching the refrigerant glide temperature and the heat sink temperature profile greatly impacts the heat transfer characteristic of gas coolers and the overall performance of HTHPs. This study conducts a pinch point analysis for the gas cooling process of supercritical R1336mzz(Z) in HTHPs. A finite volume element model of a supercritical R1336mzz(Z) counterflow gas cooler is built and solved. This study evaluates the influences of operating conditions, including inlet and outlet of coolants, refrigerant evaporation temperature and superheat degree, pinch temperature difference, and heat rejection pressure, on the performance of HTHPs. The optimal discharge pressures in gas coolers are proposed for industrial spray dryers in providing drying air at 170°C−210°C. The results of the present study help design gas coolers in transcritical R1336mzz(Z) HTHPs.

Wang, Pengtao↗

Radiative interactions in multi-dimensional chemically reacting flows using Monte Carlo simulations

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The amount and transfer of the emitted radiative energy in a finite volume element within a medium are considered in an exact manner. The spectral correlation between transmittances of two different segments of the same path in a medium makes the statistical relationship different from the conventional relationship, which only provides the non-correlated results for nongray methods is discussed. Validation of the Monte Carlo formulations is conducted by comparing results of this method of other solutions. In order to further establish the validity of the MCM, a relatively simple problem of radiative interactions in laminar parallel plate flows is considered. One-dimensional correlated Monte Carlo formulations are applied to investigate radiative heat transfer. The nongray Monte Carlo solutions are also obtained for the same problem and they also essentially match the available analytical solutions. the exact correlated and non-correlated Monte Carlo formulations are very complicated for multi-dimensional systems. However, by introducing the assumption of an infinitesimal volume element, the approximate correlated and non-correlated formulations are obtained which are much simpler than the exact formulations. Consideration of different problems and comparison of different solutions reveal that the approximate and exact correlated solutions agree very well, and so do the approximate and exact non-correlated solutions. However, the two non-correlated solutions have no physical meaning because they significantly differ from the correlated solutions. An accurate prediction of radiative heat transfer in any nongray and multi-dimensional system is possible by using the approximate correlated formulations. Radiative interactions are investigated in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The governing equations are based on the fully elliptic Navier-Stokes equations. Chemical reaction mechanisms were described by a finite rate chemistry model. The correlated Monte Carlo method developed earlier was employed to simulate multi-dimensional radiative heat transfer. Results obtained demonstrate that radiative effects on the flowfield are minimal but radiative effects on the wall heat transfer are significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and nozzle size on the radiative and conductive wall fluxes.

Liu, Jiwen↗

Assessing the Heat Transfer Modeling Capabilities of CFD Software for Involute-Shaped Plate Research Reactors

The ongoing efforts to convert High-Performance Research Reactors (HPRRs) using Highly Enriched Uranium (HEU) to Low-Enriched Uranium (LEU) fuel require reliable thermal–hydraulic assessments of modified core designs. The involute-shaped fuel plates used in several major HPRRs present unique modeling challenges due to their compact core geometries and high heat flux conditions. This study evaluates the capability of three commercial CFD tools, STAR-CCM+, COMSOL, and ANSYS CFX, to predict cladding-to-coolant heat transfer using Reynolds-Averaged Navier–Stokes (RANS) methods within the thermal–hydraulic regimes of involute-shaped plate reactors. Broad sensitivity analysis was conducted across a range of reactor-relevant parameters using two turbulence models (k−ϵ and k−ω SST) and different near-wall treatment strategies. The results were benchmarked against the Sieder–Tate correlation and experimental data from historic studies. The codes produced consistent results, showing good agreement with the empirical correlation of Sieder–Tate and the experimental measurements. The findings support the use of these commercial CFD codes as effective tools for assessing the thermal–hydraulic performance of involute-shaped plate HPRRs and guide future LEU core development.

CFD↗

A Method for Thermal Analysis of Spiral Bevel Gears

A modeling method for analyzing the three-dimensional thermal behavior of spiral bevel gears has been developed. The model surfaces are generated through application of differential geometry to the manufacturing process for face-milled spiral bevel gears. Contact on the gear surface is found by combining tooth contact analysis with three-dimensional Hertzian theory. The tooth contact analysis provides the principle curvatures and orientations of the two surfaces. This information is then used directly in the Hertzian analysis to find the contact size and maximum pressure. Heat generation during meshing is determined as a function of the applied load, sliding velocity, and coefficient of friction. Each of these factors change as the point of contact changes during meshing. A nonlinear finite element program was used to conduct the heat transfer analysis. This program permitted the time- and position-varying boundary conditions, found in operation, to be applied to a one-tooth model. An example model and analytical results are presented.

Handschuh, Robert F.↗

Numerical investigation of thermomechanical response of multiscale porous Ultra-High Temperature Ceramics

Recent advances in Ultra-High Temperature Ceramics (UHTC) manufacturing have permitted the development of multiscale porous UHTC microstructures. Within the target application of hypersonic vehicles, dense UHTCs are suitable for thermal protection on leading edges, whereas porous UHTCs may find a role in providing thermally insulated interfaces for temperature-sensitive interior components. Designing vehicles incorporating porous UHTCs requires a characterization of their thermomechanical properties across the full range of expected operating temperatures spanning –20 °C–2500 °C. This research represents a preliminary study in performing this characterization. Several numerical experiments are performed using a coupled thermomechanical implementation of the Material Point Method to determine the temperature dependence of effective material properties both with and without damage. Furthermore, complex time-dependent boundary conditions derived from known hypersonic flight profiles are simulated in order to probe the various couplings between deformation, damage and heat transfer. It is shown that the model reveals the importance of micro-buckling in determining effective material stiffness and thermal conductivity.

36 MATERIALS SCIENCE↗

Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+) (Technical Reference Manual V.4)

The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, methane, and propane 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 safety, hazards, and risk. HyRAM+ includes generic probabilities for equipment failures, probabilistic models for the impact of heat flux on humans and structures, and experimentally validated first-order models of release and flame physics. HyRAM+ integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impact on people and structures. HyRAM+ is developed at Sandia National Laboratories to support the 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 4.0. The most significant change for HyRAM+ version 4.0 from HyRAM version 3.1 is the incorporation of other alternative fuels, namely methane (as a proxy for natural gas) and propane into the toolkit. This change necessitated significant changes to the installable graphical user interface as well as changes to the back-end Python models. A second major change is the inclusion of physics models for the overpressure associated with the delayed ignition of an unconfined jet/plume of flammable gas.

02 PETROLEUM↗

Preliminary Pilot-scale Study and Techno-economic Analysis of Chemical Heat Pump with Conventional Nuclear Reactors

The energy economy is continually evolving, particularly in terms of primary energy sources, their conversions to useful forms such as electricity and heat, and their utilization in different sectors, in response to socio-political factors. Because nuclear source is clean and non-carbon-emitting energy source, it is crucial to consider its role in the evolution of the energy economy. A pilot-scale study was conducted for Chemical Heat Pump (ChHP) system at three different scales (25-, 100-, and 1000-kW thermal outputs) using steady state thermal model. To establish the viability of selling heat rather than electricity alone, techno-economic analysis was undertaken for advanced Small Modular Reactor (SMR), and SMR paired with varied thermal output ChHP systems was referred to as the nuclear hybrid energy system (NHES) in this study. Using the U.S. Energy Information Administration statistics, pricing for electricity and natural gas were anticipated for U.S. regions (California, Northwest, Midwest, Southwest, New England, and PJM). Based on thermal output from ChHP, advanced SMR with 100 MWth and four alternative NHES scenarios were explored, namely 50-, 10-, 5-, and 1-MWth. Net present value, payback period, discounted cash flow return, and levelized cost of energy were evaluated for all scenarios. Based on the economic analysis, selling heat to high-temperature industrial processes is more profitable compared to selling electricity only. Higher carbon taxes showed significant improvement in economic parameters for NHESs. Providing heat to high-temperature industries could be beneficial, helping to reduce the greenhouse gas emissions by slashing the fossil fuel consumption.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermal preparation of foods in space-vehicle environments

Convection is the primary heat transfer mechanism for most foods heated in an earth-based environment. In contrast, in the low-gravity environment of space flight, the primary heat transfer mechanism is conduction (or radiation in the absence of a conducting medium). Conduction heating is significantly slower and less efficient than convection heating. This fact poses a problem for food heating during space flight. A numerical model has been developed to evaluate this problem. This model simulates the food-heating process for Skylab. The model includes the effect of a thermally controlled on/off heat flux. Parametric studies using this model establish how the required heating time is affected by: the thermal diffusivity of the nutrient materials, the power level of the heater, the initial food temperatures, and the food container dimensions.

Bannerot, R. B.↗

On the performance of explicit and implicit algorithms for transient thermal analysis

The status of an effort to increase the efficiency of calculating transient temperature fields in complex aerospace vehicle structures is described. The advantages and disadvantages of explicit and implicit algorithms are discussed. A promising set of implicit algorithms, known as the GEAR package is described. Four test problems, used for evaluating and comparing various algorithms, have been selected and finite element models of the configurations are discribed. These problems include a space shuttle frame component, an insulated cylinder, a metallic panel for a thermal protection system and a model of the space shuttle orbiter wing. Calculations were carried out using the SPAR finite element program, the MITAS lumped parameter program and a special purpose finite element program incorporating the GEAR algorithms. Results generally indicate a preference for implicit over explicit algorithms for solution of transient structural heat transfer problems when the governing equations are stiff. Careful attention to modeling detail such as avoiding thin or short high-conducting elements can sometimes reduce the stiffness to the extent that explicit methods become advantageous.

Adelman, H. M.↗