TY - JOUR

T1 - Parallel Domain Decomposition Strategies for Stochastic Elliptic Equations Part B: Accelerated Monte Carlo Sampling with Local PC Expansions

AU - Contreras, Andres A.

AU - Mycek, Paul

AU - Le Maître, Olivier P.

AU - Rizzi, Francesco

AU - Debusschere, Bert

AU - Knio, Omar M.

N1 - KAUST Repository Item: Exported on 2020-04-23
Acknowledgements: This work was supported by the U.S. Department of Energy (DOE), Office of Science, Office of Advanced Scientific Computing Research, under award number DE-SC0010540. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA-0003525. Support of the King Abdullah University of Science and Technology is also acknowledged.

PY - 2018/1/1

Y1 - 2018/1/1

N2 - Solving stochastic partial differential equations (SPDEs) can be a computationally intensive task, particularly when the underlying parametrization of the stochastic input field involves a large number of random variables. Direct Monte Carlo (MC) sampling methods are well suited for this type of situation since their cost is independent of the input complexity. Unfortunately, MC sampling methods suffer from slow convergence. In this manuscript, we propose an acceleration framework for elliptic SPDEs that relies on domain decomposition techniques and polynomial chaos (PC) expansions of local operators to reduce the cost of solving a SPDE via MC sampling. The approach exploits the fact that, at the subdomain level, the number of variables required to accurately parametrize the input stochastic field can be significantly reduced, as covered in detail in the prequel (Part A) to this paper. This makes it feasible to construct PC expansions of the local contributions to the condensed problem (i.e., the Schur complement of the discretized operator). The approach basically consists of two main stages: (1) a preprocessing stage in which PC expansions of the condensed problem are computed and (2) a Monte Carlo sampling stage where random samples of the solution are computed. The proposed method its naturally parallelizable. Extensive numerical tests are used to validate the methodology and assess its serial and parallel performance.

AB - Solving stochastic partial differential equations (SPDEs) can be a computationally intensive task, particularly when the underlying parametrization of the stochastic input field involves a large number of random variables. Direct Monte Carlo (MC) sampling methods are well suited for this type of situation since their cost is independent of the input complexity. Unfortunately, MC sampling methods suffer from slow convergence. In this manuscript, we propose an acceleration framework for elliptic SPDEs that relies on domain decomposition techniques and polynomial chaos (PC) expansions of local operators to reduce the cost of solving a SPDE via MC sampling. The approach exploits the fact that, at the subdomain level, the number of variables required to accurately parametrize the input stochastic field can be significantly reduced, as covered in detail in the prequel (Part A) to this paper. This makes it feasible to construct PC expansions of the local contributions to the condensed problem (i.e., the Schur complement of the discretized operator). The approach basically consists of two main stages: (1) a preprocessing stage in which PC expansions of the condensed problem are computed and (2) a Monte Carlo sampling stage where random samples of the solution are computed. The proposed method its naturally parallelizable. Extensive numerical tests are used to validate the methodology and assess its serial and parallel performance.

UR - http://hdl.handle.net/10754/631514

UR - https://epubs.siam.org/doi/10.1137/17M1132197

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

U2 - 10.1137/17M1132197

DO - 10.1137/17M1132197

M3 - Article

AN - SCOPUS:85053679739

VL - 40

SP - C547-C580

JO - SIAM Journal on Scientific Computing

JF - SIAM Journal on Scientific Computing

SN - 1064-8275

IS - 4

ER -