Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal boundary conductance”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Thermal conductivity and Kapitza boundary resistance of He-4 near T(gamma)

The preliminary results of high resolution measurements of the thermal conductivity of He-4 in a cell with a nominal 1 mm gap are reported and compared with previous results obtained with a 2 mm gap. The temperature dependence of the conductivity was found to be similar for both cells, indicating that size and boundary effects are well controlled. The detailed behavior of the Kapitza resistance just below T(gamma) is also reported.

Chui, T. C. P.↗
NASA NTRS1980

Mathematical modeling of moving boundary problems in thermal energy storage

The capability for predicting the performance of thermal energy storage (RES) subsystems and components using PCM's based on mathematical and physical models is developed. Mathematical models of the dynamic thermal behavior of (TES) subsystems using PCM's based on solutions of the moving boundary thermal conduction problem and on heat and mass transfer engineering correlations are also discussed.

Solomon, A. D.↗
NASA NTRS1984

Laboratory studies of baroclinic instability at small Richardson number

As part of the support program for the Atmospheric General Circulation Experiment, laboratory studies of baroclinic and other convective instabilities were performed for a thin layer of fluid between thermally conducting horizontal discs. There were three types of modes identified. The first has a spiral-arm appearance, and exists for large enough horizontal thermal forcing, weak enough static stability, and large enough rotation. The source of this wave is shown to be the Eady mode of instability. The second mode is due to convective instability in the thermal boundary layers which exist due to the thermally conducting horizontal boundaries. Finally, for strong enough negative static stability, thermal convection of the Benard type appears. The most significant result is that the symmetric (Solberg) mode was not found, even though the infinite-plane theory predicts this mode under certain experimental conditions.

Miller, T. L.↗
NASA NTRS2013

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant↗
NASA NTRS1993

Overturn of magma ocean ilmenite cumulate layer: Implications for lunar magmatic evolution and formation of a lunar core

We explore a model for the chemical evolution of the lunar interior that explains the origin and evolution of lunar magmatism and possibly the existence of a lunar core. A magma ocean formed during accretion differentiates into the anorthositic crust and chemically stratified cumulate mantle. The cumulative mantle is gravitationally unstable with dense ilmenite cumulate layers overlying olivine-orthopyroxene cumulates with Fe/Mg that decreases with depth. The dense ilmenite layer sinks to the center of the moon forming the core. The remainder of the gravitationally unstable cumulate pile also overturns. Any remaining primitive lunar mantle rises to its level of neutral buoyancy in the cumulate pile. Perhaps melting of primitive lunar mantle due to this decompression results in early lunar Mg-rich magmatism. Because of its high concentration of incompatible heat producing elements, the ilmenite core heats the overlying orthopyroxene-bearing cumulates. As a conductively thickening thermal boundary layer becomes unstable, the resulting mantle plumes rise, decompress, and partially melt to generate the mare basalts. This model explains both the timing and chemical characteristics of lunar magmatism.

Hess, P. C.↗
NASA NTRS2012

Multiscale Modeling of Grain Boundaries in ZrB2: Structure, Energetics, and Thermal Resistance

A combination of ab initio, atomistic and finite element methods (FEM) were used to investigate the structures, energetics and lattice thermal conductance of grain boundaries for the ultra high temperature ceramic ZrB2. Atomic models of idealized boundaries were relaxed using density functional theory. Information about bonding across the interfaces was determined from the electron localization function. The Kapitza conductance of larger scale versions of the boundary models were computed using non-equilibrium molecular dynamics. The interfacial thermal parameters together with single crystal thermal conductivities were used as parameters in microstructural computations. FEM meshes were constructed on top of microstructural images. From these computations, the effective thermal conductivity of the polycrystalline structure was determined.

Lawson, John W.↗
NASA NTRS1993

A direct approach to finding unknown boundary conditions in steady heat conduction

The capability of the boundary element method (BEM) in determining thermal boundary conditions on surfaces of a conducting solid where such quantities are unknown was demonstrated. The method uses a non-iterative direct approach in solving what is usually called the inverse heat conduction problem (IHCP). Given any over-specified thermal boundary conditions such as a combination of temperature and heat flux on a surface where such data is readily available, the algorithm computes the temperature field within the object and any unknown thermal boundary conditions on surfaces where thermal boundary values are unavailable. A two-dimensional, steady-state BEM program was developed and was tested on several simple geometries where the analytic solution was known. Results obtained with the BEM were in excellent agreement with the analytic values. The algorithm is highly flexible in treating complex geometries, mixed thermal boundary conditions, and temperature-dependent material properties and is presently being extended to three-dimensional and unsteady heat conduction problems. The accuracy and reliability of this technique was very good but tended to deteriorate when the known surface conditions were only slightly over-specified and far from the inaccessible surface.

Martin, Thomas J.↗
NASA NTRS1973

Martian thermal boundary layers: Subhourly variations induced by radiative-conductive heat transfer within the dust-laden atmosphere-ground system

Radiative-conductive heat transfer has been investigated for the ground-atmosphere system of the planet Mars. The basic goal was the quantitative determination of time dependent vertical distributions of temperature and static stability for Southern-Hemispheric summer season and middle and polar latitudes, for both dust-free and dust-laden atmospheric conditions. The numerical algorithm which models at high spatial and temporal resolution the thermal energy transports in the dual ground-atmosphere system, is based on solution of the applicable heating rate equation, including radiative and molecular-conductive heat transport terms. The two subsystems are coupled by an internal thermal boundary condition applied at the ground-atmosphere interface level. Initial data and input parameters are based on Mariner 4, 6, 7, and 9 measurements and the JPL Mars Scientific Model. Numerical experiments were run for dust-free and dust-laden conditions in the midlatitudes, as well as ice-free and ice-covered polar regions. Representative results and their interpretation are presented. Finally, the theoretical framework of the generalized problem with nonconservative Mie scattering and explicit thermal-convective heat transfer is formulated, and applicable solution algorithms are outlined.

Pallmann, A. J.↗
NASA NTRS1974

Effect of wall temperature spatial variations on aerodynamic heating to an actively cooled aircraft

The design of cooled skin panels, which requires consideration of many parameters such as heat load, skin thickness, panel size, tube spacing, tube size, temperature limits, coolant flow, etc., must minimize weight and thermal stresses. A design procedure for such panels has been developed by Bell Aerospace Company (Helenbrook and Anthony, 1971). This design procedure assumes a uniform heating rate to the panel, i.e., the variation of aerodynamic heating rate due to the cyclic wall temperature distribution, which results from the discrete coolant tube locations in such a panel, has been neglected. This study was made to determine if these variations were sufficiently large to cause variations in aerodynamic heating rate which would require coupling of boundary layer and conductive thermal analysis solutions.

Jones, R. A.↗
NASA NTRS1989

Characteristic parameters of superconductor-coolant interaction including high Tc current density limits

In the area of basic mechanisms of helium heat transfer and related influence on super-conducting magnet stability, thermal boundary conditions are important constraints. Characteristic lengths are considered along with other parameters of the superconducting composite-coolant system. Based on helium temperature range developments, limiting critical current densities are assessed at low fields for high transition temperature superconductors.

Frederking, T. H. K.↗
NASA NTRS2000

Convection in Slab and Spheroidal Geometries

Three-dimensional numerical simulations of compressible turbulent thermally driven convection, in both slab and spheroidal geometries, are reviewed and analyzed in terms of velocity spectra and mixing-length theory. The same ideal gas model is used in both geometries, and resulting flows are compared. The piecewise-parabolic method (PPM), with either thermal conductivity or photospheric boundary conditions, is used to solve the fluid equations of motion. Fluid motions in both geometries exhibit a Kolmogorov-like k(sup -5/3) range in their velocity spectra. The longest wavelength modes are energetically dominant in both geometries, typically leading to one convection cell dominating the flow. In spheroidal geometry, a dipolar flow dominates the largest scale convective motions. Downflows are intensely turbulent and up drafts are relatively laminar in both geometries. In slab geometry, correlations between temperature and velocity fluctuations, which lead to the enthalpy flux, are fairly independent of depth. In spheroidal geometry this same correlation increases linearly with radius over the inner 70 percent by radius, in which the local pressure scale heights are a sizable fraction of the radius. The effects from the impenetrable boundary conditions in the slab geometry models are confused with the effects from non-local convection. In spheroidal geometry nonlocal effects, due to coherent plumes, are seen as far as several pressure scale heights from the lower boundary and are clearly distinguishable from boundary effects.

Porter, David H.↗
NASA NTRS1992

Chemical differentiation of a convecting planetary interior: Consequences for a one-plate planet such as Venus

Chemically depleted mantle forming a buoyant, refractory layer at the top of the mantle can have important implications for the evolution of the interior and surface. On Venus, the large apparent depths of compensation for surface topographic features might be explained if surface topography were supported by variations in the thickness of a 100-200 km thick chemically buoyant mantle layer or by partial melting in the mantle at the base of such a layer. Long volcanic flows seen on the surface may be explained by deep melting that generates low-viscosity MgO-rich magmas. The presence of a shallow refractory mantle layer may also explain the lack of volcanism associated with rifting. As the depleted layer thickens and cools, it becomes denser than the convecting interior and the portion of it that is hot enough to flow can mix with the convecting mantle. Time dependence of the thickness of a depleted layer may create episodic resurfacing events as needed to explain the observed distribution of impact craters on the venusian surface. We consider a planetary structure consisting of a crust, depleted mantle layer, and a thermally and chemically well-mixed convecting mantle. The thermal evolution of the convecting spherical planetary interior is calculated using energy conservation: the time rate of change of thermal energy in the interior is equated to the difference in the rate of radioactive heat production and the rate of heat transfer across the thermal boundary layer. Heat transfer across the thermal boundary layer is parameterized using a standard Nusselt number-Rayleigh number relationship. The radioactive heat production decreases with time corresponding to decay times for the U, Th, and K. The planetary interior cools by the advection of hot mantle at temperature T interior into the thermal boundary layer where it cools conductively. The crust and depleted mantle layers do not convect in our model so that a linear conductive equilibrium temperature distribution is assumed. The rate of melt production is calculated as the product of the volume flux of mantle into the thermal boundary layer and the degree of melting that this mantle undergoes. The volume flux of mantle into the thermal boundary layer is simply the heat flux divided by amount of heat lost in cooling mantle to the average temperature in the thermal boundary layer. The degree of melting is calculated as the temperature difference above the solidus, divided by the latent heat of melting. A maximum degree of melting is prescribed corresponding to the maximum amount of basaltic melt that the mantle can initially generate. As the crust thickens, the pressure at the base of the crust becomes high enough and the temperature remains low enough for basalt to transform to dense eclogite.

Parmentier, E. M.↗
NASA NTRS1972

The diurnal heat budget of the thermosphere

Detailed numerical calculations of thermospheric heat sources and sinks are presented and their relative importance is discussed in reference to the energy balance phenomena of the neutral atmosphere. It is shown that the thermal energy available from the absorption in the Schumann-Runge continuum leading to photodissociation of O2 is by far the largest energy source in the lower thermosphere. Other sources of varying importance in different altitude ranges are: (1) energy from photoelectrons, (2) energy exchange from thermal plasma, (3) chemical reactions (ion-electron dissociative recombination) energy gain, and (4) kinetic and dissipative energy associated with the neutral wind. The energy sinks of importance are thermal conduction at the lower boundary (120 km) and radiative cooling of atomic oxygen. It is shown that the combined energy from processes 2 to 4 constitute only a small fraction of the total energy available from photoelectrons and are in phase with the latter. It is suggested that a choice of a lower boundary much below 120 km, e.g. near the mesopause level (90 km), should be more appropriate for general thermospheric studies.

Chandra, S.↗
NASA NTRS1973

The diurnal heat budget of the thermosphere.

Detailed numerical calculations of thermospheric heat sources and sinks are presented and their relative importance is discussed in reference to the energy balance phenomena of the neutral atmosphere. It is shown that the thermal energy available from the absorption in the Schumann-Runge continuum leading to photo-dissociation of O2 is by far the largest energy source in the lower thermosphere. Other sources of varying importance in different altitude ranges are: (1) energy from photoelectrons; (2) energy exchange from thermal plasma; (3) chemical reaction (ion-electron dissociative recombination) energy gain; (4) kinetic and dissipative energy associated with the neutral wind. The energy sinks of importance are (1) thermal conduction at the lower boundary (120 km); and (2) radiative cooling of atomic oxygen. It is shown that the combined energy from processes 2 to 4 constitutes only a small fraction of the total energy available from photoelectrons and is in phase with the latter. These secondary sources (processes 2 to 4), therefore, do not constitute a significant energy source and their contribution can be simply incorporated into photoelectron energy (process 1) by defining an effective photoionization heating efficiency.

Chandra, S.↗
NASA NTRS1975

Solution of the equation of heat conduction with time dependent sources: Programmed application to planetary thermal history

A computer program (Program SPHERE) solving the inhomogeneous equation of heat conduction with radiation boundary condition on a thermally homogeneous sphere is described. The source terms are taken to be exponential functions of the time. Thermal properties are independent of temperature. The solutions are appropriate to studying certain classes of planetary thermal history. Special application to the moon is discussed.

Conel, J. E.↗