Improved block-Jacobi parallel algorithm for the S_N nodal method with unstructured mesh

The substantial time required for solving the neutron transport equation can be reduced through parallel calculation. The block-Jacobi algorithm is a common method for parallel neutron transport calculation in an unstructured mesh; however, iteration degradation is an inevitable problem that limits the application of this algorithm by reducing the parallel efficiency. In this study, we applied the block-Jacobi algorithm to the S_N nodal method with a triangular-z mesh and proposed an improvement to achieve high parallel efficiency. The interface prediction (IP) method was developed to prevent iteration degradation. This method is based on extrapolating the interface information instead of using the information from the preceding iteration at the interfaces in sub-domains. Meanwhile, the prediction was used recursively to accelerate the self-group scattering source iteration in the entire space; this process is referred to as the inner iteration prediction method (IIP). These two methods effectively prevent iteration degradation and reduce time required for self-group scattering source iteration. A more stable and improved parallel performance was thus achieved.