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 method of CADIS-NEE and CADIS-DXTRAN in NECP-MCX and its applications

The Study of a New Efficient Monte Carlo Method for Deep-Penetration Transport

Study on Unstructured Mesh–Based Monte Carlo/Deterministic Coupled Particle Transport Calculation Method

FW/CADIS-$Ω$: An angle-informed hybrid method for deep-penetration radiation transport

New variance reduction techniques for Monte-Carlo calculation based on first-collision source method

Preliminary study on CAD-based method of characteristics for neutron transport calculation

Improved Algorithms and Coupled Neutron-Photon Transport for Auto-Importance Sampling Method

Implementation and Verification of the DAGMC Module in Monte Carlo Code RMC

Linked calculation method based on next event estimation for directional shielding problems

Application of the First Collision Source Method to CSNS Target Station Shielding Calculation

Development of a Radiation Shielding Monte Carlo Code: RShieldMC

Study on Monte Carlo Variance-reduction Method Used for Maze Shielding Calculation

Monte Carlo Method for Calculating Particle Radiation Detector Correction Factors

Development of General Next Event Surface Crossing Estimators for Monte Carlo Method

Development of Multigroup Monte Carlo Neutron Transport Method with Regionwise Even-Parity Discontinuity Factor

An Auto-Importance Sampling Method for Deep Penetration Problems

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

Optimization of the Radiation Shielding Program QAD

Effectiveness of Radiation Transport Variance Reduction Methods for Wide-Area Environmental Contamination Assay Applications

Construction and Verification of Variance Reduction Technique System Based on Automatic Importance Sampling

Improved Adaptive Variance Reduction Algorithm Based on RMC Code for Deep Penetration Problems