Numerical simulation for a one-dimensional radiative transfer equation
This numerical simulation has been developed for the diffusion approximation of a one-dimensional stationary radiative transfer equation. Specifically, as the mean free path of the photons tends to zero this simulation depicts the formation of boundary layers and the fact that at the interior of the domain the radiation intensity becomes isotropic solving a Laplace equation. This numerical simulation supports the results obtained analyzing the radiative transfer equation via matched asymptotic expansions.
This document contains the Python code for the simulation
This document describes the analyzed problem and the method used in order to develop the simulation.