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 667 records · Page 37

Heat Loads Due to Small Penetrations in Multilayer Insulation Blankets

The main penetrations (supports and piping) through multilayer insulation systems for cryogenic tanks have been previously addressed by heat flow measurements. Smaller penetrations due to fasteners and attachments are now experimentally investigated. The use of small pins or plastic garment tag fasteners to each the handling and construction of multilayer insulation (MLI) blankets goes back many years. While it has long been understood that penetrations and other discontinuities degrade the performance of the MLI blanket, quantification of this degradation has generally been lumped into gross performance multipliers (often called degradation factors or scale factors). Small penetrations contribute both solid conduction and radiation heat transfer paths through the blanket. The conduction is down the stem of the structural element itself while the radiation is through the hole formed during installation of the pin or fastener. Analytical models were developed in conjunction with MLI perforation theory and Fouriers Law. Results of the analytical models are compared to experimental testing performed on a 10 layer MLI blanket with approximately 50 small plastic pins penetrating the test specimen. The pins were installed at 76-mm spacing inches in both directions to minimize the compounding of thermal effects due to localized compression or lateral heat transfer. The testing was performed using a liquid nitrogen boil-off calorimeter (Cryostat-100) with the standard boundary temperatures of 293 K and 78 K. Results show that the added radiation through the holes is much more significant than the conduction down the fastener. The results are shown to be in agreement with radiation theory for perforated films.

Structural attachments↗

Validation of a CFD Model Predicting the Effect of High Level Lateral Acceleration Sloshing on the Heat Transfer and Pressure Drop in a Small-Scale Tank in Normal Gravity

A two-phase CFD model is developed to study the effects of sloshing with high level lateral acceleration on the heat transfer and pressure drop in a small scale tank. Computational results are compared to the data provided by a non-isothermal sloshing experiment without phase change conducted by T. Himeno et al. at the University of Tokyo and JAXA in 2011. The results of the current model are, also, compared to CFD predictions reported by Himeno et al. A step change in lateral acceleration was applied in the experiment. Different levels of lateral acceleration amplitude, varying between 0G and 0.5G, were considered. CFD results for interface movement and tank pressure are presented and compared in this paper to the experimental data for the case in which the value of lateral acceleration was set to 0.5G. The effects of initial and boundary conditions and turbulence modeling approach on the tank pressure change during sloshing are discussed in detail. The effect of conjugate heat transfer in the tank wall is also studied to show its important role in determining the tank pressure evolution. The results of the Reynolds Averaged Navier Stokes (RANS) models are compared to the results of the Large Eddy Simulation model (LES) to underscore the importance of correctly capturing the effects of turbulence for high fidelity predictions.

multiphase flow↗

Thermal performance and radio-frequency transmissivity of candidate ablation materials for S-band antenna window application on manned spacecraft

A test program was conducted in the MSC 1.5 MW arc-heated facility to evaluate the thermal performance of ablation materials having potential application as radio frequency windows. These tests were conducted for the improvement of omnidirectional antenna operating characteristics during atmospheric reentry. Since a full scale model of the Apollo command service module was available for antenna tests, this mockup was used as a basic for the tests. Test models were subjected to heating conditions simulating the nominal lunar return trajectory (AS-501) and the design trajectories, high heat load and high heating rate. RF measurements were made before and after the arc jet tests to measure attenuation effects due to the thermal degradation of the materials under consideration. The test program demonstrated that additional development is required in materials technology to achieve an ablative system with both good RF transmission characteristics and thermal-structural integrity.

Tillian, D. J.↗

A theoretical study of the global F region for June solstice, solar maximum, and low magnetic activity

A time-dependent, three-dimensional, multi-ion model of the ionospheric F region at 120-800 km altitude is presented. Account is taken of field-aligned diffusion, cross-field electrodynamic drifts in equatorial and high latitude regions, interhemispheric flow, thermospheric winds, polar wind escape, energy-dependent chemical reactions and neutral composition changes. Attention is also given to the effects of ion production by solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow, local heating and cooling processes, offsets between the geomagnetic and geographic poles, and bending of field lines near the magnetic equator. The model incorporates all phenomena described by previous models and can be applied to tracing magnetic storm and substorm disturbances from high to low latitudes on a global scale. Sample results are provided for ionospheric features during a June solstice, the solar maximum and in a period of low geomagnetic activity. The model will eventually be used to study coupled ionosphere-thermosphere activity.

Sojka, J. J.↗

Analysis of the single-blow transient testing technique for perforated plate heat exchangers

An analysis of the single-blow transient testing technique for heat exchanger surfaces is made for perforated plate and similar discontinuous surfaces. The model assumes that there is no temperature variation across each plate, and allows for axial conduction in spacers: the resulting axial fluid temperature profile is discontinuous, rather than the usual continuous profile. Numerical solutions are obtained to the resulting set of coupled first-order differential equations for a step change in inlet fluid temperature. Results are presented in tabular form, which allow the heat transfer coefficients to be calculated from test data using the maximum slope technique.

Rodriguez, Jose I.↗

The Cool Flames Experiment

A space-based experiment is currently under development to study diffusion-controlled, gas-phase, low temperature oxidation reactions, cool flames and auto-ignition in an unstirred, static reactor. At Earth's gravity (1g), natural convection due to self-heating during the course of slow reaction dominates diffusive transport and produces spatio-temporal variations in the thermal and thus species concentration profiles via the Arrhenius temperature dependence of the reaction rates. Natural convection is important in all terrestrial cool flame and auto-ignition studies, except for select low pressure, highly dilute (small temperature excess) studies in small vessels (i.e., small Rayleigh number). On Earth, natural convection occurs when the Rayleigh number (Ra) exceeds a critical value of approximately 600. Typical values of the Ra, associated with cool flames and auto-ignitions, range from 104-105 (or larger), a regime where both natural convection and conduction heat transport are important. When natural convection occurs, it alters the temperature, hydrodynamic, and species concentration fields, thus generating a multi-dimensional field that is extremely difficult, if not impossible, to be modeled analytically. This point has been emphasized recently by Kagan and co-workers who have shown that explosion limits can shift depending on the characteristic length scale associated with the natural convection. Moreover, natural convection in unstirred reactors is never "sufficiently strong to generate a spatially uniform temperature distribution throughout the reacting gas." Thus, an unstirred, nonisothermal reaction on Earth does not reduce to that generated in a mechanically, well-stirred system. Interestingly, however, thermal ignition theories and thermokinetic models neglect natural convection and assume a heat transfer correlation of the form: q=h(S/V)(T(bar) - Tw) where q is the heat loss per unit volume, h is the heat transfer coefficient, S/V is the surface to volume ratio, and (T(bar) - Tw ) is the spatially averaged temperature excess. This Newtonian form has been validated in spatially-uniform, well-stirred reactors, provided the effective heat transfer coefficient associated with the unsteady process is properly evaluated. Unfortunately, it is not a valid assumption for spatially-nonuniform temperature distributions induced by natural convection in unstirred reactors. "This is why the analysis of such a system is so difficult." Historically, the complexities associated with natural convection were perhaps recognized as early as 1938 when thermal ignition theory was first developed. In the 1955 text "Diffusion and Heat Exchange in Chemical Kinetics", Frank-Kamenetskii recognized that "the purely conductive theory can be applied at sufficiently low pressure and small dimensions of the vessel when the influence of natural convection can be disregarded." This was reiterated by Tyler in 1966 and further emphasized by Barnard and Harwood in 1974. Specifically, they state: "It is generally assumed that heat losses are purely conductive. While this may be valid for certain low pressure slow combustion regimes, it is unlikely to be true for the cool flame and ignition regimes." While this statement is true for terrestrial experiments, the purely conductive heat transport assumption is valid at microgravity (mu-g). Specifically, buoyant complexities are suppressed at mu-g and the reaction-diffusion structure associated with low temperature oxidation reactions, cool flames and auto-ignitions can be studied. Without natural convection, the system is simpler, does not require determination of the effective heat transfer coefficient, and is a testbed for analytic and numerical models that assume pure diffusive transport. In addition, mu-g experiments will provide baseline data that will improve our understanding of the effects of natural convection on Earth.

Pearlman, Howard↗

Parametric Performance of a Spiral-Artery, Liquid-Trap-Diode Heat Pipe

A series of parametric investigations are described which determine the effect of various fluid charges on the performance of a 0.635-cm-diam spiral-artery, liquid-trap diode in both the forward and reverse modes. Specific parameters such as forward- and reverse-mode conductances, shutdown times and energies, and recovery to forward-model operation, are evaluated for ethane as a working fluid in the temperature range 170 K to 220 K. Results indicate that the heat pipe will not reliably start up in the forward mode. However, startup can be initiated when preceded by a diode reversal. Also included are data which show the susceptibility of the diode to fluid charge and tilt. The optimum fluid charge was found to be 2.67 g and transport capability at this charge was in excess of 1200 W-cm at 200 K. The diode in the reverse mode exhibited a rapid shutdown (within 9 min) with a shutdown energy of 1150 J (0.32 Wh).

Williams, R. J.↗

Modeling Heat Flow In a Calorimeter Equipped With a Textured Solar Collector

Heat engines are being considered for generating electric power for minisatellite applications, particularly for those missions in high radiation threat orbits. To achieve this objective, solar energy must be collected and transported to the hot side of the heat engine. A solar collector is needed having the combined properties of high solar absorptance, low infrared emittance, and high thermal conductivity. To test candidate solar collector concepts, a simple calorimeter was designed, manufactured, and installed in a bench top vacuum chamber to measure heat flow. In addition, a finite element analysis model of the collector/calorimeter combination was made to model this heat flow. The model was tuned based on observations from the as-manufactured collector/calorimeter combination. In addition, the model was exercised to examine other collector concepts, properties, and scale up issues.

Jaworske, Donald A.↗

Validation of A CFD Model Predicting the Effect of High Level Lateral Acceleration Sloshing on the Heat Transfer and Pressure Drop in A Small-Scale Tank in Normal Gravity

A two-phase CFD model is developed to study the effects of sloshing with high level lateral acceleration on the heat transfer and pressure drop in a small scale tank. Computational results are compared to the data provided by a non-isothermal sloshing experiment without phase change conducted by T. Himeno et al. at the University of Tokyo and JAXA in 2011 [1]. The results of the current model are, also, compared to CFD predictions reported by Himeno et al. [2]. A step change in lateral acceleration was applied in the experiment. Different levels of lateral acceleration amplitude, varying between 0G and 0.5G, were considered. CFD results for interface movement and tank pressure are presented and compared in this paper to the experimental data for the case in which the value of lateral acceleration was set to 0.5G. The effects of initial and boundary conditions and turbulence modeling approach on the tank pressure change during sloshing are discussed in detail. The effect of conjugate heat transfer in the tank wall is also studied to show its important role in determining the tank pressure evolution. The results of the Reynolds Averaged Navier Stokes (RANS) models are compared to the results of the Large Eddy Simulation model (LES) to underscore the importance of correctly capturing the effects of turbulence for high fidelity predictions.

multiphase flow↗

Thermal conductivity of heterogeneous mixtures and lunar soils

The theoretical evaluation of the effective thermal conductivity of granular materials is discussed with emphasis upon the heat transport properties of lunar soil. The following types of models are compared: probabilistic, parallel isotherm, stochastic, lunar, and a model based on nonlinear heat flow system synthesis.

Vachon, R. I.↗

Phase change paint tests on Rockwell orbiter/tank and orbiter alone configurations (OH3A/OH3B)

Wind tunnel tests were conducted on scale models of the space shuttle orbiter and external tank. The tests were designed to determine the basic heating rate and interference effects on the orbiter-tank configuration and to analyze the effectiveness of the thermal protective system on the reentry vehicle. The phase change paint techniques were used to determine areodynamic heating rates. Oil flow and schlieren photographs were used for flow visualization.

Quan, M.↗

Program documentation: Surface heating rate of thin skin models (THNSKN)

Program THNSKN computes the mean heating rate at a maximum of 100 locations on the surface of thin skin transient heating rate models. Output is printed in tabular form and consists of time history tabulation of temperatures, average temperatures, heat loss without conduction correction, mean heating rate, least squares heating rate, and the percent standard error of the least squares heating rates. The input tape used is produced by the program EHTS03.

Mcbryde, J. D.↗

Hypersonic wing test structure design, analysis, and fabrication

An investigation was conducted to provide the analyses, data, and hardware required to experimentally validate the beaded panel concept and demonstrate its usefulness as a basis for design of a hypersonic research airplane (HRA). Combinations of beaded panel structure, heat shields, channel caps, and corrugated webs for ribs and spars were analyzed for the wing of a specified HRA to operate at Mach 8 with a life span of 150 flights. Detailed analyses, conducted in accordance with established design criteria, included aerodynamic heating and load predictions, transient structural thermal calculations, extensive NASTRAN computer modeling, and structural optimization. After geometry was established for the total wing, part of the wing (85 sq ft) was designed, fabricated, and assembled into a test structure to experimentally verify the structural adequacy of the beaded panel design concept.

Plank, P. P.↗

Evaluation of non-Fourier heat waves influenced by nonlinear/linear boundary effects employing an explicit architecture and controlled stabilization

The present paper is concerned with the problem of heat waves in solids, where, the heat transport due to conduction occurs as propagating thermal disturbances which are transmitted at finite but high speeds. Starting from the general heat flux model of the Jeffrey's-type, and subsequent formulations leading to the Cattaneo-type heat flux model, an evaluation of the heat transport behavior is described for models influenced by non-Fourier effects and subjected to general nonlinear/linear boundary conditions. An explicit time-integration architecture is employed which effectively provides not only accurate representations of the relaxation effects and general boundary conditions but also seeks to provide an understanding of the representative thermal behavior and heat transport mechanisms for a variety of physical situations.

Tamma, Kumar K.↗

A mapping closure for turbulent scalar mixing using a time-evolving reference field

A general mapping-closure approach for modeling scalar mixing in homogeneous turbulence is developed. This approach is different from the previous methods in that the reference field also evolves according to the same equations as the physical scalar field. The use of a time-evolving Gaussian reference field results in a model that is similar to the mapping closure model of Pope (1991), which is based on the methodology of Chen et al. (1989). Both models yield identical relationships between the scalar variance and higher-order moments, which are in good agreement with heat conduction simulation data and can be consistent with any type of epsilon(phi) evolution. The present methodology can be extended to any reference field whose behavior is known. The possibility of a beta-pdf reference field is explored. The shortcomings of the mapping closure methods are discussed, and the limit at which the mapping becomes invalid is identified.

Girimaji, Sharath S.↗

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine↗

Numerical Simulation of Combustion and Extinction of a Solid Cylinder in Low-Speed Cross Flow

The combustion and extinction behavior of a diffusion flame around a solid fuel cylinder (PMMA) in low-speed forced flow in zero gravity was studied numerically using a quasi-steady gas phase model. This model includes two-dimensional continuity, full Navier Stokes' momentum, energy, and species equations with a one-step overall chemical reaction and second-order finite-rate Arrhenius kinetics. Surface radiation and Arrhenius pyrolysis kinetics are included on the solid fuel surface description and a parameter Phi, representing the percentage of gas-phase conductive heat flux going into the solid, is introduced into the interfacial energy balance boundary condition to complete the description for the quasi-steady gas-phase system. The model was solved numerically using a body-fitted coordinate transformation and the SIMPLE algorithm. The effects of varying freestream velocity and Phi were studied. These parameters have a significant effect on the flame structure and extinction limits. Two flame modes were identified: envelope flame and wake flame. Two kinds of flammability limits were found: quenching at low-flow speeds due to radiative loss and blow-off at high flow speeds due to insufficient gas residence time. A flammability map was constructed showing the existence of maximum Phi above which the solid is not flammable at any freestream velocity.

Tien, J. S.↗