Skip to main content

Research Repository

See what's under the surface

Advanced Search

A hierarchical finite element Monte Carlo method for stochastic two-scale elliptic equations

Brown, Donald L.; Hoang, Viet Ha

Authors

Donald L. Brown donald.brown@nottingham.ac.uk

Viet Ha Hoang



Abstract

We consider two-scale elliptic equations whose coefficients are random. In particular, we study two cases: in the first case, the coefficients are obtained from an ergodic dynamical system acting on a probability space, and in the second the case, the coefficients are periodic in the microscale but are random. We suppose that the coefficients also depend on the macroscopic slow variables. While the effective coefficient of the ergodic homogenization problem is deterministic, to approximate it, it is necessary to solve cell equations in a large but finite size “truncated" cube and compute an approximated effective coefficient from the solution of this equation. This approximated effective coefficient is, however, realization dependent; and the deterministic effective coefficient of the homogenization problem can be approximated by taking its expectation. In the periodic random setting, the effective coefficient for each realization are obtained from the solutions of cell equations which are posed in the unit cube, but to compute its average by the Monte Carlo method, we need to consider many uncorrelated realizations to accurately approximate the average. Straightforward employment of finite element approximation and the Monte Carlo method to compute this expectation with the same level of finite element resolution and the same number of Monte Carlo samples at every macroscopic point is prohibitively expensive. We develop a hierarchical finite element Monte Carlo algorithm to approximate the effective coefficients at a dense hierarchical network of macroscopic points. The method requires an optimal level of complexity that is essentially equal to that for computing the effective coefficient at one macroscopic point, and achieves essentially the same accuracy. The levels of accuracy for solving cell problems and for the Monte Carlo sampling are chosen according to the level in the hierarchy that the macroscopic points belong to. Solutions and the effective coefficients at the points where the cell problems are solved with higher accuracy and the effective coefficients are approximated with a larger number of Monte Carlo samples are employed as correctors for the effective coefficient at those points at which the cell problems are solved with lower accuracy and fewer Monte Carlo samples. The method combines the hierarchical finite element method for solving cell problems at a dense network of macroscopic points with the optimal complexity developed in D. L. Brown, Y. Efendiev and V. H. Hoang, Multiscale Model. Simul. 11 (2013), with a hierarchical Monte Carlo sampling algorithm that uses different number of samples at different macroscopic points depending on the level in the hierarchy that the macroscopic points belong to. Proof of concept numerical examples confirm the theoretical results.

Journal Article Type Article
Publication Date Oct 15, 2017
Journal Journal of Computational and Applied Mathematics
Print ISSN 0377-0427
Electronic ISSN 1879-1778
Publisher Elsevier
Peer Reviewed Peer Reviewed
Volume 323
APA6 Citation Brown, D. L., & Hoang, V. H. (2017). A hierarchical finite element Monte Carlo method for stochastic two-scale elliptic equations. Journal of Computational and Applied Mathematics, 323, doi:10.1016/j.cam.2017.04.004
DOI https://doi.org/10.1016/j.cam.2017.04.004
Keywords Multiscale; Hierarchical finite elements; Stochastic elliptic equations; Porous media
Publisher URL http://www.sciencedirect.com/science/article/pii/S0377042717301590
Copyright Statement Copyright information regarding this work can be found at the following address: http://creativecommons.org/licenses/by-nc-nd/4.0

Files

multimsmc_revised_FINAL.pdf (649 Kb)
PDF

Copyright Statement
Copyright information regarding this work can be found at the following address: http://creativecommons.org/licenses/by-nc-nd/4.0





You might also like



Downloadable Citations

;