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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effects of Galactic Irradiation on Thermal and Electronic Transport in Tungsten

The impact of irradiation on the thermal and electronic properties of materials is a persistent puzzle, particularly defect formation at the atomic and nanoscales. This work examines the nanoscale effects of low-energy irradiation on tungsten (W), focusing on defect-induced modifications to thermal and electronic transport. Using the site-projected thermal conductivity (SPTC) method [A. Gautam et al. PSS-RRL, 2400306, 2024], bulk and twin-grain boundary W with vacancy defects based on the Norgett–Robinson–Torrens displacements per atom model is analyzed. SPTC provides a detailed prediction of postcascade spatial thermal conductivity distribution. Electronic conductivity activity is estimated using the “N 2 method” [K. Nepal et al. Carbon, 119711, 2025] to explore the consequences of vacancies and grain boundaries, highlighting the defect-dependent nature of charge transport behavior. These findings offer high-resolution insights into irradiation-driven transport phenomena, with implications for space-exposed materials and nanoscale thermal/electronic management.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Vertical forces during VDEs in an ITER plasma and the role of halo currents

Vertical displacement events (VDEs) can occur in elongated tokamaks causing large currents to flow in the vessel and other adjacent metallic structures. To better understand the potential magnitude of the associated forces and the role of the so called halo currents'' on them, we have used the M3D-C1 code to simulate potential VDEs in ITER. We used actual values for the vessel resistivity and pre-quench temperatures and, unlike most of the previous studies, the halo region is naturally formed by triggering the thermal quench with an increase in the plasma thermal conductivity. We used the 2D non-linear version of the code and vary the post-thermal quench thermal conductivity profile as well as the boundary temperature in order to generate a wide range of possible cases that could occur in the experiment. We also show that, for a similar condition, increasing the halo current does not increase the total force on the wall since it is offset by a decrease in the toroidal contribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

The effect of coatings and liners on heat transfer in a dry shaft-bush tribosystem

The temperatures due to frictional heating within a solid lubricated or coated journal bearing were analyzed by using a finite element method. A solid model of the shaft-bush tribocontact was generated with an eight-node, three-dimensional, first-order isoparametric heat-transfer element and the Patran solid modeler software. The Patmar (Patran-Marc) translator was used to help develop the Marc-based finite element program for the system; this software was used on the Cray X-MP supercomputer to perform a finite element analysis of the contact. The analysis was performed for various liner materials, for thin, hard, wear-resistant coated bearings, and for different geometries and thermal cooling boundary conditions. The analyses indicated that thermal conductivity of the liner or coating material is the most vital thermal parameter that controls the interface temperature. In addition to design variations, the proximity of the cooling source to the heat-flux-generating interface is critically important to the temperature control in the system.

Ghosh, Mihir K.↗
NASA NTRS1973

Atmospheric Aerosol and Thermal Structure in the Boundary Layer Over the Los Angeles Basin

A field study using a mobile lidar was recently conducted in the L. A. Basin, California, to (1) examine the relationship between the vertical aerosol and the thermal structure, and (2) map the vertical aerosol structure in the atmospheric boundary layer over the basin. These data are needed for use in the development of a mixing-depth submodel required for photochemical air Quality simulation models. Toward these ends, a series of lidar aerosol measurements in conjunction with balloon and aircraft temperature soundings were taken at a site in El Monte, and in a mobile mode along a 90-mile freeway loop between El Monte, Santa Monica, and Long Beach. The lidar data are presented in the form of time-height and distance-height cross sections. The results indicate that, although aerosol concentrations are frequently present above the base of the marine inversion, these are generally in stratified layers in contrast to the more uniform nature of the lower convective layer, permitting the mixing depth to be distinguished on this basis. The lidar-derived mixing depths are well correlated (within 100 m) with daytime temperature inversions. Other significant features shown by the lidar data include large Basin-wide mixing-depth variations, waves with amplitudes of 200-300 m and wavelengths of 1000-1500 m on the lower aerosol layer, and apparent aerosol "chimneys" with overrunning in the vicinity of convergence zones.

Johnson, Warren B.↗
NASA NTRS1991

Prediction of thermal acoustic oscillations (TAOs) in the CLAES solid CO2/neon system

Results are presented of a study initiated to investigate the possibility that the existence of thermal acoustic oscillations (TAOs) in the Cryogenic Limb Atmospheric Etalon Spectrometer (CLAES) neon plumbing system ground configuration could be the cause of higher-than-predicted heat rates measured during thermal ground testing. Tests were conducted between warm boundary temperatures ranging from 40 to 100 K, which simulated the actual test conditions of the CLAES CO2/neon system. TAOs were observed between 6 and 106 Torr, which agreed with the analytical predictions, and verified the possible existence of TAOs in the CLAES system during ground testing. The presence of TAOs was eventually confirmed in the CLAES system during a subsequent thermal test and were determined to have caused the higher heat rates measured during the prior thermal test.

Spradley, I. E.↗
NASA NTRS1974

Basic performance of a multilayer insulation system containing 20 to 160 layers

An experimental investigation was conducted to determine the thermal effectiveness of an aluminized Mylar-silk net insulation system containing up to 160 layers. The experimentally measured heat flux was compared with results predicted by using (1) a previously developed semi-empirical equation and (2) an effective-thermal-conductivity value. All tests were conducted at a nominal hot-boundary temperature of 294 K (530 R) with liquid hydrogen as the heat sink. The experimental results show that the insulation performed as expected and that both the semi-empirical equation and effective thermal conductivity of a small number of layers were adequate in predicting the thermal performance of a large number of layers of insulation.

Stochl, R. J.↗
NASA NTRS1989

Effects of anisotropic heat conduction on solidification

Two-dimensional solidification influenced by anisotropic heat conduction has been considered. The interfacial energy balance was derived to account for the heat transfer in one direction (x or y) depending on the temperature gradient in both the x and y directions. A parametric study was made to determine the effects of the Stefan number, aspect ratio, initial superheat, and thermal conductivity ratios on the solidification rate. Because of the imposed boundary conditions, the interface became skewed and sometimes was not a straight line between the interface position at the upper and lower adiabatic walls (spatially nonlinear along the height). This skewness depends on the thermal conductivity ratio k(yy)/k(yx). The nonlinearity of the interface is influenced by the solidification rate, aspect ratio, and k(yy/k(yx).

Weaver, J. A.↗
NASA NTRS1991

A numerical simulation of magnetic reconnection and radiative cooling in line-tied current sheets

Radiative MHD equations are used for an optically thin plasma to carry out a numerical experiment related to the formation of 'postflare' loops. The numerical experiment starts with a current sheet that is in mechanical and thermal equilibrium but is unstable to both tearing-mode and thermal-condensation instabilities. The current sheet is line-tied at one end to a photospheric-like boundary and evolves asymmetrically. The effects of thermal conduction, resistivity variation, and gravity are ignored. In general, reconnection in the nonlinear stage of the tearing-mode instability can strongly affect the onset of condensations unless the radiative-cooling time scale is much smaller than the tearing-mode time scale. When the ambient plasma is less than 0.2, the reconnection enters a regime where the outflow from the reconnection region is supermagnetosonic with respect to the fast-mode wave speed. In the supermagnetosonic regime the most rapidly condensing regions occur downstream of a fast-mode shock that forms where the outflow impinges on closed loops attached to the photospheric-like boundary. A similar shock-induced condensation might occur during the formation of 'postflare' loops.

Forbes, T. G.↗
DOE OSTI2024

Heat transfer optimization of uo 2 -mo fuel using genetic algorithms

Two genetic algorithm (GA) methods were applied to thermal finite element models to optimize the heat transfer efficacy of a UO 2 -Mo composite fuel pellet with typical pressurized water reactor fuel geometry. Mo additions to UO 2 have been shown to increase the thermal conductivity, thus reducing centerline temperatures and temperature gradients. Previous studies evaluated uniformly dispersed Mo or continuous Mo internal geometries (e.g., fins, plates, discs) that were selected using engineering intuition. The current study uses two different implementations of the same GA to optimize Mo placement and minimize the fuel temperature with the only constraint being a maximum 10% Mo volume fraction. One approach superimposed Mo line elements onto the monolithic UO 2 pellet model, and the other converted entire UO 2 volume elements to Mo. The former method generated 1D heat transfer connections between nodes, whereas the latter method allowed for the formation of 3D structures. Features of the optimal fuel design produced by the GAs included dispersed Mo near the centerline that shifted the peak fuel temperature outward by 0.6 mm, Mo chains in the high-heat-flux region in the mid-to-outer radial zone, and a large continuous structure that spanned the full radius and height of the pellet and accounted for 87.7 % of the total Mo in the pellet. Analysis of this design indicates that the optimal Mo configuration is a balance between creating continuous heat transfer pathways and optimally dispersing Mo to minimize the heat transfer distance through UO 2 . This architecture ultimately produced an effective thermal conductivity of 11.3 W/m·K under the assumed boundary conditions. This result is higher than any previous values from the literature. In conclusion, potential fabrication methods and challenges are discussed in addition to the implications on fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Thermal and Electrical Properties of Isocyanate Derived Organic Aerogels for Cryogenic Insulation Applications

The use and storage of cryogens such as liquefied nitrogen, helium, hydrogen among others requires reliable and efficient thermal insulation systems. Passive insulation from high performance materials that are well-known for their inherent low thermal conductivity would reduce the overall costs involved in design, manufacture and maintenance of such systems. One such class of materials are referred to as aerogels. These materials are known for their low density, high mesoporosity, high surface areas, low thermal conductivity and high acoustic impedance. Aerogels were invented by S.S. Kistler in 1931 and the most common type are those made of silica. However, the inherent fragility of silica aerogels makes them hard to mass produce, and therefore applications have been limited. A major breakthrough was introduced by our team almost 20 years ago with the invention of polymer crosslinked silica aerogels. Those materials shifted attention to all-polymer aerogels that have overcome all fragility issues associated with their inorganic counterparts. This study focuses on such polymeric aerogels that can be mass produced as large monoliths while maintaining the low thermal conductivity of traditional silica aerogels over a wide temperature range. Manufacturing flexibility of polymeric aerogels allows fabrication of blocks and sheets that can be applied in various configurations to insulate cryogenic and superconducting devices. The thermal conductivity with 80 K and room temperature boundary are reported as well as other properties (electrical), that need to be considered when designing devices for cryogenic applications.

Churu, Gitogo↗
DOE OSTI2020

Unification of nonequilibrium molecular dynamics and the mode-resolved phonon Boltzmann equation for thermal transport simulations

Nano-size confinement induces many intriguing non-Fourier heat conduction phenomena, such as nonlinear temperature gradients, temperature jumps near the contacts, and size-dependent thermal conductivity. Over the past decades, these effects have been studied and interpreted by nonequilibrium molecular dynamics (NEMD) and phonon Boltzmann transport equation (BTE) simulations separately, but no theory that unifies these two methods has ever been established. In this work, we unify these methods using a quantitative mode-level comparison and demonstrate that they are equivalent for various thermostats. We show that different thermostats result in different non-Fourier thermal transport characteristics due to the different mode-level phonon excitations inside the thermostats, which explains the different size-dependent thermal conductivities calculated using different reservoirs, even though they give the same bulk thermal conductivity. Specifically, the Langevin thermostat behaves like a thermalizing boundary in phonon BTE and provides mode-level thermal-equilibrium phonon outlets, while the Nose-Hoover chain thermostat and velocity rescaling method behave like biased reservoirs, which provide a spatially uniform heat generation and mode-level nonequilibrium phonon outlets. These findings explain why different experimental measurement methods can yield different size-dependent thermal conductivity. They also indicate that the thermal conductivity of materials can be tuned for various applications by specifically designing thermostats. Finally, the unification of NEMD and phonon BTE will largely facilitate the study of thermal transport in complex systems in the future by, e.g., replacing computationally unaffordable first-principles NEMD simulations with computationally less expensive spectral BTE simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2011

Assessment of CFD Hypersonic Turbulent Heating Rates for Space Shuttle Orbiter

Turbulent CFD codes are assessed for the prediction of convective heat transfer rates at turbulent, hypersonic conditions. Algebraic turbulence models are used within the DPLR and LAURA CFD codes. The benchmark heat transfer rates are derived from thermocouple measurements of the Space Shuttle orbiter Discovery windward tiles during the STS-119 and STS-128 entries. The thermocouples were located underneath the reaction-cured glass coating on the thermal protection tiles. Boundary layer transition flight experiments conducted during both of those entries promoted turbulent flow at unusually high Mach numbers, with the present analysis considering Mach 10{15. Similar prior comparisons of CFD predictions directly to the flight temperature measurements were unsatisfactory, showing diverging trends between prediction and measurement for Mach numbers greater than 11. In the prior work, surface temperatures and convective heat transfer rates had been assumed to be in radiative equilibrium. The present work employs a one-dimensional time-accurate conduction analysis to relate measured temperatures to surface heat transfer rates, removing heat soak lag from the flight data, in order to better assess the predictive accuracy of the numerical models. The turbulent CFD shows good agreement for turbulent fuselage flow up to Mach 13. But on the wing in the wake of the boundary layer trip, the inclusion of tile conduction effects does not explain the prior observed discrepancy in trends between simulation and experiment; the flight heat transfer measurements are roughly constant over Mach 11-15, versus an increasing trend with Mach number from the CFD.

Wood, William A.↗
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.↗