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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Boundary formulations for shape sensitivity of temperature dependent conductivity problems

Used in concert with the Kirchhoff transformation, implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solids with temperature dependent thermal conductivity is shown to generate an accurate and economical approach for computation of shape sensitivities. For problems with specified temperature and heat flux boundary conditions, a linear problem results for both the analysis and sensitivity analysis. In problems with either convection or radiation boundary conditions, a nonlinear problem is generated. Several iterative strategies are presented for the solution of the resulting sets of nonlinear equations and the computational performances examined in detail. Multizone analysis and zone condensation strategies are demonstrated to provide substantive computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear example problems is presented that have closed form solutions. Exact analytical expressions for the shape sensitivities associated with these problems are developed and these are compared with the sensitivities computed using the boundary element formulation.

Kane, James H.↗
DOE OSTI2023

Numerical Simulations of High Prandtl Number Liquid Jets Impinging on a Flat Plate

In this work, 3D simulations of oil jets impinging on a flat, heated wall are presented. The numerical setup uses the Volume of Fluid (VoF) method to model the two-phase flow. A careful grid definition across the liquid film, along with the use of the Conjugate Heat Transfer (CHT) approach allowed local heat transfer to be solved with fine resolution at the wall. Variations of liquid flow rate, liquid temperature and surface temperature allow to cover a wide range of local Reynolds and Prandtl numbers (226 < Re < 2850, 77 < Pr < 161). Resulting surface-averaged heat transfer compares very well with experimental measurements conducted in a previous study. In-depth analysis of the flow has identified expected features from the literature. In particular, the impact of jet axial velocity profiles on the heat transfer distribution in the stagnation zone was clearly stated. The increase in heat transfer when warming the liquid film was also reproduced and explained by a decrease in oil viscosity and an increase in film velocity. All those effects were taken into account in correlations for stagnation and local values of Nusselt number. A grid sensitivity study was also conducted, showing that if the grid solving the thermal boundary layer in the stagnation zone can be coarsened without impacting local and surface-averaged predictions of heat transfer, a minimum resolution (2 to 3 cells) within the thermal boundary layer is however required for an accurate prediction of heat transfer.

conjugate heat transfer↗
NASA NTRS1974

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

Solar electric propulsion system thermal analysis

Thermal control elements applicable to the solar electric propulsion stage are discussed along with thermal control concepts. Boundary conditions are defined, and a thermal analysis was conducted with special emphasis on the power processor and equipment compartment thermal control system. Conclusions and recommendations are included.

Source record↗
NASA NTRS1993

A spectral method determination of the first critical Rayleigh number for a low-Prandtl number crystal melt in a cylindrical container

The onset of laminar axisymmetric Rayleigh-Benard convection is investigated for a low-Prandtl number liquid metal in a cylindrical container. All surfaces are considered to be solid and no-slip. Two separate cases are examined for the thermal boundary conditions at the side wall, one with conducting and the other with insulated surface. The governing Boussinesq system is first perturbed and then simplified by introducing a Stokes stream function. Subsequently, a Chebyshev Galerkin spectral model is employed to reduce the simplified system to a system of first-order nonlinear ordinary differential equations. A local stability analysis determines the two values of the first critical Rayleigh number, Ra(sub cl), for the insulated and conducting side walls. As expected, the conducting Ra(sub cl) value of 2882.5 obtained from the present approach exceeded the corresponding insulated Ra(sub cl) value of 2331.6. For the insulated case, an earlier study using a different numerical approach suggests that Ra(sub cl) = 2261.9, while an experimental study measured Ra(sub cl) = 2700.

Dietz, C. M., Jr.↗
NASA NTRS2017

High Count-Rate Study of Two TES X-Ray Microcalorimeters With Different Transition Temperatures

We have developed transition-edge sensor (TES) microcalorimeter arrays with high count-rate capability and high energy resolution to carry out x-ray imaging spectroscopy observations of various astronomical sources and the Sun. We have studied the dependence of the energy resolution and throughput (fraction of processed pulses) on the count rate for such microcalorimeters with two different transition temperatures T(sub c). Devices with both transition temperatures were fabricated within a single microcalorimeter array directly on top of a solid substrate where the thermal conductance of the microcalorimeter is dependent upon the thermal boundary resistance between the TES sensor and the dielectric substrate beneath. Because the thermal boundary resistance is highly temperature dependent, the two types of device with different T(sub c)(sup s) had very different thermal decay times, approximately one order of magnitude different. In our earlier report, we achieved energy resolutions of 1.6 and 2.eV at 6 keV from lower and higher T(sub c) devices, respectively, using a standard analysis method based on optimal filtering in the low flux limit. We have now measured the same devices at elevated x-ray fluxes ranging from 50 Hz to 1000 Hz per pixel. In the high flux limit, however, the standard optimal filtering scheme nearly breaks down because of x-ray pile-up. To achieve the highest possible energy resolution for a fixed throughput, we have developed an analysis scheme based on the socalled event grade method. Using the new analysis scheme, we achieved 5.0 eV FWHM with 96 Percent throughput for 6 keV x-rays of 1025 Hz per pixel with the higher T(sub c) (faster) device, and 5.8 eV FWHM with 97 Percent throughput with the lower T(sub c) (slower) device at 722 Hz.

Lee, Sang-Jun↗
NASA NTRS2013

Copper-Multiwall Carbon Nanotubes and Copper-Diamond Composites for Advanced Rocket Engines

This paper reports on the research effort to improve the thermal conductivity of the copper-based alloy NARloy-Z (Cu-3 wt.%Ag-0.5 wt.% Zr), the state-of-the-art alloy used to make combustion chamber liners in regeneratively-cooled liquid rocket engines, using nanotechnology. The approach was to embed high thermal conductivity multiwall carbon nanotubes (MWCNTs) and diamond (D) particles in the NARloy-Z matrix using powder metallurgy techniques. The thermal conductivity of MWCNTs and D have been reported to be 5 to 10 times that of NARloy-Z. Hence, 10 to 20 vol. % MWCNT finely dispersed in NARloy-Z matrix could nearly double the thermal conductivity, provided there is a good thermal bond between MWCNTs and copper matrix. Quantum mechanics-based modeling showed that zirconium (Zr) in NARloy-Z should form ZrC at the MWCNT-Cu interface and provide a good thermal bond. In this study, NARloy-Z powder was blended with MWCNTs in a ball mill, and the resulting mixture was consolidated under high pressure and temperature using Field Assisted Sintering Technology (FAST). Microstructural analysis showed that the MWCNTs, which were provided as tangles of MWCNTs by the manufacturer, did not detangle well during blending and formed clumps at the prior particle boundaries. The composites made form these powders showed lower thermal conductivity than the base NARloy-Z. To eliminate the observed physical agglomeration, tangled multiwall MWCNTs were separated by acid treatment and electroless plated with a thin layer of chromium to keep them separated during further processing. Separately, the thermal conductivities of MWCNTs used in this work were measured, and the results showed very low values, a major factor in the low thermal conductivity of the composite. On the other hand, D particles embedded in NARloy-Z matrix showed much improved thermal conductivity. Elemental analysis showed migration of Zr to the NARloy-Z-D interface to form ZrC, which appeared to provide a low contact thermal resistance. These results are consistent with the quantum mechanics-based model predictions. NARloy-Z-D composites have relatively high thermal conductivities and are promising for further development.

Bhat, Biliyar N.↗
NASA NTRS2013

Copper Multiwall Carbon Nanotubes and Copper-Diamond Composites for Advanced Rocket Engines

This paper reports on the research effort to improve the thermal conductivity of the copper-based alloy NARloy-Z (Cu-3 wt.%Ag-0.5 wt.% Zr), the state-of-the-art alloy used to make combustion chamber liners in regeneratively-cooled liquid rocket engines, using nanotechnology. The approach was to embed high thermal conductivity multiwall carbon nanotubes (MWCNTs) and diamond (D) particles in the NARloy-Z matrix using powder metallurgy techniques. The thermal conductivity of MWCNTs and D have been reported to be 5 to 10 times that of NARloy-Z. Hence, 10 to 20 vol. % MWCNT finely dispersed in NARloy-Z matrix could nearly double the thermal conductivity, provided there is a good thermal bond between MWCNTs and copper matrix. Quantum mechanics-based modeling showed that zirconium (Zr) in NARloy-Z should form ZrC at the MWCNT-Cu interface and provide a good thermal bond. In this study, NARloy-Z powder was blended with MWCNTs in a ball mill, and the resulting mixture was consolidated under high pressure and temperature using Field Assisted Sintering Technology (FAST). Microstructural analysis showed that the MWCNTs, which were provided as tangles of MWCNTs by the manufacturer, did not detangle well during blending and formed clumps at the prior particle boundaries. The composites made form these powders showed lower thermal conductivity than the base NARloy-Z. To eliminate the observed physical agglomeration, tangled multiwall MWCNTs were separated by acid treatment and electroless plated with a thin layer of chromium to keep them separated during further processing. Separately, the thermal conductivities of MWCNTs used in this work were measured, and the results showed very low values, a major factor in the low thermal conductivity of the composite. On the other hand, D particles embedded in NARloy-Z matrix showed much improved thermal conductivity. Elemental analysis showed migration of Zr to the NARloy-Z-D interface to form ZrC, which appeared to provide a low contact thermal resistance. These results are consistent with the quantum mechanics-based model predictions. NARloy-Z-D composites have relatively high thermal conductivities and are promising for further development.

Bhat, Biliyar N.↗
NASA NTRS1991

Thermal and solutal convection with conduction effects inside a rectangular enclosure

We numerically investigate the effects of various boundary conditions on the flow field characteristics of the physical vapor transport process. We use a prescribed temperature profile as boundary condition on the enclosure walls, and we consider parametric variations applicable to ground-based and space microgravity conditions. For ground-based applications, density gradients in the fluid phase generate buoyancy-driven convection which in turn disrupts the uniformity of the mass flux at the interface depending on the orientation. Heat conduction in the crystal can affect the fluid flow near the interface of the crystal. When considering isothermal source and sink at the interfaces, we observe a diffusive mode and three modes (i.e., thermal, solutal, and thermo-solutal). The convective modes show opposing flow field trends between thermal and solutal convection; theoretically, these trends can be used to achieve a uniform mass flux near the crystal. However, under the physical conditions chosen, the mathematical condition necessary for uniform mass flux cannot be satisfied because of thermodynamic restrictions. When a longitudinal thermal gradient is prescribed on the boundary of the crystal, a non-uniform interface temperature results, which induces a symmetrical fluid flow near the interface for the vertical case. For space microgravity applications, we show that the flow field is dominated by the Stefan wind and a uniform mass flux results at the interface.

Mennetrier, Christophe↗
NASA NTRS1994

Low-Heat-Leak Electrical Leads For Cryogenic Systems

Electrical leads offering high electrical conductivity and low thermal conductivity developed for use in connecting electronic devices inside cryogenic systems to power supplies, signal-processing circuits, and other circuitry located in nearby warmer surroundings. Strip of superconductive leads on ceramic substrate, similar to ribbon cable, connects infrared detectors at temperature of liquid helium with warmer circuitry. Electrical leads bridging thermal gradient at boundary of cryogenic system designed both to minimize conduction of heat from surroundings through leads into system and to minimize resistive heating caused by electrical currents flowing in leads.

Wise, Stephanie A.↗
NASA NTRS2024

A Thermally-Conductive Martian Core and Implications for Its Dynamo Cessation

Mars experienced a dynamo process that generated a global magnetic field ~4.3 (or earlier) to 3.6 billion years ago (Ga). The cessation of this dynamo strongly affected Mars’ history and is expected to be linked to thermochemical evolution of Mars’ iron-rich liquid core, which is strongly influenced by its thermal conductivity. Here, we directly measured thermal conductivities of solid iron-sulfur alloys to pressures relevant to the Martian core and temperatures to 1023 Kelvin. Our results show that a Martian core with 16 weight % sulfur has a thermal conductivity of ~19 to 32 Watt meter −1 Kelvin −1 from its top to the center, much higher than previously inferred from electrical resistivity measurements. Our modeled thermal conductivity profile throughout the Martian deep-mantle and core indicates a ~4- to 6-fold discontinuity across the core-mantle boundary. The core’s efficient cooling resulting from the depth-dependent, high conductivity diminishes thermal convection and forms thermal stratification, substantially contributing to cessation of Martian dynamo.

Wen-Pin Hsieh↗
DOE OSTI2020

Method of information entropy for convergence assessment of molecular dynamics simulations

The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1989

Prediction of the thermal environment and thermal response of simple panels exposed to radiant heat

A method of predicting the radiant heat flux distribution produced by a bank of tubular quartz heaters was applied to a radiant system consisting of a single unreflected lamp irradiating a flat metallic incident surface. In this manner, the method was experimentally verified for various radiant system parameter settings and used as a source of input for a finite element thermal analysis. Two finite element thermal analyses were applied to a thermal system consisting of a thin metallic panel exposed to radiant surface heating. A two-dimensional steady-state finite element thermal analysis algorithm, based on Galerkin's Method of Weighted Residuals (GFE), was formulated specifically for this problem and was used in comparison to the thermal analyzers of the Engineering Analysis Language (EAL). Both analyses allow conduction, convection, and radiation boundary conditions. Differences in the respective finite element formulation are discussed in terms of their accuracy and resulting comparison discrepancies. The thermal analyses are shown to perform well for the comparisons presented here with some important precautions about the various boundary condition models. A description of the experiment, corresponding analytical modeling, and resulting comparisons are presented.

Turner, Travis L.↗
DOE OSTI2021

Simulation of the PBX 9501 Self-Ignited Heavily Confined Cookoff Test

The PBX 9501 self-ignited thermally damaged explosive (SITX) test was carried out by Holmes et al. A 6-inch diameter sphere made of the PBX 9501 was heated inside a steel confinement vessel until it underwent a thermal explosion. This test provided measured temperature data that was used to define a boundary condition to carry out post-experiment finite element (FE) thermal simulations and comparisons of the FE thermal response with the experimental time to thermal ignition and internal temperature measurements at the center of the PBX 9501 sphere. Finite element models of this experiment were constructed using the heat transfer code Aria. A cookoff model for PBX 9501 was used to simulate the thermal decomposition of PBX 9501 and to predict the time to thermal ignition. An inverse heat conduction scheme was designed and applied to estimate a 1D temperature boundary condition.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1989

Variable thermal properties and thermal relaxation time in hyperbolic heat conduction

Numerical solutions were obtained for a finite slab with an applied surface heat flux at one boundary using both the hyperbolic (MacCormack's method) and parabolic (Crank-Nicolson method) heat conduction equations. The effects on the temperature distributions of varying density, specific heat, and thermal relaxation time were calculated. Each of these properties had an effect on the thermal front velocity (in the hyperbolic solution) as well as the temperatures in the medium. In the hyperbolic solutions, as the density or specific heat decreased with temperature, both the temperatures within the medium and the thermal front velocity increased. The value taken for the thermal relaxation time was found to determine the 'hyperbolicity' of the heat conduction model. The use of a time dependent relaxation time allowed for solutions where the thermal energy propagated as a high temperature wave initially, but approached a diffusion process more rapidly than was possible with a constant large relaxation time.

Glass, David E.↗
NASA NTRS1978

Comparison of reusable insulation systems for cryogenically-tanked earth-based space vehicles

Three reusable insulation systems concepts were developed for use with cryogenic tanks of earth-based space vehicles. Two concepts utilized double-goldized Kapton (DGK) or double-aluminized Mylar (DAM) multilayer insulation (MLI), while the third utilized a hollow-glass-microsphere, loadbearing insulation (LBI). Thermal performance measurements were made under space-hold (vacuum) conditions for insulating warm boundary temperatures of approximately 291 K. The resulting effective thermal conductivity was approximately 0.00008 W/m-K (W = weight,Kg; m = measured; K = temperature) for the MLI systems (liquid hydrogen test results) and 0.00054 W/m-K for the LBI system (liquid nitrogen test results corrected to liquid hydrogen temperature).

Sumner, I. E.↗
NASA NTRS1978

Comparison of reusable insulation systems for cryogenically-tanked earth-based space vehicles

Three reusable insulation systems concepts have been developed for use with cryogenic tanks of earth-based space vehicles. Two concepts utilized double-goldized Kapton (DGK) or double-aluminized Mylar (DAM) multilayer insulation (MLI), while the third utilized a hollow-glass-microsphere, load-bearing insulation (LBI). All three insulation systems have recently undergone experimental testing and evaluation under NASA-sponsored programs. Thermal performance measurements were made under space-hold (vacuum) conditions for insulation warm boundary temperatures of approximately 291 K. The resulting effective thermal conductivity was approximately .00008 W/m-K for the MLI systems (liquid hydrogen test results) and .00054 W/m-K for the LBI system (liquid nitrogen test results corrected to liquid hydrogen temperature). The DGK MLI system experienced a maximum thermal degradation of 38 percent, the DAM MLI system 14 percent, and the LBI system 6.7 percent due to repeated thermal cycling representing typical space flight conditions. Repeated exposure of the DAM MLI system to a high humidity environment for periods as long as 8 weeks provided a maximum degradation of only 24 percent.

Sumner, I. E.↗
NASA NTRS2004

Progress on Acoustic Measurements of the Bulk Viscosity of Near-Critical Xenon (BVX)

We plan to determine the bulk viscosity of xenon 10 times closer [in reduced temperature tau = (T-Tc)/Tc] to its liquid-vapor critical point than ever before. (Tc is the critical temperature.) To do so, we must measure the dispersion and attenuation of sound at frequencies 1/100 of those used previously. In general, sound attenuation has contributions from the bulk viscosity acting throughout the volume of the xenon as well as contributions from the thermal conductivity and the shear viscosity acting within thin thermoacoustic boundary layers at the interface between the xenon and the solid walls of the resonator. Thus, we can determine the bulk viscosity only when the boundary layer attenuation is small and well understood. We present a comparison of calculations and measurements of sound attenuation in the acoustic boundary layer of xenon near its liquid-vapor critical point.

Gillis, Keith A.↗