A boundary integral method to compute Green's functions for the radiative transport equation

The Green's function for the time-independent radiative transport equation in the whole space can be computed as an expansion in plane wave solutions. Plane wave solutions are a general class of solutions for the radiative transport equation. Because plane wave solutions are not known analytically in general, we calculate them numerically using the discrete ordinate method. We use the whole space Green's function to derive boundary integral equations. Through the solution of the boundary integral equations, we compute the Green's function for bounded domains. In particular we compute the Green's function for the half space, the slab, and the two-layered half space. The boundary conditions used here are in their most general form. Hence, this theory can be applied to boundaries with any kind of reflection and transmission law.