Abstract:In this paper, we systemically review and compare two mixed multiscale finite element methods (MMsFEM) for multiphase transport in highly heterogeneous media. In particular, we will consider the mixed multiscale finite element method using limited global information, simply denoted by MMsFEM, and the mixed generalized multiscale finite element method (MGMsFEM) with residual driven online multiscale basis functions. Both methods are under the framework of mixed multiscale finite element methods, where the pressure equation is solved in the coarse grid with carefully constructed multiscale basis functions for the velocity. The multiscale basis functions in both methods include local and global media information. In terms of MsFEM using limited global information, only one multiscale basis function is utilized in each local neighborhood while multiple basis are used in MGMsFEM. We will test and compare these two methods using the benchmark three-dimensional SPE10 model. A range of coarse grid sizes and different combinations of basis functions (offline and online) will be considered with CPU time reported for each case. In our numerical experiments, we observe good accuracy by the two above methods. Finally, we will discuss and compare the advantages and disadvantages of the two methods in terms of accuracy and computational costs.
What problem does this paper attempt to address?
This paper attempts to solve the problem of multiphase transport simulation in highly heterogeneous media. Specifically, the paper focuses on the two - phase flow model in porous media, where the pressure equation is coupled with specific transport equations. Due to the high heterogeneity and large contrast of permeability, this brings significant difficulties to traditional solution methods. Therefore, the focus of the paper is to compare the performance of two Mixed Multiscale Finite Element Methods (MMs - FEM) in dealing with such problems. These two methods are:
1. **Mixed Multiscale Finite Element Method (MMsFEM) using limited global information**: This method uses only one multiscale basis function in each local neighborhood.
2. **Mixed Generalized Multiscale Finite Element Method (MGMsFEM) with residual - driven online multiscale basis functions**: This method uses multiple basis functions in each local neighborhood.
### Main research objectives:
- **Accuracy**: Evaluate the accuracy of these two methods under different coarse grid sizes and different basis function combinations.
- **Computational cost**: Compare the computational costs of these two methods, including the time for constructing multiscale basis functions and the simulation time.
- **Flexibility**: Explore the flexibility of MGMsFEM in adjusting the dimension of the approximation space relative to MMsFEM.
### Research methods:
- **Benchmark test**: Use the three - dimensional SPE10 model for testing, which is a standard multiphase flow simulation benchmark model.
- **Numerical experiment**: Consider different coarse grid sizes and basis function combinations, and record the CPU time in each case.
- **Result analysis**: Evaluate their performance by comparing the average saturation error and the computation time of the two methods.
### Main findings:
- **Accuracy**: Both methods can provide high accuracy, but MGMsFEM performs better when using online basis functions, especially when there are oscillations in the fine - scale source terms.
- **Computational cost**: MMsFEM is more efficient in reducing errors, and at the same time, the system dimension is smaller and the computation time is shorter. Although MGMsFEM requires more computational resources in some cases, the accuracy can be significantly improved through online enrichment.
- **Flexibility**: MGMsFEM allows flexible adjustment of the dimension of the approximation space, while MMsFEM is relatively fixed.
### Conclusion:
Through systematic comparison and analysis, the paper shows the advantages and limitations of MMsFEM and MGMsFEM in dealing with multiphase transport problems in highly heterogeneous media. MMsFEM has an advantage in computational efficiency, while MGMsFEM has more advantages in flexibility and accuracy. These findings provide an important reference for choosing the appropriate multiscale method.