62 research outputs found
An adaptive GMsFEM for high-contrast flow problems
In this paper, we derive an a-posteriori error indicator for the Generalized
Multiscale Finite Element Method (GMsFEM) framework. This error indicator is
further used to develop an adaptive enrichment algorithm for the linear
elliptic equation with multiscale high-contrast coefficients. The GMsFEM, which
has recently been introduced in [12], allows solving multiscale
parameter-dependent problems at a reduced computational cost by constructing a
reduced-order representation of the solution on a coarse grid. The main idea of
the method consists of (1) the construction of snapshot space, (2) the
construction of the offline space, and (3) the construction of the online space
(the latter for parameter-dependent problems). In [12], it was shown that the
GMsFEM provides a flexible tool to solve multiscale problems with a complex
input space by generating appropriate snapshot, offline, and online spaces. In
this paper, we study an adaptive enrichment procedure and derive an
a-posteriori error indicator which gives an estimate of the local error over
coarse grid regions. We consider two kinds of error indicators where one is
based on the -norm of the local residual and the other is based on the
weighted -norm of the local residual where the weight is related to the
coefficient of the elliptic equation. We show that the use of weighted
-norm residual gives a more robust error indicator which works well for
cases with high contrast media. The convergence analysis of the method is
given. In our analysis, we do not consider the error due to the fine-grid
discretization of local problems and only study the errors due to the
enrichment. Numerical results are presented that demonstrate the robustness of
the proposed error indicators.Comment: 26 page
Generalized Multiscale Finite Element Methods for problems in perforated heterogeneous domains
Complex processes in perforated domains occur in many real-world
applications. These problems are typically characterized by physical processes
in domains with multiple scales (see Figure 1 for the illustration of a
perforated domain). Moreover, these problems are intrinsically multiscale and
their discretizations can yield very large linear or nonlinear systems. In this
paper, we investigate multiscale approaches that attempt to solve such problems
on a coarse grid by constructing multiscale basis functions in each coarse
grid, where the coarse grid can contain many perforations. In particular, we
are interested in cases when there is no scale separation and the perforations
can have different sizes. In this regard, we mention some earlier pioneering
works [14, 18, 17], where the authors develop multiscale finite element
methods. In our paper, we follow Generalized Multiscale Finite Element Method
(GMsFEM) and develop a multiscale procedure where we identify multiscale basis
functions in each coarse block using snapshot space and local spectral
problems. We show that with a few basis functions in each coarse block, one can
accurately approximate the solution, where each coarse block can contain many
small inclusions. We apply our general concept to (1) Laplace equation in
perforated domain; (2) elasticity equation in perforated domain; and (3) Stokes
equations in perforated domain. Numerical results are presented for these
problems using two types of heterogeneous perforated domains. The analysis of
the proposed methods will be presented elsewhere
- β¦