Author: Sam Pasmann
This codebase serves as a testbed for the iterative Quasi-Monte Carlo (iQMC) method for neutron transport. The theory behind iQMC is outlined in [1]. iQMC uses Quasi-Monte Carlo methods to solve successive iterations of the Source Iteration, Power Iteration and other advanced linear solvers for neutron transport.
Predefined problems may be run by editing and running any of the files in /scripts/. Adding new problems requires a new input and material file.
Currently the available generator variables are:
"random" : numpy's pseudo-random number generator
"sobol" : sobol sequence from scipy.qmc
"halton" : halton sequence from scipy.qmc
[1] Pasmann, S., Variansyah, I., and McClarren, R. G. Convergent transport source iteration calculations with quasi-monte carlo. vol. 124, American Nuclear Society, pp. 192–195.
[2] Garcia, R., Siewert, C., Radiative Transfer in Finite Inhomogeneous Plane-Parallel Atmospheres, J. Quantitative Spectroscopy & Radiative Transfer, 27, 2, pp. 141–148.