Practical Monte Carlo simulation using modified power method with preconditioning

Abstract The authors developed the modified power method (MPM) in previous publications to obtain multiple eigenmodes of an eigenvalue problem at the same time by employing a generalized eigenvalue problem (GEP) of the form of WX = VXK. Special attention has been paid to the Monte Carlo (MC) implementation of the MPM because it always suffers from the inherent statistical noises. In this paper, a preconditioning method for the GEP has been developed for the MC MPM, so that the performance is more stable and robust to the MC statistical noises. This preconditioning method is crucial for MC solving of problems with degenerated eigenmodes, which requires the accumulation of a so-called transfer matrix and the division of the system space into multiple sub-regions, the number of sub-regions being greater than the number of eigenmodes to be solved. The preconditioning method solves the issues arising from the mismatch between the number of sub-regions and the target number of eigenmodes to be calculated. The numerical results for a model cube problem and BEAVRS whole core neutron transport eigenvalue problem successfully demonstrate the validity of the preconditioning method and the extended applicability of the MC MPM for practical problems.

[1]  Peng Zhang,et al.  On the Transfer Matrix of the Modified Power Method , 2018, Comput. Phys. Commun..

[2]  Thomas E. Booth,et al.  Computing the Higher k-Eigenfunctions by Monte Carlo Power Iteration: A Conjecture , 2003 .

[3]  Peng Zhang,et al.  A general solution strategy of modified power method for higher mode solutions , 2016, J. Comput. Phys..

[4]  Peng Zhang,et al.  Extension of modified power method to two-dimensional problems , 2016, J. Comput. Phys..

[5]  Peng Zhang,et al.  Extension of the noise propagation matrix method for higher mode solutions , 2017, J. Comput. Phys..

[6]  Philip R. Page,et al.  ENDF/B-VII.0: Next Generation Evaluated Nuclear Data Library for Nuclear Science and Technology , 2006 .

[7]  T. Booth.  Power Iteration Method for the Several Largest Eigenvalues and Eigenfunctions , 2006 .

[8]  Thomas E. Booth,et al.  Multiple extremal eigenpairs by the power method , 2008, J. Comput. Phys..

[9]  Thomas E. Booth,et al.  Exact Regional Monte Carlo Weight Cancellation for Second Eigenfunction Calculations , 2010 .

[10]  Peng Zhang,et al.  Validation of UNIST Monte Carlo code MCS for criticality safety analysis of PWR spent fuel pool and storage cask , 2018 .

[11]  Peng Zhang,et al.  Calculation of Degenerated Eigenmodes with Modified Power Method , 2017 .

[12]  Thomas E. Booth,et al.  Improved criticality convergence via a modified Monte Carlo iteration method , 2009 .

[13]  T E Booth,et al.  Monte Carlo determination of multiple extremal eigenpairs. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.