A general solution strategy of modified power method for higher mode solutions. (English) Zbl 1349.65135
Summary: A general solution strategy of the modified power iteration method for calculating higher eigenmodes has been developed and applied in continuous energy Monte Carlo simulation. The new approach adopts four features: 1) the eigen decomposition of transfer matrix, 2) weight cancellation for higher modes, 3) population control with higher mode weights, and 4) stabilization technique of statistical fluctuations using multi-cycle accumulations. The numerical tests of neutron transport eigenvalue problems successfully demonstrate that the new strategy can significantly accelerate the fission source convergence with stable convergence behavior while obtaining multiple higher eigenmodes at the same time. The advantages of the new strategy can be summarized as 1) the replacement of the cumbersome solution step of high order polynomial equations required by Booth’s original method with the simple matrix eigen decomposition, 2) faster fission source convergence in inactive cycles, 3) more stable behaviors in both inactive and active cycles, and 4) smaller variances in active cycles. Advantages 3 and 4 can be attributed to the lower sensitivity of the new strategy to statistical fluctuations due to the multi-cycle accumulations. The application of the modified power method to continuous energy Monte Carlo simulation and the higher eigenmodes up to 4th order are reported for the first time in this paper.
MSC:
| 65F15 | Numerical computation of eigenvalues and eigenvectors of matrices |
| 65C05 | Monte Carlo methods |
References:
| [1] | Booth, Thomas E., Computing the higher \(k\)-eigenfunctions by Monte Carlo power iteration: a conjecture, Nucl. Sci. Eng., 143, 291-300 (2003) |
| [2] | Booth, Thomas E., Power iteration method for the several largest eigenvalues and eigenfunctions, Nucl. Sci. Eng., 154, 48-62 (2006) |
| [3] | Gubernatis, J. E.; Booth, T. E., Multiple extremal eigenpairs by the power method, J. Comput. Phys., 227, 8508-8522 (2008) · Zbl 1152.65050 |
| [4] | Booth, T. E.; Gubernatis, J. E., Monte Carlo determination of multiple extremal eigenpairs (2008), Los Alamos report LA-UR-07-7672 · Zbl 1152.65050 |
| [5] | Booth, T. E.; Gubernatis, J. E., Multiple extremal eigenpairs of very large matrices by Monte Carlo simulation (2008), Los Alamos report LA-UR-08-0043 · Zbl 1152.65050 |
| [6] | Booth, Thomas E.; Gubernatis, James E., Improved criticality convergence via a modified Monte Carlo power iteration method, (International Conference on Mathematics, Computational Methods & Reactor Physics (M&C 2009). International Conference on Mathematics, Computational Methods & Reactor Physics (M&C 2009), Saratoga Springs, New York, May 3-7, 2009 (2009), American Nuclear Society: American Nuclear Society LaGrange Park, IL) |
| [7] | Booth, T. E.; Gubernatis, J. E., Monte Carlo determination of multiple extremal eigenpairs, Phys. Rev. E, 80, Article 046704 pp. (2009) |
| [8] | Booth, Thomas E.; Gubernatis, James E., Exact regional Monte Carlo weight cancellation for second eigenfunction calculations, Nucl. Sci. Eng., 165, 283-291 (2010) |
| [9] | Shi, Bo; Petrovic, Bojan, Implementation of the modified power iteration method to two-group Monte Carlo eigenvalue problems, Ann. Nucl. Energy, 38, 781-787 (2011) |
| [10] | Yamamoto, Toshihiro, Convergence of the second eigenfunction in Monte Carlo power iteration, Ann. Nucl. Energy, 36, 7-14 (2009) |
| [11] | Shi, Bo; Petrovic, Bojan, Calculating the second eigenpair in criticality calculations using the Monte Carlo Method with source points pairing as an efficient net-weight (cancellation) algorithm, Nucl. Sci. Eng., 172, 138-150 (2012) |
| [12] | Carney, Sean; Brown, Forrest; Kiedrowski, Brian; Martin, William, Theory and application of the fission matrix method for continuous-energy Monte Carlo, Ann. Nucl. Energy, 73, 423-431 (2014) |
| [13] | Tuttelberg, Kaur; Dufek, Jan, Estimation of errors in the cumulative Monte Carlo fission source, Ann. Nucl. Energy, 72, 151-155 (2014) |
| [14] | Brown, F. B.; Carney, S. E.; Kiedrowski, B. C.; Martin, W. R., Fission matrix capability for MCNP Monte Carlo, (Joint International Conference on Supercomputing in Nuclear Applications and Monte Carlo 2013 (SNA + MC 2013). Joint International Conference on Supercomputing in Nuclear Applications and Monte Carlo 2013 (SNA + MC 2013), Paris, France, October 27-31, LA-UR-13-26962 (2013)) |
| [15] | Colin, Josey; Veit Max, D., Eigenfunction decomposition of reactor perturbations and transitions using MCNP Monte Carlo (2013), Report LA-UR-13-26449, talk LA-UR-13-26464 |
| [16] | Carney, S. E.; Brown, F. B.; Kiedrowski, B. C.; Martin, W. R., Higher-mode applications of fission matrix capability for MCNP (2013), Report LA-UR-13-27078, talk LA-UR-13-26615 |
| [17] | Brown, F. B.; Carney, S. E.; Kiedrowski, B. C.; Martin, W. R., Fission matrix capability for MCNP, Part I - Theory, (Mathematics & Computation. Mathematics & Computation, Sun Valley, Idaho, May 5-9, 2013 (2013)) |
| [18] | Brown, F. B.; Carney, S. E.; Kiedrowski, B. C.; Martin, W. R., Fission matrix capability for MCNP, Part II - Applications, (Mathematics & Computation. Mathematics & Computation, Sun Valley, Idaho, May 5-9, 2013 (2013)) |
| [19] | Carney, S. E.; Brown, F. B.; Kiedrowski, B. C.; Martin, W. R., Fission matrix capability for MCNP Monte Carlo (2012), LA-UR-12-24533 |
| [20] | Ueki, T.; Brown, F. B.; Parsons, D. K.; Kornreich, D. K., Autocorrelation and dominance ratio in Monte Carlo criticality calculations, Nucl. Sci. Eng., 145, 279-290 (2003) |
| [21] | Lee, Hyunsuk; Kong, Chidong; Lee, Deokjung, Status of Monte Carlo code development at UNIST, (PHYSOR 2014 - The Role of Reactor Physics Toward a Sustainable Future. PHYSOR 2014 - The Role of Reactor Physics Toward a Sustainable Future, The Westin Miyako, Kyoto, Japan, September 28-October 3, 2014 (2014)), on CD-ROM |
| [22] | Chadwick, M. B., ENDF/B-VII.0: next generation evaluated nuclear data library for nuclear science and technology, Nucl. Data Sheets, 107, 12, 2931-3060 (2006) |
| [24] | Yamamoto, Toshihiro; Nakamura, Takemi; Miyoshi, Yoshinori, Fission source convergence of Monte Carlo criticality calculations in weakly coupled fissile arrays, J. Nucl. Sci. Technol., 37, 41-52 (2000) |
| [25] | Report of OECD/NEA Monte Carlo source convergence benchmark problem 3 (2006) |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
