A Generalized Multiscale Finite Element Method for the Brinkman Equation

In this paper we consider the numerical upscaling of the Brinkman equation in the presence of high-contrast permeability fields. We develop and analyze a robust and efficient Generalized Multiscale Finite Element Method (GMsFEM) for the Brinkman model. In the fine grid, we use mixed finite element method with the velocity and pressure being continuous piecewise quadratic and piecewise constant finite element spaces, respectively. Using the GMsFEM framework we construct suitable coarse-scale spaces for the velocity and pressure that yield a robust mixed GMsFEM. We develop a novel approach to construct a coarse approximation for the velocity snapshot space and a robust small offline space for the velocity space. The stability of the mixed GMsFEM and a priori error estimates are derived. A variety of two-dimensional numerical examples are presented to illustrate the effectiveness of the algorithm.Comment: 22 page

Asymptotic expansions for high-contrast elliptic equations

In this paper, we present a high-order expansion for elliptic equations in high-contrast media. The background conductivity is taken to be one and we assume the medium contains high (or low) conductivity inclusions. We derive an asymptotic expansion with respect to the contrast and provide a procedure to compute the terms in the expansion. The computation of the expansion does not depend on the contrast which is important for simulations. The latter allows avoiding increased mesh resolution around high conductivity features. This work is partly motivated by our earlier work in \cite{ge09_1} where we design efficient numerical procedures for solving high-contrast problems. These multiscale approaches require local solutions and our proposed high-order expansion can be used to approximate these local solutions inexpensively. In the case of a large-number of inclusions, the proposed analysis can help to design localization techniques for computing the terms in the expansion. In the paper, we present a rigorous analysis of the proposed high-order expansion and estimate the remainder of it. We consider both high and low conductivity inclusions

Expanded mixed multiscale finite element methods and their applications for flows in porous media

We develop a family of expanded mixed Multiscale Finite Element Methods (MsFEMs) and their hybridizations for second-order elliptic equations. This formulation expands the standard mixed Multiscale Finite Element formulation in the sense that four unknowns (hybrid formulation) are solved simultaneously: pressure, gradient of pressure, velocity and Lagrange multipliers. We use multiscale basis functions for the both velocity and gradient of pressure. In the expanded mixed MsFEM framework, we consider both cases of separable-scale and non-separable spatial scales. We specifically analyze the methods in three categories: periodic separable scales, $G$- convergence separable scales, and continuum scales. When there is no scale separation, using some global information can improve accuracy for the expanded mixed MsFEMs. We present rigorous convergence analysis for expanded mixed MsFEMs. The analysis includes both conforming and nonconforming expanded mixed MsFEM. Numerical results are presented for various multiscale models and flows in porous media with shales to illustrate the efficiency of the expanded mixed MsFEMs.Comment: 33 page