Paper Info
Title
Multi-level Quasi-Monte Carlo Finite Element Methods for a Class of Elliptic PDEs with Random Coefficients
Abstract
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 $${\fancyscript{O}}(h^2)$$O(h2), where $$h$$h is the finest FE mesh width, or $${\fancyscript{O}}(N^{-1+\delta })$$O(N-1+ ) for arbitrary $$\delta >0$$ >0, where $$N$$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$$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.
Year
DOI
Venue
2015
10.1007/s10208-014-9237-5
Foundations of Computational Mathematics
Keywords
Field
DocType
Multi-level,Quasi-Monte Carlo methods,Infinite-dimensional integration,Elliptic partial differential equations with random coefficients,Finite element methods,65D30,65D32,65N30
Discretization,Mathematical optimization,Linear form,Mathematical analysis,Quasi-Monte Carlo method,Finite element method,Expected value,Parametric statistics,Elliptic partial differential equation,Mathematics,Randomness
Journal
Volume
Issue
ISSN
15
2
1615-3375
Citations 
PageRank 
References 
15
0.67
24
Authors
3
Name
Order
Citations
PageRank
Frances Y. Kuo147945.19
Christoph Schwab259558.38
Ian H. Sloan31180183.02