Abstract: A numerical method is developed for two-dimensional nonequilibrium radiation diffusion equations.Discontinuous Galerkin method is applied in spatial diseretization in which numerical flux is constructed with weighted flux averages.Implicit-explicit integration factor method for time discretization is applied to nonlinear ordinary differential equations which is obtained with discontinuous Galerkin method. Radiation diffusion equations with multiple materials are solved on unstructured grids in numerical tests.It demonstrates that the method is effective for high nonlinear and tightly coupled radiation diffusion equations.