Parallel Multilevel Monte Carlo Algorithms for Elliptic PDEs with Random Coefficients

In this work, we developed and investigated Monte Carlo algorithms for elliptic PDEs with random coefficients. We considered groundwater flow as a model problem, where a permeability field represents random coefficients. The computational complexity is the main challenge in uncertainty quantification methods. The computation contains generating of a random coefficient and solving of partial differential equations. The permeability field was generated using the circulant embedding method. Multilevel Monte Carlo (MLMC) simulation can be based on different approximations of partial differential equations. We developed three MLMC algorithms based on finite volume, finite volume with renormalization and renormalization approximation. We compared numerical simulations and parallel performance of MLMC algorithms for 2D and 3D problems.