Modern Monte Carlo variants for uncertainty quantification in neutron transport

Research output: Chapter in Book/Report/Conference proceedingChapter

2 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.

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
StatusE-pub ahead of print - 23 May 2018

ASJC Scopus subject areas

  • Mathematics(all)

Cite this

Graham, I. G., Parkinson, M. J., & Scheichl, R. (2018). Modern Monte Carlo variants for uncertainty quantification in neutron transport. In Contemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan (pp. 455-481). Springer International Publishing. https://doi.org/10.1007/978-3-319-72456-0_22

Modern Monte Carlo variants for uncertainty quantification in neutron transport. / Graham, Ivan G.; Parkinson, Matthew J.; Scheichl, Robert.

Contemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan. Springer International Publishing, 2018. p. 455-481.

Research output: Chapter in Book/Report/Conference proceedingChapter

Graham, IG, Parkinson, MJ & Scheichl, R 2018, Modern Monte Carlo variants for uncertainty quantification in neutron transport. in Contemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan. Springer International Publishing, pp. 455-481. https://doi.org/10.1007/978-3-319-72456-0_22
Graham IG, Parkinson MJ, Scheichl R. Modern Monte Carlo variants for uncertainty quantification in neutron transport. In Contemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan. Springer International Publishing. 2018. p. 455-481 https://doi.org/10.1007/978-3-319-72456-0_22
Graham, Ivan G. ; Parkinson, Matthew J. ; Scheichl, Robert. / Modern Monte Carlo variants for uncertainty quantification in neutron transport. Contemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan. Springer International Publishing, 2018. pp. 455-481
@inbook{2ab1ce42cb84498188b7db277d7f33ca,
title = "Modern Monte Carlo variants for uncertainty quantification in neutron transport",
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.",
author = "Graham, {Ivan G.} and Parkinson, {Matthew J.} and Robert Scheichl",
year = "2018",
month = "5",
day = "23",
doi = "10.1007/978-3-319-72456-0_22",
language = "English",
isbn = "9783319724553",
pages = "455--481",
booktitle = "Contemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan",
publisher = "Springer International Publishing",
address = "Switzerland",

}

TY - CHAP

T1 - Modern Monte Carlo variants for uncertainty quantification in neutron transport

AU - Graham, Ivan G.

AU - Parkinson, Matthew J.

AU - Scheichl, Robert

PY - 2018/5/23

Y1 - 2018/5/23

N2 - 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.

AB - 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.

UR - http://www.scopus.com/inward/record.url?scp=85049875184&partnerID=8YFLogxK

U2 - 10.1007/978-3-319-72456-0_22

DO - 10.1007/978-3-319-72456-0_22

M3 - Chapter

SN - 9783319724553

SP - 455

EP - 481

BT - Contemporary Computational Mathematics - A Celebration of the 80th Birthday of Ian Sloan

PB - Springer International Publishing

ER -