Stability estimates and a Lagrange-Galerkin scheme for a Navier-Stokes type model of flow in non-homogeneous porous media

TL;DR

This paper establishes the L2-stability of a Navier-Stokes type model for non-stationary flow in porous media and introduces a Lagrange-Galerkin numerical scheme with Adams-Bashforth method, validated through numerical experiments.

Contribution

It provides the first stability analysis for this model and develops a novel Lagrange-Galerkin scheme tailored for non-homogeneous porosity flows.

Findings

The scheme achieves optimal convergence rates.

Numerical results match theoretical flow profiles.

The method effectively handles complex porosity structures.

Abstract

The purposes of this work are to study the $L^{2}$-stability of a Navier-Stokes type model for non-stationary flow in porous media proposed by Hsu and Cheng in 1989 and to develop a Lagrange-Galerkin scheme with the Adams-Bashforth method to solve that model numerically.The stability estimate is obtained thanks to the presence of a nonlinear drag force term in the model which corresponds to the Forchheimer term. We derive the Lagrange-Galerkin scheme by extending the idea of the method of characteristics to overcome the difficulty which comes from the non-homogeneous porosity. Numerical experiments are conducted to investigate the experimental order of convergence of the scheme. For both simple and complex designs of porosities, our numerical simulations exhibit natural flow profiles which well describe the flow in non-homogeneous porous media.