SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
An effective temperature regulating wall consists of one layer /e.g., one of the paraffins/ relatively opaque to thermal radiation in the solid state and transparent to it in the molten state and placed between two transparent layers. A mirror coating is applied to back layer.
Phase change material heat exchangers (PCMHX) have an important role in integrating renewable energy systems. PCMHX can offer high storage density at various temperatures of interest and can be used for grid load shifting purposes. Numerical models enable engineers to estimate PCMHX performance for different design parameters and operating conditions. Modeling phase change phenomena is challenging due to the complex time-dependent nature of the process. The accuracy of models is highly sensitive to PCM thermo-physical properties. Thermal expansion coefficient (β), viscosity (μ) and melting temperature range (MR) of a PCM are important properties, especially when natural convection is not negligible. In PCMHX modeling, using less than accurate values for these properties can have significant impact on the simulation outcomes. These properties and discussions thereof are not readily available in the literature. This paper presents a brief review of the literature and a numerical study investigating the model sensitivity to the above-mentioned properties for PCMHX with and without fins. CFD is used to evaluate the charging (melting) phenomena. The study quantifies the impact of uncertainty in these properties on the melting rate and temperature profile. Constant wall temperature was considered as heat source with no heat loss to ambient. Here, the results show that β and μ has significant effect on the melting rate and evolution of the melting front. For a non-finned domain when comparing results for different published values of β and μ, the deviation in melting time can be up to 12.9% and 57.6% respectively. For high wall temperatures, change in melting range did not impact melting time. But when the wall temperature is reduced, up to 9.8% deviation in melting time is observed.
Thermal barrier coatings (TBCs) have been of interest since the 1970s for application in internal combustion (IC) engines. Thin TBCs exhibit a temperature swing phenomenon wherein wall temperatures dynamically respond to the transient working-gas temperature throughout the engine cycle, thus reducing the temperature difference driving the heat transfer. Determining these varying wall temperatures is necessary to evaluate and study the effect of coatings on wall heat transfer. This study focuses on developing a 3D computational fluid dynamics (CFD)-finite element analysis (FEA) coupled simulation, or co-simulation, routine to determine the wall temperatures of a piston coated with a thin TBC layer subject to spark ignition combustion heat flux. A CONVERGE 3D-CFD model was used to simulate the combustion process in a single-cylinder, light-duty experimental spark ignition (SI) engine. Transient piston heat transfer analysis was conducted using ABAQUS, a FEA package, under the simulated combustion heat flux load. The effect of the temperature swing phenomenon due to this TBC layer was observed in a CFD simulation by implementing the FEA results as the piston thermal boundary conditions. The boundary conditions were passed between the CFD and FEA tools until a quasi-steady state solution was achieved. Furthermore, a reduction in wall heat transfer was observed due to a reduced temperature difference between the wall and the working gas.
A simplified analytical model of leakage into a wing-elevon cove under high speed external flow conditions is presented. The model accounts for convection from the high temperature gas stream leaking into the cove, radiation exchanges between walls and with the environment, and axial conduction within the channel walls. Comparisons are made of predicted and measured wall temperatures of an experimental wing-elevon-cove model that was tested in a hypersonic wind tunnel. Analytically determined Nusselt numbers for laminar flow entrance conditions are compared with those calculated from experimentally measured gas temperatures, wall temperatures and wall heating rates. Although reasonable agreement appears to be indicated, testing times much larger than those of the present test program (5 seconds) would be desirable in evaluating the reliability and accuracy of the simplified model.
NASA has developed a Simulated Inter-Duct Research Model (SIDRM) to study ice crystal icing in aircraft engines. A numerical investigation has been completed using the NASA in-house code, GlennHT. For each flow condition, many different wall thermal boundary conditions were imposed. Six uniform wall temperature cases were generated, with wall temperature relative to the freestream total temperature of 0.85, 0.90, 0.95, 1.05, 1.10, and 1.15. In addition, five cases with uniform heat flux were generated. Heat fluxes were chosen to generate average wall temperatures in the range explored by the uniform wall temperature cases. Included in the uniform heat flux cases is the adiabatic case where heat flux equals zero. The results of these simulations are then used to quantify the degree to which the heat transfer coefficient depends on wall thermal boundary condition. It is seen that the heat transfer coefficient is fairly independent of wall thermal boundary condition, as expected. However, it would appear that, at least in some locations, it may be prudent to attempt to correlate a variation of heat transfer coefficient with wall thermal boundary condition.
Convective heat transfer characteristics and theoretical thermal stress behaviors are numerically calculated using RELAP5-3D for the helical-coiled once-through steam generator (H-OTSG) and the novel heat exchanger design known as the oval-twisted helically coiled heat exchanger (OTHCHX) under (1) fluctuating wall temperature conditions, (2) square-wave pulsating flow conditions, and (3) the combined effects of fluctuating wall temperature and square-wave pulsating flow conditions. Heat transfer coefficient models for the H-OTSG and OTHCHX were developed based on existing data and implemented into RELAP5-3D, successfully capturing the Nuavg behavior within 8% to 10% of the reported data. Under fluctuating wall temperature conditions, the OTHCHX displayed higher Nu avg behavior than the H-OTSG. As 𝑓 increased, the $𝜎^{𝑚𝑎𝑥}_{𝑡ℎ}$ decreased. The $𝜎^{𝑚𝑎𝑥}_{𝑡ℎ}$ was higher for the OTHCHX than for the H-OTSG under fluctuating wall temperature conditions. Under pulsating flow conditions, the H-OTSG and OTHCHX displayed much higher 𝑁𝑢 𝑎𝑣𝑔 than under constant flow conditions. The H-OTSG displayed a higher $𝜎^{𝑚𝑎𝑥}_{𝑡ℎ}$ over the OTHCHX. Under combined fluctuating wall temperature and pulsating flow conditions, the augmented heat transfer behavior from the pulsating flow was counteracted by the wall temperature fluctuations, producing slightly higher 𝑁𝑢 𝑎𝑣𝑔 over constant wall temperature, constant flow conditions, but much lower than only constant pulsating flow under constant wall temperature conditions. The effects of simultaneous wall temperature fluctuations and square-wave pulsating flow caused higher $𝜎^{𝑚𝑎𝑥}_{𝑡ℎ}$ than that of only wall temperature fluctuations or pulsating flow. As the Reynolds number (Re) increased, $𝜎^{𝑚𝑎𝑥}_{𝑡ℎ}$ increased. However, when 𝑓=𝑓$_{\dot{m}}$, the $𝜎^{𝑚𝑎𝑥}_{𝑡ℎ}$ showed decreasing values as Re increased. In conclusion, the results indicate that thermal-fluid resonance can help mitigate thermal stresses.
Direct numerical simulations (DNS) of Mach 6 turbulent boundary layer with nominal freestream Mach number of 6 and Reynolds number of Re(sub T) approximately 460 are conducted at two wall temperatures (Tw/Tr = 0.25, 0.76) to investigate the generated pressure fluctuations and their dependence on wall temperature. Simulations indicate that the influence of wall temperature on pressure fluctuations is largely limited to the near-wall region, with the characteristics of wall-pressure fluctuations showing a strong temperature dependence. Wall temperature has little influence on the propagation speed of the freestream pressure signal. The freestream radiation intensity compares well between wall-temperature cases when normalized by the local wall shear; the propagation speed of the freestream pressure signal and the orientation of the radiation wave front show little dependence on the wall temperature.
A calculational procedure is presented to analyze the heat-transfer and fluid-flow characteristics of a convectively air -cooled plug-nozzle operating on an afterburning turbojet engine. Anderson's method was used to predict hot-gas static pressures in the supersonic stream with fully expanded flow (high nozzle-pressure ratios); the results were excellent. For low nozzle-pressure ratios, the flow was assumed to expand one-dimensionally and isentropically to the plug back pressure. Wall temperatures predicted using this latter pressure distribution agreed well with the wall temperatures predicted using the measured hot-gas pressures (maximum deviation was about 30 K (54 deg R)). Either an in tegral boundary-layer technique or a simple pipe-flow equation may be used to calculate convective heat transfer from the hot gas to the wall. The simple pipeflow equation results in the prediction of slightly higher wall temperatures than does the integral technique. Experimental wall temperatures were generally in good agreement with the two predicted wall temperature distributions. Excellent agreement was noted b etween measured and predicted coolant static-pressure distributions. The plug-coolant temperature rise was generally overpredicted by about 22.2 K (40 deg R); possible explanations are offered. Although an an alysis of the struts, which support the plug, was purposely kept simple, reasonable results were obtained. Potential flow over an ellipse was used to calculate hot-gas static pressure; the results were satisfactory.
An analytical investigation on the effect of high aspect ratio (height/width) cooling channels, considering different coolant channel designs, on hot-gas-side wall temperature and coolant pressure drop for a liquid hydrogen cooled rocket combustion chamber, was performed. Coolant channel design elements considered were: length of combustion chamber in which high aspect ratio cooling was applied, number of coolant channels, and coolant channel shape. Seven coolant channel designs were investigated using a coupling of the Rocket Thermal Evaluation code and the Two-Dimensional Kinetics code. Initially, each coolant channel design was developed, without consideration for fabrication, to reduce the hot-gas-side wall temperature from a given conventional cooling channel baseline. These designs produced hot-gas-side wall temperature reductions up to 22 percent, with coolant pressure drop increases as low as 7.5 percent from the baseline. Fabrication constraints for milled channels were applied to the seven designs. These produced hot-gas-side wall temperature reductions of up to 20 percent, with coolant pressure drop increases as low as 2 percent. Using high aspect ratio cooling channels for the entire length of the combustion chamber had no additional benefit on hot-gas-side wall temperature over using high aspect ratio cooling channels only in the throat region, but increased coolant pressure drop 33 percent. Independent of coolant channel shape, high aspect ratio cooling was able to reduce the hot-gas-side wall temperature by at least 8 percent, with as low as a 2 percent increase in coolant pressure drop. ne design with the highest overall benefit to hot-gas-side wall temperature and minimal coolant pressure drop increase was the design which used bifurcated cooling channels and high aspect ratio cooling in the throat region. An optimized bifurcated high aspect ratio cooling channel design was developed which reduced the hot-gas-side wall temperature by 18 percent and reduced the coolant pressure drop by 4 percent. Reductions of coolant mass flow rate of up to 50 percent were possible before the hot-gas-side wall temperature reached that of the baseline. These mass flow rate reductions produced coolant pressure drops of up to 57 percent.
The structural integrity of high pressure liquid propellant rocket engine thrust chambers is typically maintained through regenerative cooling. The coolant flows through passages formed either by constructing the chamber liner from tubes or by milling channels in a solid liner. Recently, Carlile and Quentmeyer showed life extending advantages (by lowering hot gas wall temperatures) of milling channels with larger height to width aspect ratios (AR is greater than 4) than the traditional, approximately square cross section, passages. Further, the total coolant pressure drop in the thrust chamber could also be reduced, resulting in lower turbomachinery power requirements. High aspect ratio cooling channels could offer many benefits to designers developing new high performance engines, such as the European Vulcain engine (which uses an aspect ratio up to 9). With platelet manufacturing technology, channel aspect ratios up to 15 could be formed offering potentially greater benefits. Some issues still exist with the high aspect ratio coolant channels. In a coolant passage of circular or square cross section, strong secondary vortices develop as the fluid passes through the curved throat region. These vortices mix the fluid and bring lower temperature coolant to the hot wall. Typically, the circulation enhances the heat transfer at the hot gas wall by about 40 percent over a straight channel. The effect that increasing channel aspect ratio has on the curvature heat transfer enhancement has not been sufficiently studied. If the increase in aspect ratio degrades the secondary flow, the fluid mixing will be reduced. Analysis has shown that reduced coolant mixing will result in significantly higher wall temperatures, due to thermal stratification in the coolant, thus decreasing the benefits of the high aspect ratio geometry. A better understanding of the fundamental flow phenomena in high aspect ratio channels with curvature is needed to fully evaluate the benefits of this geometry. The fluid dynamic and conjugate heat transfer problem of high aspect ratio rocket engine coolant channels are being investigated numerically, but these efforts have been hampered by a lack of validating data. Wall temperature data is available for the conjugate problem for channels without curvature and aspect ratio = 5.0, and unheated fluid dynamic data are available for square and circular cross section channels with curvature at Reynold's numbers up to 40,000. But the effects of aspect ratio on secondary flow development have not been experimentally studied. To provide some insight into the effects of channel aspect ratio on secondary flow and to qualitatively provide anchoring for the numerical codes, a flow visualization experiment was initiated at the NASA Lewis Research Center.
An analytical investigation on the effect of high aspect ratio (height/width) cooling channels, considering different coolant channel designs, on hot-gas-side wall temperature and coolant pressure drop for a liquid hydrogen cooled rocket combustion chamber, was performed. Coolant channel design elements considered were: length of combustion chamber in which high aspect ratio cooling was applied, number of coolant channels, and coolant channel shape. Seven coolant channel designs were investigated using a coupling of the Rocket Thermal Evaluation code and the Two-Dimensional Kinetics code. Initially, each coolant channel design was developed, without consideration for fabrication, to reduce the hot-gas-side wall temperature from a given conventional cooling channel baseline. These designs produced hot-gas-side wall temperature reductions up to 22 percent, with coolant pressure drop increases as low as 7.5 percent from the baseline. Fabrication constraints for milled channels were applied to the seven designs. These produced hot-gas-side wall temperature reductions of up to 20 percent, with coolant pressure drop increases as low as 2 percent. Using high aspect ratio cooling channels for the entire length of the combustion chamber had no additional benefit on hot-gas-side wall temperature over using high aspect ratio cooling channels only in the throat region, but increased coolant pressure drop 33 percent. Independent of coolant channel shape, high aspect ratio cooling was able to reduce the hot-gas-side wall temperature by at least 8 percent, with as low as a 2 percent increase in coolant pressure drop. The design with the highest overall benefit to hot-gas-side wall temperature and minimal coolant pressure drop cooling can now be done in relatively short periods of time with multiple iterations.
An apparatus was built to verify an analysis of combined forced and free convection in a vertical tube with uniform wall heat flux and to determine the limits of the analysis. The test section was electrically heated by resistance heating of the tube wall and was instrumented with thermocouples in such a way that detailed thermal entrance heat-transfer coefficients could be obtained for both upflow and downflow and any asymmetry in wall temperature could be detected. The experiments showed that fully developed heat-transfer results, predicted by a previous analysis, were confirmed over the range of Rayleigh numbers investigated. The concept of "locally fully developed" heat transfer was established. This concept involves the assumption that the fully developed heat-transfer analysis can be applied locally even though the Rayleigh number is varying along the tube because of physical-property variations with temperature. Thermal entrance region data were obtained for pure forced convection and for combined forced and free convection. The analysis of laminar pure forced convection in the thermal entrance region conducted by Siegel, Sparrow, and Hallman was experimentally confirmed. A transition to an eddy motion, indicated by a fluctuation in wall temperature was found in many of the upflow runs. A stability correlation was found. The fully developed Nusselt numbers in downflow were below those for pure forced convection but fell about 10 percent above the analytical curve. Quite large circumferential variations in wall temperature were observed in downflow as compaired with those encountered in upflow, and the fully developed Nussalt numbers reported are based on average wall temperatures determined by averaging the readings of two diametrically opposite wall thermocouples at each axial position. With larger heating rates in downflow the wall temperature distributions strongly suggested a cell flow near the bottom. At still larger heating rates the wall temperatures varied in a periodic way.
This paper presents a comparison between experimental results from recent liquid hydrogen (LH 2 ) transfer line chilldown experiments at high Reynolds (Re) numbers versus liquid nitrogen (LN 2 ) experiments conducted at low Re numbers. Parasitic heat leak, inner wall temperatures, inner wall heat fluxes, and heat transfer coefficients are computed to compare between the two systems. Analysis of temperature traces and flow visualization indicates that the chilldown process evolves much more rapidly at higher Re numbers due to a quick transition from vapor flow to annular liquid flow and near immediate liquid contact along the pipe walls. The lower kinematic viscosity and surface tension of LH 2 , along with reduced parasitic heat leak and higher Re numbers relative to the LN 2 experiments, causes chilldown to proceed almost immediately into the nucleate boiling regime, in comparison to low Re flows where >75% of the chilldown is spent in vapor film boiling.
Boiling heat transfer to liquid hydrogen and nitrogen was investigated experimentally. Results are presented from a study of bulk boiling inside a cylindrical tube under vertically upward forced-flow conditions. A 0.555-inch-inside-diameter and 16 (1/8)-inch-long electrically heated stainless-steel tube was used. The range of variables studied for hydrogen were mass velocity of 2850 to 17,000 pounds per hour per square foot, local heat flux of 3600 to 40,000 Btu per hour per square foot, inlet pressure of 30 to 74 pounds per square inch absolute, and inlet subcooling of 0° to 9° R. Nitrogen test conditions were mass velocity of 15,000 to 56,000 pounds per hour per square foot, local heat flux of 2300 to 40,000 Btu per hour per square foot, inlet pressure of 47 to 56 pounds per square inch absolute, and inlet sub cooling of 1° to 6° R. The axial distribution of the tube-wall temperatures is presented. A transition in the type of boiling heat transfer was obtained. The critical heat flux corresponding to this transition was determined over a range of flow and heating rates and local qualities. At specific combinations of flow and transition location, a range of critical-heat-flux values was obtained and maximum values were determined. The maximum critical heat flux increased with increasing fluid-flow rate and decreased with increasing length of tube before transition. Similar variations of the maximum critical heat flux have been reported for water. The tube-inner-wall temperatures upstream of transition were essentially uniform and were only slightly greater than the fluid saturation temperature. The Wall-temperature profiles downstream of transition generally resembled those obtained in film-boiling studies and appeared to be strongly dependent upon local quality at the point of transition. Maximum wall temperatures of 9000 and 18000 R were obtained with hydrogen and nitrogen, respectively. Fluctuations of pressure, flow rate, and temperature occurred during some of the boiling tests. Under some conditions, maximum critical-heat-flux values were attained during steady-state operation with fluctuations. In other cases the fluctuations became uncontrolled, and critical-flux values less than the maximum values were obtained upon restabilization of the test conditions. No measurable pressure drop across the test section was obtained at any condition.
Heat transfer to a full coverage film cooled turbulent boundary layer over a flat surface was studied. The surface consisted of a discrete hole test section containing 11 rows of holes spaced 5 diameters apart in a staggered array and an instrumented recovery region. Ten diameter spacing was also studied by plugging appropriate holes. Two test sections were used, one having holes normal to the surface and the other having holes angled 30 deg to the surface in the downstream direction. Stanton number data were obtained both in the full coverage region and in the downstream recovery region for a range of blowing ratios, or mass flux ratios, from 0 to 1.3. Initial conditions at the upstream edge of the blowing region were varied from 500 to 5000 for momentum thickness Reynolds number and from 100 to 1800 for enthalpy thickness Reynolds number. The range of Reynolds numbers based on hole diameter and mainstream velocity was 6000 to 22000. Initial boundary layer thicknesses range from 0.5 to 2.0 hole diameters. Air was used as the working fluid. The data were taken for the secondary injection temperature equal to the wall temperature and also equal to the mainstream temperature. Superposition was then used to obtain Stanton number as a continuous function of the injectant temperature. The heat transfer coefficient was defined on the basis of a mainstream-to-wall temperature difference. This definition permits direct comparison of performance between film cooling and transpiration cooling.
A heat-transfer investigation was conducted with air flowing through an electrically heated silicon carbide tube with a rounded entrance, an inside diameter of 3/4 inch, and an effective heat-transfer length of 12 inches over a range of Reynolds numbers up to 300,000 and a range of average inside-tube-wall temperature up to 2500 R. The highest corresponding local outside-tube-wall temperature was 3010 R. Correlation of the heat-transfer data using the conventional Nueselt relation wherein physical properties of the fluid were evaluated at average bulk temperature resulted in a separation of data with tube-wall-temperature level. A satisfactory correlation of the heat-transfer data was obtained, however, by the use of modified correlation parameters wherein the mass velocity G (or product of average air density and velocity evaluated at bulk temperature P(sub b)V(sub b)) in the Reynolds number was replaced by the product of average air velocity evaluated at the bulk temperature and density evaluated at either the average inside-tube-wall temperature or the average film temperature; in addition, all the physical properties of air were correspondingly evaluated at either the average inside-tube-wall temperature or the average film temperature.
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.
Temperature distribution and radiation flux measurements in uranium plasma reactor with reflecting walls