An Asymptotic-Preserving Hybrid Angular Discretization for the Gray Radiative Transfer Equations

Abstract

The aim of this paper is to construct a new numerical scheme for the nonlinear gray radiative transfer (GRT) equations, namely, the asymptotic-preserving (AP) H N T -based unified gas kinetic scheme (UGKS). The constructed scheme is obtained by combing the UGKS for spatial discretization with the hybrid H N T method for angular discretization. Since the H N T is a hybrid angular discrete method of both P N and S N methods, the current H N T -based UGKS can not only mitigate the ray effects of the S N method largely, but also suppress the oscillations of the original P N method. Furthermore, we show that the current H N T -based UGKS also inherits the AP property of UGKS. A number of one-dimensional and two-dimensional numerical experiments are presented that validate the performance of the current scheme in both optically thin and thick regimes, as well as in mitigating the ray effects. Moreover, it can capture the initial layer solution without requiring additional treatments.