iQMC: Iterative Quasi-Monte Carlo for k-Eigenvalue Neutron Transport Simulations

The Iterative Quasi-Monte Carlo method, or iQMC, replaces standard quadrature techniques used in deterministic linear solvers with Quasi-Monte Carlo simulation for more accurate and efficient solutions to the neutron transport equation. This work explores employing iQMC in the Monte-Carlo Dynamic Code (MCDC) to solve k-eigenvalue problems for neutron transport with both the standard power iteration and the generalized Davidson method, a Krylov Subspace method. Results are verified with the 3-D, 2-group, Takeda-1 Benchmark problem.