Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

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 145 records · Page 8

Hypersonic engine airplane integration

A summary of the research activities is presented along with the design and analysis of an integrated scramjet engine which utilizes heat conduction flames. The scramjet engine model described was tested at a freestream Mach number of 4.0 by personnel at the General Applied Science Laboratory. The research effort was primarily concerned with the interpretation and analysis of the experimental data. The mixing tests were conducted with cold ambient air, while the combustion tests were conducted with hot vitiated air. The purpose of the cold mixing tests was to determine the degree of mixedness of the fuel at the ignitor station so as to certify if the fuel injection conditions and fuel injectors locations were properly selected to give a uniform fuel concentration profile at the ignitor station. A second objective was to assess the influence of the fuel injection on the inlet flow. Several fuel injection configurations were tested.

Source record↗

Simulation of solidification in a Bridgman cell

The Hg/Cd/Te solidification in a Bridgman cell was thermally analyzed by Continuum's VAST code. The energy equation is solved in an axisymmetric, quasi steady domain for both the molten and solid alloy regions. Alloys composition is calculated by a simplified one dimensional model to estimate its effect on melt thermal conductivity and, consequently, on the temperature field within the cell. Solidification is assumed to occur at a fixed temperature of 979 K. Simplified boundary conditions are included to model both the radiant and conductive heat exchange between the furnace walls and the alloy. Calculations are performed to show how the steady state isotherms are affected by: (1) the hot and cold furnace temperatures; (2) boundary condition parameters; and (3) the growth rate which affects the calculated alloy's composition.

Dakhoul, Y. M.↗

Thermal-capillary analysis of small-scale floating zones Steady-state calculations

Galerkin finite element analysis of a thermal-capillary model of the floating zone crystal growth process is used to predict the dependence of molten zone shape on operating conditions for the growth of small silicon boules. The model accounts for conduction-dominated heat transport in the melt, feed rod and growing crystal and for radiation between these phases, the ambient and a heater. Surface tension acting on the shape of the melt/gas meniscus counteracts gravity to set the shape of the molten zone. The maximum diameter of the growing crystal is set by the dewetting of the melt from the feed rod when the crystal radius is large. Calculations with small Bond number show the increased zone lengths possible for growth in a microgravity environment. The sensitivity of the method to the shape and intensity of the applied heating distribution is demonstrated. The calculations are compared with experimental observations.

Duranceau, J. L.↗

Crystal growth and furnace analysis

A thermal analysis of Hg/Cd/Te solidification in a Bridgman cell is made using Continuum's VAST code. The energy equation is solved in an axisymmetric, quasi-steady domain for both the molten and solid alloy regions. Alloy composition is calculated by a simplified one-dimensional model to estimate its effect on melt thermal conductivity and, consequently, on the temperature field within the cell. Solidification is assumed to occur at a fixed temperature of 979 K. Simplified boundary conditions are included to model both the radiant and conductive heat exchange between the furnace walls and the alloy. Calculations are performed to show how the steady-state isotherms are affected by: the hot and cold furnace temperatures, boundary condition parameters, and the growth rate which affects the calculated alloy's composition. The Advanced Automatic Directional Solidification Furnace (AADSF), developed by NASA, is also thermally analyzed using the CINDA code. The objective is to determine the performance and the overall power requirements for different furnace designs.

Dakhoul, Youssef M.↗

Large-scale volcano-ground ice interactions on Mars

In the present investigation of Martian volcano-ground ice interaction processes, a numerical model is developed that encompasses conductive heat transport, surface radiation, heat transfer to the atmosphere, and H2O phase-changes in an ice-rich permafrost over which lave flows erupt, followed by the intrusion of sills. An examination is made of the two large scale interaction regions formed near Aeolis Mensae and near the volcano Hadriaca Patera, northeast of Hellas. The inferred channel discharges are compared to discharge rates calculated for lava-ground ice interactions, showing that meltwater (probably accumulated under the surface) was rapidly released at a discharge rate that was limited by soil permeability. The volcano-ground ice interactions have been an important Martian geologic process, and could account for the palagonites constituting the Martian dust.

Squyres, Steven W.↗

Calorimetric determination of thermal parameters for the Li/BrCl in SOCl2 (BCX) chemistry

The heat capacity of a Li-BCX DD-cell was found to be dependent on its state of charge by drop calorimetry measurements. The method of drop calorimetry involves measuring the energy (joules) gained or lost from a sample that is transferred from a bath at temperature A to one at temperature B. The thermoneutral potential is defined as the cell potential where the cell electrochemical reactions are neither exothermic nor endothermic. A Hart scientific calorimeter system, Model No. S77XX, designed for heat conduction calorimetry and drop calorimetry was used. Calorimetric analysis yielded a thermoneutral potential of 4.14 volts and a cell heat capacity dependent on the state of charge.

Darcy, Eric C.↗

A novel pulse technique for thermal diffusivity determination of high-temperature levitated materials - Inverse analysis and experimental set-up

In order to determine the thermal diffusivity of materials at high temperatures, a two-step extended flash technique, which is applicable to levitated spherical samples, is proposed. The containerless flash method is modeled as an axisymmetric transient conduction heat transfer problem within the sphere subjected to a radiative boundary condition on its surface. The problem of nonlinearity which arises from the radiative heat transfer boundary condition is solved by replacing it with the measured time-dependent surface temperature data, thus giving rise to an 'Inverse Analysis'. Upon obtaining the analytic solution for the temperature field, the determination of the thermal diffusivity turns into a minimization problem, whereby the difference between the actual surface boundary condition and the analytic solution is minimized. In performing the proposed experiments, there is a need to undertake a cool-down experiment which should be governed by the lumped parameter analysis.

Shen, F.↗

Evaluation of Temperature Gradient in Advanced Automated Directional Solidification Furnace (AADSF) by Numerical Simulation

A numerical model of heat transfer using combined conduction, radiation and convection in AADSF was used to evaluate temperature gradients in the vicinity of the crystal/melt interface for variety of hot and cold zone set point temperatures specifically for the growth of mercury cadmium telluride (MCT). Reverse usage of hot and cold zones was simulated to aid the choice of proper orientation of crystal/melt interface regarding residual acceleration vector without actual change of furnace location on board the orbiter. It appears that an additional booster heater will be extremely helpful to ensure desired temperature gradient when hot and cold zones are reversed. Further efforts are required to investigate advantages/disadvantages of symmetrical furnace design (i.e. with similar length of hot and cold zones).

Bune, Andris V.↗

One-dimensional hybrid satellite track model for the Dynamics Explorer 2 (DE 2) satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. (1985) satellite track code but also includes elements of Roble et al. (1987b) global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the Mass Spectrometer Incoherent Scatter 1990 (MSIS-90) model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule heating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗

Thermal Capacitance (Slug) Calorimeter Theory Including Heat Losses and Other Decaying Processes

A mathematical model, termed the Slug Loss Model, has been developed for describing thermal capacitance (slug) calorimeter behavior when heat losses and other decaying processes are not negligible. This model results in the temperature time slope taking the mathematical form of exponential decay. When data is found to fit well to this model, it allows a heat flux value to be calculated that corrects for the losses and may be a better estimate of the cold wall fully catalytic heat flux, as is desired in arc jet testing. The model was applied to the data from a copper slug calorimeter inserted during a particularly severe high heating rate arc jet run to illustrate its use. The Slug Loss Model gave a cold wall heat flux 15% higher than the value of 2,250 W/sq cm obtained from the conventional approach to processing the data (where no correction is made for losses). For comparison, a Finite Element Analysis (FEA) model was created and applied to the same data, where conduction heat losses from the slug were simulated. The heat flux determined by the FEA model was found to be in close agreement with the heat flux determined by the Slug Loss Model.

Hightower, T. Mark↗

Primordial metamorphism of asteroids via electrical induction in a T Tauri-like solar wind

Possible heating mechanisms for an intense thermal cycle in the early solar system are considered, and inductive asteroid heating in a T Tauri environment is examined. Calculations of thermal models based on recent measurements of the electrical conductivity of carbonaceous meteorites are reported which are presumably more representative of actual conditions in the early solar nebula. Asteroidal thermal evolution is computed using a highly simplified code that restricts induction to the dc limit of the TM mode, approximates the complete unipolar-generation differential equation by a simple series-resistance model, and restricts the modeling of thermal conduction to a calculation of heat transfer across a thin planar shell of given thickness increasing with time, which blankets a body assumed to have an isothermal core. A grid of results is obtained for 20 combinations of asteroid radius and solar distance representing the populaton of asteroids. It is found that there is a relative heating maximum with respect to body radius, that this effect is due to the operation of two distinct heating-rate restrictions which have opposite radius dependences, and that the maximum current density is roughly proportional to the inverse of interior radius.

Herbert, F.↗

Validation Database Based Thermal Analysis of an Advanced RPS Concept

Advanced RPS concepts can be conceived, designed and assessed using high-end computational analysis tools. These predictions may provide an initial insight into the potential performance of these models, but verification and validation are necessary and required steps to gain confidence in the numerical analysis results. This paper discusses the findings from a numerical validation exercise for a small advanced RPS concept, based on a thermal analysis methodology developed at JPL and on a validation database obtained from experiments performed at Oregon State University. Both the numerical and experimental configurations utilized a single GPHS module enabled design, resembling a Mod-RTG concept. The analysis focused on operating and environmental conditions during the storage phase only. This validation exercise helped to refine key thermal analysis and modeling parameters, such as heat transfer coefficients, and conductivity and radiation heat transfer values. Improved understanding of the Mod-RTG concept through validation of the thermal model allows for future improvements to this power system concept.

numerical model↗

The effects of inlet turbulence and rotor/stator interactions on the aerodynamics and heat transfer of a large-scale rotating turbine model. Volume 2: Heat transfer data tabulation. 15 percent axial spacing

A combined experimental and analytical program was conducted to examine the effects of inlet turbulence on airfoil heat transfer. The experimental portion of the study was conducted in a large-scale (approx 5X engine), ambient temperature, rotating turbine model configured in both single stage and stage-and-a-half arrangements. Heat transfer measurements were obtained using low-conductivity airfoils with miniature thermcouples welded to a thin, electrically heated surface skin. Heat transfer data were acquired for various combinations of low or high inlet turbulence intensity, flow coefficient, first-stator/rotor axial spacing, Reynolds number and relative circumferential position of the first and second stators. Aerodynamic measurements obtained as part of the program include distributions of the mean and fluctuating velocities at the turbine inlet and, for each airfoil row, midspan airfoil surface pressures and circumferential distributions of the downstream steady state pressures and fluctuating velocities. Analytical results include airfoil heat transfer predictions produced using existing 2-D boundary layer computation schemes and an examination of solutions of the unsteady boundary layer equations. The results are reported in four separate volumes, of which this is Volume 2: Heat Transfer Data Tabulation; 15 Percent Axial Spacing.

Dring, R. P.↗

An experimental determination in Calspan Ludwieg tube of the base environment of the integrated space shuttle vehicle at simulated Mach 4.5 flight conditions (test IH5 of model 19-OTS)

A model test program was conducted to determine heat transfer and pressure distributions in the base region of the space shuttle vehicle during simulated launch trajectory conditions of Mach 4.5 and pressure altitudes between 90,000 and 210,000 feet. Model configurations with and without the solid propellant booster rockets were examined to duplicate pre- and post-staging vehicle geometries. Using short duration flow techniques, a tube wind tunnel provided supersonic flow over the model. Simultaneously, combustion generated exhaust products reproduced the gasdynamic and thermochemical structure of the main vehicle engine plumes. Heat transfer and pressure measurements were made at numerous locations on the base surfaces of the 19-OTS space shuttle model with high response instrumentation. In addition, measurements of base recovery temperature were made indirectly by using dual fine wire and resistance thermometers and by extrapolating heat transfer measurements.

Drzewiecki, R. F.↗

Solar wind heat flow

The ability of a theoretical solar wind model to predict the energy flux density due to heat conduction was investigated. In order to demonstrate this capability, the analysis used simultaneous data from the electron plasma experiment of Ogilvie and Lind, and the vector magnetic field data of Heppner and Sugiura on OGO-5. On the strength of two independent, different, and consistent experimental determinations of the conduction of heat flux by electrons at 149.6 gigameters(1 AU), it was determined that the models for the solar wind must be able to predict a value for the heat conduction in the range 0.007 to 0.008 in order to be in agreement with the observations.

Scudder, J.↗

Validating Arc-jet Condition Targets Using Inverse 1D Material Response Models

Experiments were conducted using both the Interaction Heating Facility (IHF) and Aerodynamic Heating Facility (AHF), two hypersonic arc-jet facilities at NASA Ames. These tests were run to understand the material response of ablative thermal protection system materials by exposing them to high-enthalpy flows. Using thermocouple data from the materials, test conditions were inversely determined and compared with expected results from a one-dimensional model using facility target conditions; this allowed an approximate environmental scale factor to be determined and applied to more accurately predict material response in additional tests.

material response↗

Finite elements of nonlinear continua.

The finite element method is extended to a broad class of practical nonlinear problems, treating both theory and applications from a general and unifying point of view. The thermomechanical principles of continuous media and the properties of the finite element method are outlined, and are brought together to produce discrete physical models of nonlinear continua. The mathematical properties of the models are analyzed, and the numerical solution of the equations governing the discrete models is examined. The application of the models to nonlinear problems in finite elasticity, viscoelasticity, heat conduction, and thermoviscoelasticity is discussed. Other specific topics include the topological properties of finite element models, applications to linear and nonlinear boundary value problems, convergence, continuum thermodynamics, finite elasticity, solutions to nonlinear partial differential equations, and discrete models of the nonlinear thermomechanical behavior of dissipative media.

Oden, J. T.↗

Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave

Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.

Cambier, Jean-Luc↗