Hybrid methods for radiation transport using diagonally implicit Runge-Kutta and space-time discontinuous Galerkin time integration

JOURNAL OF COMPUTATIONAL PHYSICS

卷:376

Article

Crockatt, Michael M.¹  Christlieb, Andrew J.^1,2  Garrett, C. Kristopher⁴  Hauck, Cory D.^3,5 

[1] Michigan State Univ, Dept Computat Math Sci & Engn, E Lansing, MI 48824 USA

[2] Michigan State Univ, Dept Math, E Lansing, MI 48824 USA

[3] Oak Ridge Natl Lab, Computat & Appl Math Grp, POB 2009, Oak Ridge, TN 37831 USA

[4] Los Alamos Natl Lab, Computat Phys & Methods Grp, Los Alamos, NM 87545 USA

[5] Univ Tennessee, Dept Math, Knoxville, TN 37996 USA

关键词: Hybrid methods;    Kinetic equations;    Fully implicit methods;    High-order accuracy;   

DOI  :  10.1016/j.jcp.2018.09.041

来源: Elsevier

【 摘 要 】

In this work, we describe extensions of a hybrid method for time-dependent linear, kinetic radiation transport problems to high-order time integration schemes of the diagonally-implicit Runge-Kutta (DIRK) and space-time discontinuous Galerkin (STDG) types. The hybrid methods are constructed by splitting the radiation flux into collided and uncollided components to which low- and high-resolution discrete ordinates approximations are applied. The efficiency of hybrid methods constructed using DIRK, STDG, integral deferred correction (IDC), and implicit Euler schemes is compared using a test problem in one-dimensional slab geometry containing material discontinuities. It is observed that (i) higher-order methods are more efficient than the implicit Euler method, often by an order of magnitude or more; (ii) third-order methods yield solutions of a given error in roughly half the time of the second-order method of the same type; and (iii) for a given order of accuracy it is found that the most efficient class of time integration scheme is STDG, followed by DIRK, with IDC the least efficient, for the test problem considered. Two test problems in two-dimensional xy-geometry are used to compare the computational efficiency of hybrid and standard discrete ordinates methods constructed with DIRK and STDG integrators. We observe that replacing a standard discrete ordinates method using an angular quadrature of order N with a hybrid discrete ordinates method using angular quadratures of order 2N and N/2 for the uncollided and collided fluxes, respectively, usually reduces overall solution time by a factor of 2 or more while simultaneously reducing the resulting solution error by a factor of 2 or more for the test problems considered. (C) 2018 Elsevier Inc. All rights reserved.

【 授权许可】

Free