Scalable hierarchical PDE sampler for generating spatially correlated random fields using nonmatching meshes: Scalable hierarchical PDE sampler using nonmatching meshes
Abstract
Summary This work describes a domain embedding technique between two nonmatching meshes used for generating realizations of spatially correlated random fields with applications to large‐scale sampling‐based uncertainty quantification. The goal is to apply the multilevel Monte Carlo (MLMC) method for the quantification of output uncertainties of PDEs with random input coefficients on general and unstructured computational domains. We propose a highly scalable, hierarchical sampling method to generate realizations of a Gaussian random field on a given unstructured mesh by solving a reaction–diffusion PDE with a stochastic right‐hand side. The stochastic PDE is discretized using the mixed finite element method on an embedded domain with a structured mesh, and then, the solution is projected onto the unstructured mesh. This work describes implementation details on how to efficiently transfer data from the structured and unstructured meshes at coarse levels, assuming that this can be done efficiently on the finest level. We investigate the efficiency and parallel scalability of the technique for the scalable generation of Gaussian random fields in three dimensions. An application of the MLMC method is presented for quantifying uncertainties of subsurface flow problems. We demonstrate the scalability of the sampling method with nonmatching mesh embedding, coupled with a parallelmore »
- Authors:
-
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Center for Applied Scientificc Computing
- Univ. della Svizzera Italiana, Lugano (Switzerland). Inst. of Computational Science
- Univ. of Texas, Austin, TX (United States). Inst. for Computational Engineering and Sciences
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Center for Applied Scientificc Computing; Portland State Univ., Portland, OR (United States). Fariborz Maseeh Dept. of Mathematics and Statistics
- Publication Date:
- Research Org.:
- Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Org.:
- USDOE; Swiss Commission for Technology and Innovation; Swiss National Science Foundation (SNSF); US Army Research Office (ARO)
- OSTI Identifier:
- 1438783
- Alternate Identifier(s):
- OSTI ID: 1432429
- Report Number(s):
- LLNL-JRNL-731006
Journal ID: ISSN 1070-5325
- Grant/Contract Number:
- AC52-07NA27344; DMS-1619640; W911NF-15-1-0590
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Numerical Linear Algebra with Applications
- Additional Journal Information:
- Journal Volume: 25; Journal Issue: 3; Journal ID: ISSN 1070-5325
- Publisher:
- Wiley
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; multilevel methods; PDEs with random input data; PDE sampler; non-matching meshes; H(div) problems; mixed nite elements; uncertainty quanti cation; multilevel Monte Carlo
Citation Formats
Osborn, Sarah, Zulian, Patrick, Benson, Thomas, Villa, Umberto, Krause, Rolf, and Vassilevski, Panayot S. Scalable hierarchical PDE sampler for generating spatially correlated random fields using nonmatching meshes: Scalable hierarchical PDE sampler using nonmatching meshes. United States: N. p., 2018.
Web. doi:10.1002/nla.2146.
Osborn, Sarah, Zulian, Patrick, Benson, Thomas, Villa, Umberto, Krause, Rolf, & Vassilevski, Panayot S. Scalable hierarchical PDE sampler for generating spatially correlated random fields using nonmatching meshes: Scalable hierarchical PDE sampler using nonmatching meshes. United States. https://doi.org/10.1002/nla.2146
Osborn, Sarah, Zulian, Patrick, Benson, Thomas, Villa, Umberto, Krause, Rolf, and Vassilevski, Panayot S. Tue .
"Scalable hierarchical PDE sampler for generating spatially correlated random fields using nonmatching meshes: Scalable hierarchical PDE sampler using nonmatching meshes". United States. https://doi.org/10.1002/nla.2146. https://www.osti.gov/servlets/purl/1438783.
@article{osti_1438783,
title = {Scalable hierarchical PDE sampler for generating spatially correlated random fields using nonmatching meshes: Scalable hierarchical PDE sampler using nonmatching meshes},
author = {Osborn, Sarah and Zulian, Patrick and Benson, Thomas and Villa, Umberto and Krause, Rolf and Vassilevski, Panayot S.},
abstractNote = {Summary This work describes a domain embedding technique between two nonmatching meshes used for generating realizations of spatially correlated random fields with applications to large‐scale sampling‐based uncertainty quantification. The goal is to apply the multilevel Monte Carlo (MLMC) method for the quantification of output uncertainties of PDEs with random input coefficients on general and unstructured computational domains. We propose a highly scalable, hierarchical sampling method to generate realizations of a Gaussian random field on a given unstructured mesh by solving a reaction–diffusion PDE with a stochastic right‐hand side. The stochastic PDE is discretized using the mixed finite element method on an embedded domain with a structured mesh, and then, the solution is projected onto the unstructured mesh. This work describes implementation details on how to efficiently transfer data from the structured and unstructured meshes at coarse levels, assuming that this can be done efficiently on the finest level. We investigate the efficiency and parallel scalability of the technique for the scalable generation of Gaussian random fields in three dimensions. An application of the MLMC method is presented for quantifying uncertainties of subsurface flow problems. We demonstrate the scalability of the sampling method with nonmatching mesh embedding, coupled with a parallel forward model problem solver, for large‐scale 3D MLMC simulations with up to 1.9·10 9 unknowns.},
doi = {10.1002/nla.2146},
journal = {Numerical Linear Algebra with Applications},
number = 3,
volume = 25,
place = {United States},
year = {2018},
month = {1}
}
Web of Science
Figures / Tables:

Works referenced in this record:
Tenth SPE Comparative Solution Project: A Comparison of Upscaling Techniques
journal, August 2001
- Christie, M. A.; Blunt, M. J.
- SPE Reservoir Evaluation & Engineering, Vol. 4, Issue 04
The multi-level Monte Carlo finite element method for a stochastic Brinkman Problem
journal, March 2013
- Gittelson, Claude J.; Könnö, Juho; Schwab, Christoph
- Numerische Mathematik, Vol. 125, Issue 2
Multilevel Monte Carlo methods and applications to elliptic PDEs with random coefficients
journal, January 2011
- Cliffe, K. A.; Giles, M. B.; Scheichl, R.
- Computing and Visualization in Science, Vol. 14, Issue 1
The Construction of the Coarse de Rham Complexes with Improved Approximation Properties
journal, January 2014
- Lashuk, Ilya V.; Vassilevski, Panayot S.
- Computational Methods in Applied Mathematics, Vol. 14, Issue 2
P3DFFT: A Framework for Parallel Computations of Fourier Transforms in Three Dimensions
journal, January 2012
- Pekurovsky, Dmitry
- SIAM Journal on Scientific Computing, Vol. 34, Issue 4
A Multilevel, Hierarchical Sampling Technique for Spatially Correlated Random Fields
journal, January 2017
- Osborn, Sarah; Vassilevski, Panayot S.; Villa, Umberto
- SIAM Journal on Scientific Computing, Vol. 39, Issue 5
Nonlinear multigrid solvers exploiting AMGe coarse spaces with approximation properties
journal, October 2018
- Christensen, Max la Cour; Vassilevski, Panayot S.; Villa, Umberto
- Journal of Computational and Applied Mathematics, Vol. 340
OpenCL Based Parallel Algorithm for RBF-PUM Interpolation
journal, April 2017
- Cavoretto, Roberto; Schneider, Teseo; Zulian, Patrick
- Journal of Scientific Computing, Vol. 74, Issue 1
Fast and Exact Simulation of Stationary Gaussian Processes through Circulant Embedding of the Covariance Matrix
journal, July 1997
- Dietrich, C. R.; Newsam, G. N.
- SIAM Journal on Scientific Computing, Vol. 18, Issue 4
Bayesian Spatial Modelling with R - INLA
journal, January 2015
- Lindgren, Finn; Rue, Håvard
- Journal of Statistical Software, Vol. 63, Issue 19
Element agglomeration coarse Raviart-Thomas spaces with improved approximation properties: COARSE RAVIART-THOMAS SPACES WITH IMPROVED APPROXIMATION PROPERTIES
journal, January 2012
- Lashuk, I. V.; Vassilevski, P. S.
- Numerical Linear Algebra with Applications, Vol. 19, Issue 2
A Computational Framework for Infinite-Dimensional Bayesian Inverse Problems Part I: The Linearized Case, with Application to Global Seismic Inversion
journal, January 2013
- Bui-Thanh, Tan; Ghattas, Omar; Martin, James
- SIAM Journal on Scientific Computing, Vol. 35, Issue 6
Perfect spatial hashing
journal, July 2006
- Lefebvre, Sylvain; Hoppe, Hugues
- ACM Transactions on Graphics, Vol. 25, Issue 3
Spatial variability and uncertainty in groundwater flow parameters: A geostatistical approach
journal, April 1979
- Delhomme, J. P.
- Water Resources Research, Vol. 15, Issue 2
Multilevel Techniques Lead to Accurate Numerical Upscaling and Scalable Robust Solvers for Reservoir Simulation
conference, February 2015
- Christensen, Max la Cour; Villa, Umberto; Vassilevski, Panayot
- SPE Reservoir Simulation Symposium
On Stationary Processes in the Plane
journal, January 1954
- Whittle, P.
- Biometrika, Vol. 41, Issue 3-4
Multilevel Monte Carlo Path Simulation
journal, June 2008
- Giles, Michael B.
- Operations Research, Vol. 56, Issue 3
Finite Element Error Analysis of Elliptic PDEs with Random Coefficients and Its Application to Multilevel Monte Carlo Methods
journal, January 2013
- Charrier, J.; Scheichl, R.; Teckentrup, A. L.
- SIAM Journal on Numerical Analysis, Vol. 51, Issue 1
Randomized algorithms for generalized Hermitian eigenvalue problems with application to computing Karhunen-Loève expansion: RANDOMIZED ALGORITHMS FOR GHEP
journal, November 2015
- Saibaba, Arvind K.; Lee, Jonghyun; Kitanidis, Peter K.
- Numerical Linear Algebra with Applications, Vol. 23, Issue 2
In order to make spatial statistics computationally feasible, we need to forget about the covariance function: SPDES, GMRFS, AND KERNEL METHODS
journal, October 2011
- Simpson, Daniel; Lindgren, Finn; Rue, Håvard
- Environmetrics, Vol. 23, Issue 1
PFFT: An Extension of FFTW to Massively Parallel Architectures
journal, January 2013
- Pippig, Michael
- SIAM Journal on Scientific Computing, Vol. 35, Issue 3
Further analysis of multilevel Monte Carlo methods for elliptic PDEs with random coefficients
journal, March 2013
- Teckentrup, A. L.; Scheichl, R.; Giles, M. B.
- Numerische Mathematik, Vol. 125, Issue 3
Exact de Rham Sequences of Spaces Defined on Macro-Elements in Two and Three Spatial Dimensions
journal, January 2008
- Pasciak, Joseph E.; Vassilevski, Panayot S.
- SIAM Journal on Scientific Computing, Vol. 30, Issue 5
An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach: Link between Gaussian Fields and Gaussian Markov Random Fields
journal, August 2011
- Lindgren, Finn; Rue, Håvard; Lindström, Johan
- Journal of the Royal Statistical Society: Series B (Statistical Methodology), Vol. 73, Issue 4
Multi-level Monte Carlo Finite Element method for elliptic PDEs with stochastic coefficients
journal, April 2011
- Barth, Andrea; Schwab, Christoph; Zollinger, Nathaniel
- Numerische Mathematik, Vol. 119, Issue 1
Upscaling of Mixed Finite Element Discretization Problems by the Spectral AMGe Method
journal, January 2016
- Kalchev, D. Z.; Lee, C. S.; Villa, U.
- SIAM Journal on Scientific Computing, Vol. 38, Issue 5
A Parallel Approach to the Variational Transfer of Discrete Fields between Arbitrarily Distributed Unstructured Finite Element Meshes
journal, January 2016
- Krause, Rolf; Zulian, Patrick
- SIAM Journal on Scientific Computing, Vol. 38, Issue 3
Scheduling Massively Parallel Multigrid for Multilevel Monte Carlo Methods
journal, January 2017
- Drzisga, D.; Gmeiner, B.; Rüde, U.
- SIAM Journal on Scientific Computing, Vol. 39, Issue 5
Interior penalty preconditioners for mixed finite element approximations
of elliptic problems
journal, April 1996
- Rusten, Torgeir; Vassilevski, Panayot S.; Winther, Ragnar
- Mathematics of Computation, Vol. 65, Issue 214
Numerical Multilevel Upscaling for Incompressible Flow in Reservoir Simulation: An Element-Based Algebraic Multigrid (AMGe) Approach
journal, January 2017
- la Cour Christensen, Max; Villa, Umberto; Engsig-Karup, Allan P.
- SIAM Journal on Scientific Computing, Vol. 39, Issue 1
Random numbers for large-scale distributed Monte Carlo simulations
journal, June 2007
- Bauke, Heiko; Mertens, Stephan
- Physical Review E, Vol. 75, Issue 6
A Note on Preconditioning for Indefinite Linear Systems
journal, January 2000
- Murphy, Malcolm F.; Golub, Gene H.; Wathen, Andrew J.
- SIAM Journal on Scientific Computing, Vol. 21, Issue 6
A Block Circulant Embedding Method for Simulation of Stationary Gaussian Random Fields on Block-Regular Grids
journal, January 2015
- Park, Min Ho; Tretyakov, M. V.
- International Journal for Uncertainty Quantification, Vol. 5, Issue 6
Perfect spatial hashing
conference, January 2006
- Lefebvre, Sylvain; Hoppe, Hugues
- ACM SIGGRAPH 2006 Papers on - SIGGRAPH '06
The multi-level Monte Carlo finite element method for a stochastic Brinkman Problem
text, January 2013
- Gittelson, Claude J.; Könnö, Juho; Schwab, Christoph
- ETH Zurich
On Stationary Processes in the Plane
journal, December 1954
- Whittle, P.
- Biometrika, Vol. 41, Issue 3/4
Spatial Variation.
journal, January 1988
- Kemp, C. D.; Matern, Bertil
- The Statistician, Vol. 37, Issue 1
The multi-level Monte Carlo finite element method for a stochastic Brinkman problem
text, January 2011
- Gittelson, Claude Jeffrey; Könnö, Juho; Schwab, Christoph
- ETH Zurich
Multi-Level Monte Carlo Finite Element Method for elliptic PDEs with stochastic coefficients
text, January 2011
- Barth, Andrea; Schwab, Christoph; Zollinger, Nathaniel
- ETH Zurich
Further analysis of multilevel Monte Carlo methods for elliptic PDEs with random coefficients
preprint, January 2012
- Teckentrup, A. L.; Scheichl, R.; Giles, M. B.
- arXiv
Scheduling massively parallel multigrid for multilevel Monte Carlo methods
preprint, January 2016
- Gmeiner, Björn; Drzisga, Daniel; Ruede, Ulrich
- arXiv
P3DFFT: a framework for parallel computations of Fourier transforms in three dimensions
text, January 2019
- Pekurovsky, Dmitry
- arXiv
Random numbers for large scale distributed Monte Carlo simulations
text, January 2006
- Bauke, Heiko; Mertens, Stephan
- arXiv
Works referencing / citing this record:
Multilevel approximation of Gaussian random fields: Fast simulation
journal, December 2019
- Herrmann, Lukas; Kirchner, Kristin; Schwab, Christoph
- Mathematical Models and Methods in Applied Sciences, Vol. 30, Issue 01
Multigrid Methods 2017: Multigrid Methods 2017
journal, February 2018
- Dendy, Joel E.
- Numerical Linear Algebra with Applications, Vol. 25, Issue 3
Figures / Tables found in this record: