Second-order-accurate discrete ordinates solutions of transient radiative transfer in a scattering slab with variable refractive index

The discrete ordinates method (DOM) with a second-order upwind interpolation scheme is applied to solve transient radiative transfer in a graded index slab suddenly exposed to a diffuse strong irradiation at one of its boundaries. The planar medium is absorbing and anisotropically scattering. From the comparison of the results obtained by the first-order DOM, the second-order DOM, the modified DOM and the Monte Carlo method, it can be seen that the numerical diffusion in the transient solutions obtained by the second-order DOM is less than that in the solutions obtained by the first-order DOM, but the numerical diffusion is still noticeable, especially for optically thin and moderate cases. By contrast, for optically thick cases the numerical diffusion due to the finite difference of the advection term of the transient radiative transfer equation is minor. In general, it is still necessary to adopt a DOM with a higher order scheme to capture the wave front of transient radiative transfer accurately. Besides, the influence of numerical diffusion is a little less noticeable for the case with a larger gradient of refractive index, and the distribution of direction-integrated intensity around the irradiation boundary decreases and that around the other boundary increases with the increase of the anisotropically scattering coefficient.