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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Optical Probe of the Superconducting Normal Mixed State in a Magnetic Penetration Thermometer

Using ultraviolet photon pulses, we have probed the internal behavior of a molybdenum-gold Magnetic Penetration Thermometer (MPT) that we designed for x-ray microcalorimetry. In this low-temperature detector, the diamagnetic response of a superconducting MoAu bilayer is used to sense temperature changes in response to absorbed photons. We have previously described an approximate model that explains the high responsivity of the detector to temperature changes as a consequence of a Meissner transition of the molybdenum-gold film in the magnetic field applied by the superconducting circuit used to bias the detector. We compare measurements of MPT heat capacity and thermal conductance, derived from UV photon pulse data, to our model predictions for the thermodynamic properties of the sensor and for the electron cooling obtained by quasiparticle recombination. Our data on electron cooling power is also relevant to the operation of other superconducting detectors, such as Microwave Kinetic Inductance Detectors.

superconductivity↗
NASA NTRS1993

A finite element analysis of the freeze/thaw behavior of external artery heat pipes

A two-dimensional finite element model was used to determine the freeze/thaw characteristics of an external artery heat pipe. During startup, the working fluid, which was located in the liquid channel and the circumferential wall grooves, experienced a phase transformation from a solid to a liquid state. The transient heat conduction equations with moving interfacial conditions were solved using the appropriate initial boundary conditions. The modelling results include the cross-sectional temperature distribution and the interfacial or melt front position as a function of time. A fixed grid approach was adopted in the model for the phase-change process during thawing of frozen working fluid. The interfacial position between the liquid and solid regions was found by balancing the latent heat caused by interfacial movement with the heat addition or extraction at the related grid points.

Lu, X. J.↗
NASA NTRS1992

TOPEX Microwave Radiometer - Thermal design verification test and analytical model validation

The testing of the TOPEX Microwave Radiometer (TMR) is described in terms of hardware development based on the modeling and thermal vacuum testing conducted. The TMR and the vacuum-test facility are described, and the thermal verification test includes a hot steady-state segment, a cold steady-state segment, and a cold survival mode segment totalling 65 hours. A graphic description is given of the test history which is related temperature tracking, and two multinode TMR test-chamber models are compared to the test results. Large discrepancies between the test data and the model predictions are attributed to contact conductance, effective emittance from the multilayer insulation, and heat leaks related to deviations from the flight configuration. The TMR thermal testing/modeling effort is shown to provide technical corrections for the procedure outlined, and the need for validating predictive models is underscored.

Lin, Edward I.↗
NASA NTRS1995

A Network Model for the Effective Thermal Conductivity of Rigid Fibrous Refractory Insulations

A procedure is described for computing the effective thermal conductivity of a rigid fibrous refractory insulation. The insulation is modeled as a 3-dimensional Cartesian network of thermal conductance. The values and volume distributions of the conductance are assigned to reflect the physical properties of the insulation, its constituent fibers, and any permeating gas. The effective thermal conductivity is computed by considering the simultaneous energy transport by solid conduction, gas conduction and radiation through a cubic volume of model insulation; thus the coupling between heat transfer modes is retained (within the simplifications inherent to the model), rather than suppressed by treating these heat transfer modes as independent. The model takes into account insulation composition, density and fiber anisotropy, as well as the geometric and material properties of the constituent fibers. A relatively good agreement, between calculated and experimentally derived thermal conductivity values, is obtained for a variety of rigid fibrous insulations.

Marschall, Jochen↗
NASA NTRS1996

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.↗
NASA NTRS1982

Volcanism and igneous processes in small icy satellites

Evidence for the production of endogenic processes of the small Saturnian satellites by igneous activity of low melting point NH3-H2O magma heated radiogenically is presented. An initial state of the inner satellites is modeled as a homogeneous mixture of particulate silicates and ices. Conductive transport is assumed to have either combined with heat capacity in a mixture of crystalline phases, or to have been reduced by the presence of dust, defects, vitreous or amorphous phases, and clathrates. The extent to which a eutectic melt could form is calculated in terms of the volume percentage of melt, with all heating above 175 K going toward overcoming the latent heat of the eutectic mix. Mimas and Enceladus were treated as free from radiogenic heat sources, which were significant for Tethys, Dione, Rhea, and Iapetus. The migrations of inner materials to form surface structures and seal off the interiors of the moons are described.

Stevenson, D. J.↗
NASA NTRS1989

A numerical analysis of the effects of conjugate heat transfer, vapor compressibility, and viscous dissipation in heat pipes

The effects of conjugate heat transfer, vapor compressibility, and viscous dissipation in heat pipes are discussed. The accuracy of the partially parabolic versus the elliptic presentation of the governing equations is also examined. The results show that the axial wall conduction has a tendency to make the temperature distribution more uniform for heat pipes with large ratios of pipe wall to effective liquid-wick thermal conductivity. The compressible and incompressible models show very close agreement for the total pressure drop, while the local pressure variations along the heat pipe are quite different for these two models when the radial Reynolds number at the interface is high.

Faghri, Amir↗
NASA NTRS1994

Improved Slip Casting Of Ceramic Models

Improved technique of investment slip casting developed for making precise ceramic wind-tunnel models. Needed in wind-tunnel experiments to verify predictions of aerothermodynamical computer codes. Ceramic materials used because of their low heat conductivities and ability to survive high temperatures. Present improved slip-casting technique enables casting of highly detailed models from aqueous or nonaqueous solutions. Wet shell molds peeled off models to ensure precise and undamaged details. Used at NASA Langley Research Center to form superconducting ceramic components from nonaqueous slip solutions. Technique has many more applications when ceramic materials developed further for such high-strength/ temperature components as engine parts.

Buck, Gregory M.↗
NASA NTRS1975

Thermal analytical model of a 30 cm mercury ion thruster

A lumped parameter thermal nodal network has been developed for a 30 cm engineering model mercury ion thruster. The network consists of approximately 100 nodes coded in SINDA format for use on the Univac 1106/1108 computer. This model takes into account internal dissipation, radiation, and conduction as well as environmental heating. A series of tests were performed at NASA-Lewis Research Center to simulate a wide range of thermal environments on an operating 30 cm thruster, instrumented to measure the temperature distribution within the thruster. The results of these tests were used to calibrate the analytical model. Presented in this paper are a description of the analytical model along with comparisons between analytical and experimental results for the various operating conditions.

Oglebay, J. C.↗
NASA NTRS2017

Numerical Simulation of Illumination and Thermal Conditions at the Lunar Poles Using LOLA DTMs

We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.

Glaser, P.↗
NASA NTRS1981

Glass shell manufacturing in space

Highly-uniform, hollow glass spheres (shells), which are used for inertial confinement fusion targets, were formed from metal-organic gel powder feedstock in a vertical furnace. As a result of the rapid pyrolysis caused by the furnace, the gel is transformed to a shell in five distinct stages: (a) surface closure of the porous gel; (b) generation of a closed-cell foam structure in the gel; (c) spheridization of the gel and further expansion of the foam; (d) coalescence of the closed-cell foam to a single-void shell; and (e) fining of the glass shell. The heat transfer from the furnace to the falling gel particle was modeled to determine the effective heating rate of the gel. The model predicts the temperature history for a particle as a function of mass, dimensions, specific heat, and absorptance as well as furnace temperature profile and thermal conductivity of the furnace gas. A model was developed that predicts the gravity-induced degradation of shell concentricity in falling molten shells as a function of shell characteristics and time.

Downs, R. L.↗
NASA NTRS1974

Method for determining thermo-physical properties of specimens

The square root of the product of thermophysical properties q, c and k, where p is density, c is specific heat and k is thermal conductivity, is determined directly on a test specimen such as a wind tunnel model. The test specimen and a reference specimen of known specific heat are positioned at a given distance from a heat source. The specimens are provided with a coating, such as a phase change coating, to visually indicate that a given temperature was reached. A shutter interposed between the heat source and the specimens is opened and a motion picture camera is actuated to provide a time record of the heating step. The temperature of the reference specimen is recorded as a function of time. The heat rate to which both the test and reference specimens were subjected is determined from the temperature time response of the reference specimen by the conventional thin-skin calorimeter equation.

Jones, R. A.↗
NASA NTRS1995

Influence of the vapor flux on temperature, density, and abundance distributions in a multicomponent, porous, icy body

We calculated the vapor flux of the icy components in the surface layer of a porous, short-period, Jupiter-class comet, in order to investigate the relationship of the observed relative molecular abundances in the coma with those in the nucleus. The model assumes a body containing one major ice component (H20) and up to three minor components of higher volatility (e.g., CO, CO2, CH3OH). The body's porous structure is modeled as a bundle of tubes with a given tortuosity and initially a constant pore diameter. The mass and energy equations for the different volatiles are solved simultaneously under appropiate boundary conditions. Heat is conducted by the matrix and carried by the vapors. The one-dimensional model includes radially inward and outward flowing vapor within the body, complete depletion of less volatile ices in outer layers, and recondensation of vapor in deeper, coller layers. As a result, we obtain the temperature and abundance distribution in the nucleus and the gas flux into the interior and into the coma for each of the volatiles at various positions in the orbit. The ratio of the gas flux of minor volatiles to that of H2) in the coma varies by several orders of magnitude throughout the orbit. Thus, the relative abundances of species observed in the coma are in most cases not the same as those in the nucleus. Results also indicate that it will be impossible to determine the relative abundances of ices more volatile than water from samples taken a few meters below the surface during a comet rendezvous mission. We made calculations for a wide range of different parameters, such as porosity, pore radius, and thermal conductivity of the matrix. To introduce the model we present typical results for a dust-free comet.

Benkhoff, J.↗
NASA NTRS1998

Reduction and Analysis of Phosphor Thermography Data With the IHEAT Software Package

Detailed aeroheating information is critical to the successful design of a thermal protection system (TPS) for an aerospace vehicle. This report describes NASA Langley Research Center's (LaRC) two-color relative-intensity phosphor thermography method and the IHEAT software package which is used for the efficient data reduction and analysis of the phosphor image data. Development of theory is provided for a new weighted two-color relative-intensity fluorescence theory for quantitatively determining surface temperatures on hypersonic wind tunnel models; an improved application of the one-dimensional conduction theory for use in determining global heating mappings; and extrapolation of wind tunnel data to flight surface temperatures. The phosphor methodology at LaRC is presented including descriptions of phosphor model fabrication, test facilities and phosphor video acquisition systems. A discussion of the calibration procedures, data reduction and data analysis is given. Estimates of the total uncertainties (with a 95% confidence level) associated with the phosphor technique are shown to be approximately 8 to 10 percent in the Langley's 31-Inch Mach 10 Tunnel and 7 to 10 percent in the 20-Inch Mach 6 Tunnel. A comparison with thin-film measurements using two-inch radius hemispheres shows the phosphor data to be within 7 percent of thin-film measurements and to agree even better with predictions via a LATCH computational fluid dynamics solution (CFD). Good agreement between phosphor data and LAURA CFD computations on the forebody of a vertical takeoff/vertical lander configuration at four angles of attack is also shown. In addition, a comparison is given between Mach 6 phosphor data and laminar and turbulent solutions generated using the LAURA, GASP and LATCH CFD codes. Finally, the extrapolation method developed in this report is applied to the X-34 configuration with good agreement between the phosphor extrapolation and LAURA flight surface temperature predictions. The phosphor process outlined in the paper is believed to provide the aerothermodynamic community with a valuable capability for rapidly obtaining (4 to 5 weeks) detailed heating information needed in TPS design.

Merski, N. Ronald↗
NASA NTRS1992

Properties of Graphite Fiber Reinforced Copper Matrix Composites for Space Power Applications

The thermal and mechanical properties of pitch-based graphite fiber reinforced copper matrix (Gr/Cu) composites usable for space applications such as radiator fins were investigated. Thermal conductivity was measured as a function of fiber volume fraction and architecture. Results showed for unidirectional P-100 Gr/Cu composites, the longitudinal thermal conductivity was nearly independent of fiber volume fraction. Transverse thermal conductivities (perpendicular to the fibers) were strongly affected by the fiber volume fraction with higher volume fractions resulting in lower thermal conductivities. The effect of architecture on thermal conductivity followed the cosine squared law for simple architectures. Insufficient data are available currently to model more complex architectures, but adding fibers in the direction of the heat flow increases the thermal conductivity as low conductivity plies are supplemented by high conductivity plies. Thermal expansion tests were conducted on the Gr fibers and Gr/Cu composites. The results show a considerable thermal expansion mismatch between the fibers and the Cu matrix. The longitudinal thermal expansion showed a strong dependence on the architecture of the Gr/Cu composites. The composites also show a thermal expansion hysteresis. The hysteresis was eliminated by an engineered interface. Mechanical testing concentrated on the dynamic modulus and strength of the composites. The dynamic modulus of the Gr/Cu composites was 305 GPa up to 400 C, a value equivalent to Be. The strengths of the composites were less than expected, but this is attributed to the poor bond across the interface between the Gr fibers and Cu matrix. Testing of composites with an engineered interface is expected to yield strengths nearer the values predicted by the rule of mixtures.

Ellis, David L.↗
NASA NTRS1985

Solution algorithms for nonlinear transient heat conduction analysis employing element-by-element iterative strategies

The particular problems investigated in the present study arise from nonlinear transient heat conduction. One of two types of nonlinearities considered is related to a material temperature dependence which is frequently needed to accurately model behavior over the range of temperature of engineering interest. The second nonlinearity is introduced by radiation boundary conditions. The finite element equations arising from the solution of nonlinear transient heat conduction problems are formulated. The finite element matrix equations are temporally discretized, and a nonlinear iterative solution algorithm is proposed. Algorithms for solving the linear problem are discussed, taking into account the form of the matrix equations, Gaussian elimination, cost, and iterative techniques. Attention is also given to approximate factorization, implementational aspects, and numerical results.

Winget, J. M.↗
NASA NTRS1990

LOX vaporization in high-pressure, hydrogen-rich gas

LOX droplet vaporization in high-pressure hydrogen-rich gas is analyzed, with special attention to thermodynamic effects which compel the surface to heat to the critical state and to supercritical vaporization processes on heating to criticality. Subcritical vaporization is modeled using a quasi-steady diffusion-controlled gas-phase transport formulation coupled to an effective-conductivity internal-energy-transport model accounting for circulation effects. It is demonstrated how the droplet surface might heat to the critical state, for ambient pressures slightly greater than the critical pressure of oxygen, such that the bulk of propellant within the droplet remains substantially below the critical mixing temperature.

Litchford, Ron J.↗
NASA NTRS1985

A theoretical study of non-adiabatic surface effects for a model in the NTF cryogenic wind tunnel

A theoretical analysis was made of the severity and effect of nonadiabatic surface conditions for a model in the NTF cryogenic wind tunnel. The nonadiabatic condition arises from heaters that are used to maintain a constant thermal environment for instrumentation internal to the model. The analysis was made for several axi-symmetric representations of a fuselage cavity, using a finite element heat conduction code. Potential flow and boundary layer codes were used to calculate the convection condition for the exterior surface of the model. The results of the steady state analysis show that it is possible to maintain the surface temperature very near the adiabatic value, with the judicious use of insulating material. Even for the most severe nonadiabatic condition studied, the effects on skin friction drag and displacement thickness were only marginally significant. The thermal analysis also provided an estimate of the power required to maintain a specified cavity temperature.

Macha, J. M.↗