Multi-level Quasi-Monte Carlo Finite Element Methods for a Class of Elliptic PDEs with Random Coefficients
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 repr...
Uloženo v:
| Vydáno v: | Foundations of computational mathematics Ročník 15; číslo 2; s. 411 - 449 |
|---|---|
| Hlavní autoři: | , , |
| Médium: | Journal Article |
| Jazyk: | angličtina |
| Vydáno: |
Boston
Springer US
01.04.2015
Springer Nature B.V |
| Témata: | |
| ISSN: | 1615-3375, 1615-3383 |
| On-line přístup: | Získat plný text |
| Tagy: |
Přidat tag
Žádné tagy, Buďte první, kdo vytvoří štítek k tomuto záznamu!
|
| Shrnutí: | 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. |
|---|---|
| Bibliografie: | SourceType-Scholarly Journals-1 ObjectType-Feature-1 content type line 14 ObjectType-Article-1 ObjectType-Feature-2 content type line 23 |
| ISSN: | 1615-3375 1615-3383 |
| DOI: | 10.1007/s10208-014-9237-5 |