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 145 records · Page 8

Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) Model - An Unified Concept for Earthquake Precursors Validation

The paper presents a conception of complex multidisciplinary approach to the problem of clarification the nature of short-term earthquake precursors observed in atmosphere, atmospheric electricity and in ionosphere and magnetosphere. Our approach is based on the most fundamental principles of tectonics giving understanding that earthquake is an ultimate result of relative movement of tectonic plates and blocks of different sizes. Different kind of gases: methane, helium, hydrogen, and carbon dioxide leaking from the crust can serve as carrier gases for radon including underwater seismically active faults. Radon action on atmospheric gases is similar to the cosmic rays effects in upper layers of atmosphere: it is the air ionization and formation by ions the nucleus of water condensation. Condensation of water vapor is accompanied by the latent heat exhalation is the main cause for observing atmospheric thermal anomalies. Formation of large ion clusters changes the conductivity of boundary layer of atmosphere and parameters of the global electric circuit over the active tectonic faults. Variations of atmospheric electricity are the main source of ionospheric anomalies over seismically active areas. Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model can explain most of these events as a synergy between different ground surface, atmosphere and ionosphere processes and anomalous variations which are usually named as short-term earthquake precursors. A newly developed approach of Interdisciplinary Space-Terrestrial Framework (ISTF) can provide also a verification of these precursory processes in seismically active regions. The main outcome of this paper is the unified concept for systematic validation of different types of earthquake precursors united by physical basis in one common theory.

Pulinets, S.↗

Aerodynamic heating in gaps of thermal protection system tile arrays in laminar and turbulent boundary layers

An experimental heat-transfer investigation was conducted on two staggered arrays of metallic tiles in laminar and turbulent boundary layers. This investigation was conducted for two purposes. The impingement heating distribution where flow in a longitudinal gap intersects a transverse gap and impinges on a downstream blocking tile was defined. The influence of tile and gap geometries was analyzed to develop empirical relationships for impingement heating in laminar and turbulent boundary layers. Tests were conducted in a high temperature structures tunnel at a nominal Mach number of 7, a nominal total temperature of 1800 K, and free-stream unit Reynolds numbers from 1.0 x 10 million to 4.8 x 10 million per meter. The test results were used to assess the impingement heating effects produced by parameters that include gap width, longitudinal gap length, slope of the tile forward-facing wall, boundary-layer displacement thickness, Reynolds number, and local surface pressure.

Avery, D. E.↗

A Modified Wall Matching Treatment to Account for Local Solid to Fluid Thermal Coupling

The wall-matching methodology of Wilcox is modified to include a solid-wall, thermal-conduction model. This coupled fluid-thermal-structure model is derived assuming that the wall thermal-structure behavior is locally one-dimensional and that structural deformations, due to thermally induced stresses, are not significant. The one-dimensional coupled fluid-thermal-structure model is derived such that the wall temperature is removed as an independent boundary condition variable. The one-dimensional coupled fluid-thermal-structure model is also derived for the general case of an arbitrary mixture of thermally prefect gases and a wall of arbitrary thickness and conductivity by using a compressible, streamwise-pressure-gradient-corrected, wall-matching function and Fourier's law of heat conduction. The resulting model was implemented in the VULCAN CFD code as a new boundary condition type. VULCAN was then used to simulate a two-dimensional Mach 6 wind tunnel facility nozzle flow to demonstrate/validate the one-dimensional coupled fluid-thermal-structure model. The nozzle internal-wall surface temperature and heat transfer distributions computed using the one-dimensional coupled fluid-thermal-structure model are compared to wall temperature and heat transfer distributions from an iterative multi-dimensional analysis obtained by coupling the VULCAN CFD code and the MSC/NASTRAN-thermal code. The one-dimensional coupled fluid-thermal-structure model analysis is shown to be very robust and in excellent agreement with the multi-dimensional iteratively coupled analysis. It is also shown that the one-dimensional analysis can be used as an initial guess for the multi-dimensional iteratively coupled analysis.

White, Jeffery A.↗

A Modified Wall Matching Treatment to Account for Local Solid to Fluid Thermal Coupling

The wall-matching methodology of Wilcox is modified to include a solid-wall, thermal-conduction model. This coupled fluid-thermal-structure model is derived assuming that the wall thermal-structure behavior is locally one-dimensional and that structural deformations, due to thermally induced stresses, are not significant. The one-dimensional coupled fluid-thermal-structure model is derived such that the wall temperature is removed as an independent boundary condition variable. The one-dimensional coupled fluid-thermal-structure model is also derived for the general case of an arbitrary mixture of thermally prefect gases and a wall of arbitrary thickness and conductivity by using a compressible, streamwise-pressure-gradient-corrected, wall-matching function and Fourier's law of heat conduction. The resulting model was implemented in the VULCAN CFD code as a new boundary condition type. VULCAN was then used to simulate a two-dimensional Mach 6 wind tunnel facility nozzle flow to demonstrate/validate the one-dimensional coupled fluid-thermal-structure model. The nozzle internal-wall surface temperature and heat transfer distributions computed using the one-dimensional coupled fluid-thermal-structure model are compared to wall temperature and heat transfer distributions from an iterative multi-dimensional analysis obtained by coupling the VULCAN CFD code and the MSC/NASTRAN-thermal code. The one-dimensional coupled fluid-thermal-structure model analysis is shown to be very robust and in excellent agreement with the multi-dimensional iteratively coupled analysis. It is also shown that the one-dimensional analysis can be used as an initial guess for the multi-dimensional iteratively coupled analysis.

White, Jeffery A.↗

Development and Application of Interatomic Potentials for Ultra High Temperature Ceramics (UHTC): ZrB2 and HfB2

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are characterized by high melting point, good strength, and reasonable oxidation resistance. These materials are of interest for use as sharp leading edges for hypersonic vehicles among other applications. Progress in computational modeling of UHTCs has been limited in part due to the absence of suitable interatomic potentials. We present a Tersoff style parameterization of such potentials for ZrB2 and HfB2 appropriate for atomistic simulations. Parameters are fit to data generated from ab initio calculations. The accuracy of the potentials is assessed against further ab initio data. As a first non ]trivial application, molecular dynamics simulations are performed to evaluate the thermal conductivity of single crystals and the thermal resistance of high symmetry grain boundaries.

Lawson, John W.↗

Hierarchical microstructure of Yb-filled skutterudites through ultrasonically atomized spherical powders for enhanced thermoelectric performance

Hierarchical microstructures are widely explored as a strategy to reduce lattice thermal conductivity in thermoelectric materials while preserving favorable electronic transport. Here we demonstrate that ultrasonic atomization of precursor powders provides an intrinsic pathway to generate such architectures in n-type Yb0.3Co4Sb12 skutterudites. Materials were synthesized by both conventional sealed-ampoule and ultrasonic atomization routes, each yielding essentially single-phase skutterudite after densification. However, the atomized spherical powders undergo a distinct reaction-driven microstructural evolution during thermal treatment prior to densification. Electron microscopy reveals that the skutterudite phase within the particles induces swelling, cracking, and fragmentation of the droplets before consolidation. Multimodal structural and microstructural characterizations show that densification of these fragmented powders produces a hierarchical microstructure consisting of submicron grains, micrometer-scale faceted crystallites, and powder-derived domains separated by oxide-decorated boundaries. Comprehensive thermoelectric transport measurements evidence that this multiscale architecture lowers the lattice thermal conductivity while maintaining favorable electronic transport, leading to an average ~10% enhancement in the thermoelectric figure of merit. Complementary 121Sb Mössbauer spectroscopy further indicates lattice softening in the atomized material. These results demonstrate that ultrasonic atomization provides an effective route to engineer hierarchical microstructures in skutterudite thermoelectrics and offers a promising strategy for tailoring thermal transport in energy-conversion materials.

Bouteiller, Hugo [ORNL] (ORCID:0009000421322962)↗

Effect of variable thermal conductivity on isotherms in Bridgman growth

A change in thermal conductivity associated with melting or solidification can have a profound influence on the isotherms near the solidification interface if the material is being directionally solidified in an ampoule whose walls carry a substantial portion of the heat. This analysis was prompted by a recent discovery that the thermal conductivity of Hg(1-x)CD(x)Te increased dramatically as the material is heated above the solidus curve. An illustrative example is shown in which the sample is approximated as an infinite cylinder with constant but diffferent thermal properties in the solid and melt. The boundary conditions are fixed on the surface by a conductive ampoule in a two-zone Bridgman furnace with an adiabatic region separating the two zones. The effect of the adiabatic zone in this case is to intensify the curvature of the interface rather than to lessen it.

Naumann, R. J.↗

Thermoacoustic convection of fluids in low gravity

The heat flow in a confined perfect gas in low gravity is investigated, including the effects of conduction and thermal convection. Buoyancy-driven flow is neglected, due to the low-gravity environment, but the effect of thermoacoustic motion due to fluid compressibility is included. One-dimensional mathematical models are constructed from the conservation equations for a compressible, viscous, heat-conducting fluid. A conservative, time-dependent finite-difference method is used to generate numerical solutions on a digital computer. Problems for flat plates and cylindrical segments are solved for specified thermal boundary conditions. Numerical results are given which indicate that thermoacoustic convection can significantly increase the transient heat flow over conduction model predictions for cases where a confined gas is rapidly heated.

Spradley, L. W.↗

Measuring sub-surface spatially varying thermal conductivity of silicon implanted with krypton

The thermal properties of semiconductors following exposure to ion irradiation are of great interest for the cooling of electronic devices; however, gradients in composition and structure due to irradiation often make the measurement difficult. Furthermore, the nature of spatial variations in thermal resistances due to spatially varying ion irradiation damage is not well understood. In this work, we develop an advancement in the analysis of time-domain thermoreflectance to account for spatially varying thermal conductivity in a material resulting from a spatial distribution of defects. We then use this method to measure the near-surface (≲1 μm) thermal conductivity of silicon wafers irradiated with Kr + ions, which has an approximate Gaussian distribution centered 260 nm into the sample. Our numerical analysis presented here allows for the spatial gradient of thermal conductivity to be extracted via what is fundamentally a volumetric measurement technique. We validate our findings via transmission electron microscopy, which is able to confirm the spatial variation of the sub-surface silicon structure, and provide additional insight into the local structure resulting from the effects of ion bombardment. Thermal measurements found the ion stopping region to have a nearly 50x reduction in thermal conductivity as compared to pristine silicon, while TEM showed the region was not fully amorphized. Our results suggest this drastic reduction in silicon thermal conductivity is primarily driven by structural defects in crystalline regions along with boundary scattering between amorphous and crystalline regions, with a negligible contribution being due to implanted krypton ions themselves.

36 MATERIALS SCIENCE↗

Growth of zinc selenide single crystals by physical vapor transport in microgravity

The goals of this research were the optimization of growth parameters for large (20 mm diameter and length) zinc selenide single crystals with low structural defect density, and the development of a 3-D numerical model for the transport rates to be expected in physical vapor transport under a given set of thermal and geometrical boundary conditions, in order to provide guidance for an advantageous conduct of the growth experiments. In the crystal growth studies, it was decided to exclusively apply the Effusive Ampoule PVT technique (EAPVT) to the growth of ZnSe. In this technique, the accumulation of transport-limiting gaseous components at the growing crystal is suppressed by continuous effusion to vacuum of part of the vapor contents. This is achieved through calibrated leaks in one of the ground joints of the ampoule. Regarding the PVT transport rates, a 3-D spectral code was modified. After introduction of the proper boundary conditions and subroutines for the composition-dependent transport properties, the code reproduced the experimentally determined transport rates for the two cases with strongest convective flux contributions to within the experimental and numerical error.

Rosenberger, Franz↗

CFD Analysis of Hypersonic Flowfields With Surface Thermochemistry and Ablation

In the past forty years much progress has been made in computational methods applied to the solution of problems in spacecraft hypervelocity flow and heat transfer. Although the basic thermochemical and physical modeling techniques have changed little in this time, several orders of magnitude increase in the speed of numerically solving the Navier-Stokes and associated energy equations have been achieved. The extent to which this computational power can be applied to the design of spacecraft heat shields is dependent on the proper coupling of the external flow equations to the boundary conditions and governing equations representing the thermal protection system in-depth conduction, pyrolysis and surface ablation phenomena. A discussion of the techniques used to do this in past problems as well as the current state-of-art is provided. Specific examples, including past missions such as Galileo, together with the more recent case studies of ESA/Rosetta Sample Comet Return, Mars Pathfinder and X-33 will be discussed. Modeling assumptions, design approach and computational methods and results are presented.

Henline, W. D.↗

The Adaptively-Refined MHD Solver (ARMS)

ARMS is a robust code for solving the single-fluid, fully compressible equations of magnetohydrodynamics. In this presentation I will describe how to access and run ARMS and how to use HelioSpace to analyze the results. The main features of the code are the following: 1. Equations are solved in either Cartesian or Spherical coordinates. For the latter, either a linear or a Log(r) coordinate can be used. 2. The latest flux-corrected-transport schemes developed at NRZ, are used to advance the convective equations. All waves and shocks are accurately captured. 3. A fully adaptive grid can be used. The code incorporates the PARAMESH technology for block-adaptation. Arbitrary adaptation criteria can be specified. 4. The code maintains a divergence-free magnetic field to machine accuracy throughout all refinement and de-refinement. 5. A uniform resistivity and thermal conduction can be included. Solution is either explicit or implicit. 6. Arbitrary boundary and initial conditions can be used in either Cartesian or spherical coordinates. 7. An accompanying visualization package, HelioSpace, is available for viewing and analyzing the native ARMS output, (without need for interpolating onto a structured grid).

Antichos, S. K.↗

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗

Experimental confirmation of first-principles thermal conductivity in Zirconium-doped ThO 2

The degradation of thermal conductivity in advanced nuclear fuels due to the accumulation of fission products and irradiation-induced defects is inevitable, and must be considered as part of safety and efficiency analyses of nuclear reactors. Here, this study examines the thermal conductivity of a zirconium-doped ThO 2 crystal, synthesized via the hydrothermal method using a spatial domain thermoreflectance technique. Zirconium is one of the soluble fission products in oxide fuels that can effectively scatter heat-carrying phonons in the crystalline lattice of fuel. Thus, thermal property measurements of zirconium-doped ThO 2 single crystals provide insights into the effects of substitutional zirconium doping, isolated from extrinsic factors such as grain boundary scattering. The experimental results are compared with first-principles calculations of the lattice thermal conductivity of ThO 2 , employing an iterative solution of the Peierls-Boltzmann transport equation. Additionally, the non-perturbative Green's function methodology is utilized to compute phonon-point defect scattering rates, accounting for local distortions around point defects, including mass difference changes, interatomic force constants, and structural relaxation. The congruence between the predicted results from first-principles calculations and the measured temperature-dependent thermal conductivity validates the computational methodology. Furthermore, the methodologies employed in this study enable systematic investigations of thermal conductivity reduction by fission products, potentially leading to the development of more accurate fuel performance codes.

36 - MATERIALS SCIENCE↗

Understanding Thermal Transport in Polymer –Silver Nanowire Composites

Understanding thermal transport across filler-polymer interfaces and filler-filler contacts within polymer composites is of great importance for better thermal design of the composites that are widely used in high-performance heat exchangers, energy storage devices, and flexible electronics. Over the past several decades, considerable progress has been made in improving the thermal conductivity of polymer composites, but several key questions concerning the influence of interfacial thermal resistance, or Kapitza resistance, still remain. Firstly, the thermal properties of these composites are highly dependent on thermal transport through the filler network and its contacts. For metallic nanofillers, the thermal conductivity is often estimated using the Wiedemann-Franz law based on electrical conductivity; however, it remains a question whether the Wiedemann-Franz law still holds at nanoscale contacts. Through investigation of silver nanowires of varying sizes, we were able to demonstrate that the Lorenz number for silver nanowire increases with decreasing nanowire diameter. Examination of the corresponding electrical and thermal conductivities indicate that these changes are due to that the relative contribution of phonons becomes more significant as a result of elastic stiffening. Furthermore, we show that for silver nanowires, the contact thermal resistance is ~8 times lower than that of multi-walled carbon nanotube (MWCNT) of similar diameters. Additionally, through systematic studies of electrospun polymer-silver nanowire composite nanofibers, we investigated the impact of interface morphology on the thermal conductivity enhancement of the composite system and probed the value of the Kapitza resistance for individual polymer-filler interfaces. For polymer nanofibers containing continuous, single silver nanowires, the thermal conductivity increases linearly with increasing volume fraction of silver, which is consistent with the prediction of percolation theory for samples above the percolation limit. By comparing this linear trend to the measured thermal conductivities of composite nanofibers with more complex structures, we were able to determine the resistance associated with any additional polymer-filler boundaries. In doing so we find that the thermal boundary resistance for polyvinylpyrrolidone (PVP)–silver interfaces is significantly lower than that of comparable polymer-MWCNT composite systems.Together our studies provide new insights into thermal transport in polymer nanocomposites and should help facilitate the design of high performance polymeric thermal interface material.

Thermal transport↗

Understanding Thermal Transport in Polymer - Silver Nanowire Composites

Understanding thermal transport across filler-polymer interfaces and filler-filler contacts within polymer composites is of great importance for better thermal design of the composites that are widely used in high-performance heat exchangers, energy storage devices, and flexible electronics. Over the past several decades, considerable progress has been made in improving the thermal conductivity of polymer composites, but several key questions concerning the influence of interfacial thermal resistance, or Kapitza resistance, still remain. Firstly, the thermal properties of these composites are highly dependent on thermal transport through the filler network and its contacts. For metallic nanofillers, the thermal conductivity is often estimated using the Wiedemann-Franz law based on electrical conductivity; however, it remains a question whether the Wiedemann-Franz law still holds at nanoscale contacts. Through investigation of silver nanowires of varying sizes, we were able to demonstrate that the Lorenz number for silver nanowire increases with decreasing nanowire diameter. Examination of the corresponding electrical and thermal conductivities indicate that these changes are due to that the relative contribution of phonons becomes more significant as a result of elastic stiffening. Furthermore, we show that for silver nanowires, the contact thermal resistance is ~8 times lower than that of multi-walled carbon nanotube (MWCNT) of similar diameters. Additionally, through systematic studies of electrospun polymer-silver nanowire composite nanofibers, we investigated the impact of interface morphology on the thermal conductivity enhancement of the composite system and probed the value of the Kapitza resistance for individual polymer-filler interfaces. For polymer nanofibers containing continuous, single silver nanowires, the thermal conductivity increases linearly with increasing volume fraction of silver, which is consistent with the prediction of percolation theory for samples above the percolation limit. By comparing this linear trend to the measured thermal conductivities of composite nanofibers with more complex structures, we were able to determine the resistance associated with any additional polymer-filler boundaries. In doing so we find that the thermal boundary resistance for polyvinylpyrrolidone (PVP)–silver interfaces is significantly lower than that of comparable polymer-MWCNT composite systems. Together our studies provide new insights into thermal transport in polymer nanocomposites and should help facilitate the design of high performance polymeric thermal interface material.

Matthew L. Fitzgerald↗

Developing Controlled Conductive Boundaries for JWST Cryogenic Testing

In 2017, the James Webb Space Telescope (JWST) underwent functional testing and optical metrology verification of the combined Optical Telescope Element and Integrated Science Instrument Module (OTIS) under cryogenic vacuum conditions in Chamber A at the Johnson Space Center. Maintaining flight-like thermal boundary conditions was a critical requirement for optical testing and required unique and challenging Ground Support Equipment (GSE) design solutions. Two such GSE systems, the Integrated Science Instrument Module (ISIM) Precool Straps and the Hardpoint Struts were direct conduction interfaces to the flight hardware. Hardware safety during cooldown required detailed design of their conductivity, and thermal balance testing required "zero-Q" (0-Q) heater implementation to bring the heat flow to zero, thereby cutting off these non-flight conductive links after operating temperatures were achieved. This paper describes the design considerations and approach implemented to achieve the required flight hardware cool down and return to ambient conditions, ensure flight hardware safety, and minimize the non-flight-like heat flows to or from the observatory during cryo-stable testing.

Cooke, D.↗