ライブラリ登録: Guest
Begell Digital Portal Begellデジタルライブラリー 電子書籍 ジャーナル 参考文献と会報 リサーチ集
International Journal for Uncertainty Quantification
インパクトファクター: 4.911 5年インパクトファクター: 3.179 SJR: 1.008 SNIP: 0.983 CiteScore™: 5.2

ISSN 印刷: 2152-5080
ISSN オンライン: 2152-5099

Open Access

International Journal for Uncertainty Quantification

DOI: 10.1615/Int.J.UncertaintyQuantification.2016018677
pages 57-81

A MULTIGRID MULTILEVEL MONTE CARLO METHOD USING HIGH-ORDER FINITE-VOLUME SCHEME FOR LOGNORMAL DIFFUSION PROBLEMS

Prashant Kumar
CWI-Centrum Wiskunde & Informatica, Amsterdam, The Netherlands
Cornelis W. Oosterlee
CWI-Centrum Wiskunde & Informatica, Amsterdam, The Netherlands; Delft University of Technology, Delft Institute of Applied Mathematics, Delft, The Netherlands
Richard P. Dwight
Aerodynamics Group, Faculty of Aerospace, TU Delft, P.O. Box 5058, 2600GB Delft, The Netherlands

要約

The aim of this paper is to show that a high-order discretization can be used to improve the convergence of a multilevel Monte Carlo method for elliptic partial differential equations with lognormal random coefficients in combination with the multigrid solution method. To demonstrate this, we consider a fourth-order accurate finite-volume discretization. With the help of the Matérn family of covariance functions, we simulate the coefficient field with different degrees of smoothness. The idea behind using a fourth-order scheme is to capture the additional regularity in the solution introduced due to higher smoothness of the random field. Second-order schemes previously utilized for these types of problems are not able to fully exploit this additional regularity. We also propose a practical way of combining a full multigrid solver with the multilevel Monte Carlo estimator constructed on the same mesh hierarchy. Through this integration, one full multigrid solve at any level provides a valid sample for all the preceding Monte Carlo levels. The numerical results show that the fourth-order multilevel estimator consistently outperforms the second-order variant. In addition, we observe an asymptotic gain for the standard Monte Carlo estimator.


Articles with similar content:

A MULTI-INDEX MARKOV CHAIN MONTE CARLO METHOD
International Journal for Uncertainty Quantification, Vol.8, 2018, issue 1
Ajay Jasra, Yan Zhou, Kengo Kamatani, Kody J.H. Law
Equation-Free Accelerated Simulations of the Morphological Relaxation of Crystal Surfaces
International Journal for Multiscale Computational Engineering, Vol.8, 2010, issue 4
Steven J. Plimpton, Xiaowang Zhou, Gregory J. Wagner
ON A POLYNOMIAL CHAOS METHOD FOR DIFFERENTIAL EQUATIONS WITH SINGULAR SOURCES
International Journal for Uncertainty Quantification, Vol.1, 2011, issue 1
Yunfei Song, Jae-Hun Jung
A PHOTON MONTE CARLO SOLVER UTILIZING A LOW DISCREPANCY SEQUENCE FOR THERMAL RADIATION IN COMBUSTION SYSTEMS
ICHMT DIGITAL LIBRARY ONLINE, Vol.0, 2019, issue
Somesh P. Roy, Joseph A. Farmer
A GRADIENT-BASED SAMPLING APPROACH FOR DIMENSION REDUCTION OF PARTIAL DIFFERENTIAL EQUATIONS WITH STOCHASTIC COEFFICIENTS
International Journal for Uncertainty Quantification, Vol.5, 2015, issue 1
Miroslav Stoyanov, Clayton G. Webster