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.
An extensive investigation of the changes of the local and mean heat-transfer coefficients along the pipe length, and some results.
Since 2012, researchers at the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and NASA Glenn Research Center have been collaborating on a long-term effort to study flow boiling and condensation in microgravity. The ultimate goal has been to develop the Flow Boiling and Condensation Experiment (FBCE) for the International Space Station (ISS). Based on the findings from prior flow boiling experiments both at different orientations in Earth gravity and onboard parabolic flights simulating short durations of microgravity, a final refined experiment design, construction, and operating procedure have been arrived at for long-duration microgravity flow boiling experiments onboard the ISS. This study investigates flow boiling of n-Perfluorohexane with subcooled inlet in a rectangular channel of dimensions 114.6 mm heated length, 2.5 mm width, and 5 mm height. These pre-launch experiments (Mission Sequence Testing) were conducted in vertical upflow orientation in Earth gravity using the same experimental rig that was launched to the ISS in August 2021. The various operating parameters varied are heating configuration (single- and double-sided), mass velocity (180 – 3200 kg/m2s), inlet subcooling (+0 – 32°C, encompassing both highly subcooled and near-saturated inlet conditions), and inlet pressure (119 – 191 kPa). High-speed video flow visualization images are presented to explain the two-phase interfacial physics within the channel’s heated section. Heat transfer results in terms of flow boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and averaged heat transfer coefficients are analyzed and parametric effects elucidated. Severe temporal thermodynamic equilibrium is observed for near-saturated inlet at very low velocities. Nucleate boiling degradation starts at larger heat fluxes for single-sided heating than double sided at low mass velocities with highly subcooled inlet, and conversely at high mass velocities with near-saturated inlet. Nucleate boiling degradation can be delayed to higher heat fluxes by highly subcooling the inlet and increasing mass velocity. The entire local heat transfer coefficient profiles are degraded at higher heat fluxes for near-saturated inlet, but only the downstream part for highly subcooled inlet. This study also confirmed reliability of the upcoming ISS experimental data for subcooled inlet conditions and the collected Earth-gravity data will be used for comparison against the ISS data.
Theoretical and experimental approaches to exhaust plume analysis are compared. A two-phase model is extended to include treatment of reacting gas chemistry, and thermodynamical modeling of the gaseous phase of the flow field is considered. The applicability of empirical data currently available to define particle drag coefficients, heat transfer coefficients, mean particle size, and particle size distributions is investigated. Experimental and analytical comparisons are presented for subscale solid rocket motors operating at three altitudes with attention to pitot total pressure and stagnation point heating rate measurements. The mathematical treatment input requirements are explained. The two-phase flow field solution adequately predicts gasdynamic properties in the inviscid portion of two-phase exhaust plumes. It is found that prediction of exhaust plume gas pressures requires an adequate model of flow field dynamics.
The calculation of solid propellant exhaust plume flow fields is addressed. Two major areas covered are: (1) the applicability of empirical data currently available to define particle drag coefficients, heat transfer coefficients, mean particle size and particle size distributions, and (2) thermochemical modeling of the gaseous phase of the flow field. Comparisons of experimentally measured and analytically predicted data are made. The experimental data were obtained for subscale solid propellant motors with aluminum loadings of 2, 10 and 15%. Analytical predictions were made using a fully coupled two-phase numerical solution. Data comparisons will be presented for radial distributions at plume axial stations of 5, 12, 16 and 20 diameters.
Real-time measurements of ice growth during artificial and natural icing conditions were conducted using an ultrasonic pulse-echo technique. This technique allows ice thickness to be measured with an accuracy of + or - 0.5 mm; in addition, the ultrasonic signal characteristics may be used to detect the presence of liquid on the ice surface and hence discern wet and dry ice growth behavior. Ice growth was measured on the stagnation line of a cylinder exposed to artificial icing conditions in the NASA Lewis Icing Research Tunnel (IRT), and similarly for a cylinder exposed in flight to natural icing conditions. Ice thickness was observed to increase approximately linearly with exposure time during the initial icing period. The ice accretion rate was found to vary with cloud temperature during wet ice growth, and liquid runback from the stagnation region was inferred. A steady-state energy balance model for the icing surface was used to compare heat transfer characteristics for IRT and natural icing conditions. Ultrasonic measurements of wet and dry ice growth observed in the IRT and in flight were compared with icing regimes predicted by a series of heat transfer coefficients. The heat transfer magnitude was generally inferred to be higher for the IRT than for the natural icing conditions encountered in flight. An apparent variation in the heat transfer magnitude was also observed for flights conducted through different natural icing-cloud formations.
An analysis of combined heat and mass transfer from a flat plate has been made in terms of Prandtl t s simplified physical concept of the turbulent boundary layer. The results of the analysis show that for conditions of reasonably small heat and mass transfer, the ratio of the mass-and heat-transfer coefficients is dependent on the Reynolds number of the boundary layer, the Prandtl number of the medium of diffusion, and the Schmidt number of the diffusing fluid in the medium of diffusion. For the particular case of water evaporating into air, the ratio of mass-transfer coefficient to heat-transfer coefficient is found to be slightly greater than unity.
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.
This study investigates flow boiling of n-Perfluorohexane with saturated two-phase mixture inlet in a rectangular channel of dimensions 114.6-mm heated length, 2.5-mm width, and 5 mm-height. The experiments were performed as part of the Mission Sequence Testing of the Flow Boiling and Condensation Experiment’s (FBCE) Flow Boiling Module (FBM) in the vertical upflow configuration in Earth gravity using the same experimental system that was launched to the International Space Station (ISS) in August 2021. The operating parameters varied are heating configuration (single- and double-sided), mass velocity (380 – 2400 kg/m2s), inlet quality (0.011 – 0.519), and inlet pressure (120 – 179 kPa). High-speed video photographs are presented to explain the two-phase flow patterns within the channel’s heated length. Flow patterns are constituted by low-density and high-density fronts moving along the channel, with the high-density fronts gradually reducing in length due to evaporation. Heat transfer results in terms of flow boiling curves, streamwise wall temperature profiles, streamwise heat transfer coefficient profiles, and average heat transfer coefficients are presented and trends discussed. CHF data from the present experiments are combined with prior databases to compile a consolidated FBCE-CHF database for saturated inlet to expand the ranges of operating conditions and include other flow orientations in Earth gravity. Experimental CHF trends are also discussed. The interfacial lift-off model shows a good CHF predictive accuracy evidenced by a mean absolute error of 11.97% for this consolidated database after constraining it to mass velocities greater than or equal to 500 kg/m2s. Finally, this study confirmed reliability of the upcoming ISS experiments for saturated inlet conditions and the collected Earth-gravity data will be compared to ISS microgravity data.
A numerical modeling is to examine the effects of coolant convective heat transfer coefficient and frictional heating on the local temperature characteristics of a ball element in Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump (HPOTP) bearing. The present modeling uses a control-volume based, finite-difference method to solve the non-dimensionalized heat conduction equation in spherical coordinate system. The dimensionless temperature is found as a function of Biot number, heat flux ratio between the two race contacts, and location in the ball. The current results show that, for a given cooling capability, the ball temperature generally increases almost linearly with the heat input from the race-contacts. This increase is always very high at one of the two contacts. An increase in heat transfer coefficient generally reduces the ball temperature and alleviates the temperature gradient, except for the regions very close to the race contacts. For a 10-fold increase of heat transfer coefficient, temperature decrease is 35 percent for the average over entire ball, and 10 percent at the inner-race contact. The corresponding change of temperature gradient displays opposing trends between the regions immediately adjacent to the contacts and the remaining portion of the ball. The average temperature gradient in the vicinity of both contacts increases approximately 70 to 100 percent. A higher temperature gradient produces excessive thermal stress locally which may be detrimental to the material integrity. This, however, is the only unfavorable issue for an increase of heat transfer coefficient.
A heat transfer device consisting of a heated rotating cylinder in a bath was analyzed for its effectiveness to determine heat transfer coefficient of fluids. A time dependent analysis shows that the performance is insensitive to the value of heat transfer coefficient with the given rig configuration.
The temperature distribution in spur gears operating in a state of thermal equilibrium is solved by using a finite element method. The effects of various dimensionless parameters on bulk temperature are shown. A table is provided which can be used to predict the bulk temperature on gear teeth, once the heat transfer coefficients and frictional heat input is estimated. Theoretical results for estimating heat transfer coefficients and frictional heat are also summarized.
An evaporative surface is described for heat pipes and other two-phase heat transfer applications that consists of a hybrid composition of V-grooves and capillary wicking. Characteristics of the surface include both a high heat transfer coefficient and high heat flux capability relative to conventional open faced screw thread surfaces. With a groove density of 12.6/cm and ammonia working fluid, heat transfer coefficients in the range of 1 to 2 W/sq cm K have been measured, along with maximum heat flux densities in excess of 20 W/sq cm. A peak heat transfer coefficient in excess of 2.3 W/sq cm K at 20 W/sq cm was measured with a 37.8/cm hybrid surface.
An evaporative surface is described for heat pipes and other two-phase heat transfer applications that consists of a hybrid composition of V-grooves and capillary wicking. Characteristics of the surface include both a high heat transfer coefficient and high heat flux capability relative to conventional open-faced screw thread surfaces. With a groove density of 12.6 cm/1 and ammonia working fluid, heat transfer coefficients in the range of 1 to 2 W/sq cm have been measured along with maximum heat flux densities in excess of 20 W/sq cm. A peak heat transfer coefficient in excess of 2.3 W/sq cm was measured with a 37.8 cm/1 hybrid surface.
Measurements of the velocity characteristics of the flows in two curved diffusers of rectangular cross section with C and S-shaped centerlines are presented and related to measurements of wall heat transfer coefficients along the heated flat walls of the ducts. The velocity results were obtained by laser-Doppler anemometry in a water tunnel and the heat transfer results by liquid crystal thermography in a wind tunnel. The thermographic technique allowed the rapid and inexpensive measurement of wall heat transfer coefficients along flat walls of arbitrary boundary shapes with an accuracy of about 5 percent. The results show that an increase in secondary flow velocities near the heated wall causes an increase in the local wall heat transfer coefficient, and quantify the variation for maximum secondary-flow velocities in a range from 1.5 to 17 percent of the bulk flow velocity.
This report presents results from a joint study by DRA, NASA, and ONERA for the purpose of comparing, improving, and validating the aircraft icing computer codes developed by each agency. These codes are of three kinds: (1) water droplet trajectory prediction, (2) ice accretion modeling, and (3) transient electrothermal deicer analysis. In this joint study, the agencies compared their code predictions with each other and with experimental results. These comparison exercises were published in three technical reports, each with joint authorship. DRA published and had first authorship of Part 1 - Droplet Trajectory Calculations, NASA of Part 2 - Ice Accretion Prediction, and ONERA of Part 3 - Electrothermal Deicer Analysis. The results cover work done during the period from August 1986 to late 1991. As a result, all of the information in this report is dated. Where necessary, current information is provided to show the direction of current research. In this present report on ice accretion, each agency predicted ice shapes on two dimensional airfoils under icing conditions for which experimental ice shapes were available. In general, all three codes did a reasonable job of predicting the measured ice shapes. For any given experimental condition, one of the three codes predicted the general ice features (i.e., shape, impingement limits, mass of ice) somewhat better than did the other two. However, no single code consistently did better than the other two over the full range of conditions examined, which included rime, mixed, and glaze ice conditions. In several of the cases, DRA showed that the user's knowledge of icing can significantly improve the accuracy of the code prediction. Rime ice predictions were reasonably accurate and consistent among the codes, because droplets freeze on impact and the freezing model is simple. Glaze ice predictions were less accurate and less consistent among the codes, because the freezing model is more complex and is critically dependent upon unsubstantiated heat transfer and surface roughness models. Thus, heat transfer prediction methods used in the codes became the subject for a separate study in this report to compare predicted heat transfer coefficients with a limited experimental database of heat transfer coefficients for cylinders with simulated glaze and rime ice shapes. The codes did a good job of predicting heat transfer coefficients near the stagnation region of the ice shapes. But in the region of the ice horns, all three codes predicted heat transfer coefficients considerably higher than the measured values. An important conclusion of this study is that further research is needed to understand the finer detail of of the glaze ice accretion process and to develop improved glaze ice accretion models.
Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multipass, heat transfer model with both radially inward and outward flow. Trip strips on the leading and trailing surfaces of the radial coolant passages were used to produce the rough walls. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages: coolant-to-wall temperature ratio, Rossby number, Reynolds number, and radius-to-passage hydraulic diameter ratio. The first three of these four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. Results were correlated and compared to previous results from stationary and rotating similar models with trip strips. The heat transfer coefficients on surfaces, where the heat increased with rotation and buoyancy, varied by as much as a factor of four. Maximum values of the heat transfer coefficients with high rotation were only slightly above the highest levels obtained with the smooth wall model. The heat transfer coefficients on surfaces, where the heat transfer decreased with rotation, varied by as much as a factor of three due to rotation and buoyancy. It was concluded that both Coriolis and buoyancy effects must be considered in turbine blade cooling designs with trip strips and that the effects of rotation were markedly different depending upon the flow direction.
Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multipass, heat transfer model with both radially inward and outward flow. Trip strips on the leading and trailing surfaces of the radial coolant passages were used to produce the rough walls. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages: coolant-to-wall temperature ratio, Rossby number, Reynolds number, and radius-to-passage hydraulic diameter ratio. The first three of these four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. Results were correlated and compared to previous results from stationary and rotating similar models with trip strips. The heat transfer coefficients on surfaces, where the heat increased with rotation and buoyancy, varied by as much as a factor of four. Maximum values of the heat transfer coefficients with high rotation were only slightly above the highest levels obtained with the smooth wall model. The heat transfer coefficients on surfaces, where the heat transfer decreased with rotation, varied by as much as a factor of three due to rotation and buoyancy. It was concluded that both Coriolis and buoyancy effects must be considered in turbine blade cooling designs with trip strips and that the effects of rotation were markedly different depending upon the flow direction.
The effect of the upstream wake on the blade heat transfer has been numerically examined. The geometry and the flow conditions of the first stage turbine blade of GE s E3 engine with a tip clearance equal to 2 percent of the span was utilized. Based on numerical calculations of the vane, a set of wake boundary conditions were approximated, which were subsequently imposed upon the downstream blade. This set consisted of the momentum and thermal wakes as well as the variation in modeled turbulence quantities of turbulence intensity and the length scale. Using a one-blade periodic domain, the distributions of unsteady heat transfer rate on the turbine blade and its tip, as affected by the wake, were determined. Such heat transfer coefficient distribution was computed using the wall heat flux and the adiabatic wall temperature to desensitize the heat transfer coefficient to the wall temperature. For the determination of the wall heat flux and the adiabatic wall temperatures, two sets of computations were required. The results were used in a phase-locked manner to compute the unsteady or steady heat transfer coefficients. It has been found that the unsteady wake has some effect on the distribution of the time averaged heat transfer coefficient on the blade and that this distribution is different from the distribution that is obtainable from a steady computation. This difference was found to be as large as 20 percent of the average heat transfer on the blade surface. On the tip surface, this difference is comparatively smaller and can be as large as four percent of the average.