Search
Results
-
Multi-level Monte Carlo finite volume methods for nonlinear systems of conservation laws in multi-dimensions
(2011)SAM Research ReportWe extend the Multi-Level Monte Carlo (MLMC) algorithm of [19] 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 ...Report -
Numerical solution of scalar conservation laws with random flux functions
(2012)SAM Research ReportWe consider scalar hyperbolic conservation laws in several space dimensions, with a class of random (and parametric) flux functions. We propose a Karhunen-Loève expansion on the state space of the random flux. For random flux functions which are Lipschitz continuous with respect to the state variable, we prove the existence of a unique random entropy solution. Using a Karhunen-Loève spectral decomposition of the random flux into principal ...Report -
Multi-level Monte Carlo finite difference and finite volume methods for stochastic linear hyperbolic systems
(2012)SAM Research ReportWe consider stochastic multi-dimensional linear hyperbolic systems of conservation laws. We prove existence and uniqueness of a random weak solution, provide estimates for the regularity of the solution in terms of regularities of input data, and show existence of statistical moments. Bounds for mean square error vs. expected work are proved for the Multi-Level Monte Carlo Finite Volume algorithm which is used to approximate the moments ...Report -
Sparse Tensor Multi-Level Monte Carlo Finite Volume Methods for hyperbolic conservation laws with random intitial data
(2010)SAM Research ReportWe consider scalar hyperbolic conservation laws in several $(d \ge 1)$ spatial dimensions with stochastic initial data. We prove existence and uniqueness of a random-entropy solution and show existence of statistical moments of any order k of this random entropy solution. We present a class of numerical schemes of multi-level Monte Carlo Finite Volume (MLMC-FVM) type for the approximation of random entropy solutions as well as of their ...Report -
Multi-level Monte Carlo finite volume methods for uncertainty quantification in nonlinear systems of balance laws
(2012)Research ReportReport -
Multi-Level Monte Carlo Finite Volume methods for uncertainty quantification of acoustic wave propagation in random heterogeneous layered medium
(2014)Research ReportWe consider the very challenging problem of efficient uncertainty quantification for acoustic wave propagation in a highly heterogeneous, possibly layered, random medium, characterized by possibly anisotropic, piecewise log-exponentially distributed Gaussian random fields. A multi-level Monte Carlo finite volume method is proposed, along with a novel, bias-free upscaling technique that allows to represent the input random fields, generated ...Report -
Multi-level Monte Carlo finite volume methods for shallow water equations with uncertain topography in multi-dimensions
(2011)SAM Research ReportThe initial data and bottom topography, used as inputs in shallow water models, are prone to uncertainty due to measurement errors. We model this uncertainty statistically in terms of random shallow water equations. We extend the Multi-Level Monte Carlo (MLMC) algorithm to numerically approximate the random shallow water equations efficiently. The MLMC algorithm is suitably modified to deal with uncertain (and possibly uncorrelated) data ...Report -
Static load balancing for multi-level Monte Carlo finite volume solvers
(2011)SAM Research ReportThe Multi-Level Monte Carlo finite volumes (MLMC-FVM) algorithm was shown to be a robust and fast solver for uncertainty quantification in the solutions of multi- dimensional systems of stochastic conservation laws. A novel load balancing procedure is used to ensure scalability of the MLMC algorithm on massively parallel hardware. We describe this procedure together with other arising challenges in great detail. Finally, numerical experiments ...Report