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

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

Research output: Contribution to journalArticlepeer-review

18 Scopus citations


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.
Original languageEnglish (US)
Pages (from-to)109133
JournalJournal of Computational Physics
StatePublished - Nov 25 2019

Bibliographical note

KAUST Repository Item: Exported on 2020-10-01
Acknowledged KAUST grant number(s): BAS/1/1351-01
Acknowledgements: 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.


Dive into the research topics of 'Generalized multiscale approximation of mixed finite elements with velocity elimination for subsurface flow'. Together they form a unique fingerprint.

Cite this