Non-Conforming Multiscale Finite Element Method for Stokes Flows in Heterogeneous Media. Part I: Methodologies and Numerical Experiments

Bagus Putra Muljadi,Jacek Narski,Alexei Lozinski,Pierre Degond

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

2014-04-10

Abstract:The Multiscale Finite Element Method (MsFEM) is developed in the vein of Crouzeix-Raviart element for solving viscous incompressible flows in genuine heterogeneous media. Such flows are relevant in many branches of engineering, often at multiple scales and at regions where analytical representations of the microscopic features of the flows are often unavailable. Full accounts to these problems heavily depend on the geometry of the system under consideration and are computationally expensive. Therefore, a method capable of solving multiscale features of the flow without confining itself to fine scale calculations is sought after. The approximation of boundary condition on coarse element edges when computing the multiscale basis functions critically influences the eventual accuracy of any MsFEM approaches. The weakly enforced continuity of Crouzeix - Raviart function space across element edges leads to a natural boundary condition for the multiscale basis functions which relaxes the sensitivity of our method to complex patterns of obstacles exempt from the needs of implementing any oversampling techniques. Additionally, the application of penalization method makes it possible to avoid complex unstructured domain and allows extensive use of simpler Cartesian meshes.

Numerical Analysis

What problem does this paper attempt to address?

### What problem does this paper attempt to solve?

This paper aims to solve the problem of solving viscous incompressible flow (Stokes flow) in multi - scale heterogeneous media. Specifically, the research background and problem description are as follows:

1. **Complexity of multi - scale flow problems**:
- In many engineering fields (such as petroleum engineering, micro/nano - fluid mechanics, mechanical - biological systems, etc.), flow problems often involve multiple spatio - temporal scales.
- These multi - scale problems may be caused by the highly oscillatory coefficients of the system or the heterogeneity of the domain, for example, the complex rock matrix or the random distribution of buildings, people and trees in urban canopy flow.

2. **Requirements for computational resources**:
- Considering these problems comprehensively usually depends on the geometric structure of the system and requires a large amount of computational resources.
- Although high - performance computing has developed, the scale of discrete problems is still large.

3. **Limitations of traditional methods**:
- Although traditional fine - scale computational methods are accurate, they are inefficient in dealing with multi - scale problems. Especially when facing complex obstacle patterns, it is difficult to avoid inaccuracies in boundary conditions and the generation of complex meshes.

4. **Proposed method**:
- The paper proposes a non - conforming multi - scale finite element method (MsFEM) based on Crouzeix - Raviart elements to solve the above problems.
- This method represents the micro - structure model of the flow by pre - computing multi - scale basis functions on a fine grid, thus avoiding the dependence on fine - scale calculations.

5. **Key challenges and solutions**:
- **Treatment of boundary conditions**: The boundary conditions of multi - scale basis functions have an important impact on the final accuracy. The weak continuity of the Crouzeix - Raviart space naturally provides the boundary conditions of multi - scale basis functions, making it insensitive to complex obstacle patterns and eliminating the need for oversampling techniques.
- **Simplified grid**: The introduction of a penalty method enables simulation on a simple Cartesian grid, avoiding the generation of complex unstructured domains.

### Summary

The goal of the paper is to develop an efficient computational algorithm that can solve multi - scale viscous incompressible flow problems without relying on fine - scale calculations. By improving the treatment of boundary conditions of multi - scale basis functions and introducing a penalty method, this method can provide accurate flow simulation results in complex heterogeneous media while significantly reducing the requirements for computational resources.