Improvements of the CMFD acceleration capability of OpenMOC

Due to its computational efficiency and geometrical flexibility, the Method of Characteristics (MOC) has been widely used for light water reactor lattice physics analysis. Usually acceleration methods are necessary for MOC to achieve acceptable convergence on practical reactor physics problems. Among them, Coarse Mesh Finite Difference (CMFD) is very popular and can drastically reduce the number of transport iterations. In OpenMOC, CMFD acceleration was implemented but had the limitation of supporting only a uniform CMFD mesh, which would often lead to splitting MOC source regions, thus creating an unnecessary increase in computation and memory use. In this study, CMFD acceleration with a non-uniform Cartesian mesh is implemented into OpenMOC. We also propose a quadratic fit based CMFD prolongation method in the axial direction to further improve the acceleration when multiple MOC source regions are contained in one CMFD coarse mesh. Numerical results are presented to demonstrate the improvement of the CMFD acceleration capability in OpenMOC in terms of both efficiency and stability.