Research reports

Multilevel Monte Carlo finite volume methods for random conservation laws with discontinuous flux

by J. Badwaik and C. Klingenberg and N. H. Risebro and A. M. Ruf

(Report number 2020-55)

Abstract
We consider conservation laws with discontinuous flux where the initial datum, the flux function, and the discontinuous spatial dependency coefficient are subject to randomness. We establish a notion of random adapted entropy solutions to these equations and prove well-posedness provided that the spatial dependency coefficient is piecewise constant with finitely many discontinuities. In particular, the setting under consideration allows the flux to change across finitely many points in space whose positions are uncertain. We propose a single- and multilevel Monte Carlo method based on a finite volume approximation for each sample. Our analysis includes convergence rate estimates of the resulting Monte Carlo and multilevel Monte Carlo finite volume methods as well as error versus work rates showing that the multilevel variant outperforms the single-level method in terms of efficiency. We present numerical experiments motivated by two-phase reservoir simulations for reservoirs with varying geological properties.

Keywords: uncertainty quantification, conservation laws, discontinuous flux, numerical methods

BibTeX
@Techreport{BKRR20_928,
  author = {J. Badwaik and C. Klingenberg and N. H. Risebro and A. M. Ruf},
  title = {Multilevel Monte Carlo finite volume methods for random conservation laws with discontinuous flux},
  institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
  number = {2020-55},
  address = {Switzerland},
  url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2020/2020-55.pdf },
  year = {2020}
}

Disclaimer
© Copyright for documents on this server remains with the authors. Copies of these documents made by electronic or mechanical means including information storage and retrieval systems, may only be employed for personal use. The administrators respectfully request that authors inform them when any paper is published to avoid copyright infringement. Note that unauthorised copying of copyright material is illegal and may lead to prosecution. Neither the administrators nor the Seminar for Applied Mathematics (SAM) accept any liability in this respect. The most recent version of a SAM report may differ in formatting and style from published journal version. Do reference the published version if possible (see SAM Publications).

JavaScript has been disabled in your browser