Second-order diffusion limit for the phonon transport equation-asymptotics and numerics

We investigate the numerical implementation of the limiting equation for the phonon transport equation in the small Knudsen number regime. The main contribution is that we derive the limiting equation that achieves the second order convergence, and provide a numerical recipe for computing the Robin coefficients. These coefficients are obtained by solving an auxiliary half-space equation. Numerically the half-space equation is solved by a spectral method that relies on the even-odd decomposition to eliminate corner-point singularity. Numerical evidences will be presented to justify the second order asymptotic convergence rate.