Mixed GMsFEM for second order elliptic problem in perforated domains

Abstract

We consider a class of second order elliptic problems in perforated domains with homogeneous Neumann boundary condition. It is well-known that numerically solving these problems require a very fine computational mesh, and some model reduction techniques are therefore necessary. We will develop a new model reduction technique based on the generalized multiscale finite element method (GMsFEM). The GMsFEM has been applied successfully to second order elliptic problems in perforated domains with Dirichlet boundary conditions Chung et al. (2015). However, due to the use of multiscale partition of unity functions, the same method cannot be applied to the case with Neumann boundary conditions. The aim of this paper is to develop a new mixed GMsFEM, based on a piecewise constant approximation for pressure and a multiscale approximation for velocity, giving a mass conservative method. The method can handle the Neumann boundary condition naturally. The multiscale basis functions for velocity are constructed by some carefully chosen local snapshot spaces and local spectral decompositions. The spectral convergence of the method is analyzed. Moreover, by using some local error indicators, the basis functions can be added locally and adaptively. We also consider an online procedure for the construction of new basis functions in the online stage in order to capture the distant effects. We will present some numerical examples to show the performance of the method.