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.
Nuclear rocket nozzle cooling passages, discussing heat transfer and friction correlations for single-phase hydrogen turbulent flow
Closed-form solutions of the one-dimensional heat-conduction equations for the flow of heat into a plate with a laminar boundary layer have been obtained for a configuration entering a planetary atmosphere with constant velocity and negative entry angle. The atmospheric density was assumed to obey an exponential law and the temperature was assumed constant initially. The solution is in the form of a Fourier series expansion which, for most practical applications, can be approximated by retaining only one term of the expression. The solution applies to the initial part of the entry before the maximum heating conditions are encountered.
Analytical and numerical solutions to the heat-conduction equation are obtained for the heating of absorbing media with pulsed lasers. The spatial and temporal form of the temperature is determined using several different models of the laser irradiance. Both surface and volume generation of heat are discussed. It is found that if the depth of thermal diffusion for the laser-pulse duration is large compared to the optical-attenuation depth, the surface- and volume-generation models give nearly identical results. However, if the thermal-diffusion depth for the laser-pulse duration is comparable to or less than the optical-attenuation depth, the surface-generation model can give significantly different results compared to the volume-generation model. Specific numerical results are given for a tungsten target irradiated by pulses of different temporal durations and the implications of the results are discussed with respect to the heating of metals by picosecond laser pulses.
A finite element formulation for the analysis of polymer processing is presented and its use in some typical situation including entry flow, transient Couette flow, and the Graetz (forced convection) problem is illustrated. The element formulations are constructed on the premise that momentum convection can be neglected (polymer melt flows typically have very low Reynolds' numbers), but that convective heat transfer may be significant (high Peclet numbers). Nonisothermal effects are considered important in polymer processing, due in part to the significant heating which may occur due to viscous dissipation, and also to the very strong influence of temperature on fluid viscosity. The flow is treated as Newtonian with the flow field being coupled to the heat transfer equation only through the viscous heat generation.
The vertical structure of Titan's thermosphere is calculated down to the mesopause as a function of local time based on Voyager 1 occultation data. The thermal time scales that characterize the diurnal behavior of the thermosphere are discussed, the input model atmosphere used to calculate the temperature profile is presented, and the dominant heating and cooling mechanisms in the thermosphere are examined. The temperature profiles obtained by integrating the heat transfer equation with and without electron heating are presented and discussed. The implications that derived exospheric temperatures have for the neutral hydrogen torus are investigated. The diurnal exospheric temperature is unlikely to exceed 225 K, averages between 187 and 197 K, and has a variational amplitude of 28 K or less. The vertical extent of the hydrogen cloud is too large to be explained in terms of simple thermal escape of hydrogen from the exosphere.
Partial differential equation for heat transfer in plane incompressible laminar jet is reduced by similarity transformation to ordinary differential equation
High flux evaporations from a steady meniscus formed in a 2 micron channel is modeled using the augmented Young-Laplace equation. The heat flux is found to be a function of the long range van der Waals dispersion force which represents interfacial conditions between heptane and various substrates. Heat fluxes of (1.3-1.6) x 10(exp 6) W/m(exp 2) based on the width of the channel are obtained for heptane completely wetting the substrate at 100 C. Small channels are used to obtain these large fluxes. Even though the real contact angle is 0 deg, the apparent contact angle is found to vary between 24.8 deg and 25.6 deg. The apparent contact angle, which represents viscous losses near the contact line, has a large effect on the heat flow rate because of its effect on capillary suction and the area of the meniscus. The interfacial heat flux is modeled using kinetic theory for the evaporation rate. The superheated state depends on the temperature and the pressure of the liquid phase. The liquid pressure differs from the pressure of the vapor phase due to capillarity and long range van der Waals dispersion forces which are relevant in the ultra think film formed at the leading edge of the meniscus. Important pressure gradients in the thin film cause a substantial apparent contact angle for a complete wetting system. The temperature of the liquid is related to the evaporation rate and to the substrate temperature through the steady heat conduction equation. Conduction in the liquid phase is calculated using finite element analysis except in the vicinity of the thin film. A lubrication theory solution for the thin film is combined with the finite element analysis by the method of matched asymptotic expansions.
The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.
Research on the topic of shock wave/turbulent boundary-layer interaction was carried out during the past three years at the Penn State Gas Dynamics Laboratory. This report describes the experimental research program which provides basic knowledge and establishes new data on heat transfer in swept shock wave/boundary-layer interactions. An equilibrium turbulent boundary-layer on a flat plate is subjected to impingement by swept planar shock waves generated by a sharp fin. Five different interactions with fin angle ranging from 10 deg to 20 deg at freestream Mach numbers of 3.0 and 4.0 produce a variety of interaction strengths from weak to very strong. A foil heater generates a uniform heat flux over the flat plate surface, and miniature thin-film-resistance sensors mounted on it are used to measure the local surface temperature. The heat convection equation is then solved for the heat transfer distribution within an interaction, yielding a total uncertainty of about +/- 10 percent. These experimental data are compared with the results of numerical Navier-Stokes solutions which employ a k-epsilon turbulence model. Finally, a simplified form of the peak heat transfer correlation for fin interactions is suggested.
An experimental research program providing basic knowledge and establishing new data on the heat transfer in swept shock wave/boundary-layer interactions is described. An equilibrium turbulent boundary-layer on a flat plate is subjected to impingement by swept planar shock waves generated by a sharp fin. Five different interactions with fin angles ranging from 10 to 20 deg at freestream Mach numbers of 3.0 and 4.0 produce a variety of interaction strengths from weak to very strong. A foil heater generates a uniform heat flux over the flat plate surface and miniature thin-film-resistance sensors mounted on it are used to measure the local surface temperature. The heat convection equation is then solved for the heat transfer distribution within an interaction, yielding a total uncertainty of about +/- 10 percent. These experimental data are compared with the results of numerical Navier-Stokes solutions which employ a kappa-epsilon turbulence model. Finally, a simplified form of the peak heat transfer correlation for fin interactions is suggested.
Explore the source record for details and available documents.
An experimental and computational research program providing new knowledge of the heat transfer in swept-shock-wave/boundary-layer interactions is described. An equilibrium turbulent boundary layer on a flat plate is subjected to impingement by a swept planar shock wave generated by a sharp fin. Five different interactions with fin angles ranging from 10 to 20 deg at freestream Mach numbers of 3 and 4 produce a variety of interaction strengths ranging from weak to very strong. A foil heater generates a uniform heat flux over the flat plate surface, and miniature thin-film-resistance sensors are used to measure the local surface temperature. The heat convection equation is then solved for the heat transfer distribution within an interaction, yielding an uncertainty of about +/- 10%. These data are compared with numerical Navier-Stokes solutions that employ a k-epsilon turbulence model. A simple peak heat transfer correlation for fin interactions is suggested.
The equations for superfluid helium flowing through a straight tube with heated walls are integrated. The model equations are based on those of Landau as modified for superfluid turbulence by Gorter and Mellink (1949). The model is implemented by the program SUPERFLOW which runs on a personal computer. The effect of the heating level on the mass flux, the energy flux, and the pressure and temperature profiles is investigated. The four types of profiles which occur without sidewall heating are also found with heated walls. The progression through these four types is shown to depend primarily on the dimensionless parameter, defined previously to characterize the unheated profiles, and a dimensionless ratio of heat fluxes. The pressure and temperature maxima which rise well above the boundary values increase significantly with sidewall heating. Approximate design equations for estimating the mass flux and the profiles are presented. The physical basis of the results is discussed. These results are useful in the design of optimal cooling loops and other superfluid flow systems.
Microstructure observations in the Pacific cold tongue reveal that strong turbulence often penetrates into the thermocline producing hundreds of W/m2 of downward heat transport during nighttime and early morning. However, virtually all observations of this deep-cycle turbulence (DCT) are from 0◦N,140◦W. Here, a hierarchy of ocean process simulations including submesoscale-permitting regional models and turbulence-permitting large eddy simulations (LES) embedded in a regional model provide insight into DCT at and beyond 0◦N,140◦W. A regional hindcast quantifies the spatio-temporal variability of subsurface turbulent heat fluxes throughout the cold tongue from 1999-2016. Climatologically, the mean heat flux is strongest (∼100W/m2) within2◦ of the equator. The heat flux is slightly (∼10 W/m2) stronger in the northern than southern hemisphere throughout the cold tongue. The seasonal cycle of the heat flux, which is not everywhere aligned with the surface flux, ranges from 150 W/m2 near the equator to 30 W/m2 and 10 W/m2 at 4◦N and S respectively. Aseasonal variability of turbulent heat fluxes is logarithmically distributed and temporally uncorrelated with surface heat fluxes. The aseasonal variability of turbulence is highlighted in 34-day LES of Boreal autumn at 0◦N and 3◦N,140◦W. Intense DCT occurs frequently above the undercurrent at 0◦N and intermittently at 3◦N. Daily-mean heat fluxes scale with the bulk vertical shear and the wind stress, which together explain ∼90% of the daily variance across both LES. Observational validation of the scaling at 0◦N,140◦W is encouraging, but observations beyond 0◦N,140◦W can facilitate needed refinement of mixing parameterization in ocean models.
The thermal electron heating rate is an important heat source term in the ionospheric electron energy balance equation, representing heating by photoelectrons or by precipitating higher energy electrons. A formula for the thermal electron heating rate is derived from the kinetic equation using the electron-electron collision operator as given by the unified theory of Kihara and Aono. This collision operator includes collective interactions to produce a finite collision operator with an exact Coulomb logarithm term. The derived heating rate O(e) is the sum of three terms, O(e) = O(p) + S + O(int), which are respectively: (1) primary electron production term giving the heating from newly created electrons that have not yet suffered collisions with the ambient electrons; (2) a heating term evaluated on the energy surface m(e)/2 = E(T) at the transition between Maxwellian and tail electrons at E(T); and (3) the integral term representing heating of Maxwellian electrons by energetic tail electrons at energies ET. Published ionospheric electron temperature studies used only the integral term O(int) with differing lower integration limits. Use of the incomplete heating rate could lead to erroneous conclusions regarding electron heat balance, since O(e) is greater than O(int) by as much as a factor of two.
The thermal electron heating rate is an important heat source term in the ionospheric electron energy balance equation, representing heating by photoelectrons or by precipitating higher energy electrons. A formula for the thermal electron heating rate is derived from the kinetic equation using the electron-electron collision operator as given by the unified theory of Kihara and Aono. This collision operator includes collective interactions to produce a finite collision operator with an exact Coulomb logarithm term. The derived heating rate O(e) is the sum of three terms, O(e) = O(P) + S + O(int), which are respectively: (1) primary electron production term giving the heating from newly created electrons that have not yet suffered collisions with the ambient electrons; (2) a heating term evaluated on the energy surface m(e)/2 = E(T) at the transition between Maxwellian and tail electrons at E(T); and (3) the integral term representing heating of Maxwellian electrons by eneegetic tail electrons at energies ET. Published ionospheric electron temperature studies used only the integral term O(int) with differing lower integration limits. Use of the incomplete heating rate could lead to erroneous conclusions regarding electron heat balance, since O(e) is greater than O(int) by as much as a factor of two. Previously announced in STAR as N84-15941
Radiative-conductive heat transfer has been investigated for the ground-atmosphere system of the planet Mars. The basic goal was the quantitative determination of time dependent vertical distributions of temperature and static stability for Southern-Hemispheric summer season and middle and polar latitudes, for both dust-free and dust-laden atmospheric conditions. The numerical algorithm which models at high spatial and temporal resolution the thermal energy transports in the dual ground-atmosphere system, is based on solution of the applicable heating rate equation, including radiative and molecular-conductive heat transport terms. The two subsystems are coupled by an internal thermal boundary condition applied at the ground-atmosphere interface level. Initial data and input parameters are based on Mariner 4, 6, 7, and 9 measurements and the JPL Mars Scientific Model. Numerical experiments were run for dust-free and dust-laden conditions in the midlatitudes, as well as ice-free and ice-covered polar regions. Representative results and their interpretation are presented. Finally, the theoretical framework of the generalized problem with nonconservative Mie scattering and explicit thermal-convective heat transfer is formulated, and applicable solution algorithms are outlined.
Control of both the magnitude and the shape of tokamak profiles will be necessary to achieve stable, high-performance plasmas. In order to reject disturbances in real time, feedback-control algorithms rely on accurate real-time knowledge of the plasma state. When diagnostics alone are insufficient, either because they are limited in number or their measurements are too noisy, observers can be used to combine diagnostic data with a response model to provide a better estimation of different plasma properties. An observer has been developed to estimate the electron temperature profile in real time using both diagnostic data from the Thomson scattering system and a model based on the electron heat transport equation describing the evolution of the electron temperature profile. Here, neural network surrogate models are leveraged to help improve the overall model prediction while staying within computation time constraints for real-time use. The observer algorithm is shown in offline tests to produce smooth profiles that are consistent with both the diagnostic data and the electron heat transport equation. When implemented into the real-time plasma control system, this observer will provide valuable information on the electron temperature profile to many potential feedback-control applications.