High-Order Variable Eddington Factor Methods for Thermal Radiative Transfer

We derive an arbitrary-order Discontinuous Galerkin Variable Eddington Factor method that is effective on a linearized thermal radiative transfer problem in both transport outer iterations and inner preconditioned BiCGStab iterations while supporting the use of a negative flux fixup inside the iteration. The inner iterations are shown to be uniform in the mesh size, polynomial order, and penalty parameter.