Locally energy-stable finite element schemes for incompressible flow problems: Design and analysis for equal-order interpolations

TL;DR

This paper introduces finite element schemes for incompressible flows that ensure local energy stability using equal-order interpolations, combining energy-preserving discretizations with pressure stabilization to achieve convergence and non-oscillatory solutions.

Contribution

It develops a novel finite element method that guarantees local energy stability and convergence for incompressible flow problems with equal-order interpolations.

Findings

Method is non-oscillatory in numerical tests

Achieves optimal convergence behavior

Ensures semi-discrete energy inequality and stability

Abstract

We show that finite element discretizations of incompressible flow problems can be designed to ensure preservation/dissipation of kinetic energy not only globally but also locally. In the context of equal-order (piecewise-linear) interpolations, we prove the validity of a semi-discrete energy inequality for a quadrature-based approximation to the nonlinear convective term, which we combine with the Becker--Hansbo pressure stabilization. An analogy with entropy-stable algebraic flux correction schemes for the compressible Euler equations and the shallow water equations yields a weak `bounded variation' estimate from which we deduce the semi-discrete Lax--Wendroff consistency and convergence towards dissipative weak solutions. The results of our numerical experiments for standard test problems confirm that the method under investigation is non-oscillatory and exhibits optimal convergence behavior.