Hostname: page-component-cd9895bd7-hc48f Total loading time: 0 Render date: 2024-12-29T04:21:35.640Z Has data issue: false hasContentIssue false
  • English
  • Français

Multilevel Markov Chain Monte Carlo Method for High-Contrast Single-Phase Flow Problems

Published online by Cambridge University Press:  19 December 2014

Yalchin Efendiev ,
Bangti Jin ,
Presho Michael  and
Xiaosi Tan
Yalchin Efendiev
Affiliation:
Department of Mathematics, Texas A&M University, College Station, TX 77843-3368, USA SRI Center for Numerical Porous Media (NumPor), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Kingdom of Saudi Arabia
Bangti Jin
Affiliation:
Department of Computer Science, University College London, Gower Street, London WC1E 6BT, UK
Presho Michael*
Affiliation:
Department of Mathematics, Texas A&M University, College Station, TX 77843-3368, USA
Xiaosi Tan
Affiliation:
Department of Mathematics, Texas A&M University, College Station, TX 77843-3368, USA
*
*Email addresses: [email protected] (Y. Efendiev), [email protected] (B. Jin), [email protected] (M. Presho), [email protected] (X. Tan)
Get access

Abstract

In this paper we propose a general framework for the uncertainty quantification of quantities of interest for high-contrast single-phase flow problems. It is based on the generalized multiscale finite element method (GMsFEM) and multilevel Monte Carlo (MLMC) methods. The former provides a hierarchy of approximations of different resolution, whereas the latter gives an efficient way to estimate quantities of interest using samples on different levels. The number of basis functions in the online GMsFEM stage can be varied to determine the solution resolution and the computational cost, and to efficiently generate samples at different levels. In particular, it is cheap to generate samples on coarse grids but with low resolution, and it is expensive to generate samples on fine grids with high accuracy. By suitably choosing the number of samples at different levels, one can leverage the expensive computation in larger fine-grid spaces toward smaller coarse-grid spaces, while retaining the accuracy of the final Monte Carlo estimate. Further, we describe a multilevel Markov chain Monte Carlo method, which sequentially screens the proposal with different levels of approximations and reduces the number of evaluations required on fine grids, while combining the samples at different levels to arrive at an accurate estimate. The framework seamlessly integrates the multiscale features of the GMsFEM with the multilevel feature of the MLMC methods following the work in [26], and our numerical experiments illustrate its efficiency and accuracy in comparison with standard Monte Carlo estimates.

Keywords

65N3076M10Generalized multiscale finite element methodmultilevel Monte Carlo method, multilevel Markov chain Monte Carlouncertainty quantification
Type
Research Article
Information
Communications in Computational Physics , Volume 17 , Issue 1 , January 2015 , pp. 259 - 286
Copyright
Copyright © Global Science Press Limited 2015 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Aarnes, J.E., Krogstad, S., and Lie, K.-A.. A hierarchicalmultiscalemethod for two-phase flow based upon mixed finite elements and nonuniform coarse grids. Multiscale Model. Simul., 5(2):337–363, 2006.CrossRefGoogle Scholar
[2]Arbogast, T., Pencheva, G., Wheeler, M.F., and Yotov, I.. A multiscale mortar mixed finite element method. Multiscale Model. Simul., 6(1):319–346, 2007.CrossRefGoogle Scholar
[3]Babuska, I. and Lipton, R.. Optimal local approximation spaces for generalized finite element methods with application to multiscale problems. Multiscale Model. Simul., 9(1):373–406, 2011.CrossRefGoogle Scholar
[4]Babušska, I., Tempone, R., and Zouraris, G.E.. Galerkin finite element approximations of stochastic elliptic partial differential equations. SIAM J. Numer. Anal., 42(2):800–825, 2004.CrossRefGoogle Scholar
[5]Bal, G., Langmore, I., and Marzouk, Y.. Bayesian inverse problems with Monte Carlo forward models. Inverse Probl. Imaging, 7(1):81–105, 2013.CrossRefGoogle Scholar
[6]Barrault, M., Maday, Y., Nguyen, N.C., and Patera, A.T.. An “empirical interpolation” method: application to efficient reduced-basis discretization of partial differential equations. C. R. Math. Acad. Sci. Paris, 339(9):667–672, 2004.CrossRefGoogle Scholar
[7]Barth, A., Schwab, C., and Zollinger, N.. Multi-level Monte Carlo finite element method for elliptic PDEs with stochastic coefficients. Numer. Math., 119(1):123–161, 2011.CrossRefGoogle Scholar
[8]Boyaval, S.. Reduced-basis approach for homogenization beyond the periodic setting. Multiscale Model. Simul., 7(1):466–494, 2008.CrossRefGoogle Scholar
[9]Christen, J.A. and Fox, C.. Markov chain Monte Carlo using an approximation. J. Comput. Graph. Statist., 14(4):795–810, 2005.CrossRefGoogle Scholar
[10]Cliffe, K.A., Giles, M.B., Scheichl, R., and Teckentrup, A.L.. Multilevel Monte Carlo methods and applications to elliptic PDEs with random coefficients. Comput. Vis. Sci., 14(1):3–15, 2011.CrossRefGoogle Scholar
[11]Efendiev, Y., Galvis, J., and Thomines, F.. A systematic coarse-scale model reduction technique for parameter-dependent flows in highly heterogeneous media and its applications. Multiscale Model. Simul., 10(4):1317–1343, 2012.CrossRefGoogle Scholar
[12]Efendiev, Y., Galvis, J., and Wu, X.-H.. Multiscale finite element methods for high-contrast problems using local spectral basis functions. J. Comput. Phys., 230(4):937–955, 2011.CrossRefGoogle Scholar
[13]Efendiev, Y., Ginting, V., Hou, T., and Ewing, R.. Accurate multiscale finite element methods for two-phase flow simulations. J. Comput. Phys., 220(1):155–174, 2006.CrossRefGoogle Scholar
[14]Efendiev, Y., Hou, T., and Ginting, V.. Multiscale finite element methods for nonlinear problems and their applications. Commun. Math. Sci., 2(4):553–589, 2004.CrossRefGoogle Scholar
[15]Efendiev, Y., Hou, T., and Luo, W.. Preconditioning Markov chain Monte Carlo simulations using coarse-scale models. SIAMJ. Sci. Comput., 28(2):776–803, 2006.CrossRefGoogle Scholar
[16]Efendiev, Y. and Hou, T.Y.. Multiscale Finite Element Methods. Springer, New York, 2009. Theory and applications.Google Scholar
[17]Ghanem, R.G. and Spanos, P.D.. Stochastic Finite Elements: A Spectral Approach. Springer-Verlag, New York, 1991.CrossRefGoogle Scholar
[18]Giles, M.. Improved multilevel Monte Carlo convergence using the Milstein scheme. In Monte Carlo and quasi-Monte Carlo methods 2006, pages 343–358. Springer, Berlin, 2008.CrossRefGoogle Scholar
[19]Giles, M.B.. Multilevel Monte Carlo path simulation. Oper. Res., 56(3):607–617, 2008.CrossRefGoogle Scholar
[20]Ginting, V., Pereira, F., Presho, M., and Wo, S.. Application of the two-stage Markov chain Monte Carlo method for characterization of fractured reservoirs using a surrogate flow model. Comput. Geosci., 15(4):691–707, 2011.CrossRefGoogle Scholar
[21]Heinrich, S.. Multilevel Monte Carlo methods. In Margenov, S., Waśniewski, J., and Yalamov, P., editors, LectureNotes in Computer Science, volume 2179, pages 58–67. Springer-Verlag, Berlin, 2001.Google Scholar
[22]Hou, T.Y. and Wu, X.-H.. A multiscale finite element method for elliptic problems in composite materials and porous media. J. Comput. Phys., 134(1):169–189, 1997.CrossRefGoogle Scholar
[23]Hughes, T.J.R., Feijóo, G.R., Mazzei, L., and Quincy, J.-B.. The variational multiscale method—a paradigm for computational mechanics. Comput. Methods Appl.Mech. Engrg., 166(1-2):3–24, 1998.CrossRefGoogle Scholar
[24]Ito, K. and Jin, B.. Inverse Problems: Tikhonov Theory and Algorithms. World Scientific, Singapore, 2014.Google Scholar
[25]Jenny, P., Lee, S.H., and Tchelepi, H.A.. Multi-scale finite-volumemethod for elliptic problems in subsurface flow simulation. J. Comput. Phys., 187(1):47–67, 2003.CrossRefGoogle Scholar
[26]Ketelsen, C., Scheichl, R., and Teckentrup, A.. A hierarchical multilevel markov chain monte carlo algorithm with applications to uncertainty quantification in subsurface flow. Submitted, arXiv:1303.7343.Google Scholar
[27]Loève, M.. Probability theory. II. Springer-Verlag, New York, fourth edition, 1978.CrossRefGoogle Scholar
[28]Nguyen, N.C.. A multiscale reduced-basismethod for parametrized elliptic partial differential equations with multiple scales. J. Comput. Phys., 227(23):9807–9822, 2008.CrossRefGoogle Scholar
[29]Robert, C.P. and Casella, G.. Monte Carlo Statistical Methods. Springer-Verlag, New York, second edition, 2004.CrossRefGoogle Scholar
[30]Rozza, G., Huynh, D.B.P., and Patera, A.T.. Reduced basis approximation and a posteriori error estimation for affinely parametrized elliptic coercive partial differential equations: application to transport and continuum mechanics. Arch. Comput. Methods Eng., 15(3):229–275, 2008.CrossRefGoogle Scholar
[31]Speight, A.. A multilevel approach to control variates. J. Comput. Finance, 12(4):3–27, 2009.CrossRefGoogle Scholar
[32]Xiu, D. and Karniadakis, G.E.. TheWiener-Askey polynomial chaos for stochastic differential equations. SIAM J. Sci. Comput., 24(2):619–644, 2002.CrossRefGoogle Scholar
[33]Xiu, D. and Karniadakis, G.E.. Modeling uncertainty in flow simulations via generalized polynomial chaos. J. Comput. Phys., 187(1):137–167, 2003.CrossRefGoogle Scholar