Solution of Multimaterial Equilibrium Radiation Diffusion Problems by using the Discontinuous Galerkin Method

The discontinuous Galerkin method is used for solving the two-dimensional equilibrium radiation diffusion equation. We construct the weighted interior penalty method based on the geometric average weight. The semi-implicit integration factor method is applied to the nonlinear ordinary differential equations obtained by the discontinuous Galerkin spatial discretization. Numerical results are presented to demonstrate the validity and reliability of using the discontinuous Galerkin method for solving the highly nonlinear radiation diffusion equation.