Multi-level Monte Carlo Finite Element method for elliptic PDEs with stochastic coefficients

被引:207
作者
Barth, Andrea [1 ]
Schwab, Christoph [1 ]
Zollinger, Nathaniel [1 ]
机构
[1] ETH Zentrum, Seminar Angew Math, CH-8092 Zurich, Switzerland
基金
欧洲研究理事会;
关键词
PARTIAL-DIFFERENTIAL-EQUATIONS; RANDOM INPUT DATA; COLLOCATION METHOD; CONSERVATIVE TRANSPORT; ADDITIVE NOISE; SIMULATION; APPROXIMATION; SPDES; FLOW;
D O I
10.1007/s00211-011-0377-0
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
In Monte Carlo methods quadrupling the sample size halves the error. In simulations of stochastic partial differential equations (SPDEs), the total work is the sample size times the solution cost of an instance of the partial differential equation. A Multi-level Monte Carlo method is introduced which allows, in certain cases, to reduce the overall work to that of the discretization of one instance of the deterministic PDE. The model problem is an elliptic equation with stochastic coefficients. Multi-level Monte Carlo errors and work estimates are given both for the mean of the solutions and for higher moments. The overall complexity of computing mean fields as well as k-point correlations of the random solution is proved to be of log-linear complexity in the number of unknowns of a single Multi-level solve of the deterministic elliptic problem. Numerical examples complete the theoretical analysis.
引用
收藏
页码:123 / 161
页数:39
相关论文
共 50 条
[21]   MULTILEVEL QUASI-MONTE CARLO INTEGRATION WITH PRODUCT WEIGHTS FOR ELLIPTIC PDES WITH LOGNORMAL COEFFICIENTS [J].
Herrmann, Lukas ;
Schwab, Christoph .
ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS-MODELISATION MATHEMATIQUE ET ANALYSE NUMERIQUE, 2019, 53 (05) :1507-1552
[22]   COMPUTING MEAN FIRST EXIT TIMES FOR STOCHASTIC PROCESSES USING MULTI-LEVEL MONTE CARLO [J].
Higham, Desmond J. ;
Roj, Mikolaj .
2012 WINTER SIMULATION CONFERENCE (WSC), 2012,
[23]   Robust Optimization of PDEs with Random Coefficients Using a Multilevel Monte Carlo Method [J].
Van Barel, Andreas ;
Vandewalle, Stefan .
SIAM-ASA JOURNAL ON UNCERTAINTY QUANTIFICATION, 2019, 7 (01) :174-202
[24]   A MULTISCALE DATA-DRIVEN STOCHASTIC METHOD FOR ELLIPTIC PDEs WITH RANDOM COEFFICIENTS [J].
Zhang, Zhiwen ;
Ci, Maolin ;
Hou, Thomas Y. .
MULTISCALE MODELING & SIMULATION, 2015, 13 (01) :173-204
[25]   A multilevel sparse kernel-based stochastic collocation finite element method for elliptic problems with random coefficients [J].
Dong, Zhaonan ;
Georgoulis, Emmanuil H. ;
Levesley, Jeremy ;
Usta, Fuat .
COMPUTERS & MATHEMATICS WITH APPLICATIONS, 2018, 76 (08) :1950-1965
[26]   A multilevel Monte Carlo finite element method for the stochastic Cahn-Hilliard-Cook equation [J].
Khodadadian, Amirreza ;
Parvizi, Maryam ;
Abbaszadeh, Mostafa ;
Dehghan, Mehdi ;
Heitzinger, Clemens .
COMPUTATIONAL MECHANICS, 2019, 64 (04) :937-949
[27]   Analysis of quasi-Monte Carlo methods for elliptic eigenvalue problems with stochastic coefficients [J].
Gilbert, A. D. ;
Graham, I. G. ;
Kuo, F. Y. ;
Scheichl, R. ;
Sloan, I. H. .
NUMERISCHE MATHEMATIK, 2019, 142 (04) :863-915
[28]   Estimation of Distribution Systems Expected Energy Not Supplied Index by Multi-level Monte Carlo Method [J].
Huda, A. S. Nazmul ;
Zivanovic, Rastko .
ELECTRIC POWER COMPONENTS AND SYSTEMS, 2019, 47 (9-10) :810-822
[29]   A finite element method for elliptic problems with stochastic input data [J].
Harbrecht, Helmut .
APPLIED NUMERICAL MATHEMATICS, 2010, 60 (03) :227-244
[30]   Reduced basis method for the adapted mesh and Monte Carlo methods applied to an elliptic stochastic problem [J].
Morcos, Noura ;
Sayah, Toni .
COMPUTATIONAL & APPLIED MATHEMATICS, 2019, 38 (02)