Shubin Fu

Yingjie Zhou, Shubin Fu, Eric Tsz Shun Chung

We present a two-grid algebraic multiscale preconditioner for large sparse symmetric positive definite systems arising from elliptic problems with highly heterogeneous coefficients. The coarse space is constructed directly from the system matrix by graph partitioning and local generalized eigenvalue solvers, yielding basis functions that capture the low-energy modes responsible for slow convergence. The method requires no geometric information, making it suitable for unstructured and matrix-only settings, and its construction is naturally parallelizable. Numerical results for heterogeneous Darcy flow problems show robustness with respect to coefficient contrast and problem size, better performance than standard algebraic multigrid on challenging large-scale cases, and good parallel scalability.

Eric T. Chung, Yalchin Efendiev, Shubin Fu

In this paper, we discuss the application of Generalized Multiscale Finite Element Method (GMsFEM) to elasticity equation in heterogeneous media. Our applications are motivated by elastic wave propagation in subsurface where the subsurface properties can be highly heterogeneous and have high contrast. We present the construction of main ingredients for GMsFEM such as the snapshot space and offline spaces. The latter is constructed using local spectral decomposition in the snapshot space. The spectral decomposition is based on the analysis which is provided in the paper. We consider both continuous Galerkin and discontinuous Galerkin coupling of basis functions. Both approaches have their cons and pros. Continuous Galerkin methods allow avoiding penalty parameters though they involve partition of unity functions which can alter the properties of multiscale basis functions. On the other hand, discontinuous Galerkin techniques allow gluing multiscale basis functions without any modifications. Because basis functions are constructed independently from each other, this approach provides an advantage. We discuss the use of oversampling techniques that use snapshots in larger regions to construct the offline space. We provide numerical results to show that one can accurately approximate the solution using reduced number of degrees of freedom.

Shubin Fu, Eric Chung, Guanglian Li

In this paper, we proposed two new types of edge multiscale methods motivated by \cite{GL18} to solve Partial Differential Equations (PDEs) with high-contrast heterogeneous coefficients: Edge spectral multiscale Finte Element method (ESMsFEM) and Wavelet-based edge multiscale Finite Element method (WEMsFEM). Their convergence rates for elliptic problems with high-contrast heterogeneous coefficients are demonstrated in terms of the coarse mesh size $H$, the number of spectral basis functions and the level of the wavelet space $\ell$, which are verified by extensive numerical tests.

Shubin Fu, Robert Altmann, Eric T. Chung et al.

In this work, we employ the Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM) to solve the problem of linear heterogeneous poroelasticity with coefficients of high contrast. The proposed method makes use of the idea of energy minimization with suitable constraints in order to generate efficient basis functions for the displacement and the pressure. These basis functions are constructed by solving a class of local auxiliary optimization problems based on eigenfunctions containing local information on the heterogeneity. Techniques of oversampling are adapted to enhance the computational performance. Convergence of first order is shown and illustrated by a number of numerical tests.

Yanfang Yang, Shubin Fu, Eric T Chung

In this paper, we consider an online basis enrichment mixed generalized multiscale method with oversampling, for solving flow problems in highly heterogeneous porous media. This is an exten- sion of the online mixed generalized multiscale method [6]. The multiscale online basis functions are computed by solving a Neumann problem in an over-sampled domain, instead of a standard neighborhood of a coarse face. We are motivated by the restricted domain decomposition method. Extensive numerical experiments are presented to demonstrate the performance of our methods for both steady-state flow, and two-phase flow and transport problems.

Yanfang Yang, Shubin Fu, Eric T. Chung

In this paper, we consider flow simulation in highly heterogeneous media that has many practical applications in industry. To enhance mass conservation, we write the elliptic problem in a mixed formulation and introduce a robust two-grid preconditioner to seek the solution. We first need to transform the indefinite saddle problem to a positive definite problem by preprocessing steps. The preconditioner consists of a local smoother and a coarse preconditioner. For the coarse preconditioner, we design an adaptive spectral coarse space motivated by the GMsFEM (Generalized Multiscale Finite Element Method). We test our preconditioner for both Darcy flow and two phase flow and transport simulation in highly heterogeneous porous media. Numerical results show that the proposed preconditioner is highly robust and efficient.

Shubin Fu, Eric Chung

In this paper, we propose a local-global multiscale mortar mixed finite element method (MMMFEM) for multiphase transport in heterogeneous media. We consider the two-phase flow system, the pressure equation is solved via the multiscale mortar mixed finite element method, a mass conservation velocity field can be obtained, then we use explicit finite volume method to solve the saturation equation. We use polynomials and multiscale basis to form the coarse mortar space. The multiscale basis is the restriction of global pressure field obtained at previous time step on the coarse interface. We solve the pressure equation on the fine grid to initialize the simulation. Numerical experiments on some benchmark 2D and 3D heterogeneous models are provided to validate the performance of our method.

Yanfang Yang, Eric T. Chung, Shubin Fu

In this paper, we develop an online basis enrichment method with the mortar mixed finite element method, using the oversampling technique, to solve for flow problems in highly heterogeneous media. We first compute a coarse grid solution with a certain number of offline basis functions per edge, which are chosen as standard polynomials basis functions. We then iteratively enrich the multiscale solution space with online multiscale basis functions computed by using residuals. The iterative solution converges to the fine scale solution rapidly. We also propose an oversampling online method to achieve faster convergence speed. The oversampling refers to using larger local regions in computing the online multiscale basis functions. We present extensive numerical experiments(including both 2D and 3D) to demonstrate the performance of our methods for both steady state flow, and two-phase flow and transport problems. In particular, for the time dependent two-phase flow and transport problems, we apply the online method to the initial model, without updating basis along the time evolution. Our numerical results demonstrate that by using a few number of online basis functions, one can achieve a fast convergence.

Peiqi Li, Jie Chen, Shubin Fu

In this paper, we study the efficient numerical solution of Darcy equations in strongly heterogeneous media with high-contrast permeability and propose a hybrid framework that combines learning with multiscale numerical methods. The learning component is used for the prediction of multiscale basis functions in the mixed generalized multiscale finite element method (mixed GMsFEM), with the goal of reducing the repeated local computations required in the offline stage. Once these basis functions are predicted, the global system is assembled and the pressure field is computed by a two-grid preconditioned solver. The resulting method accelerates the costly local basis-construction stage while retaining the multiscale discretization and preconditioned iterative structure of the underlying solver. Numerical experiments on two-dimensional heterogeneous Darcy problems show that the proposed framework yields more accurate final pressure reconstruction than several representative learning-based methods and remains stable under strong heterogeneity and high-contrast coefficients. In comparison with the traditional mixed GMsFEM, its main advantage lies in the efficiency of the basis-generation stage, while the quality of the global solve is still ensured by the two-grid preconditioner. These results indicate that accelerating multiscale basis construction through learning, while preserving a mature numerical solver for the global problem, provides a viable approach for high-resolution Darcy-type simulations.

Shubin Fu, Eric T. Chung

In this paper, we consider the offline and online Constraint Energy Minimizing Generalized Mul- tiscale Finite Element Method (CEM-GMsFEM) for high-contrast linear elasticity problem. Offline basis construction starts with an auxiliary multiscale space by solving local spectral problems. We select eigenfunctions that correspond to a few small eigenvalues to form the auxiliary space. Using the auxiliary space, we solve a constraint energy minimization problem to construct offline multiscale spaces. The minimization problem is defined in the oversampling domain, which is larger than the target coarse block. To get a good approximation space, the oversampling domain should be large enough. We also propose a relaxed minimization problem to construct multiscale basis functions, which will yield more accurate and robust solution. To take into account the influence of input pa- rameters, such as source terms, we propose the construction of online multiscale basis and an adaptive enrichment algorithm. We provide extensive numerical experiments on 2D and 3D models to show the performance of the proposed method.

Eric T. Chung, Shubin Fu, Yanfang Yang

Mortar methods are widely used techniques for discretizations of partial differential equations and preconditioners for the algebraic systems resulting from the discretizations. For problems with high contrast and multiple scales, the standard mortar spaces are not robust, and some enrichments are necessary in order to obtain an efficient and robust mortar space. In this paper, we consider a class of flow problems in high contrast heterogeneous media, and develop a systematic approach to obtain an enriched multiscale mortar space. Our approach is based on the constructions of local multiscale basis functions. The multiscale basis functions are constructed from local problems by following the framework of the Generalized Multiscale Finite Element Method (GMsFEM). In particular, we first create a local snapshot space. Then we select the dominated modes within the snapshot space using an appropriate Proper Orthogonal Decomposition (POD) technique. These multiscale basis functions show better accuracy than polynomial basis for multiscale problems. Using the proposed multiscale mortar space, we will construct a multiscale finite element method to solve the flow problem on a coarse grid and a preconditioning technique for the fine scale discretization of the flow problem. In particular, we develop a multiscale mortar mixed finite element method using the mortar space. In addition, we will design a two-level additive preconditioner and a two-level hybrid preconditioner based on the proposed mortar space for the iterative method applied to the fine scale discretization of the flow problem. We present several numerical examples to demonstrate the efficiency and robustness of our proposed mortar space with respect to both the coarse multiscale solver and the preconditioners.

Kai Gao, Shubin Fu, Richard L. Gibson et al.

It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.