Discontinuous Galerkin finite element method for the continuum radiative transfer problem inside axis-symmetric circumstellar envelopes

Abstract: The study of the continuum radiative transfer problem inside circumstellar envelopes is both a theoretical and numerical challenge, especially in the frequency-dependent and multi-dimensional case. While approximate methods are easier to handle numerically, they often fail to accurately describe the radiation field inside complex geometries. For these cases, it is necessary to directly solve numerically the radiative transfer equation. We investigate the accuracy of a discontinuous Galerkin finite element method (DGFEM hereafter) applied to the frequency-dependent two dimensional radiative transfer equation, coupled with the radiative equilibrium equation, inside axis-symmetric circumstellar envelopes. The DGFEM is a variant of finite element methods. It employs discontinuous elements and flux integrals along their boundaries, ensuring local conservation. However, as opposed to the classical finite-element methods, the solution is discontinuous across element edges. We implemented the method in a code and tested its accuracy by comparing our results with the benchmarks from the literature. For all the tested cases, the temperature profile agrees within one percent. Additionally, the emerging spectral energy distributions (SEDs) and images, obtained subsequently by ray-tracing techniques from the DGFEM solution, agree on average within $5~\mathrm{\%}$ and $10~\mathrm{\%}$, respectively. We show that the DGFEM can accurately describe the temperature profile inside axis-symmetric circumstellar envelopes. Consecutively the emerging SEDs and images are also well reproduced. The discontinuous Galerkin finite element method provides an alternative method (other than Monte Carlo methods for instance) for solving the radiative transfer equation, and could be used in cases that are more difficult to handle with the other methods.