Iterative and FEM methods to solve the 2?D Radiative Transfer Equation with specular reflexion

Eurotherm Conference 105: Computational Thermal Radiation in Participating Media V

物理学;计算机科学

Le Hardy, David^1 ; Favennec, Yann^1 ; Rousseau, Benoît^1

LTN UMR CNRS 6607, University of Nantes, France^1

关键词: Discontinuous Galerkin methods;    Discrete ordinate methods;    Heterogenous materials;    Iterative algorithm;    Radiative Transfer Equation(RTE);    Radiative transfer equations;    Streamlineupwind / petrov-galerkin methods (SUPG);    Successive over relaxation;   

Others  :  https://iopscience.iop.org/article/10.1088/1742-6596/676/1/012013/pdf
DOI  :  10.1088/1742-6596/676/1/012013

学科分类：计算机科学（综合）

来源: IOP

【 摘 要 】

The present paper deals with iterative algorithms coupled with finite element methods (FEM) to solve the Radiative Transfer Equation (RTE) within semi-transparent heterogenous materials where specular reflexions occur on their boundaries. As our intention is to use such solution for inversion, the forward model should be solved as fastly as possible. This communication compares, in terms of both accuracy and CPU, the Discontinuous Galerkin (DG) method with the Streamline Upwind Petrov-Galerkin (SUPG) method, both being coupled with the Discrete Ordinate Method. Next, several iteratives methods used to accelerate the convergence are compared. These methods are the Gauss-Siedel (GS), the Source-Iteration (SI) and the Successive Over-Relaxation (SOR) methods.