Search NASA⌕ Search

SEARCH · Search NASA

Results for “Heat flux calculation”

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 307 records · Page 17

A Streamlined Approach to Spatial Mapping of Complex 3D Thermal Boundary Condition Data

Spatial mapping of 1D, 2D, and particularly 3D thermal boundary condition data is essential for accurate thermal simulation of spacecraft and launch vehicles. This information typically is available in simple ASCII file format but can be problematic to use for thermal simulations. Difficulties that can be encountered include the ability to directly use the available file format, transformation of the data into the desired reference coordinate system in 3D space, and accurate interpolation of the data onto the thermal calculation domain. This paper presents a streamlined approach to spatial mapping of complex unsymmetrical 3D cold wall heat flux data using the fields’ capability in Simcenter3D. A recent example of 3D thermal plume simulation at MSFC will be employed to illustrate the use of text, csv, and Excel files, specifying coordinate systems for transformation of data in 3D space, and use of the various interpolation schemes available for spatial mapping. Spatial boundary condition verification is also very important, and the spacecraft thermal analyst needs effective visualization tools to develop confidence in the boundary condition definitions. Visualization of the spacecraft plume raw field data, the calculated field data, and the interpolated field data onto the thermal calculation domain will be discussed, with a particular focus on real-time viewing of the interpolated data for available interpolation schemes. The objective of this paper is to outline the modeling techniques used with the integrated thermal analysis of a terrestrial lunar lander demonstrator. Specifically, the techniques outlined herein focus on spatial mapping of convective and radiative exhaust plume results from a computational fluid dynamics (CFD) analysis. Two cases were investigated, however, the modeling techniques that were utilized translate to both cases seamlessly. Those two cases were as follows: • Lunar lander demonstrator firing on-pad (worst case). • Lunar lander demonstrator hovering at an elevated position above the launch pad. For the purposes of information control, this paper will focus on the modeling techniques only, using altered or dimensionless data where necessary. BACKGROUND The thermal model examined in this paper is of the lunar lander demonstrator known as XL-1T (terrestrial), born from a collaborative effort between Masten Space Systems (MSS) and NASA, under NASA’s Lunar CATALYST (also known as Lunar Cargo Transportation and Landing by Soft Touchdown) initiative. The lander is a reusable vertical takeoff/vertical landing (VTVL) test bed which is controlled by four throttleable main engines utilizing green hypergolic propellants.

Hawkins, Robert F.↗

Hypersonic rarefied flow over sharp slender cones

Drag, heat transfer, and number flux were measured on sharp cones in the near free, molecule flow regime, and the results were compared with available Monte Carlo calculations. In general, the calculations predicted the magnitude of the data; however, the heat transfer and drag increased with increasing Knudsen number at a faster rate than predicted. Also the drag coefficients measured for the slender cones at high Knudsen number were higher than predicted for free molecule flow. These disagreements between theory and experiment could possibly be attributed to the simplicity of the surface interaction laws assumed in the theory. Reynolds analogy factors obtained from the experimental measurements agreed with free-molecule values and also with that obtained by the Monte Carlo technique.

Kussoy, M. I.↗

Experimental and Numerical Characterization of High Heat Fluxes During Transient Blackbody Calibrations

High heat fluxes are encountered in numerous applications, such as hypersonic vehicles in flight, fires, and engines, Calibration of heat flux gages may be performed in a dual cavity cylindrical blackbody resulting in a transient calibration environment. To characterize the transient heat fluxes. experiments were performed on a dual cavity cylindrical blackbody at nominal temperatures varying from 800 C to 1900 C in increments of 100 C. Based on experiments, the optimum heat flux sensor insertion location as measured from the center partition was determined. The pre-insertion steady state axial temperature profile is compared experimentally, numerically, and analytically. The effect of convection in the blackbody cavity during the insertion is calculated and found to be less than 2 per cent. Also, an empirical correlation for predicting the emissivity of the blackbody is included. Detailed transient thermal models have been developed to simulate the heat flux calibration process at two extreme fluxes. The high (1MW/sq m) and relatively low (70 kw/sq m) fluxes are reported in this article. The transient models show the effect of inserting a heat flux gage at room temperature on the thermal equilibrium of the blackbody at 1800 C and 800 C nominal temperatures, respectively. Also, heat flux sensor outputs are derived from computed sensor temperature distributions and compared to experimental results.

Abdelmessih, Amanie N.↗

Development of Solar Wind Model Driven by Empirical Heat Flux and Pressure Terms

We are developing a time stationary self-consistent 2D MHD model of the solar corona and solar wind as suggested by Sittler et al. (2003). Sittler & Guhathakurta (1999) developed a semiempirical steady state model (SG model) of the solar wind in a multipole 3-streamer structure, with the model constrained by Skylab observations. Guhathakurta et al. (2006) presented a more recent version of their initial work. Sittler et al. (2003) modified the SG model by investigating time dependent MHD, ad hoc heating term with heat conduction and empirical heating solutions. Next step of development of 2D MHD models was performed by Sittler & Ofman (2006). They derived effective temperature and effective heat flux from the data-driven SG model and fit smooth analytical functions to be used in MHD calculations. Improvements of the Sittler & Ofman (2006) results now show a convergence of the 3-streamer topology into a single equatorial streamer at altitudes > 2 R(sub S). This is a new result and shows we are now able to reproduce observations of an equatorially confined streamer belt. In order to allow our solutions to be applied to more general applications, we extend that model by using magnetogram data and PFSS model as a boundary condition. Initial results were presented by Selwa et al. (2008). We choose solar minimum magnetogram data since during solar maximum the boundary conditions are more complex and the coronal magnetic field may not be described correctly by PFSS model. As the first step we studied the simplest 2D MHD case with variable heat conduction, and with empirical heat input combined with empirical momentum addition for the fast solar wind. We use realistic magnetic field data based on NSO/GONG data, and plan to extend the study to 3D. This study represents the first attempt of fully self-consistent realistic model based on real data and including semi-empirical heat flux and semi-empirical effective pressure terms.

Sittler, Edward C., Jr.↗

Interior radiances in optically deep absorbing media. I - Exact solutions for one-dimensional model.

An exact analytic solution to the one-dimensional scattering problem with arbitrary single scattering albedo and arbitrary surface albedo is presented. Expressions are given for the emergent flux from a homogeneous layer, the internal flux within the layer, and the radiative heating. A comparison of these results with the values calculated from the matrix operator theory indicates an exceedingly high accuracy. A detailed study is made of the error in the matrix operator results and its dependence on the accuracy of the starting value.

Kattawar, G. W.↗

The computation of radiation from nonequilibrium hypersonic flows

The results of the solution of the equations that describe a hypersonic ionized flow about an elliptically blunted cone are presented. The flow conditions correspond to those of the proposed Aeroassist Flight Experiment (AFE) vehicle at altitudes between the perigee at 78 km and the approximate limit of the continuum regime at 90 km. For the free-stream velocities of interest, about 9 km/sec, the flowfield is out of thermo-chemical equilibrium, electronically excited, ionized and radiating. The gas consists of eight-chemical species including free electrons. The thermal state of the gas is modeled with a translational-rotational temperature, four vibrational temperatures for the diatomic species and an electron-electronic temperature. The electronic excitation of molecules is included. The nonequilibrium air radiation from each fluid element is computed and the radiative heat flux at the body surface is determined. The stagnation point radiative heating result agrees with previous calculations.

Candler, Graham↗

Turbine rotor/stator flowfield analysis

A numerical study of the unsteady aerodynamic and thermal environment associated with axial turbine stages is presented. Computations are performed using a modification of the ROTOR1 rotor/stator interaction code. Two different turbine states are analyzed: the first state of the United Technologies Research Center large scale rotating rig and the first state of the Space Shuttle main engine (SSME) high pressure fuel turbopump. Time-averaged blade midspan pressure and heat transfer profiles are calculated using the following different surface boundary conditions: adiabatic wall, prescribed wall temperature, and prescribed heat flux. Numerical solutions for the large scale rotating rig are compared with experimental data. Unsteady pressure envelopes are also presented for each geometry. In addition, instantaneous contours are plotted for the SSME configuration.

Griffin, Lisa↗

Rates of fluvio-thermal erosion on Mars

In order to get some ideas about the order of magnitude of fluvio-thermal recession rate, a mathematical thermal model for Martian outflow channels is proposed. This model corresponds to a system undergoing a permanent thermal regime where the surface temperature is constant and equal to the phase change temperature (due to the immediate removal of melted materials). This is an ablation model. For its application, estimations of the heat transfer coefficient h and thermal flux q are necessary. Determination of these coefficients needs the calculation of dimensionless numbers (Reynolds, Prandtl, and Nusselt), and the consideration of turbulent regime of the flow.

Aguirre-Puente, J.↗

Evaluation of a Thermodynamic Ice Crystal Icing Model Using Experimental Ice Accretion Data

This paper presents the evaluation of a thermodynamic ice crystal icing model, previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. It has been proposed that there are two types of distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. In the former, ice accretion occurs where a freeze fraction (0 to 1) of melted ice crystals freezes on a surface, along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized by having strong adhesion to the surface. In the latter, icing occurs from accumulated unmelted ice on a surface, where a melt fraction (0 to 1) dictates the amount of unmelted impinged ice. This melt-dominated icing is characterized by weakly bonded surface adhesion. The experimentally observed ice growth rates suggest that only a small fraction of the impinging ice remains on the surface, implying a mass loss mechanism such as splash, runback, bounce, or erosion. This mass loss parameter must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface. This loss parameter, however, along with the freeze and melt fraction, are the only experimental parameters that are currently not measured directly. Using reported icing growth rates from published ice crystal icing experiments, a methodology is proposed to determine these unknown parameters. This work takes reported ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) in 2016 and tests NASA collaborated on with the National Research Council (NRC) of Canada in 2012 that examined the fundamental physics of ice crystal icing. Those research efforts sought to generate icing conditions representative of those that occur inside a jet engine when ingesting ice crystals. This paper presents the fundamental equations of the thermodynamic model, the methodology used to determine the aforementioned unknown icing parameters, and results from model evaluation using experimental data. In addition, this paper builds on the previously proposed model by adding a transient conduction term to explain ice growth behavior at the onset of experimental tests that was observed to be different from steady-state ice growth that occurred later in the test run.With the addition of this energy term, this becomes a quasi-steady model. A key finding from this work suggests that mass loss fractions can exceed 0.90 for steady ice growth periods. In addition, due to conductive heat fluxes when using a warmer-than-freezing airfoil, lower mass loss fraction values were calculated during the initial transient period.

Ice-Crystal Ice Accretion↗

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.↗

Nonisothermal Pluto atmosphere models

The present thermal profile calculation for a Pluto atmosphere model characterized by a high number fraction of CH4 molecules encompasses atmospheric heating by solar UV flux absorption and conductive transport cooling to the surface of Pluto. The stellar occultation curve predicted for an atmosphere of several-microbar surface pressures (which entail the existence of a substantial temperature gradient close to the surface) agrees with observations and implies that the normal and tangential optical depth of the atmosphere is almost negligible. The minimum period for atmospheric methane depletion is calculated to be 30 years.

Hubbard, W. B.↗

Oceanic Fluxes of Mass, Heat and Freshwater: A Global Estimate and Perspective

Data from fifteen globally distributed, modern, high resolution, hydrographic oceanic transects are combined in an inverse calculation using large scale box models. The models provide estimates of the global meridional heat and freshwater budgets and are used to examine the sensitivity of the global circulation, both inter and intra-basin exchange rates, to a variety of external constraints provided by estimates of Ekman, boundary current and throughflow transports. A solution is found which is consistent with both the model physics and the global data set, despite a twenty five year time span and a lack of seasonal consistency among the data. The overall pattern of the global circulation suggested by the models is similar to that proposed in previously published local studies and regional reviews. However, significant qualitative and quantitative differences exist. These differences are due both to the model definition and to the global nature of the data set.

MacDonald, Alison Marguerite↗

Photoelectron fluxes in the Martian ionosphere

Calculations are presented of the steady-state photoelectron distribution in the upper atmosphere of Mars, consistent with the neutral upper atmosphere and ionosphere particle concentrations and temperatures measured by Viking 1. Uncertainties in the calculations affect the thermal electron gas heating rate. Major conclusions are that (1) over most of the altitude range of the Martian ionosphere, the steady-state photoelectron flux amplitude is larger than that in the earth's, so that photoelectron-impact-excited airglow on Mars is generally more significant than it is on earth; (2) the steady-state photoelectron energy distribution in the Martian ionosphere is softer and more structured than that in the terrestrial ionosphere; (3) photoelectron impact ionization contributes about 30% to the total ionization rate in the Martian ionosphere; and (4) photoelectron impact excitation contributes 20-30% of the CO2(+) and CO zenith airglow emissions on Mars.

Mantas, G. P.↗

Experimental investigation of two-phase heat transfer in a porous matrix.

One-dimensional two-phase flow transpiration cooling through porous metal is studied experimentally. The experimental data is compared with a previous one-dimensional analysis. Good agreement with calculated temperature distribution is obtained as long as the basic assumptions of the analytical model are satisfied. Deviations from the basic assumptions are caused by nonhomogeneous and oscillating flow conditions. Preliminary derivation of nondimensional parameters which characterize the stable and unstable flow conditions is given. Superheated liquid droplets observed sputtering from the heated surface indicated incomplete evaporation at heat fluxes well in access of the latent energy transport. A parameter is developed to account for the nonequilibrium thermodynamic effects. Measured and calculated pressure drops show contradicting trends which are attributed to capillary forces.

Von Reth, R.↗

Calculation procedure for transient heat transfer to a cooled plate in a heated stream whose temperature varies arbitrarily with time

Solutions for the surface temperature and surface heat flux are found for laminar, constant property, slug flow over a plate convectively cooled from below, when the temperature of the fluid over the plate varies arbitrarily with time at the plate leading edge. A simple technique is presented for handling arbitrary fluid temperature variation with time by approximating it by a sequence of ramps or steps for which exact analytical solutions are available.

Sucec, J.↗

Linear modes of convection in the solar envelope

The structure of the solar convective envelope calculated according to mixing-length theory is not consistent with the transport of the convective heat flux by any superposition of the fundamental linear convective modes. If the amplitudes of these modes can be estimated by consideration of the most important nonlinear terms in the equations of motion, it becomes possible to find solar models which are consistent with mode transport of the heat through the convective zone. These models may be characterized by a mixing length which varies from point to point through the envelope. The linear modes in such models show peaks in velocity at certain preferred length scales, but these scales do not necessarily correspond to those of the solar motions, such as the supergranulation.

Bogart, R. S.↗

Radiation-dominated Snowmelt on Mars

A set of microclimatic conditions conducive to the generation of liquid water is being sought to see whether their requirements are any less stringent than those suggested by the global-climatic models. In particular, the melting of snow and ice in the equatorial zone (where the valley networks are preferentially situated) is investigated with a model that includes detailed infrared radiation transfer between the atmosphere and the snow, penetration of sunlight into the snow, sensible heat transfer to the atmosphere, and water vapor diffusion within the surface boundary layer. IR fluxes are calculated by using pressure and temperature dependent forms for the band absorptances of CO2 and H2O in the spectral range 50 to 2000/cm.

Clow, G. D.↗

Observations and theory of the solar semidiurnal tide in the mesosphere of Venus

The 1978 Pioneer Venus mission, in addition to producing strong IR radiometer experiment evidence for the existence of a semidiurnal solar tide in the Venus mesosphere, has yielded radiative flux and cloud particle distribution data about the structure of solar heating and IR cooling which facilitate calculations of tidal forcing functions and estimates of radiative dissipation rates. It is presently shown that a physically reasonable mean flow can be found which leads to theoretically predicted tides that are in excellent agreement with those observed. The zonal mean flow obtained exhibits a maximum velocity of about 130 m/sec at 70 km altitude, decreasing to about zero at 88 km, together with a jet-like meridional structure that strongly deviates from solid body rotation.

Fels, S. B.↗