Document Type: Original Article
Faculty of Engineering, Shahid Beheshti University, G.C, P.O. Box: 1983963113, Tehran, Iran
The Jacobian-Free Newton-Krylov (JFNK) method has been widely used in solving nonlinear equations arising in many applications. In this paper, the JFNK solver is examined as an alternative to the traditional power iteration method for calculation of the fundamental eigenmode in reactor analysis based on even-parity neutron transport theory. Since the Jacobian is not formed the only extra storage required is associated with the workspace of the Krylov solver used at every Newton step. A new nonlinear function is developed for the even-parity neutron transport equation utilized to solve the eigenvalue problem using the JFNK. This Newton-based method is compared with the standard iterative power method for a number of multi-groups, one and two dimensional neutron transport benchmarks. The results show that the proposed algorithm generally ends with fewer iterations and shorter run times than those of the traditional power method.
- A new nonlinear function has been developed for solving the eigenvalue problem.
- The JFNK approximation casts as an accelerated iterative approach.
- Numerical results were generated using FEMPT code.
- The results indicate that the JFNK method can converge faster than the standard procedure.