Discontinuous finite element transport solutions in thick diffusive problems

The performance of discontinuous finite element methods (DFEMs) on problems that contain optically thick diffusive regions is analyzed and tested. The asymptotic analysis is quite general; it holds for an entire family of DFEMs in slab, XY, and XYZ geometries on arbitrarily connected polygonal or polyhedral spatial grids. The main contribution of the work is a theory that predicts and explains how DFEMs behave when applied to thick diffusive regions. It is well known that in the interior of such a region, the exact transport solution satisfies (to leading order) a diffusion equation,,vith boundary conditions that are known. Thus, in the interiors of such regions, the ideal discretized transport solution would satisfy (to leading order) an accurate discretization of the same diffusion equation and boundary conditions. The theory predicts that one class of DFEMs, which we call "zero-resolution" methods, fails dramatically in thick diffusive regions, yielding solutions that are completely meaningless. Another class-full-resolution methods-has leading-order solutions that satisfy discretizations of the correct diffusion equation. Full-resolution DFEMs are classified according to several categories of performance: continuity, robustness, accuracy and boundary condition. Certain kinds of lumping, some of which are believed to be new improve DFEM behavior in the continuity, robustness, and boundary-condition categories Theoretical results are illustrated using different variations of linear and bilinear DFEMs on several test problems in XY geometry. In every case, numerical results agree precisely with the predictions of the asymptotic theory.