SEARCH · Search NASA
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.
Generalization of laws governing heat transfer and pressure drop during transverse flow of gases in contact with finned tube banks
No-vent fill / no-vent top-off (NVF/NVTO) is a method to enable the trans-fer of cryogenic propellants without the need to vent. Transferring propellants inmicrogravity environments is necessary to support the long-duration space missionsplanned for the coming decades. A key developmental milestone to developing thiscapability fully is the simulation and experimentation of NVF/NVTO in 1-g settings.This work seeks to advance the former by assessing the effectiveness of different heattransfer correlations used in a 1-g non-equilibrium model. Heat transfer correlationsfor natural convection and two-phase boiling heat transfer are compared against ex-perimental results for 34 different tests. Two distinct experiments, comprised of 21tests, are presented in detail here. These tests cover a variety of different injectionmethods, receiver tank geometries, receiver tank materials, cryogenic propellants, andinitial tank states. Overall for the cases tested in this work, the model was able topredict the pressure response within the receiver tank to a mean absolute percentageerror of 25.9%. The temperature response error rate was 29.5% and 24.9% for wall-gasand two-phase nodes, respectively. Compared to the worst case correlation set tested,the pressure prediction error rate represents a 31% improvement. For the cases pre-sented in this work, the mean pressure prediction error rate was 15.5%. Future workcould evaluate the effectiveness of this model against experiments conducted withdifferent propellants or different initial conditions. In its current form, however, themodel can still be used to help reduce design and testing time for 1-g experiments,enabling the quicker iteration that is necessary to meet the technological demands offuture space missions.
A 2200-w capacity spacecraft heat rejection system designed for the SEPS and utilizing heat pipe radiator panels has been investigated. The total thermal control system consists of two radiator panels connected to the heat source by variable conductance heat pipes (VCHP's). The system was designed to operate in the 223 to 333 temperature range. The radiators have an emittance of 0.88 at their operational temperature and a fin efficiency of approximately 80 percent. The radiators are thermally isolated from the SEPS and environment by multilayer insulation and thermal shields. Butane was selected as the working fluid for the VCHP because of its low freezing point (135), which is necessary to prevent diffusion freezeout of the liquid during the cold outbond missions. Helium was selected for the control gas. This paper describes the VCHP system, discusses the system design parameters and presents the results of the analyses.
We are living through a “heatshed” moment, a critical turning point in which heat must be rapidly addressed to protect health and livelihoods worldwide. The year 2024 was the hottest on record since global records began in 1850, and the past 10 years were the hottest decade in recorded history. Although extreme heat is the top weather-related killer in the United States, affecting livelihoods and economies, cities and nations dedicate more attention and resources to other climate hazards, such as hurricanes, floods, droughts, and wildfires. In this commentary, we argue that heat resilience must be rapidly advanced through heat mitigation, heat management, and heat governance efforts to ensure thriving communities in the face of an overheating world.
The elemental equation governing heat transfer in aerodynamic flows is the internal energy equation. For a boundary layer flow, a double integration of the Reynolds-averaged form of this equation provides an expression of the wall heat flux in terms of the integrated effects, over the boundary layer, of various physical processes: turbulent dissipation, mean dissipation, turbulent heat flux, etc. Recently available direct numerical simulation data for a Mach 11 cold-wall turbulent boundary layer allows a comparison of the exact contributions of these terms in the energy equation to the wall heat flux with their counterparts modeled in the Reynolds-averaged Navier-Stokes (RANS) framework. Various approximations involved in RANS, both closure models as well as approximations involved in adapting incompressible RANS models to a compressible form, are assessed through examination of the internal energy balance. There are a number of potentially problematic assumptions and terms identified through this analysis. Here, the effect of compressibility corrections of the dilatational dissipation type is explored, as is the role of the modeled turbulent dissipation, in the context of wall heat flux predictions. The results indicate several potential avenues for RANS model improvement for hypersonic cold-wall boundary-layer flows.
Analytical models were developed to model the heat transfer through high-temperature fibrous insulation used in metallic thermal protection systems on reusable launch vehicles. The optically thick approximation was used to simulate radiation heat transfer through the insulation. Different models for gaseous conduction and solid conduction in the fibers, and for combining the various modes of heat transfer into a local, volume-averaged, thermal conductivity were considered. The governing heat transfer equations were solved numerically, and effective thermal conductivities were calculated from the steady-state results. An experimental apparatus was developed to measure the apparent thermal conductivity of insulation subjected to pressures, temperatures and temperature gradients representative of re-entry conditions for launch vehicles. The apparent thermal conductivity of an alumina fiber insulation was measured at nominal densities of 24, 48 and 96 kg/cu m. Data were obtained at environmental pressures from 10(exp 4) to 760 torr, with the insulation cold side maintained at room temperature and its hot side temperature varying up to 1000 C. The experimental results were used to evaluate the analytical models. The best analytical model resulted in effective thermal conductivity predictions that were within 8% of experimental results.
The objective was to develop and implement a numerical method capable of solving the nonlinear partial differential equations governing heat, mass, and momentum transfer in a 3-D cylindrical geometry in order to examine the character of convection in an asymmetrically heated cylindrical ampoule. The details of the numerical method, including verification tests involving comparison with results obtained from other methods, are presented. The results of the study of 3-D convection in an asymmetrically heated cylinder are described.
The pursuit of higher energy density in lithium-ion batteries has made high-nickel (Ni) layered oxides leading cathode candidates for next-generation electric vehicles. However, their poor thermal stability, particularly at Ni contents ≥ 90%, increases the risk of cathode-initiated thermal runaway. Furthermore, we present a data-driven framework combining linear and nonlinear machine learning models to predict key thermal runaway descriptors from a high-throughput differential scanning calorimetry database. With cathode composition and state of charge (SOC) as input features, the ensemble model accurately predicts peak temperature, heat release, and peak heat flow. SHAP analysis identifies Ni content and SOC as the dominant factors controlling thermal runaway temperature, while SOC primarily governs heat release and peak heat flow. Al, Mg, and Mn improve thermal stability by strengthening metal–oxygen bonding and delaying structural transformation, whereas B mainly reduces heat release through surface passivation. Validation with a new cathode composition confirms accurate prediction of SOC-dependent thermal runaway behavior and critical SOC.
The DEMO-FTES project investigated the feasibility of Fracture Thermal Energy Storage (FTES) as a seasonal energy storage solution in crystalline rock formations. FTES leverages hydraulically induced fractures to exchange heat between circulating fluids and the surrounding rock mass, enabling long-term thermal energy retention due to the high specific heat and low thermal conductivity of rock. This approach has the potential to reduce heating and cooling energy demands and enhance building energy resilience. The project combined dimensional analysis, numerical modeling, laboratory experiments, and meso-scale field tests to evaluate FTES performance and advance its technology readiness level from 3 to 5. Scaling analysis identified key dimensionless parameters governing heat transfer and fluid flow, ensuring laboratory and field tests were representative of larger-scale systems. Numerical simulations using TOUGH and iTOUGH2 frameworks supported experiment design and interpretation, modeling fracture geometry, thermal-hydraulic behavior, and thermo-mechanical coupling. Laboratory tests at EPFL involved creating single and multiple fractures in 25 cm cubic samples of Gabbro and Granite under true triaxial stress.
The precise control and understanding of heat flow in heterostructures is pivotal for advancements in thermoelectric energy conversion, thermal barrier coatings, and efficient heat management in electronic and optoelectronic devices. In this study, we employ high-angular-resolution time-resolved X-ray diffraction to structurally measure thermal resistance in a laser-excited AlGaAs/GaAs semiconductor heterostructure. Our methodology offers femtometer-scale spatial sensitivity and nanosecond time resolution, enabling us to directly observe heat transport across a buried interface. We corroborate established Thermal Boundary Resistance (TBR) values for AlGaAs/GaAs heterostructures and demonstrate that TBR arises from material property discrepancies on either side of a nearly flawless atomic interface. This work not only sheds light on the fundamental mechanisms governing heat flow across buried interfaces but also presents a robust experimental framework that can be extended to other heterostructure systems, paving the way for optimized thermal management in next-generation devices.
The thermal resistance associated with the interface between a solid and a liquid is analysed from an atomistic point of view. Partial evaluation of the associated Green–Kubo integral elucidates the various factors governing heat transport across the interface and leads to a quantitative model for the thermal resistance in terms of atomistic-level system parameters. The model is validated using molecular dynamics simulations.
Grain growth in the wake of the melt pool formed during alloy-based additive manufacturing (AM) is complex and multifaceted, depending on parameters governing heat transport, fluid flow, and solidification itself. Cellular automata (CA) models have proven effective in providing computationally efficient and physically sound predictions of grain structure for several AM problems, but their efficiency is tied to the performance of heat transport models. CA models use only a small portion of the problem's temperature data (near the moving melt pool boundary), and much of the CA calculations do not affect the final result due to re-melting of material. Coupling of and communication between heat transport and solidification models, and eliminating operations irrelevant towards final grain structure prediction, will be necessary for using these methods for efficient simulation of large parts. In this work, we introduce a procedure of decoupling the CA from temperature field simulation, using files of relevant temperature data written by the heat transport model. This approach is validated against the standard coupling approach using data obtained through the computational fluid dynamics software OpenFOAM. Negligible differences are seen in grain size, volume, and texture distributions for multilayer simulation of test problems, while the quantity of temperature data for these test problems was reduced by four orders of magnitude (from 100s of GB to 10s of MB) and the code performance sped up by a factor of around 50. Furthermore, variability in microstructure as a function of cell size, substrate, time step, and nucleation parameters is studied, and it is found that cell sizes less than or equal to 1.67 μm and sufficiently small time steps yield Αstatistically equivalent microstructures. Finally, a potential use case for this CA approach—the layer-wise convergence in grain structure starting from extremes in initial grain size—is examined. This approach's ability to simulate expected trends in nucleation and epitaxial grain growth for large regions of microstructure, simulated independently of heat transport models themselves, should prove useful for investigation of various microstructure uncertainties and prediction of part-scale experimental results.
Lignocellulosic biomass has a complex, species-specific microstructure that governs heat and mass transport during conversion processes. A quantitative understanding of the evolution of pore size and structure is critical to optimize conversion processes for biofuel and bio-based chemical production. Further, improving our understanding of the microstructure of biochar coproduct will accelerate development of its myriad applications. This work quantitatively compares the microstructural features and the anisotropic permeabilities of two woody feedstocks, red oak and Douglas fir, using X-ray computed tomography (XCT) before and after the feedstocks are subjected to pyrolysis. Quantitative analysis of the three-dimensional (3D) reconstructions allows for direct calculations of void fractions, pore size distributions and tortuosity factors. Next, 3D images are imported into an immersed boundary based finite volume solver to simulate gas flow through the porous structure and to directly calculate the principal permeabilities along longitudinal, radial, and tangential directions. The permeabilities of native biomass are seen to differ by three to four orders of magnitude in the different principal directions, but we find that this anisotropy is substantially reduced in the biochar formed during pyrolysis. The quantitative transport properties reported here enhance the ability of pyrolysis simulations to account for feedstock-specific effects and thereby provide a useful touchstone for the biorefining community.
Powered-lift aircraft may produce severe high-temperature environments which are potentially damaging to a landing surface or the aircraft. The interaction betweean the high temperature flow field and a nonadiabatic landing surface is analyzed with a coupled computational fluid dynamics/solid thermal conduction computer code, HOTJET. The HOTJET code couples time-accurate, implicit, factored solution schemes for the governing fluid dynamics equations (Reynolds-averaged Navier-Stokes equations) to the unsteady thermal conduction equation, which governs heat flux within a solid. HOTJET is validated against exact solutions to the thermal conduction and Navier-Stokes equations. First-of-a-kind results are included which show the impact of surface material properties on the fluid physics and the coupled fluid/material thermal fields.
The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.
Abstract General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks—where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin—commonly display standing shocks within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable to low-luminosity sources like Sgr A* and M87*, the shocks have transrelativistic speed, moderate plasma beta (the ratio of ion thermal pressure to magnetic pressure is β pi1 ∼ 1–8), and low sonic Mach number (the ratio of shock speed to sound speed is M s ∼ 1–6). We study such shocks with 2D particle-in-cell simulations, and we quantify the efficiency and mechanisms of electron heating for the special case of preshock magnetic fields perpendicular to the shock direction of propagation. We find that the postshock electron temperature T e2 exceeds the adiabatic expectation T e2,ad by an amount T e 2 / T e 2 , ad − 1 ≃ 0.0016 M s 3.6 , nearly independent of the plasma beta and of the preshock electron-to-ion temperature ratio T e1 / T i1 , which we vary from 0.1 to unity. We investigate the heating physics for M s ∼ 5–6 and find that electron superadiabatic heating is governed by magnetic pumping at T e1 / T i1 = 1, whereas heating by B -parallel electric fields (i.e., parallel to the local magnetic field) dominates at T e1 / T i1 = 0.1. Our results provide physically motivated subgrid prescriptions for electron heating at the collisionless shocks seen in GRMHD simulations of black hole accretion flows.
The technological problems in which heat transfer research is required are briefly reviewed, discussing the role of the government laboratory in this research. The Hot Section Technology Program at NASA Lewis is addressed as an example. A building-block approach to turbine aerothermal research is shown, as is an interactive computational/experimental synergism and an interdisciplinary research methodology for life prediction of turbine engine hot-section components. The main issues for government laboratories as organizers and focusers in this area of research are listed, and challenges and opportunities for the future are indicated.
Numerical schemes are employed to investigate heat transfer in the turbulent flow of He II. FEM is used to solve a set of equations governing the heat transfer and hydrodynamics of He II in the turbulent regime. Numerical results are compared with available experimental data and interpreted in terms of conventional heat transfer parameters such as the Prandtl number, the Peclet number, and the Nusselt number. Within the prescribed Reynolds number domain, the Gorter-Mellink thermal counterflow mechanism becomes less significant, and He II acts like an ordinary fluid. The convection heat transfer characteristics of He II in the highly turbulent regime can be successfully described by using the conventional turbulence and heat transfer theories.