Multiscale steady discrete unified gas kinetic scheme with macroscopic coarse mesh acceleration using preconditioned Krylov subspace method for multigroup neutron Boltzmann transport equation

A multiscale steady discrete unified gas kinetic scheme with macroscopic coarse mesh acceleration [accelerated steady discrete unified gas kinetic scheme (SDUGKS)] is proposed to improve the convergence of the original SDUGKS for an optically thick system in solving the multigroup neutron Boltzmann transport equation (NBTE) to analyze the distribution of fission energy in the reactor core. In the accelerated SDUGKS, by solving the coarse mesh macroscopic governing equations (MGEs) derived from the moment equations of the NBTE, the numerical solutions of the NBTE on fine meshes at the mesoscopic level can be rapidly obtained from the prolongation of the coarse mesh solutions of the MGE. Furthermore, the use of the coarse mesh can greatly reduce the computational variables and improve the computational efficiency of the MGE. The biconjugate gradient stabilized Krylov subspace method with the modified incomplete LU preconditioner and the lower-upper symmetric-Gauss-Seidel sweeping method are implemented to solve the discrete systems of the macroscopic coarse mesh acceleration model and mesoscopic SDUGKS to further improve the numerical efficiency. Numerical solutions validate good numerical accuracy and high acceleration efficiency of the proposed accelerated SDUGKS for the complicated multiscale neutron transport problems.