Multi-level Monte Carlo finite volume methods for nonlinear systems of conservation laws in multi-dimensions

JOURNAL OF COMPUTATIONAL PHYSICS

卷:231

Article

Mishra, S.¹  Schwab, Ch.¹  Sukys, J.¹ 

[1] ETH, Seminar Appl Math, Zurich, Switzerland

关键词: Conservation laws;    Euler;    MHD;    Uncertainty quantification;    Multi-level Monte Carlo;    Parallelization;   

DOI  :  10.1016/j.jcp.2012.01.011

来源: Elsevier

【 摘 要 】

We extend the multi-level Monte Carlo (MLMC) in order to quantify uncertainty in the solutions of multi-dimensional hyperbolic systems of conservation laws with uncertain initial data. The algorithm is presented and several issues arising in the massively parallel numerical implementation are addressed. In particular, we present a novel load balancing procedure that ensures scalability of the MLMC algorithm on massively parallel hardware. A new code is described and applied to simulate uncertain solutions of the Euler equations and ideal magnetohydrodynamics (MHD) equations. Numerical experiments showing the robustness, efficiency and scalability of the proposed algorithm are presented. (C) 2012 Elsevier Inc. All rights reserved.

【 授权许可】

Free