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

Pasmann, Samuel; Variansyah, Ilham; Kelley, C.T.; McClarren, Ryan G.

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

Jun 20, 2023

e-Print:

2306.11600 [physics.comp-ph]

View in:

ADS Abstract Service

pdf

reference search0 citations

Citations per year

0 Citations

Abstract: (arXiv)

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.

References(14)

Figures(6)

[1]

Monte Carlo and Quasi-Monte Carlo Sampling

P. Bickel
,
P. Diggle
,
S. Feinberg
,
U. Gather
,
I. Olkin

et al.

[2]

Quasi-Monte Carlo Methods for Particle Transport Problems.

J. Spanier

[3]

A Quasi-Monte Carlo Method With Krylov Linear Solvers for Multigroup Neutron Transport Simulations.

S. Pasmann
,
I. Variansyah
,
C.T. Kelley
,
R. McClarren

- https://doiorg/101080/0029563920222143704

[3]

pp. 1-15 . URL

S. Pasmann
,
I. Variansyah
,
C.T. Kelley
,
R. McClarren

[4]

The Davidson Method as an Alternative to Power Iterations for Criticality Calculations.

C. Subramanian
,
S.V. Criekingen
,
V. Heuveline
,
F. Nataf
,
P. Hav

[5]

3-D Neutron Transport Benchmarks.

T. Takeda
,
H. Ikeda

[6]

Development of MC/DC: a performant, scalable, and portable Python-based Monte Carlo neutron transport code.

I. Variansyah
,
J.P. Morgan
,
J. Northrop
,
K.E. Niemeyer
,
R.G. McClarren

[7]

Computational Methods of Neutron Transport. John

E. Miller
,
W.J. Lewis

[8]

Numerical Solution of the k-Eigenvalue Problem. Ph.D. thesis, Emory University

S.P. Hamilton

[9]

Krylov Subspace Iterations for Deterministic k-Eigenvalue Calculations.

J.S. Warsa
,
T.A. Wareing
,
J.E. Morel
,
J.M. McGhee
,
R.B. Lehoucq

[10]

The Iterative Calculation of a Few of the Lowest Eigenvalues and Corresponding Eigenvectors of Large Real-Symmetric Matrices.

E.R. Davidson

[11]

Iterative Methods for Linear and Nonlinear Equations

C.T. Kelley

[12]

Numba: A llvm-based python jit compiler.

S.K. Lam
,
A. Pitrou
,
S. Seibert

[13]

A 2-D/3-D Cartesian Geometry Non-Conforming Spherical Harmonic Neutron Transport Solver.

S.V. Criekingen