8,398 research outputs found
Numerical approximation of statistical solutions of scalar conservation laws
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 multi-level 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 multi-level Monte Carlo
method over the standard Monte Carlo method. Numerical experiments illustrating
the ability of both methods to accurately compute multi-point statistical
quantities of interest are also presented
Numerical approximation of statistical solutions of scalar conservation laws
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 multi-level 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 multi-level Monte Carlo
method over the standard Monte Carlo method. Numerical experiments illustrating
the ability of both methods to accurately compute multi-point statistical
quantities of interest are also presented
Multilevel Monte Carlo for Random Degenerate Scalar Convection Diffusion Equation
We consider the numerical solution of scalar, nonlinear degenerate
convection-diffusion problems with random diffusion coefficient and with random
flux functions. Building on recent results on the existence, uniqueness and
continuous dependence of weak solutions on data in the deterministic case, we
develop a definition of random entropy solution. We establish existence,
uniqueness, measurability and integrability results for these random entropy
solutions, generalizing \cite{Mishr478,MishSch10a} to possibly degenerate
hyperbolic-parabolic problems with random data. We next address the numerical
approximation of random entropy solutions, specifically the approximation of
the deterministic first and second order statistics. To this end, we consider
explicit and implicit time discretization and Finite Difference methods in
space, and single as well as Multi-Level Monte-Carlo methods to sample the
statistics. We establish convergence rate estimates with respect to the
discretization parameters, as well as with respect to the overall work,
indicating substantial gains in efficiency are afforded under realistic
regularity assumptions by the use of the Multi-Level Monte-Carlo method.
Numerical experiments are presented which confirm the theoretical convergence
estimates.Comment: 24 Page
Multi-scale control variate methods for uncertainty quantification in kinetic equations
Kinetic equations play a major rule in modeling large systems of interacting
particles. Uncertainties may be due to various reasons, like lack of knowledge
on the microscopic interaction details or incomplete informations at the
boundaries. These uncertainties, however, contribute to the curse of
dimensionality and the development of efficient numerical methods is a
challenge. In this paper we consider the construction of novel multi-scale
methods for such problems which, thanks to a control variate approach, are
capable to reduce the variance of standard Monte Carlo techniques
A posteriori error analysis and adaptive non-intrusive numerical schemes for systems of random conservation laws
In this article we consider one-dimensional random systems of hyperbolic
conservation laws. We first establish existence and uniqueness of random
entropy admissible solutions for initial value problems of conservation laws
which involve random initial data and random flux functions. Based on these
results we present an a posteriori error analysis for a numerical approximation
of the random entropy admissible solution. For the stochastic discretization,
we consider a non-intrusive approach, the Stochastic Collocation method. The
spatio-temporal discretization relies on the Runge--Kutta Discontinuous
Galerkin method. We derive the a posteriori estimator using continuous
reconstructions of the discrete solution. Combined with the relative entropy
stability framework this yields computable error bounds for the entire
space-stochastic discretization error. The estimator admits a splitting into a
stochastic and a deterministic (space-time) part, allowing for a novel
residual-based space-stochastic adaptive mesh refinement algorithm. We conclude
with various numerical examples investigating the scaling properties of the
residuals and illustrating the efficiency of the proposed adaptive algorithm
- …