Generalized multiscale approximation of mixed finite elements with velocity elimination for subsurface flow

Embargo End Date


Chen, Jie
Chung, Eric T.
He, Zhengkang
Sun, Shuyu

KAUST Department
Computational Transport Phenomena Lab
Earth Science and Engineering Program
Physical Science and Engineering (PSE) Division

KAUST Grant Number

Online Publication Date

Print Publication Date


A frame work of the mixed generalized multiscale finite element method (GMsFEM) for solving Darcy's law in heterogeneous media is studied in this paper. Our approach approximates pressure in multiscale function space that is between fine-grid space and coarse-grid space and solves velocity directly in the fine-grid space. To construct multiscale basis functions for each coarse-grid element, three types of snapshot space are raised. The first one is taken as the fine-grid space for pressure and the other two cases need to solve a local problem on each coarse-grid element. We describe a spectral decomposition in the snapshot space motivated by the analysis to further reduce the dimension of the space that is used to approximate the pressure. Since the velocity is directly solved in the fine-grid space, in the linear system for the mixed finite elements, the velocity matrix can be approximated by a diagonal matrix without losing any accuracy. Thus it can be inverted easily. This reduces computational cost greatly and makes our scheme simple and easy for application. Comparing to our previous work of mixed generalized multiscale finite element method [E. T. Chung, Y. Efendiev, and C. S. Lee. Mixed generalized multiscale finite element methods and applications. Multiscale Modeling and Simulation, 13(1):338-366, 2015.], both the pressure and velocity space in this approach are bigger. As a consequence, this method enjoys better accuracy. While the computational cost does not increase because of the good property of velocity matrix. Moreover, the proposed method preserves the local mass conservation property that is important for subsurface problems. Numerical examples are presented to illustrate the good properties of the proposed approach. If offline spaces are appropriately selected, one can achieve good accuracy with only a few basis functions per coarse element according to the numerical results.

Chen, J., Chung, E. T., He, Z., & Sun, S. (2019). Generalized multiscale approximation of mixed finite elements with velocity elimination for subsurface flow. Journal of Computational Physics, 109133. doi:10.1016/

The work is supported by the National Natural Science Foundation of China (No.11401467 and No.51874262), Introduction plan of senior foreign experts (G20190027012), and XJTLU Key Programme Special Fund (KSF-P-02). The work is also supported in part by funding from King Abdullah University of Science and Technology (KAUST) through the grant BAS/1/1351-01. The research of Eric Chung is supported by Hong Kong RGC General Research Fund (project 14317516) and CUHK Direct Grant for Research 2017-18.

Elsevier BV

Journal of Computational Physics


Additional Links

Permanent link to this record