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 505 records · Page 28

Modeling the seasonal cycle of CO2 on Mars: A fit to the Viking lander pressure curves

We have constructed a more accurate Mars thermal model, similar to the one used by Leighton and Murray in 1966, which solves radiative, conductive, and latent heat balance at the surface as well as the one-dimensional heat conduction equation for 40 layers to a depth of 15 meters every 1/36 of a Martian day. The planet is divided into 42 latitude bands with a resolution of two degrees near the poles and five degrees at lower latitudes, with elevations relative to the 6.1 mbar reference areoid. This estimate of the Martian zonally averaged topography was derived primarily from radio occultations. We show that a realistic one-dimensional thermal model is able to reproduce the VL1 pressure curve reasonably well without having to invoke complicated atmospheric effects such as dust storms and polar hoods. Although these factors may cause our deduced values for each model parameter to differ from its true value, we believe that this simple model can be used as a platform to study many aspects of the Martian CO2 cycle over seasonal, interannual, and long-term climate timescales.

Wood, S. E.↗

Thermal math model analysis of FRSI test article subjected to cold soak and entry environments

A multi-objective test program was conducted at the NASA/JSC Radiant Heat Test Facility in which an aluminum skin/stringer test panel insulated with FRSI (Flexible Reusable Surface Insulation) was subjected to 24 simulated Space Shuttle Orbiter ascent/entry heating cycles with a cold soak in between in the 10th and 20th cycles. A two-dimensional thermal math model was developed and utilized to predict the thermal performance of the FRSI. Results are presented which indicate that the modeling techniques and property values have been proven adequate in predicting peak structure temperatures and entry thermal responses from both an ambient and cold soak condition of an FRSI covered aluminum structure.

Gallegos, J. J.↗

Effective thermal conductivity of composites with fibre-matrix debonding

Debonding of the fiber-matrix interface is a major cause for the degradation of the mechanical properties and the loss of thermal conductivity of fiber-reinforced composites. This paper discusses two analytical approaches for modeling the thermal conduction problem of composites. One is based on the concept of modeling the thermal barrier by an equivalent heat transfer coefficient at the fiber-matrix interface, as described by Hasselman and Johnson (1987) and Benveniste and Miloh (1986). The other approach, suggested by Hatta and Taya (1986), is by treating a composite with debonded interface as a coated-fiber composite. The major advantage of the latter aproach is that the thickness of the fiber coating can be realistically modeled depending upon the extent of degradation of the composite with the thermal conductivity of the coating as that of air.

Fadale, T. D.↗

Temperature and density structure of the corona and inner solar wind

A combined theoretical and observational approach is used to construct models of the inner corona applicable to quiet open-field regions and coronal holes. These models, extending from the top of the transition region at about 1.003 solar radii out to 3 solar radii, are calculated by employing parameters that describe the coronal geometry and heating. Thermal conductivity, heating, radiative losses, and solar-wind convection are included in the full energy-balance equation, which is integrated to obtain density and temperature models; the unknown heating term is parameterized in terms of the mechanical flux incident at the base and the characteristic scale height for the dissipation of this flux; and boundary conditions at the base are fixed by EUV-derived models of the transition region. Two representative models, one for a coronal hole and one for a quiet region, are examined in detail as illustrative examples of the technique used. It is noted that the models obtained predict most of the pertinent observations reasonably well but significantly overestimate the radio emission at 80 and 160 MHz and significantly underestimate the intensity of the Fe XV line at 284 A, especially in coronal holes.

Kopp, R. A.↗

Controls on the CO2 seasonal cycle

Surface pressure measurement performed by the Viking landers show substantial variations in pressure on seasonal timescales that are characterized by two local minima and two local maxima. These variations have widely been attributed to the seasonal condensation and sublimation of CO2 in the two polar regions. It has been somewhat of a surprise that the amplitude of the minimum and maximum that is dominated by the CO2 cycle in the north was much weaker than the corresponding amplitude of the south-dominated extrema. Another surprise was that the seasonal pressure cycle during years 2 and 3 of the Viking mission was so similar to that for year 1, despite the occurrence of two global dust storms during year 1 and none during years 2 and 3. An energy balance model that incorporates dynamical factors from general circulation model (GCM) runs in which the atmospheric dust opacity and seasonal date were systematically varied was used to model the observed seasonal pressure variations. The energy balance takes account of the following processes in determining the rates of CO2 condensation and sublimation at each longitudinal and latitudinal grid point: solar radiation, infrared radiation from the atmosphere and surface, subsurface heat conduction, and atmospheric heat advection. Condensation rates are calculated both at the surface and in the atmosphere. In addition, the energy balance model also incorporates information from the GCM runs on seasonal redistribution of surface pressure across the globe. Estimates of surface temperature of the seasonal CO2 caps were used to define the infrared radiative losses from the seasonal polar caps. The seasonal pressure variations measured at the Viking lander sites were closely reproduced.

Pollack, J. B.↗

Combustion hot section technology

The overall objective of the Turbine Engine Hot Section Technology Combustion Project is to develop and verify improved and more accurate analysis methods for increasing the ability to design with confidence the combustion system for advanced aircraft turbine engines. The analysis methods developed will be generically applicable to combustion systems and not restricted to one specific engine or manufacturer. This project's approach was to first assess and evaluate existing combustor aerothermal analysis models by means of a contracted effort initiated during FY 1982. This evaluation effort has assessed and quantified known models' strengths and deficiencies. During FY 1984 the Aerothermal Modeling Program, Phase 2 will be initiated, which is expected to have contracted model development efforts in the areas of improved numerical methods for turbulent viscous flows, flow interactions, and fuel spray flow foekd interactions. A Phase 3 effort is planned to address remaining model deficiencies. The primary inhouse effort in this area will be the determination of high pressure flame radiation characteristics in a full annular combustor. This experiment will be conducted in the NASA LeRC High Pressure Facility with the results compiled into a comprehensive flame radiation and liner heat flux model.

Ercegovic, D. B.↗

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of non-thermal electrons and x-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick HI, CaII and MgII transitions using the non-LTE prescription in the program MULTI. Hydrogen ionization rates from x-ray photo-ionization and collisional ionization by non-thermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring.

Hawley, Suzanne L.↗

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane-parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of nonthermal electrons and X-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick H I, Ca II and Mg II transitions using the non-local thermodynamic equilibrium (non-LTE) prescription in the program MULTI. Hydrogen ionization rates from X-ray photoionization and collisional ionization by nonthermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in the pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring. We have computed the chromospheric radiation that results from a range of coronal heating rates, with particular emphasis on the widely observed diagnostic H(alpha). Our conclusion is that the H(alpha) fluxes and profiles actually observed in flares can only be produced under conditions of a low-pressure corona with strong beam heating. Therefore we suggest that H(alpha) in flares is produced primarily at the footprints of newly heated loops where significant evaporation has not yet occurred. As a single loop evolves in time, no matter how strong the heating rate may become, the H(alpha) flux will diminish as the corona becomes denser and hence more effective at stopping the beam. This prediction leads to several observable consequences regarding the spatial and temporal signatures of the X-ray and H(alpha) radiation during flares.

Hawley, Suzanne L.↗

Apparatus for determining thermophysical properties of test specimens

Apparatus is described for directly measuring the quantity square root of pck of a test specimen such as a wind tunnel model where p is density, c is the specific heat and k is the thermal conductivity of the specimen. The test specimen and a reference specimen are simultaneously subjected to the heat from a heat source. A thermocouple is attached to the reference specimen for producing a first electrical analog signal proportional to the heat rate Q that the test specimen is subjected to and an infrared radiometer that is aimed at the test specimen produces a second electrical analog signal proportional to the surface temperature T of the test specimen. An analog-to-digital converter converts the first and second electrical analog signals to digital signals. These digital signals are applied to a computer for determining the quantity.

Creel, T. R., Jr.↗

Theoretical aspects of tidal and planetary wave propagation at thermospheric heights

A simple semiquantitative model is presented which allows analytic solutions of tidal and planetary wave propagation at thermospheric heights. This model is based on perturbation approximation and mode separation. The effects of viscosity and heat conduction are parameterized by Rayleigh friction and Newtonian cooling. Because of this simplicity, one gains a clear physical insight into basic features of atmospheric wave propagation. In particular, we discuss the meridional structures of pressure and horizontal wind (the solutions of Laplace's equation) and their modification due to dissipative effects at thermospheric heights. Furthermore, we solve the equations governing the height structure of the wave modes and arrive at a very simple asymptotic solution valid in the upper part of the thermosphere. That 'system transfer function' of the thermosphere allows one to estimate immediately the reaction of the thermospheric wave mode parameters such as pressure, temperature, and winds to an external heat source of arbitrary temporal and spatial distribution. Finally, the diffusion effects of the minor constituents due to the global wind circulation are discussed, and some results of numerical calculations are presented.

Volland, H.↗

A Model of the Soft X-ray Background as a Blast Wave Viewed from Inside

The suggestion that the soft X-ray background arises in part from the Sun which is inside a large supernova blastwave was examined by models of spherical blastwaves. The models can produce quantitative fits to both surface brightnesses and energy band ratios when t = 10 to the 5th power E sub o = 5 x 10 to the 50th power ergs, and n sub approx. 0.004 cm to the -3 power. The models are generalized by varying the relative importance of factors such as thermal conduction, Coulomb heating of electrons, and external pressure; and to allow the explosions to occur in preexisting cavities with steep density gradients, or by examination of the effects of large obstructions or other anisotrophies in the ambient medium.

Edgar, R. J.↗

The loss and depth of CO2 ice in comet nuclei

An analytical model was developed to simulate the material differentiation of a cometary nucleus composed of water ice, putative unclathrated CO2 ice and silicate dust in specified proportions. Selective sublimation of any free CO2 ice present in a comet would produce a surface layer of water ice and dust overlying the original CO2 rich material. This surface layer reduces the temperature of buried CO2 ice and restricts the outflow of gaseous CO2. On each orbit, water sublimation at smaller heliocentric distances temporarily reduces the thickness of the water ice and dust layer and liberates dust. The model includes the effects of nucleus rotation, arbitrary orientation of the rotation axis, latitude, heat conduction into the deep interior of the nucleus and restriction of CO2 gas outflow by the water ice and dust layer. The effects of the permeability of the surface water ice layer, the nucleus rotation rate, and the latitude were investigated. Comparison of these and similar results with observations could yield information regarding the permeability and chemical composition of cometary material and suggest sampling strategies to minimize fractionation effects.

Fanale, F. P.↗

Rate limits in silicon sheet growth - The connections between vertical and horizontal methods

Meniscus-defined techniques for the growth of thin silicon sheets fall into two categories: vertical and horizontal growth. The interactions of the temperature field and the crystal shape are analyzed for both methods using two-dimensional finite-element models which include heat transfer and capillarity. Heat transfer in vertical growth systems is dominated by conduction in the melt and the crystal, with almost flat melt/crystal interfaces that are perpendicular to the direction of growth. The high axial temperature gradients characteristic of vertical growth lead to high thermal stresses. The maximum growth rate is also limited by capillarity which can restrict the conduction of heat from the melt into the crystal. In horizontal growth the melt/crystal interface stretches across the surface of the melt pool many times the crystal thickness, and low growth rates are achievable with careful temperature control. With a moderate axial temperature gradient in the sheet a substantial portion of the latent heat conducts along the sheet and the surface of the melt pool becomes supercooled, leading to dendritic growth. The thermal supercooling is surpressed by lowering the axial gradient in the crystal; this configuration is the most desirable for the growth of high quality crystals. An expression derived from scaling analysis relating the growth rate and the crucible temperature is shown to be reliable for horizontal growth.

Thomas, Paul D.↗

Microgravity ignition experiment

The purpose of this project is to develop a flight ready apparatus of the microgravity ignition experiment for the GASCan 2 program. This involved redesigning, testing, and making final modifications to the existing apparatus. The microgravity ignition experiment is intended to test the effect of microgravity on the time to ignition of a sample of alpha-cellulose paper. An infrared heat lamp is used to heat the paper sample within a sealed canister. The interior of the canister was redesigned to increase stability and minimize conductive heat transfer to the sample. This design was fabricated and tested and a heat transfer model of the paper sample was developed.

Motevalli, Vahid↗

Quantitative Global Heat Transfer in a Mach-6 Quiet Tunnel

This project developed quantitative methods for obtaining heat transfer from temperature sensitive paint (TSP) measurements in the Mach-6 quiet tunnel at Purdue, which is a Ludwieg tube with a downstream valve, moderately-short flow duration and low levels of heat transfer. Previous difficulties with inferring heat transfer from TSP in the Mach-6 quiet tunnel were traced to (1) the large transient heat transfer that occurs during the unusually long tunnel startup and shutdown, (2) the non-uniform thickness of the insulating coating, (3) inconsistencies and imperfections in the painting process and (4) the low levels of heat transfer observed on slender models at typical stagnation temperatures near 430K. Repeated measurements were conducted on 7 degree-half-angle sharp circular cones at zero angle of attack in order to evaluate the techniques, isolate the problems and identify solutions. An attempt at developing a two-color TSP method is also summarized.

Sullivan, John P.↗

MHD Modeling of Coronal Loops: the Transition Region Throat

Context. The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross-sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. Aims. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in terms of steady heating models. Methods. We study the area response with a time-dependent 2D magnetohydrodynamic (MHD) loop model, including the description of the expanding magnetic field, coronal heating and losses by thermal conduction, and radiation from optically thin plasma. We run a simulation for a loop 50 Mm long and quasi-statically heated to about 4 millikelvin. Results. We find that the area can change substantially with the quasi-steady heating rate, e.g., by approx. 40% at 0.5 millikelvin as the loop temperature varies between 1 millikelvin and 4 millikelvin, and, therefore, affects the interpretation of the differential emission measure vs. temperature (DEM(T)) curves.

magnetohydrodynamics (MHD)↗