A Constraint energy minimizing generalized multiscale finite element method for parabolic equations

Mengnan Li,Eric Chung,Lijian Jiang

DOI: https://doi.org/10.48550/arXiv.1806.04816

2018-06-13

Abstract:In this paper, we present a Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM) for parabolic equations with multiscale coefficients, arising from applications in porous media. We will present the construction of CEM-GMsFEM and rigorously analyze its convergence for the parabolic equations. The convergence rate is characterized by the coarse grid size and the eigenvalue decay of local spectral problems, but is independent of the scale length and contrast of the media. The analysis shows that the method has a first order convergence rate with respect to coarse grid size in the energy norm and second order convergence rate with respect to coarse grid size in $L^2$ norm under some appropriate assumptions. For the temporal discretization, finite difference techniques are used and the convergence analysis of full discrete scheme is given. Moreover, a posteriori error estimator is derived and analyzed. A few numerical results for porous media applications are presented to confirm the theoretical findings and demonstrate the performance of the approach.

Numerical Analysis

What problem does this paper attempt to address?

This paper attempts to solve the numerical simulation problem of parabolic equations with multi - scale coefficients in the application of porous media. Specifically, the author proposes a Constrained Energy Minimization - based Generalized Multiscale Finite Element Method (CEM - GMsFEM) to address the challenges encountered by traditional numerical methods in dealing with high - contrast porous media. Such problems are very important in practical applications, such as subsurface fluid flow modeling.

### Problem Background

In porous media, the multi - scale characteristics of geological features pose significant challenges to subsurface fluid modeling. Usually, the geological features that capture these effects are developed on a very fine scale, which makes it impractical to directly conduct fluid simulations. Therefore, simplified models are required to calculate the solutions of fluid problems.

### Existing Methods

There are currently some methods for solving multi - scale problems:
- **Upscaling Method**: Based on the homogenization theory, form coarse - scale equations through pre - calculated effective coefficients.
- **Multiscale Finite Element Method (MsFEM)**: Obtain coarse - scale equations through variational formulations and multiscale basis functions containing fine - scale information.
- **Variational Multiscale Method** and **Multiscale Finite Volume Method**: Use different frameworks to solve multi - scale problems.

### Solution Proposed in the Paper

To improve existing methods, this paper proposes the CEM - GMsFEM method, which is constructed through the following steps:
1. **Construction of Auxiliary Basis Functions**: Solve the local spectral problem in each coarse - grid block to obtain auxiliary basis functions.
2. **Constrained Energy Minimization Problem**: Solve the constrained energy minimization problem in the oversampling domain to obtain the required multiscale basis functions.

### Key Contributions

- **Convergence Analysis**: It is proved that the convergence rate of this method is first - order in the energy norm and second - order in the \(L^2\) norm, and is independent of the ratio of the medium.
- **A Posteriori Error Estimation**: Derive and analyze the a posteriori error estimator, providing a basis for the adaptive enrichment scheme.
- **Numerical Results**: Verify the theoretical results through multiple numerical experiments, demonstrating the effectiveness and robustness of the method.

### Summary

The main objective of this paper is to solve parabolic equations with high - contrast coefficients through the CEM - GMsFEM method, and conduct detailed convergence analysis and numerical verification. This method not only improves computational efficiency but also can maintain high precision in complex multi - scale environments.