Radiation diffusion equations arise in a wide range of physical applications such as astrophysics and inertial confinement fusion. Two key ingredients in constructing their numerical schemes are both positivity-preserving property and good accuracy for discontinuous and/or anisotropic problems on arbitrary distorted meshes. This talk will introduce our recent research progress including positivity-preserving finite volume schemes, theoretical analysis, and numerical applications.