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 181 records · Page 10

Computational analyses of crystal growth

Two important aspects of Hg/Cd/Te crystal growth processes are discussed. First, the thermal field and second, the fluid movement in the melt zone. The thermal analysis includes numerical calculation of axisymmetric heat conduction within the sample. It also includes a three-dimensional radiation model to calculate the radiative heat exchange between the furnace and the crystal as determined by the complex geometry of the furnace and the adiabatic shield. The thermal analysis also includes a crystal conductivity which is dependent on temperature and composition. To tackle the fluid flow aspect of the problem, an attempt was made to use a newly developed incompressible flow code based on the slight compressibility, and hence the finite sound speed, of all real fluids.

Dakhoul, Youssef M.↗

Unsteady temperature distribution in volume reflectors

The development of unsteady temperature distributions in semitransparent scattering materials under the influence of arbitrary incident radiative flux with other than specified temperature boundary conditions is considered. Specifically, such a problem is considered for an approximate radiative transfer model which allows analytic solutions. The model represents the unsteady, one-dimensional, radiative heating of a weakly absorbing, semitransparent medium with specified conductive heat flux at one boundary. The radiation field is modeled using an approximation to the radiative heat flux obtained by solution of the Kubelka-Munk differential equations.

Weston, K. C.↗

Experimental determination of convective heat transfer coefficients in the separated flow region of the Space Shuttle Solid Rocket Motor

A series of cold flow heat transfer tests was conducted with a 7.5-percent scale model of the Space Shuttle Rocket Motor (SRM) to measure the heat transfer coefficients in the separated flow region around the nose of the submerged nozzle. Modifications were made to an existing 7.5 percent scale model of the internal geometry of the aft end of the SRM, including the gimballed nozzle in order to accomplish the measurements. The model nozzle nose was fitted with a stainless steel shell with numerous thermocouples welded to the backside of the thin wall. A transient 'thin skin' experimental technique was used to measure the local heat transfer coefficients. The effects of Reynolds number, nozzle gimbal angle, and model location were correlated with a Stanton number versus Reynolds number correlation which may be used to determine the convective heating rates for the full scale Space Shuttle Solid Rocket Motor nozzle.

Whitesides, R. Harold↗

Particle simulation of rarefied aeropass maneuvers of the Magellan spacecraft

Rarefied and free molecular flowfields corresponding to proposed aeropass maneuvers through the atmosphere of Venus by the Magellan spacecraft are computed with a vectorized particle simulation method. The significance of rarefaction and reaction effects are assessed and surface heat flux, drag coefficients, and flowfield properties are computed. A simple surface heat transfer model is coupled directly to the simulation and assumes that each surface element is in radiative equilibrium with deep space. This allows direct computation of surface temperature distributions rather than requiring prescribed isothermal boundary conditions. Uncoupled dual-node heat transfer models, which account for heat capacity and thermal conductivity within each structural component, permit more accurate determination of surface temperatures. Such temperatures restrict the allowable altitudes and entry speeds of the aeropass maneuvers, and must therefore be estimated accurately. Simulations with an entry velocity of 8600 m/s at altitudes between 125 and 140 km reveal that allowable surface temperatures occur only at the highest altitudes where the flowfield is appropriately modeled as free-molecular. Excessive surface heat flux occurs at lower altitudes where molecular collisions are significant.

Haas, Brian L.↗

Thermal Testing and Analysis of an Efficient High-Temperature Multi-Screen Internal Insulation

Conventional multi-layer insulations exhibit excellent insulation performance but they are limited to the temperature range to which their components reflective foils and spacer materials are compatible. For high temperature applications, the internal multi-screen insulation IMI has been developed that utilizes unique ceramic material technology to produce reflective screens with high temperature stability. For analytical insulation sizing a parametric material model is developed that includes the main contributors for heat flow which are radiation and conduction. The adaptation of model-parameters based on effective steady-state thermal conductivity measurements performed at NASA Langley Research Center (LaRC) allows for extrapolation to arbitrary stack configurations and temperature ranges beyond the ones that were covered in the conductivity measurements. Experimental validation of the parametric material model was performed during the thermal qualification test of the X-38 Chin-panel, where test results and predictions showed a good agreement.

Weiland, Stefan↗

Shock Formation of Slow Magnetosonic Waves in Coronal Plumes

We investigate the height of shock formation in coronal plumes for slow magnetosonic waves. The models take into account plume geometric spreading, heat conduction, and radiative damping. The wave parameters as well as the spreading functions of the plumes and the base magnetic field strengths are given by empirical constraints mostly from Solar and Heliospheric Observatory (SOHO)/ Ultraviolet Coronograph Spectrometer (UVCS), Extreme Ultraviolet Imaging Telescope (EIT), Michelson Doppler Imager (MDI), and Large Angle Spectrometric Coronagraph (LASCO). Our models show that shock formation occurs at relatively low coronal heights, typically within 1.2 RsuN, depending on the model parameters. The shock formation is calculated using the well-established wave breaking criterion given by the intersection of C+ characteristics in the space-time plane. Our models show that shock heating by slow magnetosonic waves is expected to be relevant at most heights in solar coronal plumes, although such waves are probably not the main energy supply mechanism.

Cuntz, Manfred↗

The rebound shock model for solar spicules - Dynamics at long times

The spicule model due to Hollweg is extended and developed. The dynamics is emphasized here; radiative and ionization losses, heat conduction, and nonshock heat input, are not included. In the model, a series of rebound shocks results in chromospheric material with spicule-like properties below a raised transition region. The shocks result from a single quasi-impulsive source in the photosphere. It is found that at long times, the model approaches a new hydrostatic equilibrium with the transition region remaining raised, and with a region of shock-heated chromosphere below it. Attention is given to the variation of the properties of the model in response to different values for the magnitude and location of the source, and to different initial transition region heights. It is concluded that the model is capable of generating structures with properties consistent with observations of spicules (with the exception of temperature) when only the dynamics is considered.

Sterling, Alphonse C.↗

Gap heating with pressure gradients

The heating rate distribution and temperature response on the gap walls of insulating tiles is analyzed to determine significant phenomena and parameters in flows where there is an external surface pressure gradient. Convective heating due to gap flow, modeled as fully developed pipe flow, is coupled with a two-dimensional thermal model of the tiles that includes conduction and radiative heat transfer. To account for geometry and important environmental parameters, scale factors are obtained by curve-fitting measured temperatures to analytical solutions. These scale factors are then used to predict the time-dependent gap heat flux and temperature response of tile gaps on the Space Shuttle Orbiter during entry.

Scott, C. D.↗

Flow instability in particle-bed nuclear reactors

A three-dimensional model of the stability of the particle-bed reactor is presented, in which the fluid has mobility in three dimensions. The model accurately represents the stability at low Re numbers as well as the effects of the cold and hot frits and of the heat conduction and radiation in the particle bed. The model can be easily extended to apply to the cylindrical geometry of particle-bed reactors. Exemplary calculations are carried out, showing that a particle bed without a cold frit would be subject to instability if operated at the high-temperature ratios used for nuclear rockets and at power densities below about 4 MW/l; since the desired power density for such a reactor is about 40 MW/l, the operation at design exit temperature but at reduced power could be hazardous. Calculations show however that it might be possible to remove the instability problem by appropriate combinations of cold and hot frits.

Kerrebrock, J. L.↗

Finite Volume Numerical Methods for Aeroheating Rate Calculations from Infrared Thermographic Data

The use of multi-dimensional finite volume heat conduction techniques for calculating aeroheating rates from measured global surface temperatures on hypersonic wind tunnel models was investigated. Both direct and inverse finite volume techniques were investigated and compared with the standard one-dimensional semi-infinite technique. Global transient surface temperatures were measured using an infrared thermographic technique on a 0.333-scale model of the Hyper-X forebody in the NASA Langley Research Center 20-Inch Mach 6 Air tunnel. In these tests the effectiveness of vortices generated via gas injection for initiating hypersonic transition on the Hyper-X forebody was investigated. An array of streamwise-orientated heating striations was generated and visualized downstream of the gas injection sites. In regions without significant spatial temperature gradients, one-dimensional techniques provided accurate aeroheating rates. In regions with sharp temperature gradients caused by striation patterns multi-dimensional heat transfer techniques were necessary to obtain more accurate heating rates. The use of the one-dimensional technique resulted in differences of 20% in the calculated heating rates compared to 2-D analysis because it did not account for lateral heat conduction in the model.

Daryabeigi, Kamran↗

A finite element model of conduction, convection, and phase change near a solid/melt interface

Detailed understanding of heat transfer and fluid flow is required for many aerospace thermal systems. These systems often include phase change and operate over a range of accelerations or effective gravitational fields. An approach to analyzing such systems is presented which requires the simultaneous solution of the conservation laws of energy, momentum, and mass, as well as an equation of state. The variable property form of the governing equations are developed in two-dimensional Cartesian coordinates for a Newtonian fluid. A numerical procedure for solving the governing equations is presented and implemented in a computer program. The Galerkin form of the finite element method is used to solve the spatial variation of the field variables, along with the implicit Crank-Nicolson time marching algorithm. Quadratic Langrangian elements are used for the internal energy and the two components of velocity. Linear Lagrangian elements are used for the pressure. The location of the solid/liquid interface as well as the temperatures are determined form the calculated internal energy and pressure. This approach is quite general in that it can describe heat transfer without phase change, phase change with a sharp interface, and phase change without an interface. Analytical results from this model are compared to those of other researchers studying transient conduction, convection, and phase change and are found to be in good agreement. The numerical procedure presented requires significant computer resources, but this is not unusual when compared to similar studies by other researchers. Several methods are suggested to reduce the computational times.

Viterna, Larry A.↗

A numerical study of a three dimensional spherical thermospheric density and wind model

Numerical calculations of the generation and propagation of the two important fundamental symmetric tidal wave modes - the diurnal mode (1, 1, 1,) and the semidiurnal mode (2, 2, 2) - were performed applying a realistic model thermosphere and taking into account heat conduction and the temporally and spatially varying ion-neutral collision number. Both wave modes are predominantly generated by the solar EUV heat input. It is shown that the latitude structure of the (1, 1, 1)-mode which is identical with the Hough function(1, -1) within the lower non-dissipative atmosphere degenerates into the spherical function P sub 1, 1 at thermospheric heights. The pressure field of this mode constitutes the observed pressure bulge of the thermosphere, the diurnal component of which peaks at 15 h L. T. The electric polarization field of the geomagnetic Sq current generates a significant fraction of this wave mode at F layer heights. This wave component shifts the total horizontal wind system to earlier times by about 1 hour in agreement with ionospheric observations. The latitude structure of the (2, 2, 2) mode is identical with the Hough function (2, 2) within the lower non-dissipative atmosphere. It degenerates to the spherical function P sub 2, 2 at thermospheric heights.

Volland, H.↗

A mixed pseudospectral/finite difference method for a thermally driven fluid in a nonuniform gravitational field

A numerical study of the steady, axisymmetric flow in a heated, rotating spherical shell is conducted to model the Atmospheric General Circulation Experiment (AGCE) proposed to run aboard a later shuttle mission. The AGCE will consist of concentric rotating spheres confining a dielectric fluid. By imposing a dielectric field across the fluid a radial body force will be created. The numerical solution technique is based on the incompressible Navier-Stokes equations. In the method a pseudospectral technique is based on the incompressible Navier-Stokes equations. In the method a pseudospectral technique is used in the latitudinal direction, and a second-order accurate finite difference scheme discretizes time and radial derivatives. This paper discusses the development and performance of this numerical scheme for the AGCE which has been modelled in the past only by pure FD formulations. In addition, previous models have not investigated the effect of using a dielectric force to simulate terrestrial gravity. The effect of this dielectric force on the flow field is investigated as well as a parameter study of varying rotation rates and boundary temperatures. Among the effects noted are the production of larger velocities and enhanced reversals of radial temperature gradients for a body force generated by the electric field.

Macaraeg, M. G.↗

A mixed pseudospectral/finite difference method for a thermally driven fluid in a nonuniform gravitational field

A numerical study of the steady, axisymmetric flow in a heated, rotating spherical shell is conducted to model the Atmospheric General Circulation Experiment (AGCE) proposed to run aboard a later Shuttle mission. The AGCE will consist of concentric rotating spheres confining a dielectric fluid. By imposing a dielectric field across the fluid a radial body force will be created. The numerical solution technique is based on the incompressible Navier-Stokes equations. In the method a pseudospectral technique is used in the latitudinal direction, and a second-order accurate finite difference scheme discretizes time and radial derivatives. This paper discusses the development and performance of this numerical scheme for the AGCE which has been modeled in the past only by pure FD formulations. In addition, previous models have not investigated the effect of using a dielectric force to simulate terrestrial gravity. The effect of this dielectric force on the flow field is investigated as well as a parameter study of varying rotation rates and boundary temperatures. Among the effects noted are the production of larger velocities and enhanced reversals of radial temperature gradients for a body force generated by the electric field.

Macaraeg, M. G.↗

Modeling of combustion processes of stick propellants via combined Eulerian-Lagrangian approach

This research is motivated by the improved ballistic performance of large-caliber guns using stick propellant charges. A comprehensive theoretical model for predicting the flame spreading, combustion, and grain deformation phenomena of long, unslotted stick propellants is presented. The formulation is based upon a combined Eulerian-Lagrangian approach to simulate special characteristics of the two phase combustion process in a cartridge loaded with a bundle of sticks. The model considers five separate regions consisting of the internal perforation, the solid phase, the external interstitial gas phase, and two lumped parameter regions at either end of the stick bundle. For the external gas phase region, a set of transient one-dimensional fluid-dynamic equations using the Eulerian approach is obtained; governing equations for the stick propellants are formulated using the Lagrangian approach. The motion of a representative stick is derived by considering the forces acting on the entire propellant stick. The instantaneous temperature and stress fields in the stick propellant are modeled by considering the transient axisymmetric heat conduction equation and dynamic structural analysis.

Kuo, K. K.↗

Experimental determination of in situ utilization of lunar regolith for thermal energy storage

A Lunar Thermal Energy from Regolith (LUTHER) experiment has been designed and fabricated at the NASA Lewis Research Center to determine the feasibility of using lunar soil as thermal energy storage media. The experimental apparatus includes an alumina ceramic canister which contains simulated lunar regolith, a heater, nine heat shields, a heat transfer cold jacket, and 19 type-B platinum rhodium thermocouples. The simulated lunar regolith is a basalt that closely resembles the lunar basalt returned to earth by the Apollo missions. The experiment will test the effects of vacuum, particle size, and density on the thermophysical properties of the regolith, which include melt temperature, specific heat thermal conductivity, and latent heat of storage. Two separate tests, using two different heaters, will be performed to study the effect of heating the system using radiative and conductive heat transfer. A finite differencing SINDA model was developed at NASA Lewis Research Center to predict the performance of the LUTHER experiment. The code will predict the effects of vacuum, particle size, and density has on the heat transfer to the simulated regolith.

Richter, Scott W.↗

Experimental and Computational Investigations of Phase Change Thermal Energy Storage Canisters

Two sets of experimental data are examined in this paper, ground and space experiments, for cylindrical canisters with thermal energy storage applications. A 2-D computational model was developed for unsteady heat transfer (conduction and radiation) with phase-change. The radiation heat transfer employed a finite volume method. The following was found in this study: (1) Ground Experiments: the convection heat transfer is equally important to that of the radiation heat transfer; radiation heat transfer in the liquid is found to be more significant than that in the void; including the radiation heat transfer in the liquid resulted in lower temperatures (about 15 K) and increased the melting time (about 10 min.); generally, most of the heat flow takes place in the radial direction. (2) Space Experiments: radiation heat transfer in the void is found to be more significant than that in the liquid (exactly the opposite to the Ground Experiments); accordingly, the location and size of the void affects the performance considerably; including the radiation heat transfer in the void resulted in lower temperatures (about 40 K).

Ibrahim, Mounir↗

Comprehensive 3D Multiphysics Model on Electrochemical Recovery of O 2 from Metabolic CO 2 at the International Space Station (ISS)

The International Space Station (ISS) is presently equipped with an elaborate, heavy, and high-power consuming system that recovers approximately 50% of O 2 from metabolic CO 2 as part of the atmospheric revitalization (AR) at the ISS habitat. Future long-duration missions will require a sustainable and efficient system capable of yielding a minimum of 75% O 2 recovery to reach the self-sufficiency required for long space missions beyond earth’s low orbit. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to not only increase significantly current O 2 recovery efficiency, improving self-sufficiency on AR at the ISS habitat and future long missions, but also reduce the complexity of the system. The authors have developed and deployed a comprehensive 3D multiphysics model that thoroughly replicates the actual configuration and fluid/material domains of the MFECR. The coupled physics in this multiphysics model include multicomponent-multiphase electrochemical-driven reactions, non-ideal mass transport mechanism, free and porous flow, heat transfer, CO 2 solubility on alkaline electrolyte, water condensation on porous medium, and DC electrical current generation along with Joule heating effect. This model is aimed to conduct quantitive benchmark on three different MFECR’s layouts, one without serpentine paths (plain) and two with serpentines leading to four and twelve paths respectively. Once experimental data is generated via a test matrix of 200 tests, the model will be validated to conduct MFECR’s process optimization and revalidate the quantitive benchmark on three different MFECR’s layouts.

Jesus A. Dominguez↗