Stability analysis of nonlinear two-grid method for multigroup neutron diffusion problems

Abstract We present theoretical analysis of a nonlinear acceleration method for solving multigroup neutron diffusion problems. This method is formulated with two energy grids that are defined by (i) fine-energy groups structure and (ii) coarse grid with just a single energy group. The coarse-grid equations are derived by averaging of the multigroup diffusion equations over energy. The method uses a nonlinear prolongation operator. We perform stability analysis of iteration algorithms for inhomogeneous (fixed-source) and eigenvalue neutron diffusion problems. To apply Fourier analysis the equations of the method are linearized about solutions of infinite-medium problems. The developed analysis enables us to predict convergence properties of this two-grid method in different types of problems. Numerical results of problems in 2D Cartesian geometry are presented to confirm theoretical predictions.

[1]  Edward W. Larsen,et al.  A grey transport acceleration method for time-dependent radiative transfer pro lems , 1988 .

[2]  Edward W. Larsen,et al.  The Relationship between the Coarse-Mesh Finite Difference and the Coarse-Mesh Diffusion Synthetic Acceleration Methods , 2014 .

[3]  Dmitriy Y. Anistratov,et al.  Multilevel Quasidiffusion Methods for Solving Multigroup Neutron Transport k-Eigenvalue Problems in One-Dimensional Slab Geometry , 2011 .

[4]  Edward W. Larsen,et al.  A synthetic acceleration scheme for radiative diffusion calculations , 1985 .

[5]  D. A. Knoll,et al.  Acceleration of k-Eigenvalue/Criticality Calculations Using the Jacobian-Free Newton-Krylov Method , 2011 .

[6]  Dmitriy Y. Anistratov,et al.  Nonlinear Diffusion Acceleration Method with Multigrid in Energy for k-Eigenvalue Neutron Transport Problems , 2016 .

[7]  E. Lewis,et al.  Benchmark specification for Deterministic 2-D/3-D MOX fuel assembly transport calculations without spatial homogenisation (C5G7 MOX) , 2001 .

[9]  Ser Gi Hong,et al.  Fourier Convergence Analysis of the Rebalance Methods for Discrete Ordinates Transport Equations in Eigenvalue Problems , 2010 .

[10]  Jim E. Morel,et al.  A Two-Grid Acceleration Scheme for the Multigroup Sn Equations with Neutron Upscattering , 1993 .

[11]  Jim E. Morel,et al.  Linear multifrequency-grey acceleration recast for preconditioned Krylov iterations , 2007, J. Comput. Phys..

[12]  A. M. Winslow.  Multifrequency-gray method for radiation diffusion with Compton scattering , 1995 .

[13]  Dmitriy Anistratov,et al.  A Multilevel Projective Method for Solving the Space-Time Multigroup Neutron Kinetics Equations Coupled with the Heat Transfer Equation , 2014 .

[14]  Yousry Y. Azmy,et al.  Newton’s Method for Solving k-Eigenvalue Problems in Neutron Diffusion Theory , 2011 .

[15]  E. Aristova.  Simulation of Radiation Transport in a Channel Based on the Quasi-Diffusion Method , 2008 .

[16]  Dmitriy Y. Anistratov,et al.  Multilevel NDA Methods for Solving Multigroup Eigenvalue Neutron Transport Problems , 2013 .

[17]  Edward W. Larsen,et al.  Stability of Monte Carlo k-eigenvalue simulations with CMFD feedback , 2016, J. Comput. Phys..

[18]  Dmitriy Anistratov,et al.  Nonlinear methods for solving particle transport problems , 1993 .

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

[20]  正忠 小笠原.  J. J. Duderstadt and G. A. Moses: Inertial Confinement Fusion, John Wiley & Sons, New York, 1982, x+347ページ, 24×17cm, 16,000円. , 1983 .

[21]  Markus Berndt,et al.  A Nonlinear Krylov Accelerator for the Boltzmann k-Eigenvalue Problem , 2012 .