Diffusion Synthetic Acceleration for High-Order Discontinuous Finite Element SN Transport Schemes and Application to Locally Refined Unstructured Meshes

Abstract Diffusion synthetic acceleration (DSA) schemes compatible with adaptive mesh refinement (AMR) grids are derived for the SN transport equations discretized using high-order discontinuous finite elements. These schemes are directly obtained from the discretized transport equations by assuming a linear dependence in angle of the angular flux along with an exact Fick’s law and, therefore, are categorized as partially consistent. These schemes are akin to the symmetric interior penalty technique applied to elliptic problems and are all based on a second-order discontinuous finite element discretization of a diffusion equation (as opposed to a mixed or P1 formulation). Therefore, they only have the scalar flux as unknowns. A Fourier analysis has been carried out to determine the convergence properties of the three proposed DSA schemes for various cell optical thicknesses and aspect ratios. Out of the three DSA schemes derived, the modified interior penalty (MIP) scheme is stable and effective for realistic problems, even with distorted elements, but loses effectiveness for some highly heterogeneous configurations. The MIP scheme is also symmetric positive definite and can be solved efficiently with a preconditioned conjugate gradient method. Its implementation in an AMR SN transport code has been performed for both source iteration and GMRes-based transport solves, with polynomial orders up to 4. Numerical results are provided and show good agreement with the Fourier analysis results. Results on AMR grids demonstrate that the cost of DSA can be kept low on locally refined meshes.

[1]  Marvin L. Adams,et al.  Diffusion-synthetic acceleration given anisotropic scattering, general quadratures, and multidimensions , 1993 .

[2]  Edward W. Larsen,et al.  Diffusion-synthetic acceleration methods for discrete-ordinates problems , 1984 .

[3]  Jim E. Morel,et al.  Krylov Iterative Methods and the Degraded Effectiveness of Diffusion Synthetic Acceleration for Multidimensional SN Calculations in Problems with Material Discontinuities , 2004 .

[4]  P. Colella,et al.  An Adaptive Mesh Refinement Algorithm for the Radiative Transport Equation , 1998 .

[5]  Marvin L. Adams,et al.  Diffusion Synthetic Acceleration of Discontinuous Finite Element Transport Iterations , 1992 .

[6]  Roland W. Freund,et al.  QMRPACK: a package of QMR algorithms , 1996, TOMS.

[7]  Jean C. Ragusa,et al.  A two-mesh adaptive mesh refinement technique for SN neutral-particle transport using a higher-order DGFEM , 2010, J. Comput. Appl. Math..

[8]  Jean C. Ragusa,et al.  A high-order discontinuous Galerkin method for the SN transport equations on 2D unstructured triangular meshes , 2009 .

[9]  Edward W. Larsen,et al.  UNCONDITIONALLY STABLE DIFFUSION-SYNTHETIC ACCELERATION METHODS FOR THE SLAB GEOMETRY DISCRETE ORDINATES EQUATIONS. PART II: NUMERICAL RESULTS. , 1982 .

[10]  Y. Azmy.  Unconditionally stable and robust adjacent-cell diffusive preconditioning of weighted-difference particle transport methods is impossible , 2002 .

[11]  D. Arnold.  An Interior Penalty Finite Element Method with Discontinuous Elements , 1982 .

[12]  Shawn D. Pautz,et al.  Discontinuous Finite Element SN Methods on Three-Dimensional Unstructured Grids , 2001 .

[13]  B. Stamm.  Spectral Discontinuous Galerkin Method for Hyperbolic Problems , 2005 .

[14]  J. Nitsche.  Über ein Variationsprinzip zur Lösung von Dirichlet-Problemen bei Verwendung von Teilräumen, die keinen Randbedingungen unterworfen sind , 1971 .

[15]  Ludmil Zikatanov,et al.  A posteriori error estimator and AMR for discrete ordinates nodal transport methods , 2009 .

[16]  Khosro Shahbazi,et al.  An explicit expression for the penalty parameter of the interior penalty method , 2022 .

[17]  Randal S. Baker,et al.  A Diffusion Synthetic Acceleration Method for Block Adaptive Mesh Refinement , 2006 .

[18]  Christian Aussourd,et al.  Styx: A Multidimensional AMR SN Scheme , 2003 .

[19]  Jim E. Morel,et al.  Fully Consistent Diffusion Synthetic Acceleration of Linear Discontinuous SN Transport Discretizations on Unstructured Tetrahedral Meshes , 2002 .

[20]  J. Douglas,et al.  Interior Penalty Procedures for Elliptic and Parabolic Galerkin Methods , 1976 .

[21]  R. E. Alcouffe,et al.  Diffusion synthetic acceleration methods for the diamond-differenced discrete-ordinates equations , 1977 .

[22]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[23]  R. Baker.  A Block Adaptive Mesh Refinement Algorithm for the Neutral Particle Transport Equation , 2002 .

[24]  Edward W. Larsen,et al.  UNCONDITIONALLY STABLE DIFFUSION-SYNTHETIC ACCELERATION METHODS FOR THE SLAB GEOMETRY DISCRETE ORDINATES EQUATIONS. PART I: THEORY. , 1982 .

[25]  Edward W. Larsen,et al.  Fast iterative methods for discrete-ordinates particle transport calculations , 2002 .

[26]  Jean C. Ragusa,et al.  On the Convergence of DGFEM Applied to the Discrete Ordinates Transport Equation for Structured and Unstructured Triangular Meshes , 2009 .

[27]  Yaqi Wang,et al.  Standard and goal-oriented adaptive mesh refinement applied to radiation transport on 2D unstructured triangular meshes , 2011, J. Comput. Phys..