A fast-converging scheme for the phonon Boltzmann equation with dual relaxation times

Abstract

The phonon Boltzmann transport equation with dual relaxation times is often used to describe the heat conduction in semiconductor materials when the classical Fourier's law is no longer valid. For practical engineering designs, accurate and efficient numerical methods are highly demanded to solve the equation. At a large Knudsen number (i.e., the ratio of the phonon mean free path to a characteristic system length), steady-state solutions can be obtained via the conventional iterative scheme (CIS) within a few iterations. However, when the Knudsen number becomes small, i.e., when the phonon transport occurs in the diffusive or hydrodynamic regime, thousands of iterations are required to obtain converged results. In this work, a general synthetic iterative scheme (GSIS) is proposed to tackle the inefficiency of CIS. The key ingredient of the GSIS is that a set of macroscopic synthetic equations, which is exactly derived from the Boltzmann transport equation, is simultaneously solved with the kinetic equation to obtain the temperature and heat flux. During the iteration, the macroscopic quantities are used to evaluate the equilibria in the scattering terms of the kinetic equation, thus guiding the evolution of the phonon distribution function, while the distribution function, in turn, provides closures to the synthetic equations. The Fourier stability analysis is conducted to reveal the superiority of the GSIS over the CIS in terms of fast convergence in periodic systems. It is shown that the convergence rate of the GSIS can always be maintained under 0.2 so that only two iterations are required to reduce the iterative error by one order of magnitude. Numerical results in wall-bounded systems are presented to demonstrate further the efficiency of GSIS, where the CPU time is reduced by up to three orders of magnitude, especially in both the diffusive and hydrodynamic regimes where the Knudsen number is small.