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 433 records · Page 24

Physical effects on compact high-velocity clouds in the circumgalactic medium

ABSTRACT We numerically investigate the evolution of compact high-velocity clouds (CHVCs) passing through a hot, tenuous gas representing the highly ionized circumgalactic medium (CGM) by applying the adaptive-mesh refinement code flash. The model clouds start from both hydrostatic and thermal equilibrium and are in pressure balance with the CGM. Here, we present 14 models, divided into two mass categories and two metallicities each and different velocities. We allow for self-gravity and thermal conduction or not. All models experience mass diffusion, radiative cooling, and external heating leading to dissociation and ionization. Our main findings are (1) self-gravity stabilizes clouds against Rayleigh–Taylor instability, which is disrupted within 10 sound-crossing times without; (2) clouds can develop Jeans-instable regions internally even though they are initially below Jeans mass; (3) all clouds lose mass by ram pressure and Kelvin–Helmholtz instability; (4) thermal conduction substantially lowers mass-loss rates, by this, extending the clouds’ lifetimes, particularly, more than doubling the lifetime of low-mass clouds; (5) thermal conduction leads to continuous, filamentary stripping, while the removed gas is heated up quickly and mixes efficiently with the ambient CGM; (6) without thermal conduction the removed gas consists of dense, cool, clumpy fragments; (7) thermal conduction might prevent CHVCs from forming stars; and (8) clouds decelerated by means of drag from the ambient CGM form head-tail shapes and collapse after they reach velocities characteristic for intermediate-velocity clouds. Conclusively, only sophisticated modelling of CHVCs as non-homogeneous and non-isothermal clouds with thermal conduction and self-gravity explains observed morphologies and naturally leads to the suppression of star formation.

Sander, Bastian↗

High-Temperature Linear Receiver Enabled by Multicomponent Aerogels

Concentrating solar thermal (CST) technology has significant potential mainly due to its dispatchability and low cost of storage. However, to compete with other sources, including utility-scale solar PV, its final cost (cents/kWh) still needs to be lowered. Cost reduction can be achieved by improving the system level efficiency of the CST plants through the deployment of advanced power cycles, which operate at high temperatures of ~700°C. However, optical, and thermal losses pose a major challenge to the efficiency of such CST systems. The overall aim of this project is to investigate and de-risk a linear solar receiver concept called an Aerogel Insulated Receiver (AIR) that generates high temperatures (up to 700°C) at a low solar concentration ratio (<100) and with a high collection efficiency (optical × receiver). Our prior work has demonstrated the thermal stability1 and optical and heat-insulating properties2,3 of transparent aerogel insulation at a one-inch scale. The focus of this work is on (1) co-optimization of the geometry of the aerogel tile and receiver enclosure to fit a standard parabolic collector (PTC), (2) scale-up of aerogels into 4-inch tiles while preserving key properties, (3) experimental measurement of receiver heat loss (W/m) in a >70-cm test stand and validation of anticipated receiver performance at high temperatures. Regarding (1), appropriate optical and thermal models for a parabolic trough receiver (PTR) are developed and validated. The geometry of the aerogels and the receiver enclosure are co-optimized to maximize the collection efficiency. The model predicts a 54% collection efficiency at 700°C for an AIR design based on flat aerogels. By combining the collection efficiency with the power block efficiency of supercritical CO 2 cycles, we predict >10% improvements in peak plant efficiency relative to existing line-focusing CST systems. The application of curved plasmonic aerogels is predicted to further increase the collection efficiency to 64%. Regarding (2), we demonstrate the successful development of 6-inch-long refractory aerogel tiles with optical, thermal, and stability characteristics consistent with our prior work. This scale-up requires a transition to a larger ALD station and modifying the ALD process variables such as exposure time and the number of precursor doses. Regarding (3), we design and develop an AIR test stand measuring 3 feet in length. Heat loss performance analysis is carried out using the test stand. The results indicate that aerogel insulation can significantly reduce receiver thermal losses at the high operating temperatures required for next-generation PTRs. The experimental results agree with the heat loss performance predicted by our receiver model. Lastly, we conducted preliminary failure mode and effects (FMEA) and techno-economic (TEA) analyses to identify failure mitigation strategies and commercial opportunities, respectively. Overall, this project identifies key opportunities and challenges in deploying aerogel insulating receivers in next-generation line-focusing CST technologies.

14 SOLAR ENERGY↗

Spark Channel Dynamics of Electrostatic Discharges (Rev. 1) [Slides]

To model spark discharges, this work implemented, verified, and validated: 1-D radial conservative Lagrangian hydrodynamic scheme (with heat conduction); Eddington/P1 radiative transfer approximation, time-implicit scheme; RLC circuit solver, time-implicit scheme (adaptable to other circuit types). Novel analytics EOS fit for air up to temperatures of ~150,000 K. Novel "seed" electron concept to initialize finite conductivity in spark without pre-heating the air in the channel, as typically used in spark modeling literature. Reasonable, if conservative (from an engineering viewpoint), agreement with experimental data for indirect (circuit) and direct (interferometry) measurements of spark. Many potential avenues of future work: Two-temperature dryodynamics (requires two-temperature EOS development); 2-D axisymmetric geometry for axial (z) variation, modeling delectric electrodes; Kinetic equations for time-dependent charged species populations.

42 ENGINEERING↗

Dew Point Simulation and Analysis of the Heat Melt Compactor Gas Effluent

The effluent gases and water vapor of the Trash Compaction and Processing System (TCPS) can be vented through the Vacuum Exhaust System (VES) if requirements are met. One of the constraints in this venting scenario is the dew point of the exhaust gases disposed into the VES must be less than 15.5 °C. With simulations using Aspen Plus and the effluent gas sample results of the Heat Melt Compactor (HMC), flash calculations were conducted in the modeling study to calculate the feed temperature and dew point at fixed pressures. Saturated vapor curves were produced and provided a preliminary result on optimal feed conditions that satisfy the dew point and vapor-phase only requirements after flashing across the throttle valve towards the VES.

trash↗

Damage Accumulation and Failure of Plasma-Sprayed Thermal Barrier Coatings under Thermal Gradient Cyclic Conditions

Thermal barrier coatings will be more aggressively designed to protect gas turbine engine hot-section components in order to meet future engine higher fuel efficiency and lower emission goals. A fundamental understanding of the sintering and thermal cycling induced delamination of thermal barrier coating systems under engine-like heat flux conditions will potentially help to improve the coating temperature capability. In this study, a test approach is established to emphasize the real-time monitoring and assessment of the coating thermal conductivity, which can initially increase under the steady-state high temperature thermal gradient test due to coating sintering, and later decrease under the thermal gradient cyclic test due to coating cracking and delamination. Thermal conductivity prediction models have been established for a ZrO2-(7- 8wt%)Y2O3 model coating system in terms of heat flux, time, and testing temperatures. The coating delamination accumulation is then assessed based on the observed thermal conductivity response under the combined steady-state and cyclic thermal gradient tests. The coating thermal gradient cycling associated delaminations and failure mechanisms under simulated engine heat-flux conditions will be discussed in conjunction with the coating sintering and fracture testing results.

Zhu, Dongming↗

Alfvenic pulses in the solar atmosphere

Some nonlinear aspects of Alfvenic pulses propagating in coronal loops and the underlying chromosphere are numerically investigated. Heat conduction and radiation are included. The Alfvenic pulses are modeled as axisymmetric twists on a vertical cylindrical flux tube. They nonlinearly couple into acoustic-gravity waves propagating along the flux tube. A single Alfvenic pulse is found to leave two acoustic-gravity pulses in its wake. These pulses can result in significant motions of the transition region and underlying chromosphere. These motions do not resemble spicules, but they may correspond to a variety of observations indicating that the solar atmosphere is in a continual dynamic state. It is suggested that a dynamic chromosphere and transition region may be the inevitable consequence of the coronal heating process itself.

Mariska, J. T.↗

Entry heat transfer tests of the 0.006-scale space shuttle orbiter model (50-0) in Langley Research Center freon tunnel at Mach 6 (OH45)

Results are presented of heat transfer tests of a 147B configuration orbiter model (50-0) conducted in the NASA Langley Research Center Freon Tunnel (LRC/CF4). These tests were conducted at a nominal Mach number of 6, and at Reynolds numbers of 0.3 and 0.5 x 1,000,000 per foot. The objectives of the tests were to determine the effects of the low freon specific heat ratio, gamma, on the heating distributions and to determine the impingement of the orbiter bow shock on the wing. The data presented include thin skin heat transfer data (tabulated data and plotted data).

Foust, J. W.↗

Passive heat removal in horizontally oriented micro-HTGRs

There are novel, horizontally oriented high-temperature gas-cooled micro-reactor (HTGR) designs under consideration because of their portability, simplicity, and reliable operation. These features of the micro-HTGR design can meet the fission battery attributes, such as economic, standardize, installed, unattended and reliable. Although there have been HTGR related safety studies, micro-HTGRs when horizontally oriented are expected to exhibit significantly different thermal physics. During the unavailability of helium circulation, the internal reactor core is designed to cool by block-to-block conduction and radiation, and the reactor vessel surface is cooled by the ambient air. This scenario is anticipated during the transport of the micro-HTGR in a shipping container. The conduction and radiation between the prismatic micro-HTGR blocks in the core can be influenced by variances in the thermal contacts. This work investigated the conduction within a simulated horizontal HTGR core. An experimental setup was used to validate a numerical model. The experimental setup consisted of a hexagonal assembly with scaled prismatic blocks placed within a high-temperature vacuum environment. The gaps between the blocks were well controlled and monitored. The experimental setup was designed in such a way that the temperature variation in the axial direction was minimal, such that the experiment could be observed as a 2D (r, ) heat transfer problem. The pressurized tube was equipped with Infra-Red (IR) transparent windows, and an IR camera was used to capture the spatio-temporal temperature evolution in the hexagonal block assembly during the cooldown experiments. The experimental scenario was computationally modeled with a finite element analysis (FEA) program. Once validated, the computational model was used to investigate the impact of gap conductance on overall decay heat removal. Using a conservative estimate for gap conductance value (100 W/m 2 -K) between the prismatic blocks, there is a negligible increase in temperature observed during decay heat generation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermal ripples in a resistive and radiative instability

The development of the resistive tearing instability in the case of sheared magnetic fields is considered, taking into account also the occurrence of a radiatively driven thermal instability. It is pointed out that thermal conduction has generally been neglected in theories similar to those discussed. The present investigation is concerned with a consideration of both parallel and perpendicular thermal conduction, in addition to finite resistivity and radiative loss. Attention is given to the equations and the model, the spatial singularity which arises with consideration of only the parallel heat conduction, the removal of this singularity and the formation of temperature oscillations (thermal ripples) by inclusion of the perpendicular heat-flux component, and details regarding the numerical procedure. A brief explanation is provided of the conditions required for the oscillations, and potential implications of the results with respect to the solar flare are discussed.

Steinolfson, R. S.↗

Steady-state temperature distribution in living tissue modeled as cylindrical shells.

Closed form, analytical solutions to the problem of steady-state heat transfer in living tissue modeled as cylindrical shells are presented and discussed. These solutions are particularly useful for the study of temperature distributions in the extremities. Metabolic heat generation, conduction, and heat transported by the blood perfusing the tissue are considered in the model. The results demonstrate the important role that the blood stream plays in the transfer of heat inside living tissue. Solutions are also presented for the limiting cases of diminishing blood flow that would occur during vasoconstriction or occlusion of blood by external means.

Shitzer, A.↗

An analysis of parameters affecting ampacity in aircraft bipolar MVDC power cables via coupled electrical, thermal, and computational fluid dynamic modelling

Abstract The next generation of aircraft, including more electric aircraft and all‐electric aircraft (AEA), requires electric power systems with high power density and low system mass specifications. Increasing the voltage of the system to the range of a few kV, medium voltage (MV), is a reasonable approach to achieving high‐power‐density and low‐system‐mass EPSs for aircraft applications. Higher voltages, however, pose many challenges for aviation MV power cables such as arcs and arc tracking, partial discharges (PDs), and thermal management. In this regard, thermal management is more challenging since heat transfer by convection is greatly reduced at wide‐body aircraft's cruising altitudes due to the reduced air pressure. In this paper, a finite element method (FEM) model is developed in COMSOL Multiphysics for an aircraft bipolar MVDC (±5 kV) power cable. Using the model, the maximum permissible cable current at a low pressure of 18.8 kPa (at an altitude of 12.2 km from sea level, the usual cruising altitude for wide‐body aircraft) is calculated. Also, an analytical model is developed based on analytical and proven empirical correlations governing conductive, radiative, and convective heat transfers at the steady state to estimate the ampacity of the bipolar cable system at reduced pressure. It was shown that the proposed analytical model can be used for atmospheric pressure and systems with a larger number of poles, expanding its range of applications. The results of the FEM and analytical models correlate at wide ranges of parameters such as ambient temperature, duct size, distance between the positive and negative pole cables, and the overall diameter of the cables. The influence of horizontal and vertical arrangement of poles is included in the analytical model. The results of this study can be used to design bipolar MVDC power cable systems for the envisaged wide‐body AEA.

Azizi, Arian↗

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Lunar temperature and global heat flux from laboratory electrical conductivity and lunar magnetometer data

Three-layer monotonic electrical conductivity models for the lunar interior to a depth of 600 km are used in conjunction with laboratory measurements of the electrical conductivity of olivine and pyroxene to estimate a temperature-depth profile. The temperatures calculated for depths of 400-600 km are consistent with attenuation of the seismic shear wave. The temperature calculated at a depth of 100-250 km yields a heat flow that is in good agreement with the directly measured lunar heat flow. The temperature, however, is sufficiently close to melting that mascon anisostasy would not be maintained. Thus a better conductor is required at this depth.

Sonett, C. P.↗

Hybrid Modeling Study on Grain Evolution in the Metal Welding Process and Its Potential Lunar Application

Metal is most commonly used structural material in a wide range of spacecraft, and welding is the principal method for joining metal components into functional systems. However, conducting welding experiments under extreme environments—such as microgravity or vacuum conditions in space—is prohibitively expensive and experimentally challenging. To overcome these limitations, multi-physics computational welding models provide a cost-effective and versatile alternative. In this work, the authors have developed a coupled thermal (fluid) microstructure simulation framework to model metal welding under varying gravity conditions. The framework integrates a mixed-mode heat transfer formulation (conduction, convection, and radiation) with molten pool fluid dynamics, enabling accurate prediction of temperature fields and weld-pool geometry. A grain growth model is further incorporated to capture the spatial and temporal evolution of microstructure, including grain size distribution and morphological transitions during solidification. This approach provides detailed insight into molten pool evolution and grain-level microstructure development throughout the welding process. By explicitly parameterizing environmental conditions, the model supports extrapolation to off-Earth manufacturing scenarios such as welding on the lunar surface. Tantalum—chosen in this study due to its high melting point, oxidation resistance, and mechanical stability at elevated temperatures—serves as the material system for model demonstration. Beyond Tantalum, the integrated multi-physics framework offers broad applicability for predictive welding simulations of various structural and refractory metals or alloys used in extreme terrestrial or extraterrestrial environments.

kinetic Monte Carlo (SPPARKS)↗

MHD Simulations of the Plasma Flow in the Magnetic Nozzle

The magnetohydrodynamic (MHD) flow of plasma through a magnetic nozzle is simulated by solving the governing equations for the plasma flow in the presence of an static magnetic field representing the applied nozzle. This work will numerically investigate the flow and behavior of the plasma as the inlet plasma conditions and magnetic nozzle field strength are varied. The MHD simulations are useful for addressing issues such as plasma detachment and to can be used to gain insight into the physical processes present in plasma flows found in thrusters that use magnetic nozzles. In the model, the MHD equations for a plasma, with separate temperatures calculated for the electrons and ions, are integrated over a finite cell volume with flux through each face computed for each of the conserved variables (mass, momentum, magnetic flux, energy) [1]. Stokes theorem is used to convert the area integrals over the faces of each cell into line integrals around the boundaries of each face. The state of the plasma is described using models of the ionization level, ratio of specific heats, thermal conductivity, and plasma resistivity. Anisotropies in current conduction due to Hall effect are included, and the system is closed using a real-gas equation of state to describe the relationship between the plasma density, temperature, and pressure.A separate magnetostatic solver is used to calculate the applied magnetic field, which is assumed constant for these calculations. The total magnetic field is obtained through superposition of the solution for the applied magnetic field and the self-consistently computed induced magnetic fields that arise as the flowing plasma reacts to the presence of the applied field. A solution for the applied magnetic field is represented in Fig. 1 (from Ref. [2]), exhibiting the classic converging-diverging field pattern. Previous research was able to demonstrate effects such as back-emf at a super-Alfvenic flow, which significantly alters the shape of the magnetic field in both the near- and far-field regions. However, in that work the downstream domain was constrained to a channel of constant cross-sectional area. In the present work we seek to address this issue by modeling the downstream region with a domain that permits free expansion of the plasma, permitting a better evaluation of the downstream effects the applied field has on the plasma. The inlet boundary conditions and applied magnetic field values will also be varied to determine the effect the initial plasma energy content and applied magnetic field energy density have on the near- and far-field plasma properties on the MHD code. This will determine the effect of inlet boundary conditions on the results downstream and address issues related to the restrictive numerical domain previously used.

Smith, T. E. R.↗

Mercury - Internal structure and thermal evolution

Mercury's thermal evolution and internal structure are modeled based on the planet's gross physical properties (which imply a high metallic iron content) and predictions for its chemistry made from the Lewis-Cameron model of condensation of the primitive solar nebula (which implies that Mercury may be composed only of those materials that condensed at temperatures near that of metallic iron condensation in the cooling nebula). Various heat sources, initial temperatures, and thermal conductivities are considered for a homogeneous model and a differentiated two-layer model. Density distributions are calculated from the mean density and estimates of the present-day temperature. The moment of inertia and the hydrostatic value of the second degree harmonic coefficient of Mercury's gravity field are found for the differentiated and undifferentiated models. These results should be useful for preliminary interpretation of the Mariner 10 measurements of Mercury's gravitational field.

Siegfried, R. W., II↗

A mathematical model of the middle and high latitude ionosphere

A time-dependent three-dimensional model of the middle and high latitude ionosphere is described. The density distributions of six ion species NO(+), O(2+), N(2+), O(+), N(+), He(+), and the electron and ion temperatures are obtained from a numerical solution of the appropriate continuity, momentum, and energy equations. The equations are solved as a function of height for an inclined magnetic field at E and F region altitudes. The three-dimensional nature of the model is obtained by following flux tubes of plasma as they convect or corotate through a moving neutral atmosphere. The model takes account of field-aligned diffusion, cross-field electrodynamic drifts, thermospheric winds, polar wind escape, energy-dependent chemical reactions, neutral composition changes, ion production due to solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow and local heating and cooling processes. The model also takes account of the offset between the geomagnetic and geographic poles.

Schunk, R. W.↗

Wind-US Code Physical Modeling Improvements to Complement Hypersonic Testing and Evaluation

This report gives an overview of physical modeling enhancements to the Wind-US flow solver which were made to improve the capabilities for simulation of hypersonic flows and the reliability of computations to complement hypersonic testing. The improvements include advanced turbulence models, a bypass transition model, a conjugate (or closely coupled to vehicle structure) conduction-convection heat transfer capability, and an upgraded high-speed combustion solver. A Mach 5 shock-wave boundary layer interaction problem is used to investigate the benefits of k- s and k-w based explicit algebraic stress turbulence models relative to linear two-equation models. The bypass transition model is validated using data from experiments for incompressible boundary layers and a Mach 7.9 cone flow. The conjugate heat transfer method is validated for a test case involving reacting H2-O2 rocket exhaust over cooled calorimeter panels. A dual-mode scramjet configuration is investigated using both a simplified 1-step kinetics mechanism and an 8-step mechanism. Additionally, variations in the turbulent Prandtl and Schmidt numbers are considered for this scramjet configuration.

Georgiadis, Nicholas J.↗