Fully consistent lowest-order finite element methods for generalised Stokes flows with variable viscosity

Computers & Mathematics with Applications, vol. 188, pp. 40–49

Abstract

In finite element methods for incompressible flows, the most popular approach to allow equal-order velocity-pressure pairs are residual-based stabilisations. When using first-order elements, however, the viscous part of the residual cannot be approximated, which often degrades accuracy. For constant viscosity, or by assuming a Lipschitz condition on the viscosity field, we can construct stabilisation methods that fully approximate the residual, regardless of the polynomial order of the finite element spaces. This work analyses and tests two variants of such a fully consistent approach, with the generalised Stokes system as a model problem. We prove unique solvability and derive expressions for the stabilisation parameter, generalising some classical results for constant viscosity. Numerical results illustrate how our method completely eliminates the spurious pressure boundary layers typically induced by low-order PSPG-like stabilisations.

Authors 2

1. Felipe Galarce

   Pontificia Universidad Católica de Valparaíso

   Affiliation as printed

   School of Civil Engineering, Pontificia Universidad Católica de Valparaíso, Chile

2. Douglas R. Q. Pacheco corresponding Aachen Center for Simulation and Data Science Chair for Computational Analysis of Technical Systems (CATS) Chair of Methods for Model‐based Development in Computational Engineering

   RWTH Aachen University

   Affiliation as printed

   Center for Simulation and Data Science (JARA-CSD), RWTH Aachen University, Germany

   Chair for Computational Analysis of Technical Systems, RWTH Aachen University, Germany

   Chair of Methods for Model-based Development in Computational Engineering, RWTH Aachen University, Germany