Two-dimensional equations of the surface harmonics method for solving problems of spatial neutron kinetics in square-lattice reactors

Two-dimensional time-dependent finite-difference equations of the surface harmonics method (SHM) for the description of the neutron transport are derived for square-lattice reactors. These equations are implemented in the SUHAM-TD code. Verification of the derived equations and the developed code are performed by the example of known test problems, and the potential and efficiency of the SHM as applied to the solution of the time-dependent neutron transport equation in the diffusion approximation in two-dimensional geometry are demonstrated. These results show the substantial advantage of SHM over direct finite-difference modeling in computational costs.