TY - GEN

T1 - Converged accelerated finite difference scheme for the multigroup neutron diffusion equation

AU - Terranova, Nicholas

AU - Mostacci, Domiziano

AU - Ganapol, Barry D.

N1 - Copyright:
Copyright 2013 Elsevier B.V., All rights reserved.

PY - 2013

Y1 - 2013

N2 - Computer codes involving neutron transport theory for nuclear engineering applications always require verification to assess improvement. Generally, analytical and semi-analytical benchmarks are desirable, since they are capable of high precision solutions to provide accurate standards of comparison. However, these benchmarks often involve relatively simple problems, usually assuming a certain degree of abstract modeling. In the present work, we show how semi-analytical equivalent benchmarks can be numerically generated using convergence acceleration. Specifically, we investigate the error behavior of a 1D spatial finite difference scheme for the multigroup (MG) steady-state neutron diffusion equation in plane geometry. Since solutions depending on subsequent discretization can be envisioned as terms of an infinite sequence converging to the true solution, extrapolation methods can accelerate an iterative process to obtain the limit before numerical instability sets in. The obtained results have been compared to the analytical solution to the 1D multigroup diffusion equation when available, using FORTRAN as the computational language. Finally, a slowing down problem has been solved using a cascading source update, showing how a finite difference scheme performs for ultra-fine groups (104 groups) in a reasonable computational time using convergence acceleration.

AB - Computer codes involving neutron transport theory for nuclear engineering applications always require verification to assess improvement. Generally, analytical and semi-analytical benchmarks are desirable, since they are capable of high precision solutions to provide accurate standards of comparison. However, these benchmarks often involve relatively simple problems, usually assuming a certain degree of abstract modeling. In the present work, we show how semi-analytical equivalent benchmarks can be numerically generated using convergence acceleration. Specifically, we investigate the error behavior of a 1D spatial finite difference scheme for the multigroup (MG) steady-state neutron diffusion equation in plane geometry. Since solutions depending on subsequent discretization can be envisioned as terms of an infinite sequence converging to the true solution, extrapolation methods can accelerate an iterative process to obtain the limit before numerical instability sets in. The obtained results have been compared to the analytical solution to the 1D multigroup diffusion equation when available, using FORTRAN as the computational language. Finally, a slowing down problem has been solved using a cascading source update, showing how a finite difference scheme performs for ultra-fine groups (104 groups) in a reasonable computational time using convergence acceleration.

KW - Convergence acceleration

KW - Extrapolation methods

KW - Multigroup

KW - Neutron diffusion equation

KW - Ultra-fine group

UR - http://www.scopus.com/inward/record.url?scp=84883340979&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84883340979&partnerID=8YFLogxK

M3 - Conference contribution

AN - SCOPUS:84883340979

SN - 9781627486439

T3 - International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, M and C 2013

SP - 2088

EP - 2102

BT - International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, M and C 2013

T2 - International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering, M and C 2013

Y2 - 5 May 2013 through 9 May 2013

ER -