Abstract:A spatial second-order scheme for the nonlinear radiative transfer equations is introduced in this paper. The discretization scheme is based on the filtered spherical harmonics ($FP_N$) method for the angular variable and the unified gas kinetic scheme (UGKS) framework for the spatial and temporal variables respectively. In order to keep the scheme positive and second-order accuracy, firstly, we use the implicit Monte Carlo linearization method [6] in the construction of the UGKS numerical boundary fluxes. Then, by carefully analyzing the constructed second-order fluxes involved in the macro-micro decomposition, which is induced by the $FP_N$ angular discretization, we establish the sufficient conditions that guarantee the positivity of the radiative energy density and material temperature. Finally, we employ linear scaling limiters for the angular variable in the $P_N$ reconstruction and for the spatial variable in the piecewise linear slopes reconstruction respectively, which are shown to be realizable and reasonable to enforce the sufficient conditions holding. Thus, the desired scheme, called the $PPFP_N$-based UGKS, is obtained. Furthermore, in the regime $\epsilon\ll 1$ and the regime $\epsilon=O(1)$, a simplified spatial second-order scheme, called the $PPFP_N$-based SUGKS, is presented, which possesses all the properties of the non-simplified one. Inheriting the merit of UGKS, the proposed schemes are asymptotic preserving. By employing the $FP_N$ method for the angular variable, the proposed schemes are almost free of ray effects. To our best knowledge, this is the first time that spatial second-order, positive, asymptotic preserving and almost free of ray effects schemes are constructed for the nonlinear radiative transfer equations without operator splitting. Various numerical experiments are included to validate the properties of the proposed schemes.

A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations

An Improved Cylindrical FDTD Method and Its Application to Field-Tissue Interaction Study in MRI

A Novel High‐order Time‐domain Scheme for Three‐dimensional Maxwell's Equations

A Finite Element Method for Angular Discretization of the Radiation Transport Equation on Spherical Geodesic Grids

Using Spherical Harmonics to solve the Boltzmann equation: an operator based approach

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

An efficient solver for volumetric scattering based on fast spherical harmonics transforms

3D surface filtering using spherical harmonics

An Implicit Finite Volume Scheme to Solve the Time Dependent Radiation Transport Equation Based on Discrete Ordinates

Angular adaptivity with spherical harmonics for Boltzmann transport

Spherical Harmonics for the 1D Radiative Transfer Equation II: Thermal Emission

Spatial second-order positive and asymptotic preserving filtered $P_N$ schemes for nonlinear radiative transfer equations

Analysis of a Combined Spherical Harmonics and Discontinuous Galerkin Discretization for the Boltzmann Transport Equation

A simplified spherical harmonics approximated-radiation transport model for three-dimensional (3D) photon migration in small volume tissues

Spherical Harmonics for the 1D Radiative Transfer Equation. II. Thermal Emission

Finite Element Formulation of Exact Dirichlet-to-Neumann Radiation Conditions on Elliptic and Spheroidal Boundaries

Spherical Harmonics for the 1D Radiative Transfer Equation I: Reflected Light

Solving the complete pseudo-impulsive radiation and diffraction problem using a spectral element method

A 3D radiative transfer framework: V. Homologous Flows