Original version
SIAM Journal on Numerical Analysis. 2018, 56 (5), 2989-3009, DOI: https://doi.org/10.1137/17M1154874
Abstract
We propose efficient numerical algorithms for approximating statistical solutions of scalar conservation laws. The proposed algorithms combine finite volume spatio-temporal approximations with Monte Carlo and multilevel Monte Carlo discretizations of the probability space. Both sets of methods are proved to converge to the entropy statistical solution. We also prove that there is a considerable gain in efficiency resulting from the multilevel Monte Carlo method over the standard Monte Carlo method. Numerical experiments illustrating the ability of both methods to accurately compute multipoint statistical quantities of interest are also presented.