Truncated Fourier-series approximation of the time-domain radiative transfer equation using finite elements

2013 ◽  
Vol 30 (3) ◽  
pp. 470 ◽  
Author(s):  
Aki Pulkkinen ◽  
Tanja Tarvainen
Author(s):  
M.A. Badri ◽  
P. Jolivet ◽  
B. Rousseau ◽  
S. Le Corre ◽  
H. Digonnet ◽  
...  

2009 ◽  
Vol 132 (2) ◽  
Author(s):  
Gisela Widmer

The stationary monochromatic radiative transfer equation is posed in five dimensions, with the intensity depending on both a position in a three-dimensional domain as well as a direction. For nonscattering radiative transfer, sparse finite elements [2007, “Sparse Finite Elements for Non-Scattering Radiative Transfer in Diffuse Regimes,” ICHMT Fifth International Symposium of Radiative Transfer, Bodrum, Turkey; 2008, “Sparse Adaptive Finite Elements for Radiative Transfer,” J. Comput. Phys., 227(12), pp. 6071–6105] have been shown to be an efficient discretization strategy if the intensity function is sufficiently smooth. Compared with the discrete ordinates method, they make it possible to significantly reduce the number of degrees of freedom N in the discretization with almost no loss of accuracy. However, using a direct solver to solve the resulting linear system requires O(N3) operations. In this paper, an efficient solver based on the conjugate gradient method with a subspace correction preconditioner is presented. Numerical experiments show that the linear system can be solved at computational costs that are nearly proportional to the number of degrees of freedom N in the discretization.


Sign in / Sign up

Export Citation Format

Share Document