The Boltzmann Transport equation (BTE) was solved numerically in cylindrical coordinates and in time domain to simulate a Frequency Domain Thermo-Reflectance (FDTR) experiment. First, a parallel phonon BTE solver that accounts for all phonon modes, frequencies, and polarizations was developed and tested. The solver employs the finite-volume method (FVM) for discretization of physical space, and the finite-angle method (FAM) for discretization of angular space. The solution was advanced in time using explicit time marching. The simulations were carried out in time domain and band-based parallelization of the BTE solver was implemented. The phase lag between the temperature averaged over the probed region of the transducer and the modulated laser pump signal was extracted for a pump laser modulation frequency ranging from 20–200 MHz. It was found that with the relaxation time scales used in the present study, the computed phase lag is underpredicted when compared to experimental data, especially at smaller modulation frequencies. The challenges in solving the BTE for such applications are highlighted.
This content will become publicly available on August 1, 2025
- Award ID(s):
- 2003747
- PAR ID:
- 10525069
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- International Journal of Heat and Mass Transfer
- Volume:
- 228
- Issue:
- C
- ISSN:
- 0017-9310
- Page Range / eLocation ID:
- 125698
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
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
-
more » « less
Abstract The Fourier heat conduction and the hyperbolic heat conduction equations were solved numerically to simulate a frequency-domain thermoreflectance (FDTR) experimental setup. Numerical solutions enable use of realistic boundary conditions, such as convective cooling from the various surfaces of the substrate and transducer. The equations were solved in time domain and the phase lag between the temperature at the center of the transducer and the modulated pump laser signal were computed for a modulation frequency range of 200 kHz to 200 MHz. It was found that the numerical predictions fit the experimentally measured phase lag better than analytical frequency-domain solutions of the Fourier heat equation based on Hankel transforms. The effects of boundary conditions were investigated and it was found that if the substrate (computational domain) is sufficiently large, the far-field boundary conditions have no effect on the computed phase lag. The interface conductance between the transducer and the substrate was also treated as a parameter, and was found to have some effect on the predicted thermal conductivity, but only in certain regimes. The hyperbolic heat conduction equation yielded identical results as the Fourier heat conduction equation for the particular case studied. The thermal conductivity value (best fit) for the silicon substrate considered in this study was found to be 108 W/m/K, which is slightly different from previously reported values for the same experimental data.
-
Zintl phase Mg 3 Sb 2 , which has ultra-low thermal conductivity, is a promising anisotropic thermoelectric material. It is worth noting that the prediction and experiment value of lattice thermal conductivity ( κ ) maintain a remarkable difference, troubling the development and application. Thus, we firstly included the four-phonon scattering processes effect and performed the Peierls–Boltzmann transport equation (PBTE) combined with the first-principles lattice dynamics to study the lattice thermal transport in Mg 3 Sb 2 . The results showed that our theoretically predicted κ is consistent with the experimentally measured, breaking through the limitations of the traditional calculation methods. The prominent four-phonon scatterings decreased phonon lifetime, leading to the κ of Mg 3 Sb 2 at 300 K from 2.45 (2.58) W m −1 K −1 to 1.94 (2.19) W m −1 K −1 along the in (cross)-plane directions, respectively, and calculation accuracy increased by 20%. This study successfully explains the lattice thermal transport behind mechanism in Mg 3 Sb 2 and implies guidance to advance the prediction accuracy of thermoelectric materials.more » « less
-
more » « less
Nanocrystalline silicon can have unique thermal transport and mechanical properties governed by its constituent grain microstructure. Here, we use phonon ray-tracing and molecular dynamics simulations to demonstrate the largely tunable thermomechanical behaviors with varying grain sizes (a0) and aspect ratios (ξ). Our work shows that, by selectively increasing the grain size along the heat transfer direction while keeping the grain area constant, the in-plane lattice thermal conductivity (kx) increases more significantly than the cross-plane lattice thermal conductivity (ky) due to anisotropic phonon–grain boundary scattering. While kx generally increases with increasing ξ, a critical value exists for ξ at which kx reaches its maximum. Beyond this transition point, further increases in ξ result in a decrease in kx due to substantial scattering of low-frequency phonons with anisotropic grain boundaries. Moreover, we observe reductions in the elastic and shear modulus with decreasing grain size, and this lattice softening leads to significant reductions in phonon group velocity and thermal conductivity. By considering both thermal and mechanical size effects, we identify two distinct regimes of thermal transport, in which anisotropic phonon–grain boundary scattering becomes more appreciable at low temperatures and lattice softening becomes more pronounced at high temperatures. Through phonon spectral analysis, we attribute the significant thermal conductivity anisotropy in nanograined silicon to grain boundary scattering of low-frequency phonons and the softening-driven thermal conductivity reduction to Umklapp scattering of high-frequency phonons. These findings offer insights into the manipulation of thermomechanical properties of nanocrystalline silicon via microstructure engineering, carrying profound implications for the development of future nanomaterials.
