Abstract The Fourier and the hyperbolic heat conduction equations were solved numerically to simulate a frequency-domain thermoreflectance (FDTR) experiment. Numerical solutions enable isolation of pump and probe laser spot size effects and use of realistic boundary conditions. The equations were solved in time domain and the phase lag between the temperature of the transducer (averaged over the probe laser spot) and the modulated pump laser signal was computed for a modulation frequency range of 200 kHz–200 MHz. Numerical calculations showed that extracted values of the thermal conductivity are sensitive to both the pump and probe laser spot sizes, while analytical solutions (based on Hankel transform) cannot isolate the two effects. However, for the same effective (combined) spot size, the two solutions are found to be in excellent agreement. If the substrate (computational domain) is sufficiently large, the far-field boundary conditions were found to have no effect on the computed phase lag. The interface conductance between the transducer and the substrate was found to have some effect on the extracted thermal conductivity. The hyperbolic heat conduction equation yielded almost the same results as the Fourier heat conduction equation for the particular case studied. The numerically extracted thermal conductivity value (best fit) for the silicon substrate considered in this study was found to be about 82–108 W/m/K, depending on the pump and probe laser spot sizes used.
more »
« less
Frequency-domain probe beam deflection method for measurement of thermal conductivity of materials on micron length scale
Time-domain thermoreflectance and frequency-domain thermoreflectance (FDTR) have been widely used for non-contact measurement of anisotropic thermal conductivity of materials with high spatial resolution. However, the requirement of a high thermoreflectance coefficient restricts the choice of metal coating and laser wavelength. The accuracy of the measurement is often limited by the high sensitivity to the radii of the laser beams. We describe an alternative frequency-domain pump-probe technique based on probe beam deflection. The beam deflection is primarily caused by thermoelastic deformation of the sample surface, with a magnitude determined by the thermal expansion coefficient of the bulk material to measure. We derive an analytical solution to the coupled elasticity and heat diffusion equations for periodic heating of a multilayer sample with anisotropic elastic constants, thermal conductivity, and thermal expansion coefficients. In most cases, a simplified model can reliably describe the frequency dependence of the beam deflection signal without knowledge of the elastic constants and thermal expansion coefficients of the material. The magnitude of the probe beam deflection signal is larger than the maximum magnitude achievable by thermoreflectance detection of surface temperatures if the thermal expansion coefficient is greater than 5 × 10 −6 K −1 . The uncertainty propagated from laser beam radii is smaller than that in FDTR when using a large beam offset. We find a nearly perfect matching of the measured signal and model prediction, and measure thermal conductivities within 6% of accepted values for materials spanning the range of polymers to gold, 0.1–300 W/(m K).
more »
« less
- Award ID(s):
- 1720633
- PAR ID:
- 10412205
- Date Published:
- Journal Name:
- Review of Scientific Instruments
- Volume:
- 94
- Issue:
- 1
- ISSN:
- 0034-6748
- Page Range / eLocation ID:
- 014903
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The anisotropic Fourier Heat Conduction Equation (FHCE) and the multidimensional phonon Boltzmann Transport Equation (BTE) were solved numerically in cylindrical coordinates and in time domain to simulate a Time Domain Thermo-Reflectance (TDTR) experimental silicon/aluminum substrate/transducer setup. The out-of-phase response of the probe laser was predicted at various beam offset distances for a pump laser pulse frequency of 80 MHz and modulation frequency of 10 MHz and compared against experimental measurements for a silicon substrate. The isotropic FHCE was also solved for comparison. Results show that the isotropic FHCE with bulk thermal conductivity of 145 W/m/K significantly underpredicts the out-of-phase temperature difference, particularly at smaller beam offsets. With an isotropic thermal conductivity of 105 W/m/K, the computed results match experimental data at smaller beam offsets well, but overpredicts the experimental data at larger beam offsets. An almost-perfect match is obtained by using an anisotropic thermal conductivity wherein the radial (in-plane) thermal conductivity is set to 85 W/m/K and the axial (through-plane) conductivity is set to 130 W/m/K. The multidimensional frequency and polarization dependent phonon BTE is next solved. The BTE results for the out-of-phase temperature difference match experimental observations well at small and intermediate beam offsets, but overpredicts the experimental data at larger beam offsets. FHCE results are fitted to the BTE predictions, and the extracted (best fit) thermal conductivity is found to be 110 W/m/K.more » « less
-
Purpose AlSi10Mg alloy is commonly used in laser powder bed fusion due to its printability, relatively high thermal conductivity, low density and good mechanical properties. However, the thermal conductivity of as-built materials as a function of processing (energy density, laser power, laser scanning speed, support structure) and build orientation, are not well explored in the literature. This study aims to elucidate the relationship between processing, microstructure, and thermal conductivity. Design/methodology/approach The thermal conductivity of laser powder bed fusion (L-PBF) AlSi10Mg samples are investigated by the flash diffusivity and frequency domain thermoreflectance (FDTR) techniques. Thermal conductivities are linked to the microstructure of L-PBF AlSi10Mg, which changes with processing conditions. The through-plane exceeded the in-plane thermal conductivity for all energy densities. A co-located thermal conductivity map by frequency domain thermoreflectance (FDTR) and crystallographic grain orientation map by electron backscattered diffraction (EBSD) was used to investigate the effect of microstructure on thermal conductivity. Findings The highest through-plane thermal conductivity (136 ± 2 W/m-K) was achieved at 59 J/mm 3 and exceeded the values reported previously. The in-plane thermal conductivity peaked at 117 ± 2 W/m-K at 50 J/mm 3 . The trend of thermal conductivity reducing with energy density at similar porosity was primarily due to the reduced grain size producing more Al-Si interfaces that pose thermal resistance. At these interfaces, thermal energy must convert from electrons in the aluminum to phonons in the silicon. The co-located thermal conductivity and crystallographic grain orientation maps confirmed that larger colonies of columnar grains have higher thermal conductivity compared to smaller columnar grains. Practical implications The thermal properties of AlSi10Mg are crucial to heat transfer applications including additively manufactured heatsinks, cold plates, vapor chambers, heat pipes, enclosures and heat exchangers. Additionally, thermal-based nondestructive testing methods require these properties for applications such as defect detection and simulation of L-PBF processes. Industrial standards for L-PBF processes and components can use the data for thermal applications. Originality/value To the best of the authors’ knowledge, this paper is the first to make coupled thermal conductivity maps that were matched to microstructure for L-PBF AlSi10Mg aluminum alloy. This was achieved by a unique in-house thermal conductivity mapping setup and relating the data to local SEM EBSD maps. This provides the first conclusive proof that larger grain sizes can achieve higher thermal conductivity for this processing method and material system. This study also shows that control of the solidification can result in higher thermal conductivity. It was also the first to find that the build substrate (with or without support) has a large effect on thermal conductivity.more » « less
-
In this study, we report the length dependence of thermal conductivity ( k ) of zinc blende-structured Zinc Selenide (ZnSe) and Zinc Telluride (ZnTe) for length scales between 10 nm and 10 μm using first-principles computations, based on density-functional theory. The k value of ZnSe is computed to decrease significantly from 22.9 W m −1 K −1 to 1.8 W m −1 K −1 as the length scale is diminished from 10 μm to 10 nm. The k value of ZnTe is also observed to decrease from 12.6 W m −1 K −1 to 1.2 W m −1 K −1 for the same decrease in length. We also measured the k of bulk ZnSe and ZnTe using the Frequency Domain Thermoreflectance (FDTR) technique and observed a good agreement between the FDTR measurements and first principles calculations for bulk ZnSe and ZnTe. Understanding the thermal conductivity reduction at the nanometer length scale provides an avenue to incorporate nanostructured ZnSe and ZnTe for thermoelectric applications.more » « less
-
Gallium nitride (GaN) high electron mobility transistors (HEMTs) are key components of modern radio frequency (RF) power amplifiers. However, device self-heating negatively impacts both the performance and reliability of GaN HEMTs. Accordingly, laser-based pump-probe methods have been used to characterize the thermal resistance network of epitaxial material stacks that are used to fabricate HEMT structures. However, validation studies of these measurement results at the device level are lacking. In the present work, a GaN-on-SiC wafer was characterized using frequency-domain thermoreflectance and steady-state thermoreflectance techniques. The thermal conductivity of the GaN channel/buffer layer, SiC substrate, and the interfacial thermal boundary resistance at the GaN/SiC interface were determined. Results were validated by performing thermal imaging and modeling of a transmission line measurement (TLM) structure fabricated on the GaN-on-SiC wafer.more » « less