56 research outputs found
On quantifying uncertainties for the linearized BGK kinetic equation
We consider the linearized BGK equation and want to quantify uncertainties in
the case of modelling errors. More specifically, we want to quantify the error
produced if the pre-determined equilibrium function is chosen inaccurately. In
this paper we consider perturbations in the velocity and in the temperature of
the equilibrium function and consider how much the error is amplified in the
solution
Computation of Electromagnetic Fields Scattered From Objects With Uncertain Shapes Using Multilevel Monte Carlo Method
Computational tools for characterizing electromagnetic scattering from
objects with uncertain shapes are needed in various applications ranging from
remote sensing at microwave frequencies to Raman spectroscopy at optical
frequencies. Often, such computational tools use the Monte Carlo (MC) method to
sample a parametric space describing geometric uncertainties. For each sample,
which corresponds to a realization of the geometry, a deterministic
electromagnetic solver computes the scattered fields. However, for an accurate
statistical characterization the number of MC samples has to be large. In this
work, to address this challenge, the continuation multilevel Monte Carlo
(CMLMC) method is used together with a surface integral equation solver. The
CMLMC method optimally balances statistical errors due to sampling of the
parametric space, and numerical errors due to the discretization of the
geometry using a hierarchy of discretizations, from coarse to fine. The number
of realizations of finer discretizations can be kept low, with most samples
computed on coarser discretizations to minimize computational cost.
Consequently, the total execution time is significantly reduced, in comparison
to the standard MC scheme.Comment: 25 pages, 10 Figure
Multilevel Double Loop Monte Carlo and Stochastic Collocation Methods with Importance Sampling for Bayesian Optimal Experimental Design
An optimal experimental set-up maximizes the value of data for statistical
inferences and predictions. The efficiency of strategies for finding optimal
experimental set-ups is particularly important for experiments that are
time-consuming or expensive to perform. For instance, in the situation when the
experiments are modeled by Partial Differential Equations (PDEs), multilevel
methods have been proven to dramatically reduce the computational complexity of
their single-level counterparts when estimating expected values. For a setting
where PDEs can model experiments, we propose two multilevel methods for
estimating a popular design criterion known as the expected information gain in
simulation-based Bayesian optimal experimental design. The expected information
gain criterion is of a nested expectation form, and only a handful of
multilevel methods have been proposed for problems of such form. We propose a
Multilevel Double Loop Monte Carlo (MLDLMC), which is a multilevel strategy
with Double Loop Monte Carlo (DLMC), and a Multilevel Double Loop Stochastic
Collocation (MLDLSC), which performs a high-dimensional integration by
deterministic quadrature on sparse grids. For both methods, the Laplace
approximation is used for importance sampling that significantly reduces the
computational work of estimating inner expectations. The optimal values of the
method parameters are determined by minimizing the average computational work,
subject to satisfying the desired error tolerance. The computational
efficiencies of the methods are demonstrated by estimating the expected
information gain for Bayesian inference of the fiber orientation in composite
laminate materials from an electrical impedance tomography experiment. MLDLSC
performs better than MLDLMC when the regularity of the quantity of interest,
with respect to the additive noise and the unknown parameters, can be
exploited
Fast r-adaptivity for multiple queries of heterogeneous stochastic material fields
We present an r-adaptivity approach for boundary value problems with randomly fluctuating material parameters solved through the Monte Carlo or stochastic collocation methods. This approach tailors a specific mesh for each sample of the problem. It only requires the computation of the solution of a single deterministic problem with the same geometry and the average parameter, whose numerical cost becomes marginal for large number of samples. Starting from the mesh used to solve that deterministic problem, the nodes are moved depending on the particular sample of mechanical parameter field. The reduction in the error is small for each sample but sums up to reduce the overall bias on the statistics estimated through the Monte Carlo scheme. Several numerical examples in 2D are presented.Peer ReviewedPostprint (author's final draft
Multilevel Monte Carlo methods
The author's presentation of multilevel Monte Carlo path simulation at the
MCQMC 2006 conference stimulated a lot of research into multilevel Monte Carlo
methods. This paper reviews the progress since then, emphasising the
simplicity, flexibility and generality of the multilevel Monte Carlo approach.
It also offers a few original ideas and suggests areas for future research
- …