Regularity of the diffusion-dispersion tensor and error analysis of Galerkin FEMs for a porous media flow

TL;DR

This paper introduces a new error analysis method for Galerkin finite element methods applied to porous media flow, overcoming regularity limitations of the diffusion-dispersion tensor by using a parabolic projection approach.

Contribution

A novel parabolic projection approach is developed for error analysis, requiring only Lipschitz continuity of the diffusion tensor, unlike traditional methods needing higher regularity.

Findings

Established optimal $L^p((0,T);L^q)$ error estimates.

Achieved almost optimal $L^((0,T);L^)$ error estimate.

Provided a new framework for error analysis under weaker regularity assumptions.

Abstract

We study Galerkin finite element methods for an incompressible miscible flow in porous media with the commonly-used Bear-Scheidegger diffusion-dispersion tensor $D({\bf u}) = \Phi d_m I + |{\bf u}| \big ( \alpha_T I + (\alpha_L - \alpha_T) \frac{{\bf u} \otimes {\bf u}}{|{\bf u}|^2}\big)$. The traditional approach to optimal $L^\infty((0,T);L^2)$ error estimates is based on an elliptic Ritz projection, which usually requires the regularity of $\nabla_x\partial_tD({\bf u}(x,t)) \in L^p(\Omega_T)$. However, the Bear-Scheidegger diffusion-dispersion tensor may not satisfy the regularity condition even for a smooth velocity field ${\bf u}$. A new approach is presented in this paper, in terms of a parabolic projection, which only requires the Lipschitz continuity of $D({\bf u})$. With the new approach, we establish optimal $L^p((0,T);L^q)$ error estimates and an almost optimal $L^\infty((0,T);L^\infty)$ error estimate.