1,926 research outputs found
Strict bounding of quantities of interest in computations based on domain decomposition
This paper deals with bounding the error on the estimation of quantities of
interest obtained by finite element and domain decomposition methods. The
proposed bounds are written in order to separate the two errors involved in the
resolution of reference and adjoint problems : on the one hand the
discretization error due to the finite element method and on the other hand the
algebraic error due to the use of the iterative solver. Beside practical
considerations on the parallel computation of the bounds, it is shown that the
interface conformity can be slightly relaxed so that local enrichment or
refinement are possible in the subdomains bearing singularities or quantities
of interest which simplifies the improvement of the estimation. Academic
assessments are given on 2D static linear mechanic problems.Comment: Computer Methods in Applied Mechanics and Engineering, Elsevier,
2015, online previe
Polyharmonic homogenization, rough polyharmonic splines and sparse super-localization
We introduce a new variational method for the numerical homogenization of
divergence form elliptic, parabolic and hyperbolic equations with arbitrary
rough () coefficients. Our method does not rely on concepts of
ergodicity or scale-separation but on compactness properties of the solution
space and a new variational approach to homogenization. The approximation space
is generated by an interpolation basis (over scattered points forming a mesh of
resolution ) minimizing the norm of the source terms; its
(pre-)computation involves minimizing quadratic (cell)
problems on (super-)localized sub-domains of size .
The resulting localized linear systems remain sparse and banded. The resulting
interpolation basis functions are biharmonic for , and polyharmonic
for , for the operator -\diiv(a\nabla \cdot) and can be seen as a
generalization of polyharmonic splines to differential operators with arbitrary
rough coefficients. The accuracy of the method ( in energy norm
and independent from aspect ratios of the mesh formed by the scattered points)
is established via the introduction of a new class of higher-order Poincar\'{e}
inequalities. The method bypasses (pre-)computations on the full domain and
naturally generalizes to time dependent problems, it also provides a natural
solution to the inverse problem of recovering the solution of a divergence form
elliptic equation from a finite number of point measurements.Comment: ESAIM: Mathematical Modelling and Numerical Analysis. Special issue
(2013
A multiscale method for heterogeneous bulk-surface coupling
In this paper, we construct and analyze a multiscale (finite element) method
for parabolic problems with heterogeneous dynamic boundary conditions. As
origin, we consider a reformulation of the system in order to decouple the
discretization of bulk and surface dynamics. This allows us to combine
multiscale methods on the boundary with standard Lagrangian schemes in the
interior. We prove convergence and quantify explicit rates for low-regularity
solutions, independent of the oscillatory behavior of the heterogeneities. As a
result, coarse discretization parameters, which do not resolve the fine scales,
can be considered. The theoretical findings are justified by a number of
numerical experiments including dynamic boundary conditions with random
diffusion coefficients
Multilevel Markov Chain Monte Carlo Method for High-Contrast Single-Phase Flow Problems
In this paper we propose a general framework for the uncertainty
quantification of quantities of interest for high-contrast single-phase flow
problems. It is based on the generalized multiscale finite element method
(GMsFEM) and multilevel Monte Carlo (MLMC) methods. The former provides a
hierarchy of approximations of different resolution, whereas the latter gives
an efficient way to estimate quantities of interest using samples on different
levels. The number of basis functions in the online GMsFEM stage can be varied
to determine the solution resolution and the computational cost, and to
efficiently generate samples at different levels. In particular, it is cheap to
generate samples on coarse grids but with low resolution, and it is expensive
to generate samples on fine grids with high accuracy. By suitably choosing the
number of samples at different levels, one can leverage the expensive
computation in larger fine-grid spaces toward smaller coarse-grid spaces, while
retaining the accuracy of the final Monte Carlo estimate. Further, we describe
a multilevel Markov chain Monte Carlo method, which sequentially screens the
proposal with different levels of approximations and reduces the number of
evaluations required on fine grids, while combining the samples at different
levels to arrive at an accurate estimate. The framework seamlessly integrates
the multiscale features of the GMsFEM with the multilevel feature of the MLMC
methods following the work in \cite{ketelson2013}, and our numerical
experiments illustrate its efficiency and accuracy in comparison with standard
Monte Carlo estimates.Comment: 29 pages, 6 figure
- …