A Fast Hybrid Fourier–Boltzmann Transport Equation Solver for Nongray Phonon Transport

Nongray phonon transport solvers based on the Boltzmann transport equation (BTE) are being increasingly employed to simulate submicron thermal transport in semiconductors and dielectrics. Typical sequential solution schemes encounter numerical difficulties because of the large spread in scattering rates. For frequency bands with very low Knudsen numbers, strong coupling between other BTE bands result in slow convergence of sequential solution procedures. This is due to the explicit treatment of the scattering kernel. In this paper, we present a hybrid BTE-Fourier model which addresses this issue. By establishing a phonon group cutoff Kn_c, phonon bands with low Knudsen numbers are solved using a modified Fourier equation which includes a scattering term as well as corrections to account for boundary temperature slip. Phonon bands with high Knudsen numbers are solved using the BTE. A low-memory iterative solution procedure employing a block-coupled solution of the modified Fourier equations and a sequential solution of BTEs is developed. The hybrid solver is shown to produce solutions well within 1% of an all-BTE solver (using Kn_c = 0.1), but with far less computational effort. Speedup factors between 2 and 200 are obtained for a range of steady-state heat transfer problems. The hybrid solver enables efficient and accurate simulation of thermal transport in semiconductors and dielectrics across the range of length scales from submicron to the macroscale.