Modern Monte Carlo variants for uncertainty quantification in neutron transport

Ivan G. Graham, Matthew J. Parkinson, Robert Scheichl

Research output: Chapter in Book/Report/Conference proceedingChapter

4 Citations (Scopus)

Abstract

We describe modern variants of Monte Carlo methods for Uncertainty Quantification (UQ) of the Neutron Transport Equation, when it is approximated by the discrete ordinates method with diamond differencing. We focus on the mono-energetic 1D slab geometry problem, with isotropic scattering, where the cross-sections are log-normal correlated random fields of possibly low regularity. The paper includes an outline of novel theoretical results on the convergence of the discrete scheme, in the cases of both spatially variable and random crosssections. We also describe the theory and practice of algorithms for quantifying the uncertainty of a functional of the scalar flux, using Monte Carlo and quasi- Monte Carlo methods, and their multilevel variants. A hybrid iterative/direct solver for computing each realisation of the functional is also presented. Numerical experiments show the effectiveness of the hybrid solver and the gains that are possible through quasi-Monte Carlo sampling and multilevel variance reduction. For the multilevel quasi-Monte Carlo method, we observe gains in the computational ε- cost of up to two orders of magnitude over the standard Monte Carlo method, and we explain this theoretically. Experiments on problems with up to several thousand stochastic dimensions are included.

Original languageEnglish
Title of host publicationContemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan
PublisherSpringer International Publishing
Pages455-481
Number of pages27
ISBN (Electronic)9783319724560
ISBN (Print)9783319724553
DOIs
Publication statusE-pub ahead of print - 23 May 2018

ASJC Scopus subject areas

  • Mathematics(all)

Fingerprint Dive into the research topics of 'Modern Monte Carlo variants for uncertainty quantification in neutron transport'. Together they form a unique fingerprint.

Cite this