Fourier analysis of an equal-order incompressible flow solver stabilized by pressure gradient projection

Abstract

Fourier analysis techniques are applied to the stabilized finite element method (FEM) recently proposed by Codina and Blasco for the approximation of the incompressible Navier-Stokes equations, here denoted by pressure gradient projection (SPGP) method. The analysis is motivated by spurious waves that pollute the computed pressure in start-up flow simulation. An example of this spurious phenomenon is reported. It is shown that Fourier techniques can predict the numerical behaviour of stabilized methods with remarkable accuracy, even though the original Navier-Stokes setting must be significantly simplified to apply them. In the steady state case, good estimates for the stabilization parameters are obtained. In the transient case, spurious long waves are shown to be persistent when the element Reynolds number is large and the Courant number is small. This can be avoided by treating the pressure gradient projection implicitly, though this implies additional computing effort. Standard extrapolation variants are unfortunately unstable. Comparisons with Galerkin-least-squares (GLS) method and Chorin's projection method are also addressed. Copyright (C) 2000 John Wiley and Sons, Ltd.