skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A Multilevel, Hierarchical Sampling Technique for Spatially Correlated Random Fields

Abstract

In this paper, we propose an alternative method to generate samples of a spatially correlated random field with applications to large-scale problems for forward propagation of uncertainty. A classical approach for generating these samples is the Karhunen--Loève (KL) decomposition. However, the KL expansion requires solving a dense eigenvalue problem and is therefore computationally infeasible for large-scale problems. Sampling methods based on stochastic partial differential equations provide a highly scalable way to sample Gaussian fields, but the resulting parametrization is mesh dependent. We propose a multilevel decomposition of the stochastic field to allow for scalable, hierarchical sampling based on solving a mixed finite element formulation of a stochastic reaction-diffusion equation with a random, white noise source function. Lastly, numerical experiments are presented to demonstrate the scalability of the sampling method as well as numerical results of multilevel Monte Carlo simulations for a subsurface porous media flow application using the proposed sampling method.

Authors:
 [1];  [1];  [2]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Center for Applied Scienti c Computing
  2. Univ. of Texas, Austin, TX (United States). Institute for Computational Engineering and Sciences
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1438756
Report Number(s):
LLNL-JRNL-696879; LLNL-JRNL-695979
Journal ID: ISSN 1064-8275
Grant/Contract Number:  
AC52-07NA27344
Resource Type:
Accepted Manuscript
Journal Name:
SIAM Journal on Scientific Computing
Additional Journal Information:
Journal Volume: 39; Journal Issue: 5; Journal ID: ISSN 1064-8275
Publisher:
SIAM
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; multilevel methods; PDEs with random input data; mixed nite elements; uncertainty quanti cation; multilevel Monte Carlo

Citation Formats

Osborn, Sarah, Vassilevski, Panayot S., and Villa, Umberto. A Multilevel, Hierarchical Sampling Technique for Spatially Correlated Random Fields. United States: N. p., 2017. Web. https://doi.org/10.1137/16M1082688.
Osborn, Sarah, Vassilevski, Panayot S., & Villa, Umberto. A Multilevel, Hierarchical Sampling Technique for Spatially Correlated Random Fields. United States. https://doi.org/10.1137/16M1082688
Osborn, Sarah, Vassilevski, Panayot S., and Villa, Umberto. Thu . "A Multilevel, Hierarchical Sampling Technique for Spatially Correlated Random Fields". United States. https://doi.org/10.1137/16M1082688. https://www.osti.gov/servlets/purl/1438756.
@article{osti_1438756,
title = {A Multilevel, Hierarchical Sampling Technique for Spatially Correlated Random Fields},
author = {Osborn, Sarah and Vassilevski, Panayot S. and Villa, Umberto},
abstractNote = {In this paper, we propose an alternative method to generate samples of a spatially correlated random field with applications to large-scale problems for forward propagation of uncertainty. A classical approach for generating these samples is the Karhunen--Loève (KL) decomposition. However, the KL expansion requires solving a dense eigenvalue problem and is therefore computationally infeasible for large-scale problems. Sampling methods based on stochastic partial differential equations provide a highly scalable way to sample Gaussian fields, but the resulting parametrization is mesh dependent. We propose a multilevel decomposition of the stochastic field to allow for scalable, hierarchical sampling based on solving a mixed finite element formulation of a stochastic reaction-diffusion equation with a random, white noise source function. Lastly, numerical experiments are presented to demonstrate the scalability of the sampling method as well as numerical results of multilevel Monte Carlo simulations for a subsurface porous media flow application using the proposed sampling method.},
doi = {10.1137/16M1082688},
journal = {SIAM Journal on Scientific Computing},
number = 5,
volume = 39,
place = {United States},
year = {2017},
month = {10}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 2 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Finite element exterior calculus, homological techniques, and applications
journal, May 2006


Random numbers for large-scale distributed Monte Carlo simulations
journal, June 2007


Spatial models generated by nested stochastic partial differential equations, with an application to global ozone mapping
journal, March 2011

  • Bolin, David; Lindgren, Finn
  • The Annals of Applied Statistics, Vol. 5, Issue 1
  • DOI: 10.1214/10-AOAS383

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
  • DOI: 10.1137/12089586X

Strong and Weak Error Estimates for Elliptic Partial Differential Equations with Random Coefficients
journal, January 2012

  • Charrier, Julia
  • SIAM Journal on Numerical Analysis, Vol. 50, Issue 1
  • DOI: 10.1137/100800531

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
  • DOI: 10.1007/s00791-011-0160-x

Multilevel Monte Carlo Path Simulation
journal, June 2008


A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs
journal, January 1998


A note on the optimality of the Karhunen-Loeve expansion
journal, May 1983


Parallel Auxiliary Space AMG Solver for $H(div)$ Problems
journal, January 2012

  • Kolev, Tzanio V.; Vassilevski, Panayot S.
  • SIAM Journal on Scientific Computing, Vol. 34, Issue 6
  • DOI: 10.1137/110859361

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
  • DOI: 10.1515/cmam-2014-0004

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
  • DOI: 10.1002/nla.1819

Bayesian Spatial Modelling with R - INLA
journal, January 2015

  • Lindgren, Finn; Rue, Håvard
  • Journal of Statistical Software, Vol. 63, Issue 19
  • DOI: 10.18637/jss.v063.i19

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
  • DOI: 10.1111/j.1467-9868.2011.00777.x

Preconditioning discretizations of systems of partial differential equations
journal, April 2010

  • Mardal, Kent-Andre; Winther, Ragnar
  • Numerical Linear Algebra with Applications, Vol. 18, Issue 1
  • DOI: 10.1002/nla.716

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
  • DOI: 10.1137/S1064827599355153

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
  • DOI: 10.1137/070698178

Interior penalty preconditioners for mixed finite element approximations of elliptic problems
journal, April 1996


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
  • DOI: 10.1002/nla.2026

Inverse problems: A Bayesian perspective
journal, May 2010


On Stationary Processes in the Plane
journal, January 1954


    Works referencing / citing this record:

    Scalable hierarchical PDE sampler for generating spatially correlated random fields using nonmatching meshes: Scalable hierarchical PDE sampler using nonmatching meshes
    journal, January 2018

    • Osborn, Sarah; Zulian, Patrick; Benson, Thomas
    • Numerical Linear Algebra with Applications, Vol. 25, Issue 3
    • DOI: 10.1002/nla.2146

    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
    • DOI: 10.1142/s0218202520500050

    Second-order control of complex systems with correlated synthetic data
    journal, November 2019


    Eigenvalue Problems for Exponential-Type Kernels
    journal, January 2020

    • Cai, Difeng; Vassilevski, Panayot S.
    • Computational Methods in Applied Mathematics, Vol. 20, Issue 1
    • DOI: 10.1515/cmam-2018-0186