Multi-level Quasi-Monte Carlo Finite Element Methods for a Class of Elliptic PDEs with Random Coefficients

• Published in: Foundations of computational mathematics Vol. 15; no. 2; pp. 411 - 449
• Main Authors: Kuo, Frances Y., Schwab, Christoph, Sloan, Ian H.
• Format: Journal Article
• Language: English
• Published: Boston Springer US 01.04.2015; Springer; Springer Nature B.V
• Subjects: Algorithms; Analysis; Applications of Mathematics; Computational mathematics; Computer Science; Discretization; Economics; Errors; Estimates; Finite element analysis; Finite element method; Linear and Multilinear Algebras; Math Applications in Computer Science; Mathematical analysis; Mathematical models; Mathematics; Mathematics and Statistics; Matrix Theory; Monte Carlo methods; Monte Carlo simulation; Multilevel analysis; Numerical Analysis; Partial differential equations; Sampling; United States; 65D30; Infinite-dimensional integration; Elliptic partial differential equations with random coefficients; 65D32; Multi-level; Quasi-Monte Carlo methods; Finite element methods; 65N30
• ISSN: 1615-3375; 1615-3383
• DOI: 10.1007/s10208-014-9237-5

This paper is a sequel to our previous work (Kuo et al. in SIAM J Numer Anal, 2012) where quasi-Monte Carlo (QMC) methods (specifically, randomly shifted lattice rules) are applied to finite element (FE) discretizations of elliptic partial differential equations (PDEs) with a random coefficient represented by a countably infinite number of terms. We estimate the expected value of some linear functional of the solution, as an infinite-dimensional integral in the parameter space. Here, the (single-level) error analysis of our previous work is generalized to a multi-level scheme, with the number of QMC points depending on the discretization level and with a level-dependent dimension truncation strategy. In some scenarios, it is shown that the overall error (i.e., the root-mean-square error averaged over all shifts) is of order O ( h 2 ) , where h is the finest FE mesh width, or O ( N - 1 + δ ) for arbitrary δ > 0 , where N denotes the maximal number of QMC sampling points in the parameter space. For these scenarios, the total work for all PDE solves in the multi-level QMC FE method is essentially of the order of one single PDE solve at the finest FE discretization level , for spatial dimension d = 2 with linear elements. The analysis exploits regularity of the parametric solution with respect to both the physical variables (the variables in the physical domain) and the parametric variables (the parameters corresponding to randomness). As in our previous work, families of QMC rules with “POD weights” (“product and order dependent weights”) which quantify the relative importance of subsets of the variables are found to be natural for proving convergence rates of QMC errors that are independent of the number of parametric variables.