A nodal method based on CMFD for pin-by-pin SP3 calculation

• Published in: Annals of nuclear energy Vol. 167; p. 108849
• Main Authors: Zhao, Wenbo, Chai, Xiaoming, Zhang, Bin, Liu, Jiayi, Peng, Xingjie, Li, Tianya, Yu, Junchong
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 01.03.2022
• Subjects: Coarse mesh finite difference; Jacobi-Davidson; Pin-by-pin; SP3; Coarse mesh finite difference; Jacobi-Davidson; SP3; Pin-by-pin
• ISSN: 0306-4549; 1873-2100
• DOI: 10.1016/j.anucene.2021.108849

•CMFD-like nodal method is developed to solve pin-wise SP3 and diffusion problems.•Correction factor of 2nd moment flux is neglected to ensure stability.•Jacobi-Davidson method is more efficient than other Krylov subspace methods.•CMFD with n × n mesh per cell is more accurate than FMFD with 2n × 2n.•1 × 1 mesh per cell is a suitable mesh division for pin-by-pin calculation. A nodal method based on Coarse Mesh Finite Difference (CMFD) is proposed for pin-by-pin core simulation with SP3 approximation. The Laplace operators of the 0th moment flux of SP3 equations are treated using CMFD, while the Laplace operators of the 2nd moment flux are treated using fine mesh finite difference (FMFD). Correction factor is determined by solving local two-node problems. Transverse 0th moment flux is expanded to second-order Legendre polynomials. SP3 equations are then transformed to a generalized eigenvalue problem which is solved using Krylov subspace methods, including Jacobi-Davidson method, Generalized Davidson method, and Krylov-Schur method. Standard nonlinear iterative strategy is carried out to obtain converged correction factor. A prototype code CORCA-PIN is developed. Numerical results show that solver option of CMFD using 1 × 1 radial mesh per cell and Jacobi Davidson iteration method is suitable for pin-by-pin whole core calculations with SP3 approximation.