Abstract: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.

The numerical solution of space-dependent neutron kinetics equations in hexagonal-z geometry using Diagonally Implicit Runge Kutta method

The Numerical Solution of Space-Dependent Neutron Kinetics Equations in Hexagonal-Z Geometry Using Backward Differentiation Formula with Adaptive Step Size

Discontinuous Finite Element Method for Neutron Transport Equations on No Matching Mesh

A Three-Dimensional Flux Expansion Nodal Method for Hexagonal Geometry Application

Simulation of S transport equation for hexagonal-z reactor using the COMSOL Multiphysics software

Development of Three-Dimensional Time-Spatial Kinetics Code Basing on Hexagonal Nodal Code

Development of a neutron space–time kinetics solver with improved quasi-static method based on OpenFOAM

A Benchmark for Neutron Kinetics Methods in Hexagonal Geometry

Generalized Integral Method for Solving the Neutron Transport Equation with the Hybridized Discontinuous Galerkin Method

Development of 3D Multi-group Neutron Diffusion Code for Hexagonal Geometry

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

Application of Monte Carlo Methods to Neutron Space-Time Kinetics

Investigation of the second-order nodal expansion method for the forward and adjoint neutron flux in the hexagonal geometry

Research of 3-D Hexagonal Nodal Transport Method for Fast Reactor

Time Implicit Unified Gas Kinetic Scheme for 3D Multi-Group Neutron Transport Simulation

A symmetric Gauss-Seidel method for the steady-state Boltzmann equation

Improved quasi-static nodal Green's function method for cylindrical geometry

Delta-form-based method of solving high order spatial discretization schemes for neutron transport

An Accurate S$_N$ Method for Solving Static Multigroup Neutron Transport Equations in Slab Geometry

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

A stochastic method for transient, three-dimensional neutron transport