Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat Flux”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Heat transfer to throat tubes in a square-chambered rocket engine at the NASA Lewis Research Center

A gaseous H2/O2 rocket engine was constructed at the NASA-Lewis to provide a high heat flux source representative of the heat flux to the blades in the high pressure fuel turbopump (HPFTP) during startup of the space shuttle main engines. The high heat flux source was required to evaluate the durability of thermal barrier coatings being investigated for use on these blades. The heat transfer, and specifically, the heat flux to tubes located at the throat of the test rocket engine was evaluated and compared to the heat flux to the blades in the HPFTP during engine startup. Gas temperatures, pressures and heat transfer coefficients in the test rocket engine were measured. Near surface metal temperatures below thin thermal barrier coatings were also measured at various angular orientations around the throat tube to indicate the angular dependence of the heat transfer coefficients. A finite difference model for a throat tube was developed and a thermal analysis was performed using the measured gas temperatures and the derived heat transfer coefficients to predict metal temperatures in the tube. Near surface metal temperatures of an uncoated throat tube were measured at the stagnation point and showed good agreement with temperatures predicted by the thermal model. The maximum heat flux to the throat tube was calculated and compared to that predicted for the leading edge of an HPFTP blade. It is shown that the heat flux to an uncooled throat tube is slightly greater than the heat flux to an HPFTP blade during engine startup.

Nesbitt, James A.↗
NASA NTRS2023

Universal Two-Phase Convection Heat Transfer Correlations for Cryogenic Pipe Chilldown

This paper presents a set of universal two-phase convection heat transfer correlations for modeling boiling heat transfer during pipe chilldown fit over the widest available range of cryogenic fluids and thermodynamic conditions. The correlations improve upon prior correlations that were developed separately for liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching datasets. The new correlations include equations to calculate the single-phase vapor heat transfer, film boiling heat transfer, transition boiling heat transfer, nucleate boiling heat transfer, single-phase liquid heat transfer, bulk vapor temperature during high quality film boiling, Leidenfrost temperature, critical heat flux, critical heat flux temperature, and the onset of nucleate boiling temperature. The correlations were validated against LH2, LN2, liquid methane, liquid oxygen, and liquid argon pipe quenching datasets. The eight datasets cover the following parameter ranges: pipe lengths of 0.1 to 6.5m; outer pipe diameters of 12.7 to 25.4 mm; pipe wall thicknesses from 0.51 to 1.64 mm; flow directions of upward, downward, and horizontal; gravity levels of 1g and 0g±0.01g. A numerical model assumes homogeneous mixing between the vapor and liquid and implicitly integrates the coupled energy equations for the pipe and fluid, as well as continuity for the fluid. The model was used to estimate the equilibrium quality and fluid mass accumulation along the pipes during chilldown from estimates of the fluid-pipe heat transfer extracted from temperature measurements. The correlations can be implemented in lumped parameter codes such as SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP) to improve the accuracy of chilldown time and chilldown boiloff mass predictions. Such predictions are useful for designing ground or in-space cryogenic liquid transfer systems.

Universal Correlations↗
DOE OSTI2025

Gyrokinetic profile prediction and validation of a negative triangularity plasma in ASDEX Upgrade

In this work, gyrokinetic simulations are performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73–93) for a negative triangularity H-mode plasma in ASDEX Upgrade, and compared with experimental measurements. The PORTALS framework (Rodriguez-Fernandez et al 2024 Nucl. Fusion 64 076034) is used to accelerate the prediction of kinetic profiles for this plasma, using surrogate modeling and Bayesian optimization. Ion heat flux, electron heat flux, and electron particle flux are simultaneously matched across the simulated radial regime of the plasma (normalized radius $r/a = 0.35-0.90$), and the resulting ion temperature, electron temperature, and electron density profiles match well with the experimental profile data within this radial range. A synthetic Correlation Electron Cyclotron Emission diagnostic is applied to find well-matched electron temperature fluctuation properties between simulation and experiment. The flux-matched profiles provide a basis for investigation of the turbulence nature across the plasma radius, revealing the dominance of Trapped Electron Mode turbulence at $r/a = 0.35$, the dominance of Ion Temperature Gradient turbulence at $r/a = 0.55$, 0.75, and 0.83, and an instability boundary at $r/a = 0.90$.

gyrokinetic simulation↗
NASA NTRS1987

Experimental evaluation of heat transfer on a 1030:1 area ratio rocket nozzle

A 1030:1 carbon steel, heat-sink nozzle was tested. The test conditions included a nominal chamber pressure of 2413 kN/sq m and a mixture ratio range of 2.78 to 5.49. The propellants were gaseous oxygen and gaseous hydrogen. Outer wall temperature measurements were used to calculate the inner wall temperature and the heat flux and heat rate to the nozzle at specified axial locations. The experimental heat fluxes were compared to those predicted by the Two-Dimensional Kinetics (TDK) computer model analysis program. When laminar boundary layer flow was assumed in the analysis, the predicted values were within 15% of the experimental values for the area ratios of 20 to 975. However, when turbulent boundary layer conditions were assumed, the predicted values were approximately 120% higher than the experimental values. A study was performed to determine if the conditions within the nozzle could sustain a laminar boundary layer. Using the flow properties predicted by TDK, the momentum-thickness Reynolds number was calculated, and the point of transition to turbulent flow was predicted. The predicted transition point was within 0.5 inches of the nozzle throat. Calculations of the acceleration parameter were then made to determine if the flow conditions could produce relaminarization of the boundary layer. It was determined that if the boundary layer flow was inclined to transition to turbulent, the acceleration conditions within the nozzle would tend to suppress turbulence and keep the flow laminar-like.

Kacynski, Kenneth J.↗
NASA NTRS1987

Experimental evaluation of heat transfer on a 1030:1 area ratio rocket nozzle

A 1030:1 carbon steel, heat-sink nozzle was tested. The test conditions included a nominal chamber pressure of 2413 kN/sq m and a mixture ratio range of 2.78 to 5.49. The propellants were gaseous oxygen and gaseous hydrogen. Outer wall temperature measurements were used to calculate the inner wall temperature and the heat flux and heat rate to the nozzle at specified axial locations. The experimental heat fluxes were compared to those predicted by the Two-Dimensional Kinetics (TDK) computer model analysis program. When laminar boundary layer flow was assumed in the analysis, the predicted values were within 15 percent of the experimental values for the area ratios of 20 to 975. However, when turbulent boundary layer conditions were assumed, the predicted values were approximately 120 percent higher than the experimental values. A study was performed to determine if the conditions within the nozzle could sustain a laminar boundary layer. Using the flow properties predicted by TDK, the momentum-thickness Reynolds number was calculated, and the point of transition to turbulent flow was predicted. The predicted transition point was within 0.5 inches of the nozzle throat. Calculations of the acceleration parameter were then made to determine if the flow conditions could produce relaminarization of the boundary layer. It was determined that if the boundary layer flow was inclined to transition to turbulent, the acceleration conditions within the nozzle would tend to suppress turbulence and keep the flow laminar-like.

Kacynski, Kenneth J.↗
NASA NTRS1988

Methods for heat transfer and temperature field analysis of the insulated diesel, phase 3

Work during Phase 3 of a program aimed at developing a comprehensive heat transfer and thermal analysis methodology for design analysis of insulated diesel engines is described. The overall program addresses all the key heat transfer issues: (1) spatially and time-resolved convective and radiative in-cylinder heat transfer, (2) steady-state conduction in the overall structure, and (3) cyclical and load/speed temperature transients in the engine structure. These are all accounted for in a coupled way together with cycle thermodynamics. This methodology was developed during Phases 1 and 2. During Phase 3, an experimental program was carried out to obtain data on heat transfer under cooled and insulated engine conditions and also to generate a database to validate the developed methodology. A single cylinder Cummins diesel engine was instrumented for instantaneous total heat flux and heat radiation measurements. Data were acquired over a wide range of operating conditions in two engine configurations. One was a cooled baseline. The other included ceramic coated components (0.050 inches plasma sprayed zirconia)-piston, head and valves. The experiments showed that the insulated engine has a smaller heat flux than the cooled one. The model predictions were found to be in very good agreement with the data.

Morel, Thomas↗
NASA NTRS1993

Verification of thermal analysis codes for modeling solid rocket nozzles

One of the objectives of the Solid Propulsion Integrity Program (SPIP) at Marshall Space Flight Center (MSFC) is development of thermal analysis codes capable of accurately predicting the temperature field, pore pressure field and the surface recession experienced by decomposing polymers which are used as thermal barriers in solid rocket nozzles. The objective of this study is to provide means for verifications of thermal analysis codes developed for modeling of flow and heat transfer in solid rocket nozzles. In order to meet the stated objective, a test facility was designed and constructed for measurement of the transient temperature field in a sample composite subjected to a constant heat flux boundary condition. The heating was provided via a steel thin-foil with a thickness of 0.025 mm. The designed electrical circuit can provide a heating rate of 1800 W. The heater was sandwiched between two identical samples, and thus ensure equal power distribution between them. The samples were fitted with Type K thermocouples, and the exact location of the thermocouples were determined via X-rays. The experiments were modeled via a one-dimensional code (UT1D) as a conduction and phase change heat transfer process. Since the pyrolysis gas flow was in the direction normal to the heat flow, the numerical model could not account for the convection cooling effect of the pyrolysis gas flow. Therefore, the predicted values in the decomposition zone are considered to be an upper estimate of the temperature. From the analysis of the experimental and the numerical results the following are concluded: (1) The virgin and char specific heat data for FM 5055 as reported by SoRI can not be used to obtain any reasonable agreement between the measured temperatures and the predictions. However, use of virgin and char specific heat data given in Acurex report produced good agreement for most of the measured temperatures. (2) Constant heat flux heating process can produce a much higher heating rate than the radiative heating process. The results show that heating rates of 125 C/s (225 F/s) can be achieved. (3) The electrical resistance of the FM 5055 samples were about 150 Omega in the virgin state, and decreased to about 35 Omega when charred. A reliable scheme must be developed to electrically insulate the composite from the heater foil in order to prevent any current leakage into the sample which can result in volumetric heating of the char zone of the composite.

Keyhani, M.↗
DOE OSTI2023

Experimental investigation on nucleate boiling heat transfer of low GWP refrigerants over metal-foam enhanced tube bundles

The objective of this work is to develop an enhanced tube bundle to improve the nucleate boiling heat transfer for the application of flooded evaporators. The enhanced tube bundle was made by deploying metal foam around the outer surface of aluminum tubes, which provide higher heat transfer area and nucleation site density. The pool boiling experiments for horizontal tube bundles were conducted in a pressure vessel at difference heat fluxes. The heat transfer performance of bare aluminum tube bundle was compared against the metal-foam enhanced tube bundle using low GWP refrigerant, R-1234yf. The experimental results showed that the heat transfer coefficient of the metal-foam enhanced tube bundle was at least two times higher than that of the bare tube bundle, and the enhancement ratio depended on the heat flux condition.

Yang, Cheng-Min↗
NASA NTRS2024

Microgravity Flow Boiling Experiments With Liquid-Vapor Mixture Inlet Onboard the International Space Station

This article is part of a series of studies culminating from the multi-objective Flow Boiling and Condensation Experiment (FBCE) onboard the International Space Station, which utilized the Flow Boiling Module (FBM) for experiments between February and July 2022. This study investigates microgravity flow boiling of n-Perfluorohexane with liquid-vapor mixture (two-phase) inlet conditions to the FBM with either one or two opposite walls heated. The FBM’s channel has a rectangular cross-sectional area of 5.0 × 2.5 mm 2 and a heated length of 114.6 mm. Key parameters of interest include mass velocity (180 – 2400 kg/m 2 s), inlet quality (-0.01 – 0.87), inlet pressure (120 – 200 kPa), and heat flux (1.8 W/cm 2 to critical heat flux), and a large database is amassed. The flow is visualized via a high-speed video camera and photographs are recorded at each heating increment to assess the periodic flow patterns within the channel and the near-wall interfacial behavior. Flow patterns are complex and mainly characterized by high- and low-density fronts alternately traversing the channel to yield high- and low-density-dominant periods of boiling. At all operating conditions, high-density fronts are faster during high-density-dominant periods. At low inlet qualities, the flow is annular near the channel inlet with a central vapor core surrounded by an annular liquid layer. Each high-density front having a high liquid fraction leaves a thin liquid layer sheared onto the heated walls. Boiling occurs within the liquid layer and a vapor layer is formed next to the heated wall. Inlet quality and mass velocity most dictate the overall flow patterns followed by heating configuration, and to a much lesser extent, heat flux and inlet pressure. Heat transfer characteristics are assessed via averaged boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and parametric curves of local and averaged heat transfer coefficients. Inlet pressure has an insignificant effect on heat transfer. At similar operating conditions, both the heating configurations yield similar trends and values of heat transfer coefficient and critical heat flux (CHF, slightly higher for single-sided) even though double-sided heating adds twice the heat to the fluid and doubly raises local quality. The heat fluxes required for both onset of nucleate boiling degradation and CHF are larger at high mass velocities and low inlet qualities. For a fixed inlet quality, high mass velocities yield higher average heat transfer coefficients at both lower and higher heat fluxes, while the nucleate boiling regime at intermediate heat fluxes is unaffected. For a fixed mass velocity, higher inlet qualities yield higher and lower average heat transfer coefficients at lower and higher heat fluxes, respectively.

Microgravity↗
DOE OSTI2025

Validation and moisture content sensitivity analysis of cross-laminated timber wall assemblies in EnergyPlus

Cross-laminated timber buildings are becoming more common in North America, with many numerical studies showing potential energy savings. However, no studies have validated any EnergyPlus heat transfer algorithms or quantified their accuracy in simulating CLT in building envelopes. This study empirically validates the heat flux predictions for each of EnergyPlus's heat transfer algorithms (Conduction Transfer Functions (CTF), Effective Moisture Penetration Depth (EMPD), Conduction Finite Difference (CondFD), and Heat and Moisture Transfer (HAMT)) for two different CLT ply thicknesses with both summer and winter boundary conditions measured in controlled lab experiments. It also evaluates the model sensitivity of heat flux and heating and cooling loads to moisture content. The 1D validation shows that the HAMT model is the most accurate among all algorithms. All EnergyPlus's heat flux predictions are accurate independent of CLT plate thickness for summer conditions. However, the three constant property algorithms (CTF, EMPD, and CondFD) underpredict heat flux throughout the whole day during winter conditions. The 1D sensitivity analysis indicates that elevated moisture content can increase peak heat fluxes through the material by up to 20 %. Finally, the whole building model sensitivity analysis shows increased heating load and slight cooling load variation due to increased moisture content when using constant property models. The analysis shows significantly lower peak thermal demand (7 % lower heating and 6 % lower cooling) and monthly thermal load (8 % less cooling and 6 % less heating) predictions when using HAMT vs a constant property model.

42 ENGINEERING↗
NASA NTRS1994

Heat transfer measurements for Stirling machine cylinders

The primary purpose of this study was to measure the effects of inflow-produced heat turbulence on heat transfer in Stirling machine cylinders. A secondary purpose was to provide new experimental information on heat transfer in gas springs without inflow. The apparatus for the experiment consisted of a varying-volume piston-cylinder space connected to a fixed volume space by an orifice. The orifice size could be varied to adjust the level of inflow-produced turbulence, or the orifice plate could be removed completely so as to merge the two spaces into a single gas spring space. Speed, cycle mean pressure, overall volume ratio, and varying volume space clearance ratio could also be adjusted. Volume, pressure in both spaces, and local heat flux at two locations were measured. The pressure and volume measurements were used to calculate area averaged heat flux, heat transfer hysteresis loss, and other heat transfer-related effects. Experiments in the one space arrangement extended the range of previous gas spring tests to lower volume ratio and higher nondimensional speed. The tests corroborated previous results and showed that analytic models for heat transfer and loss based on volume ratio approaching 1 were valid for volume ratios ranging from 1 to 2, a range covering most gas springs in Stirling machines. Data from experiments in the two space arrangement were first analyzed based on lumping the two spaces together and examining total loss and averaged heat transfer as a function of overall nondimensional parameter. Heat transfer and loss were found to be significantly increased by inflow-produced turbulence. These increases could be modeled by appropriate adjustment of empirical coefficients in an existing semi-analytic model. An attempt was made to use an inverse, parameter optimization procedure to find the heat transfer in each of the two spaces. This procedure was successful in retrieving this information from simulated pressure-volume data with artificially generated noise, but it failed with the actual experimental data. This is evidence that the models used in the parameter optimization procedure (and to generate the simulated data) were not correct. Data from the surface heat flux sensors indicated that the primary shortcoming of these models was that they assumed turbulence levels to be constant over the cycle. Sensor data in the varying volume space showed a large increase in heat flux, probably due to turbulence, during the expansion stroke.

Kornhauser, Alan A.↗
DOE OSTI2021

Energetics of surface melt in West Antarctica

We use reanalysis data and satellite remote sensing of cloud properties to examine how meteorological conditions alter the surface energy balance to cause surface melt that is detectable in satellite passive microwave imagery over West Antarctica. This analysis can detect each of the three primary mechanisms for inducing surface melt at a specific location: thermal blanketing involving sensible heat flux and/or longwave heating by optically thick cloud cover, all-wave radiative enhancement by optically thin cloud cover, and föhn winds. We examine case studies over Pine Island and Thwaites glaciers, which are of interest for ice shelf and ice sheet stability, and over Siple Dome, which is more readily accessible for field work. During January 2015 over Siple Dome we identified a melt event whose origin is an all-wave radiative enhancement by optically thin clouds. During December 2011 over Pine Island and Thwaites glaciers, we identified a melt event caused mainly by thermal blanketing from optically thick clouds. Over Siple Dome, those same 2011 synoptic conditions yielded a thermal-blanketing-driven melt event that was initiated by an impulse of sensible heat flux and then prolonged by cloud longwave heating. The December 2011 synoptic conditions also generated föhn winds at a location on the Ross Ice Shelf adjacent to the Transantarctic Mountains, and we analyze this case with additional support from automatic weather station data. In contrast, a late-summer thermal blanketing period over Pine Island and Thwaites glaciers during February 2013 showed surface melt initiated by cloud longwave heating and then prolonged by enhanced sensible heat flux. One limitation thus far with this type of analysis involves uncertainties in the cloud optical properties. Nevertheless, with improvements this type of analysis can enable quantitative prediction of atmospheric stress on the vulnerable Antarctic ice shelves in a steadily warming climate.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Regional surface fluxes from satellite-derived surface temperatures (AVHRR) and radiosonde profiles

Radiometric surface temperatures, derived from measurements by the AVHRR instrument aboard the NOAA-9 and the NOAA-11 polar orbiting satellites, were used in combination with wind velocity and temperature profiles measured by radiosondes, to calculate surface fluxes of sensible heat. The measurements were made during FIFE, the First ISLSCP (International Satellite Land Surface Climatology Project) Field Experiment, in a hilly tall grass prairie area of northeastern Kansas. The method of calculation was based on turbulent similarity formulations for the atmospheric boundary layer. Good agreement (r = 0.7) was obtained with reference values of sensible heat flux, taken as arithmetric means of measurements with the Bowen ratio method at six ground stations. The values of evaporation (latent heat fluxes), derived from these sensible heat fluxes by means of the energy budget, were also in good agreement (r = 0.94) with the corresponding reference values from the ground stations.

Brutsaert, Wilfried↗
NASA NTRS2008

Simulated Reentry Heating by Torching

The two first order reentry heating parameters are peak heating flux (W/cm2) and peak heat load (kJ/cm2). Peak heating flux (and deceleration, gs) is higher for a ballistic reentry and peak heat load is higher for a lifting reentry. Manned vehicle reentries are generally lifting reentries at nominal 1-5 gs so that personnel will not be crushed by high deceleration force. A few off-nominal manned reentries have experienced 8 or more gs with corresponding high heating flux (but below nominal heat load). The Shuttle Orbiter reentries provide about an order of magnitude difference in peak heating flux at mid-bottom (TPS tiles, approximately 6 W/cm2 or 5 BTU/ft2- sec) and leading edge (RCC, approximately 60 W/cm2 or 50 BTU/ft2- sec). Orion lunar return and Mars sample lander are of the same order of magnitude as orbiter leading edge peak heat loads. Flight temperature measurements are available for some orbiter TPS tile and RCC locations. Return-to-Flight on-orbit tile-repair-candidate-material-heating performance was evaluated by matching propane torch heating of candidate-materials temperatures at several depths to orbiter TPS tile flight-temperatures. Char and ash characteristics, heat expansion, and temperature histories at several depths of the cure-in-place ablator were some of the TPS repair material performance characteristics measured. The final char surface was above the initial surface for the primary candidate (silicone based) material, in contrast to a receded surface for the Apollo-type ablative heat shield material. Candidate TPS materials for Orion CEV (LEO and lunar return), and for Mars sample lander are now being evaluated. Torching of a candidate ablator material, PICA, was performed to match the ablation experienced by the STARDUST PICA heat shield. Torching showed that the carbon fiberform skeleton in a sample of PICA was inhomogeneous in that sample, and allowed measurements (of the clumps and voids) of the inhomogeneity. Additional reentry heating-performance characterizations of high temperature insulation materials were performed.

Harvey, Gale A.↗
NASA NTRS1992

Measuring surface fluxes in CAPE

Two stations (site 1612 and site 2008) were operated by the University of Georgia group from 6 July 1991 to 18 August 1991. The following data were collected continuously: surface energy fluxes (i.e., net radiation, soil heat fluxes, sensible heat flux and latent heat flux), air temperature, vapor pressure, soil temperature (at 1 cm depth), and precipitation. Canopy reflectance and light interception data were taken three times at each site between 6 July and 18 August. Soil moisture content was measured twice at each site.

Kanemasu, E. T.↗
DOE OSTI2024

Q1 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q1: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $/circ$ Initially run CESOL with SOLPS and map heat flux using simple HEAT method: ▪ Generate SOLPS grid for CAT reference case, ▪ Implement and apply ‘lower’ fidelity method of using HEAT-like analytic method to map SOLPS heat flux from charged particles to first wall, and ▪ Provide initial heat flux + neutron loading for further engineering analysis. 2. Demonstrate multiphysics magnet analysis: $/circ$ Demonstrate validation of the Elmer workflow by comparing induced stresses due to EM forces generated by TF coils using multiple codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Use of ground-based remotely sensed data for surface energy balance calculations during Monsoon '90

Surface energy balance was evaluated at a semiarid watershed using direct and indirect measurements of the turbulent fluxes, a remote technique based on measurements of surface reflectance and temperature, and conventional meteorological information. Comparison of remote estimates of net radiant flux and soil heat flux densities with measured values showed errors on the order of +/-40 W/sq m. To account for the effects of sparse vegetation, semi-empirical adjustments to aerodynamic resistance were required for evaluation of sensible heat flux density (H). However, a significant scatter in estimated versus measured latent heat flux density (LE) was still observed, +/-75 W/sq m over a range from 100-400 W/sq m. The errors of H and LE estimates were reduced to +/-50 W/sq m when observations were restricted to clear sky conditions.

Moran, M. S.↗