SEARCH · Engineering Papers
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.
This report's contents focus on making use of NASA Glenn on-site computational facilities,to develop, validate, and apply models for use in advanced 3D Navier-Stokes Computational Fluid Dynamics (CFD) codes to enhance the capability to compute heat transfer and losses in turbomachiney.
Heat pipe reactors have been considered by the Space Nuclear Propulsion program for Nuclear Electric Propulsion (NEP) power conversion systems and will require the use of heat exchangers to transfer heat via heat pipes to the Brayton working fluid from the reactor. Sodium (Na) and lithium (Li) were considered as viable working fluids inside the heat pipes which were assumed to have the same geometry based on studies and information from the Los Alamos National Laboratory. The heat exchanger was assumed to be a rectangular duct with heat pipes serving as tubes from previous NEP work and recommendations. Based on this geometry, Zukauskas correlations were used to model the convective heat transfer and pressure losses. Parametric sizing of the reactor component involved operational limits-based heat pipe thermal hydraulic modeling in cohesion with required user input geometry for the in-core lattice and various subcomponents. This work considered various power conversion inlet temperatures (PCIT) of 1100 K, 1150 K, and 1200 K for Na heat pipes and 1100 K, 1150 K, 1200 K, and 1400 K for Li heat pipes based on recommendations from prior work. Using these different PCITs, the subsystem masses and pressure losses were determined and analyzed. Na showed a lower overall operating temperature and about a fifth of the maximum heat throughput capability than that of Li for the same geometry. Due to this, the entire Na-based subsystem ended up being three times more massive than the Li-based subsystem given five times the required number of heat pipes. At the low PCIT of 1100 K, the Na-based subsystem exhibited the lowest pressure losses given the large overall cross sectional flow area and relatively low frictional pressure losses. However, as the PCIT increased, the frictional pressure losses increased resulting in higher pressure losses at the 1200 K PCIT than Li-based subsystem. However, the Li-based subsystem exhibited the largest pressure losses of all analyzed cases at the 1400 K PCIT due to the low density of the Brayton working fluid at this temperature.
Capillary-pumped loops (CPL's) are devices that are used to transport heat from one location to another--specifically to transfer heat away from something. In low-gravity applications, such as satellites (and possibly the International Space Station), CPL's are used to transfer heat from electrical devices to space radiators. This is accomplished by evaporating one liquid surface on the hot side of the CPL and condensing the vapor produced onto another liquid surface on the cold side. Capillary action, the phenomenon that causes paper towels to absorb spilled liquids, is used to "pump" the liquid back to the evaporating liquid surface (hot side) to complete the "loop." CPL's require no power to operate and can transfer heat over distances as large as 30 ft or more. Their reliance upon evaporation and condensation to transfer heat makes them much more economical in terms of weight than conventional heat transfer systems. Unfortunately, they have proven to be unreliable in space operations, and the explanation for this unreliability has been elusive. The Capillary-Driven Heat Transfer (CHT) experiment is investigating the fundamental fluid physics phenomena thought to be responsible for the failure of CPL's in low-gravity operations. If the failure mechanism can be identified, then appropriate design modifications can be developed to make capillary phase-change heat-transport devices a more viable option in space applications. CHT was conducted onboard the Space Shuttle Columbia during the first Microgravity Science Laboratory (MSL-1) mission, STS-94, which flew from July 1 to 17, 1997. The CHT glovebox investigation, which was conceived by Dr. Kevin Hallinan and Jeffrey Allen of the University of Dayton, focused on studying the dynamics associated with the heating and cooling at the evaporating meniscus within a capillary phase-change device in a low-gravity environment. The CHT experimental hardware was designed by a small team of engineers from Aerospace Design & Fabrication (ADF), the NASA Lewis Research Center, and the University of Dayton. The hardware consisted of two experiment modules that each contained an instrumented test loop (idealized capillary-pumped loop), a base unit for power conversion and backlighting, a display unit with 15 LED's (light-emitting diodes) to display temperatures, pressure, heater power, and time, a control unit to select heaters and heater settings, a cooling fan, and associated cables.
A combined numerical and experimental investigation of radiation-induced convection is presented to show that the convective stability of the top-heated enclosure is disrupted by heat transfer conditions at the wall. When the enclosure is not insulated the thermal stratification of the fluid is modified by convective and radiative losses to the surrounding environment. This results in a double annular cell flow which, when cut by the laser sheet, shows a four-vortex pattern with a weak annular cell at the bottom and a large counter-rotating annular cell at the top. When the enclosure is insulated the convective stability of the fluid is again disrupted - this time as a result of radiative heat transfer between the enclosing surfaces which drives two annular flow cells of relatively equal size. Comparison between model and experiment shows that radiation effects are important even at temperature levels as low as 300 C and, if these effects are not included, numerical predictions can be highly erroneous.
A combined numerical and experimental investigation of radiation-induced convection is presented to show that the convective stability of the top-heated enclosure is disrupted by heat transfer conditions at the wall. When the enclosure is not insulated the thermal stratification of the fluid is modified by convective and radiative losses to the surrounding environment. This results in a double annular cell flow which, when cut by the laser sheet, shows a four-vortex pattern with a weak annular cell at the bottom and a large counter-rotating annular cell at the top. When the enclosure is insulated the convective stability of the fluid is again disrupted - this time as a result of radiative heat transfer between the enclosing surfaces which drives two annular flow cells of relatively equal size. Comparison between model and experiment shows that radiation effects are important even at temperature levels as low as 300 C and, if these effects are not included, numerical predictions can be highly erroneous.
Detailed heat-transfer rate distributions measured laterally over the windward surface of an orbiter-like configuration using thin-film resistance heat-transfer gauges and globally using the newly developed relative intensity, two-color thermographic phosphor technique are presented for Mach 6 and 10 in air. The angle of attack was varied from 0 to 40 deg, and the freestream Reynolds number based on the model length was varied from 4 x 10(exp 5) to 6 x 10(exp 6) at Mach 6, corresponding to laminar, transitional, and turbulent boundary layers; the Reynolds number at Mach 10 was 4 x 10(exp 5), corresponding to laminar flow. The primary objective of the present study was to provide detailed benchmark heat-transfer data for the calibration of computational fluid-dynamics codes. Predictions from a Navier-Stokes solver referred to as the Langley aerothermodynamic upwind relaxation algorithm and an approximate boundary-layer solving method known as the axisymmetric analog three-dimensional boundary layer code are compared with measurement. In general, predicted laminar heat-transfer rates are in good agreement with measurements.
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.
Here, this work aims to improve capabilities for modeling localized effects in porous media models of Pebble Bed Reactors. The wall-channeling effect is the primary local phenomenon of interest in a PBR, where the presence of the reflector wall disrupts the pebble packing, causing the pebbles near the wall to pack less efficiently and creating large void regions. Accurate modeling of the near-wall region is important as it will affect core bypass flow and temperature predictions. Porous media models are commonly used for design scoping and plant-level simulations of PBRs. Although these models have some capabilities to model the near-wall region, the correlations that are available in porous media codes are often inaccurate when a multi-region model is used to discretize the near-wall region. This work employs a high-to-low analysis to study the accuracy of available interphase heat transfer closures. NekRS, a spectral element computational fluid dynamics code, is used to perform Large Eddy Simulations. These LES simulation results are compared to porous media model results from the Pronghorn porous media code. The friction term of the KTA drag closure is first improved, reducing the error in the prediction of the near-wall velocity from over 50% to less than 5%. This is combined with improvements to the form term from previous works to produce a drag closure that is capable of accurately modeling the wall-channeling effect across a variety of flow conditions. The Nusselt number predictions of several heat transfer correlations are compared to the high-fidelity results where it is found that the KTA heat transfer correlation is capable of accurately predicting the local Nusselt numbers that were determined in the high-fidelity simulation. Comparison of the radial solid temperature profiles, however, reveal discrepancies between NekRS and Pronghorn. It is discovered that the implementation of the interphase heat transfer coefficient that exists in many current porous media codes is not valid when local porosities are modeled. Instead, it is suggested that the interphase heat transfer coefficient should be dependent on the local porosity, the Nusselt number, and the local solid surface-to-volume ratio. Implementation of this change produces improvement in the agreement between the results obtained by NekRS and Pronghorn while using the KTA heat transfer correlation.
Radiator heat transfer and geometric characteristics of Brayton space power cycle using gas-solid suspension working fluid
Research and development programs in variable conductance heat pipe technology were conducted. The treatment has been comprehensive, involving theoretical and/or experimental studies in hydrostatics, hydrodynamics, heat transfer into and out of the pipe, fluid selection, and materials compatibility, in addition to the principal subject of variable conductance control techniques. Efforts were not limited to analytical work and laboratory experimentation, but extended to the development, fabrication and test of spacecraft hardware, culminating in the successful flight of the Ames Heat Pipe Experiment on the OAO-C spacecraft.
An experimental study was conducted (1) to experimentally measure, assess and analyze the heat transfer within the internal cooling configuration of a radial turbine rotor blade and (2) to obtain heat transfer data to evaluate and improve computational fluid dynamics (CFD) procedures and turbulent transport models of internal coolant flows. A 1.15 times scale model of the coolant passages within the NASA LERC High Temperature Radial Turbine was designed, fabricated of Lucite and instrumented for transient beat transfer tests using thin film surface thermocouples and liquid crystals to indicate temperatures. Transient heat transfer tests were conducted for Reynolds numbers of one-fourth, one-half, and equal to the operating Reynolds number for the NASA Turbine. Tests were conducted for stationary and rotating conditions with rotation numbers in the range occurring in the NASA Turbine. Results from the experiments showed the heat transfer characteristics within the coolant passage were affected by rotation. In general, the heat transfer increased and decreased on the sides of the straight radial passages with rotation as previously reported from NASA-HOST-sponsored experiments. The heat transfer in the tri-passage axial flow region adjacent to the blade exit was relatively unaffected by rotation. However, the heat transfer on one surface, in the transitional region between the radial inflow passage and axial, constant radius passages, decreased to approximately 20 percent of the values without rotation. Comparisons with previous 3-D numerical studies indicated regions where the heat transfer characteristics agreed and disagreed with the present experiment.
High temperature solar thermal facilities are looking to increase operating temperatures through novel heat transfer media, one such being solid particles. These particles operating at high temperatures will require transferring their thermal energy into another working fluid like supercritical carbon dioxide which can be used in advanced power cycles. Achieving high heat transfer between the particles and supercritical carbon dioxide is essential to high efficiency and low-cost operation. Therefore, optimizing the thermal conductivity of these particles is one potential way to ensure high performance. Traditionally, unimodal particle distributions have been employed in high temperature particle solar power plants. However, ambient temperature testing of bimodal particle distributions has revealed a superior thermal conductivity when compared to its unimodal counterpart at the same temperature. This data was obtained by certified, off-the-shelf instruments that can effectively simulate the conditions a particle would be exposed to in a high temperature solar thermal system. Data obtained in this way suggests that the increased thermal conductively imputed by a bimodal particle distribution is significant at working temperatures in solar facilities. Furthermore, the thermal conductivity of these bimodal particle distributions peaks when the best combination of large and small particles is applied. At high temperatures, binary particle distributions are compared to monodispersed distributions of larger particles where heat transfer is more prolific due to the increased surface radiation. Various thermal conductivity, porosity and heat exchanger models are explored in conjunction with data acquired up to 700 C.
A decade ago several important fundamental heat transfer phenomena were identified which were considered basic to the ability to predict heat transfer loads in aircraft gas turbines. The progress in addressing these fundamentals over the past ten years is assessed. Much reseach effort has been devoted to their study in university, industry and government labs and significant progress has been achieved. Advances in computer technology have enabled the modeling of complex 3-D fluid flow in gas turbines so necessary for heat transfer calculations. Advances in instrumentation plus improved data acquisition have brought about more reliable data sets. While much has advanced in the 1980's, much challenging research remains to be done. Several of these areas are suggested.
This paper presents the effects of downstream film cooling, with and without leading edge showerhead film cooling, on turbine-vane external heat transfer. Steady-state experimental measurements were made in a three-vane linear two-dimensional cascade. The principal independent parameters were maintained over ranges consistent with actual engine conditions. The test matrix was structured to provide an assessment of the independent influence of parameters of interest, namely, exit Mach number, exit Reynolds number, coolant-to-gas temperature ratio, and coolant-to-gas pressure ratio. The data obtained indicate that considerable cooling benefits can be achieved by utilizing downstream film cooling. The downstream film cooling process was shown to be a complex interaction of two competing mechanisms. The thermal dilution effect, associated with the injection of relatively cold fluid, results in a decrease in the heat transfer to the airfoil. Conversely, the turbulence augmentation, produced by the injection process, results in increased heat transfer to the airfoil.
Investigation of a hot-film wedge probe in viscoelastic fluids by imposing on the probe a sinusoidal vibration of known amplitude and frequency. Root-mean-square (rms) velocities calculated from the displacement of the probe are compared to rms velocities obtained with a constant temperature anemometer. The tests were performed under turbulent flow conditions and also at flow rates where viscoelastic effects, i.e., decrease of heat transfer rates from the probe to the fluid and drag reduction, were observed. The experimental technique was established by measuring the response of the probe in mineral oil.
A research and development program in variable conductance heat pipe technology is reported. The project involved: (1) theoretical and/or experimental studies in hydrostatics, (2) hydrodynamics, (3) heat transfer into and out of the pipe, (4) fluid selection, and (5) materials compatibility. The development, fabrication, and test of the space hardware resulted in a successful flight of the heat pipe experiment on the OAO-3 satellite. A summary of the program is provided and a guide to the location of publications on the project is included.
The objective of this paper is to define empirical parameters (or closwre models) for an initial thermai non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two CFD codes currently being used at Glenn Research Center (GRC) for Stirling engine modeling are Fluent and CFD-ACE. The porous-media models available in each of these codes are equilibrium models, which assmne that the solid matrix and the fluid are in thermal equilibrium at each spatial location within the porous medium. This is believed to be a poor assumption for the oscillating-flow environment within Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, we non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location end time during the cycle. A NASA regenerator research grant has been providing experimental and computational results to support definition of various empirical coefficients needed in defining a noa-equilibrium, macroscopic, porous-media model (i.e., to define "closure" relations). The grant effort is being led by Cleveland State University, with subcontractor assistance from the University of Minnesota, Gedeon Associates, and Sunpower, Inc. Friction-factor and heat-transfer correlations based on data taken with the NASAlSunpower oscillating-flow test rig also provide experimentally based correlations that are useful in defining parameters for the porous-media model; these correlations are documented in Gedeon Associates' Sage Stirling-Code Manuals. These sources of experimentally based information were used to define the following terms and parameters needed in the non-equilibrium porous-media model: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity (including themal dispersion and estimate of tortuosity effects}, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity (including the effect of tortuosity) was also estimated. Determination of the porous-media model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Convertor (TDC), which uses a random-fiber regenerator matrix. The non-equilibrium porous-media model presented is considered to be an initial, or "draft," model for possible incorporation in commercial CFD codes, with the expectation that the empirical parameters will likely need to be updated once resulting Stirling CFD model regenerator and engine results have been analyzed. The emphasis of the paper is on use of available data to define empirical parameters (and closure models) needed in a thermal non-equilibrium porous-media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates. However, it is anticipated that a thermal non-equilibrium model such as that presented here, when iacorporated in the CFD codes, will improve our ability to accurately model Stirling regenerators with CFD relative to current thermal-equilibrium porous-media models.
A multiphysics analysis of a film cooling flow was conducted using a computational fluid dynamics (CFD) solver in conjunction with a conjugate heat transfer (CHT) solver. The film cooling problem consisted of a square nozzle producing a heated subsonic jet over a perforated plate with three sections of 45 cooling holes each. Results from the multiphysics analysis were compared to a CFD-only approach and experimental results. Reynolds-averaged Navier-Stokes (RANS) computations were employed first, with and without the CHT solver. The multiphysics approach provided much closer agreement with experimentally measured plate temperatures than the CFD-only approach. Comparisons of computed velocities and temperatures were also made with experimental data. Hybrid RANS-large eddy simulation calculations employing the Delayed Detached Eddy Simulation (DDES) approach were also conducted within the multiphysics framework. The DDES approach was unable to resolve turbulent structures in the near wall regions. Hence, the RANS part of the DDES model was used in near wall regions, and produced surface temperatures along with velocity and temperature profiles close to the wall that were very similar to the RANS results. The DDES approach did resolve structures in the outer portion of the boundary layer especially near the end of the plate where the boundary layer was thicker. The key result was the improved surface temperatures produced by the multiphysics approach employing CHT over the fluid-only approach.