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 253 records · Page 14

Heat transfer model for predicting squib ignition times

A squib ignition model based on transient heat condition from the hot bridgewire to the pyrotechnic is described. No Arrhenius-type chemical reaction is included. Instead, a thermal contact resistance is postulated to exist between the hot bridgewire and the pyrotechnic. Ignition is assumed to occur when a 2.5 micron layer of pyrotechnic next to the bridgewire reaches a characteristic ignition temperature for that pyrotechnic. This model was applied to the JPL squib, which uses a 50 micron (0.002-in.) diameter Tophet A bridgewire to ignite a boron, potassium perchlorate mix. A computer program was utilized that solves the transient heat condition problem with the boundary conditions stipulated by the model. The thermal contact conductance at the interface was determined by trial and error so that the experimentally determined ignition time for one firing condition would be properly predicted by the model. The agreement was quite good for tests run between -129 C and +93.3 C at current levels of 3.5 and 5 A. Axial heat conduction along the bridgewire is shown to be negligible.

Sernas, V.↗

Inhibition of conductive heat flow by magnetic constriction in the corona and transition region - Dependence on the shape of the constriction

Attention is given to the process by which the shape of the constriction of field lines inhibiting the conduction of heat down from the corona acts in conjunction with the degree of constriction to inhibit the heat flow. The analytical model for the heat flow involves tapered flux tubes in which the plasma properties are constant on cross sections, while the plasma is static and the only energy transfer is by thermal conduction. It is determined that only two specific solutions to the model are applicable to the solar atmosphere: the steady state, which is appropriate for quiet regions and for active ones that are not flaring, and the time-dependent case, in which no heat enters the hot end, which is appropriate for conductive cooling of flare loops.

Moore, R. L.↗

Model of hot-film sensor with substrate

A detailed model is constructed to investigate parameters influencing a hot-film anemometer system used to study cross-flow phenomena on swept wings. This anemometer is designed to detect vortices set up by cross-flow. The nature of the vortices is not well known, and a well-verified instrumentation model is considered essential to allow accurate interpretation of experimental results. Previous investigators have indicated that heat conduction to the film's substrate is significant. Two variations on the inclusion of substrate effects in the model are discussed. In the first, the substrate is assumed to have a constant temperature throughput: in the second more elaborate model, the temperature of the substrate is allowed to vary. This variation is related to the convective heat removal at the substrate/gas interface.

Mcree, Griffith J.↗

Finite Element A Posteriori Error Estimation for Heat Conduction

This research investigates residual-based a posteriori error estimates for finite element approximations of heat conduction in single-layer and multi-layered materials. The finite element approximation, based upon hierarchical modelling combined with p-version finite elements, is described with specific application to a two-dimensional, steady state, heat-conduction problem. Element error indicators are determined by solving an element equation for the error with the element residual as a source, and a global error estimate in the energy norm is computed by collecting the element contributions. Numerical results of the performance of the error estimate are presented by comparisons to the actual error. Two methods are discussed and compared for approximating the element boundary flux. The equilibrated flux method provides more accurate results for estimating the error than the average flux method. The error estimation is applied to multi-layered materials with a modification to the equilibrated flux method to approximate the discontinuous flux along a boundary at the material interfaces. A directional error indicator is developed which distinguishes between the hierarchical modeling error and the finite element error. Numerical results are presented for single-layered materials which show that the directional indicators accurately determine which contribution to the total error dominates.

Lang, Christapher G.↗

Assessment of a One-Dimensional Finite Element Charring Ablation Material Response Model for Phenolic-Impregnated Carbon Ablator

In this study, mathematical formulations to model charring ablation problems were numerically implemented using finite element analysis (FEA) with ABAQUS, which account for the material decomposition and progressive surface removal in the heat conduction and the surface energy balance equations. FEA was performed for a one-dimensional model to predict the temperature and ablation histories of a phenolic-impregnated carbon ablator sample (i.e., a common heat shield material for hypersonic vehicles and spacecraft) subjected to oxy-acetylene torch flame (i.e., 0.8 SLPM acetylene gas to 2.7 SLPM oxygen gas). The recovery enthalpy and convective heat transfer coefficient for the ablation model were calculated based on gas compositions and two assumed surface conditions (i.e., equilibrium and frozen). Simulations using the calculated recovery enthalpy and convective heat transfer coefficient resulted in a recession rate of 6.38 times (equilibrium) and 14.08 times (frozen) higher than the experimental data, despite fair agreement of the surface temperature. In addition, the effect of the heat transfer coefficient was investigated through a steady-state ablation analysis. The results of the analysis indicate that there is not one single value for the heat transfer coefficient that would allow the prediction to match both measured recession rate and surface temperature. Possible reasons for such an inconsistency are provided and discussed.

Thermal Protection System↗

Phase change paint tests to investigate effects of TPS tiles on heating rates of the Rockwell space shuttle orbiter (test OH4C, model 21-0)

Information and data from wind tunnel tests conducted on 0.0175-scale models of the space shuttle orbiter are presented. The primary objective of the tests was to evaluate aerodynamic heating effects of the tiles in the thermal protection system (TPS). Tile gap depth and flow orientation effects on the TPS were investigated. Tile patterns were cut into the undersides of the orbiter models to simulate the gaps. One model was left smooth for comparison.

Quan, M.↗

Thermal stress tectonics on the satellites of Saturn and Uranus

Thermal stress histories of the Saturnian and Uranian satellites are investigated. To this end, the thermal evolution of an icy satellite subjected to accretional and radiogenic heating, thermal conduction, and solid-state convection is modeled, and changes in the internal stress that occur during satellite evolution are examined. Results show that internal temperature changes that occur during normal evolution of many of the satellites of Saturn and Uranus can be expected to generate large extensional stresses in the satellites' outer regions. These stresses arise from three sources: (1) radiogenic warming, causing thermal expansion of materials in the satellite's deep interior; (2) radiogenic warming in larger satellites that can induce a phase transition from ice II to ice I and to produce a volume increase in the deep interior; and (3) accretional heating depositing heat in the satellite'e outer regions.

Hillier, John↗

The nature of the residual Martian polar caps

A model of the behavior of the Martian polar caps is described which incorporates the heating effects of the atmosphere, as well as insolation and conduction. This model is used to try to match the observed regression curves of the polar caps, and it predicts that all the seasonally condensed CO2 will be lost by around the summer solstice. The implication is that the residual caps are composed of water ice which, it is found by further modeling, should be stable during the Martian summers. However, it is also argued that this model may be too simplistic, and that the effects of wind in redistributing the seasonal condensate may lead to sufficient thicknesses of CO2 in the central polar region to allow the year-long existence of CO2 without significantly changing the retreat characteristics of the cap, and it is, therefore, concluded that at the present, the nature of the residual caps cannot be reliably determined.

Briggs, G. A.↗

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗

Aerodynamic heating in hypersonic flows

Aerodynamic heating in hypersonic space vehicles is an important factor to be considered in their design. Therefore the designers of such vehicles need reliable heat transfer data in this respect for a successful design. Such data is usually produced by testing the models of hypersonic surfaces in wind tunnels. Most of the hypersonic test facilities at present are conventional blow-down tunnels whose run times are of the order of several seconds. The surface temperatures on such models are obtained using standard techniques such as thin-film resistance gages, thin-skin transient calorimeter gages and coaxial thermocouple or video acquisition systems such as phosphor thermography and infrared thermography. The data are usually reduced assuming that the model behaves like a semi-infinite solid (SIS) with constant properties and that heat transfer is by one-dimensional conduction only. This simplifying assumption may be valid in cases where models are thick, run-times short, and thermal diffusivities small. In many instances, however, when these conditions are not met, the assumption may lead to significant errors in the heat transfer results. The purpose of the present paper is to investigate this aspect. Specifically, the objectives are as follows: (1) to determine the limiting conditions under which a model can be considered a semi-infinite body; (2) to estimate the extent of errors involved in the reduction of the data if the models violate the assumption; and (3) to come up with correlation factors which when multiplied by the results obtained under the SIS assumption will provide the results under the actual conditions.

Reddy, C. Subba↗

Temporal behaviour of the thermal model of hard X-ray bursts

A simple, analytic model is presented of a hot, thermal hard X-ray source, continuously heated, bounded by ion-acoustic conduction fronts, and expanding in a loop. The model is used to investigate the assumption that the 'rise time' of the X-ray emission is approximately given by the loop length divided by the ion-sound speed appropriate to the peak temperature. It is found that a freely-expanding source does not behave in this way; instead, the rise time is symptomatic of the timescale for primary energy release. If the energy release rate does not fall significantly before the source fills the loop, however, then this assumption may be approximately satisfied, if a condition on the temporal behavior of the energy release is satisfied. Finally, some remarks on the relative timing of temperature and emission measure peaks are made, and possible further applications mentioned of the results presented herein.

Mackinnon, A. L.↗

Decoupled Method for Reconstruction of Surface Conditions From Internal Temperatures On Ablative Materials With Uncertain Recession Model

Obtaining measurements of flight environments on ablative heat shields is both critical for spacecraft development and extremely challenging due to the harsh heating environment and surface recession. Thermocouples installed several millimeters below the surface are commonly used to measure the heat shield temperature response, but an ill-posed inverse heat conduction problem must be solved to reconstruct the surface heating environment from these measurements. Ablation can contribute substantially to the measurement response making solutions to the inverse problem strongly dependent on the recession model, which is often poorly characterized. To enable efficient surface reconstruction for recession model sensitivity analysis, a method for decoupling the surface recession evaluation from the inverse heat conduction problem is presented. The decoupled method is shown to provide reconstructions of equivalent accuracy to the traditional coupled method but with substantially reduced computational effort. These methods are applied to reconstruct the environments on the Mars Science Laboratory heat shield using diffusion limit and kinetically limited recession models.

Oliver, A. Brandon↗

Electron Energetics in the Martian Dayside Ionosphere: Model Comparisons with MAVEN Data

This paper presents a study of the energetics of the dayside ionosphere of Mars using models and data from several instruments on board the Mars Atmosphere and Volatile EvolutioN spacecraft. In particular, calculated photoelectron fluxes are compared with suprathermal electron fluxes measured by the Solar Wind Electron Analyzer, and calculated electron temperatures are compared with temperatures measured by the Langmuir Probe and Waves experiment. The major heat source for the thermal electrons is Coulomb heating from the suprathermal electron population, and cooling due to collisional rotational and vibrational CO2 dominates the energy loss. The models used in this study were largely able to reproduce the observed high topside ionosphere electron temperatures (e.g., 3000 K at 300 km altitude) without using a topside heat flux when magnetic field topologies consistent with the measured magnetic field were adopted. Magnetic topology affects both suprathermal electron transport and thermal electron heat conduction. The effects of using two different solar irradiance models were also investigated. In particular, photoelectron fluxes and electron temperatures found using the Heliospheric Environment Solar Spectrum Radiation irradiance were higher than those with the Flare Irradiance Spectrum Model-Mars. The electron temperature is shown to affect the O2(+) dissociative recombination rate coefficient, which in turn affects photochemical escape of oxygen from Mars.

magnetic topology↗

Spatial and Temporal Patterns of Tidal Dissipation in Synchronous Satellites

Tidal heating is an important energy source for several solar system bodies, and there is a wide-spread perception that the pattern of surface heat flow is diagnostic of internal structure. We wish to clarify that situation. Our analysis depends upon two important assumptions: First, that heat transport is dominated by conduction. Second, that the body can be modeled by a sequence of spherically symmetric layers, each with a linear visco-elastic rheology. Under these assumptions, surface heat flow patterns in tidally dominated satellites will reflect radially integrated dissipation patterns. For synchronously rotating satellites with zero obliquity, this pattern depends quite strongly on orbital eccentricity but relatively little on purely radial variations in internal structure. The total amount of heat generated within the body does depend sensitively on internal structure, but the spatial pattern is rather insensitive to structure, especially at low orbital eccentricities.

Bills, Bruce G.↗

Thermal evolution of the earth

The earth's heat budget and models of the earth's thermal evolution are discussed. Sources of the planetary heat are considered and modes of heat transport are addressed, including conduction, convection, and chemical convection. Thermal and convectional models of the earth are covered, and models of thermal evolution are discussed in detail, including changes in the core, the influence of layered mantle convection on the thermal evolution, and the effect of chemical differentiation on the continents.

Spohn, T.↗

The prediction of sea-surface temperature variations by means of an advective mixed-layer ocean model

An advective mixed layer ocean model was developed by eliminating the assumption of horizontal homogeneity in an already existing mixed layer model, and then superimposing a mean and anomalous wind driven current field. This model is based on the principle of conservation of heat and mechanical energy and utilizes a box grid for the advective part of the calculation. Three phases of experiments were conducted: evaluation of the model's ability to account for climatological sea surface temperature (SST) variations in the cooling and heating seasons, sensitivity tests in which the effect of hypothetical anomalous winds was evaluated, and a thirty-day synoptic calculation using the model. For the case studied, the accuracy of the predictions was improved by the inclusion of advection, although nonadvective effects appear to have dominated.

Atlas, R. M.↗