Hybrid multi-level solvers for discontinuous Galerkin finite element discrete ordinate (DG-FEM-SN) diffusion synthetic acceleration (DSA) of radiation transport algorithms
File(s)1-s2.0-S0306454916305047-main.pdf (1.11 MB)
Published version
Author(s)
O'Malley, B
Kophazi, J
Smedley-Stevenson, RP
Eaton, MD
Type
Journal Article
Abstract
his paper examines two established preconditioners which were developed
to accelerate the solution of discontinuous Galerkin nite element method (DG-
FEM) discretisations of the elliptic neutron di usion equation. They are each
presented here as a potential way to accelerate the solution of the Modi ed In-
terior Penalty (MIP) form of the discontinuous di usion equation, for use as a
di usion synthetic acceleration (DSA) of DG-FEM discretisations of the neutron
transport equation. The preconditioners are both two-level schemes, di ering
in the low-level space utilised. Once projected to the low-level space a selection
of algebraic multigrid (AMG) preconditioners are utilised to obtain a further
correction step, these are therefore \hybrid" preconditioners. The rst precon-
ditioning scheme utilises a continuous piece-wise linear nite element method
(FEM) space, while the second uses a discontinuous piece-wise constant space.
Both projections are used alongside an element-wise block Jacobi smoother in
order to create a symmetric preconditioning scheme which may be used along-
side a conjugate gradient algorithm. An eigenvalue analysis reveals that both
should aid convergence but the piece-wise constant based method struggles with
some of the smoother error modes. Both are applied to a range of problems in-
cluding some which are strongly heterogeneous. In terms of conjugate gradient
(CG) iterations needed to reach convergence and computational time required,
both methods perform well. However, the piece-wise linear continuous scheme
appears to be the more e ective of the two. An analysis of computer memory
usage found that that the discontinuous piece-wise constant method had the
lowest memory requirements.
to accelerate the solution of discontinuous Galerkin nite element method (DG-
FEM) discretisations of the elliptic neutron di usion equation. They are each
presented here as a potential way to accelerate the solution of the Modi ed In-
terior Penalty (MIP) form of the discontinuous di usion equation, for use as a
di usion synthetic acceleration (DSA) of DG-FEM discretisations of the neutron
transport equation. The preconditioners are both two-level schemes, di ering
in the low-level space utilised. Once projected to the low-level space a selection
of algebraic multigrid (AMG) preconditioners are utilised to obtain a further
correction step, these are therefore \hybrid" preconditioners. The rst precon-
ditioning scheme utilises a continuous piece-wise linear nite element method
(FEM) space, while the second uses a discontinuous piece-wise constant space.
Both projections are used alongside an element-wise block Jacobi smoother in
order to create a symmetric preconditioning scheme which may be used along-
side a conjugate gradient algorithm. An eigenvalue analysis reveals that both
should aid convergence but the piece-wise constant based method struggles with
some of the smoother error modes. Both are applied to a range of problems in-
cluding some which are strongly heterogeneous. In terms of conjugate gradient
(CG) iterations needed to reach convergence and computational time required,
both methods perform well. However, the piece-wise linear continuous scheme
appears to be the more e ective of the two. An analysis of computer memory
usage found that that the discontinuous piece-wise constant method had the
lowest memory requirements.
Date Issued
2016-12-28
Date Acceptance
2016-11-29
Citation
Annals of Nuclear Energy, 2016, 102, pp.134-147
ISSN
1873-2100
Publisher
Elsevier Masson
Start Page
134
End Page
147
Journal / Book Title
Annals of Nuclear Energy
Volume
102
Copyright Statement
© 2016 The Authors. Published by Elsevier Ltd. This is an open access article under the CC-BY license (http://
creativecommons.org/licenses/by/4.0/
)
creativecommons.org/licenses/by/4.0/
)
Subjects
0915 Interdisciplinary Engineering
0299 Other Physical Sciences
Publication Status
Published