Discrete spherical harmonics method for radiative transfer in scalar planar inhomogeneous atmosphere

The radiative transfer problems in a participating inhomogeneous scalar planar atmosphere, subjected to diffuse or collimated incidence, are investigated using the discrete spherical harmonics method. In developing the method, the radiative intensity is expanded in a finite series of Legendre polynomials and the resulting first-order coupled differential equations of radiance moments are expressed in a set of discrete polar directions. The method is applied to homogeneous/inhomogeneous atmospheres of various anisotropic scattering degrees and thicknesses, and reflective boundary conditions. The discrete spherical harmonics method albedo, transmittance, and radiative intensity predictions agree well with benchmark literature results. Additionally, numerical predictions show that the discrete spherical harmonics method using Mark boundary conditions are more efficient than using Marshak boundary conditions.