A fully coupled two-level Schwarz preconditioner based on smoothed aggregation for the transient multigroup neutron diffusion equations
The neutron multigroup diffusion equations is an important approximation of the neutron transport equation, and has been widely used for studying the motion of neutrons and their interactions with the stationary background materials. Solving the multigroup diffusion equations faces challenging because many variables are coupled together and a lot of computing resources are required. In this paper, we focus on the development of a scalable parallel algorithm framework for the system of equations arising from the discretization of the multigroup diffusion equations using a finite element method in space and using an implicit scheme in time. The parallel agorithm framweork consists of an inexact Jacobian-free Newton for the nonlinear equations, a Krylov subspace method for the Jacobian system and a parallel Schwarz preconditioner for the acceleration of the algorithm convergence. A coarse space based on smoothed aggregation is introduced to construct a two-level method to improve the parallel performance of the Schwarz preconditioner. We numerically shown that compared with the one-level method, the proposed two-level method works significantly better in terms of the number of linear iterations and the compute time for a system of equations with millions of unknowns and a supercomputer with thousands of processors.
- Research Organization:
- Idaho National Lab. (INL), Idaho Falls, ID (United States)
- Sponsoring Organization:
- USDOE Office of Nuclear Energy (NE)
- DOE Contract Number:
- AC07-05ID14517
- OSTI ID:
- 1478629
- Report Number(s):
- INL/CON-17-40867-Rev002
- Resource Relation:
- Journal Volume: 25; Journal Issue: 3; Conference: 18th Copper Mountain Conference on Multigrid Methods, 2017, Copper Mountain, Colorado, 03/26/2017 - 03/30/2017
- Country of Publication:
- United States
- Language:
- English
Similar Records
Scalability of preconditioners as a strategy for parallel computation of compressible fluid flow
Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation