Mitigating Spatial Error in the iterative-Quasi-Monte Carlo (iQMC) Method for Neutron Transport Simulations with Linear Discontinuous Source Tilting and Effective Scattering and Fission Rate Tallies
The iterative Quasi-Monte Carlo (iQMC) method is a recently proposed method for multigroup neutron transport simulations. iQMC can be viewed as a hybrid between deterministic iterative techniques, Monte Carlo simulation, and Quasi-Monte Carlo techniques. iQMC holds several algorithmic characteristics that make it desirable for high performance computing environments including a $O(N^{-1})$ convergence scheme, ray tracing transport sweep, and highly parallelizable nature similar to analog Monte Carlo. While there are many potential advantages of using iQMC there are also inherent disadvantages, namely the spatial discretization error introduced from the use of a mesh across the domain. This work introduces two significant modifications to iQMC to help reduce the spatial discretization error. The first is an effective source transport sweep, whereby the source strength is updated on-the-fly via an additional tally. This version of the transport sweep is essentially agnostic to the mesh, material, and geometry. The second is the addition of a history-based linear discontinuous source tilting method. Traditionally, iQMC utilizes a piecewise-constant source in each cell of the mesh. However, through the proposed source tilting technique iQMC can utilize a piecewise-linear source in each cell and reduce spatial error without refining the mesh. Numerical results are presented from the 2D C5G7 and Takeda-1 k-eigenvalue benchmark problems. Results show that the history-based source tilting significantly reduces error in global tallies and the eigenvalue solution in both benchmarks. Through the effective source transport sweep and linear source tilting iQMC was able to converge the eigenvalue from the 2D C5G7 problem to less than $0.04\%$ error on a uniform Cartesian mesh with only $204\times204$ cells.