Multiscale method for Oseen problem in porous media with non-periodic grain patterns

TL;DR

This paper introduces a multiscale finite element method to accurately simulate inertial flow in porous media with complex, non-periodic grain patterns, addressing challenges in modeling flow at multiple scales.

Contribution

The paper develops a novel multiscale finite element approach based on Crouzeix-Raviart elements for non-periodic porous media, enabling accurate flow predictions without oversampling.

Findings

Demonstrates convergence and accuracy of the method

Shows effective modeling of complex grain patterns

Maintains solution accuracy near grains without oversampling

Abstract

Accurate prediction of the macroscopic flow parameters needed to describe flow in porous media relies on a good knowledge of flow field distribution at a much smaller scale---in the pore spaces. The extent of the inertial effect in the pore spaces can not be underestimated yet is often ignored in large-scale simulations of fluid flow. We present a multiscale method for solving Oseen's approximation of incompressible flow in the pore spaces amid non-periodic grain patterns. The method is based on the multiscale finite element method (MsFEM [Hou and Wu, 1997]) and is built in the vein of Crouzeix-Raviart elements [Crouzeix and Raviart, 1973]. Simulations of inertial flow in highly non-periodic settings are conducted and presented. Convergence studies in terms of numerical errors relative to the reference solution are given to demonstrate the accuracy of our method. The weakly enforced continuity across coarse element edges is shown to maintain accurate solutions in the vicinity of the grains without the need for any oversampling methods. The penalization method is employed to allow a complicated grain pattern to be modeled using a simple Cartesian mesh. This work is a stepping stone towards solving the more complicated Navier-Stokes equations with a non-linear inertial term.