Modified spherical harmonics method for solving the radiative transport equation

A new effective approach to solving the three-dimensional radiative transport equation with an arbitrary phase function is proposed. The solution depends on eigenvectors and eigenvalues of several symmetrical tridiagonal matrices of infinite size. The matrices must be truncated and diagonalized numerically. Then, given eigenvectors and eigenvalues of these matrices, the dependence of the solution on position and direction is found analytically. The approach is based on expanding the angular part of the specific intensity in q-dependent spherical functions for each spatial Fourier component characterized by the vector q. Apart from the truncation of the matrices, no other approximations are made.