The second-order formulation of the $P_N$ equations with Marshak boundary conditions

1 Nov 2019  ·  Matthias Andres, Florian Schneider ·

We consider a reformulation of the classical $P_N$ method with Marshak boundary conditions for the approximation of the monoenergetic stationary linear transport equation as a system of second-order PDEs. Our derivation allows the automatic generation of a model hierarchy which can then be handed to standard PDE tools. This method allows for heterogeneous coefficients, irregular grids, anisotropic boundary sources and anisotropic scattering. The wide applicability is demonstrated in several numerical test cases. We make our implementation available online, which allows for fast prototyping.

PDF Abstract

Categories


Numerical Analysis Numerical Analysis 35L40, 35Q70, 65M60, 65M70