Numerical codes for the time-dependent radiative transfer equation by FDM-DOM

The numerical codes provide the calculation of the time-dependent radiative transfer equation (RTE) using the Galerkin quadrature method with the 6th order level-symmetric even quadrature set based on the finite difference and discrete ordinates methods for highly forward scattering media under the non-reentry boundary condition.

【Numerical schemes for the RTE-calculations】
・3rd order weighted essentially non oscillatory scheme (spatial discretization)
・3rd order total variation diminishing-Runge-Kutta method (temporal discretization)
・Galerkin quadrature method (angular discretization)

【The medium geometry and optical properties (default)】
・Rectangular medium with the source-detector distance of 0.4 cm
・Spatial grid size: 0.02 cm and timestep: 0.5 ps
・Absorption coefficient: 0.2 cm-1
・Scattering coefficient: 100 cm-1
・Anisotropic factor: 0.9
・Refractive index: 1.4