Abstract
We present a new bilinear discontinuous (Galerkin) finite element discretization of the P1 (spherical harmonics) equations, a first order system of equations used for describing neutral particle radiation transport or modeling radiative transfer problems. The discrete equations are described for two-dimensional rectangular meshes; we solve the linear system with Krylov iterative methods. We have developed a novel, two-level preconditioner to improve convergence of the Krylov solvers that is based on a linear continuous finite element discretization of the diffusion equation, solved with a conjugate gradient iteration, preceded and followed by one of several different smoothing relaxations. A Fourier analysis shows that our approach is very effective over a wide range of problems. Numerical experiments confirm the results of the Fourier analysis. Computations for a realistic problem show that the preconditioner is effective and the solution method is efficient in practice.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.