The partial differential equation describing equilibrium radiation diffusion is strongly nonlinear, which has been widely utilized in various fields such as astrophysics and others. The equilibrium radiation diffusion equation usually appears over multiple complicated domains, and the material characteristics vary between each domain. The diffusion coefficient near the interface is discontinuous. In this paper, the equilibrium radiation diffusion equation with discontinuous diffusion coefficient will be solved numerically by the unsymmetric radial basis function collocation method. The energy term T4 is linearized by utilizing the Picard-Newton and Richtmyer linearization methods on the basis of the fully implicit scheme discretization. And the successive permutation iteration and direct linearization methods are applied to linearize the diffusion terms. The accuracy of the proposed methods is proved by numerical experiments for regular and irregular domains with different types of interfaces.
Read full abstract