Performance of several stabilized finite element methods for the Stokes equations based on the lowest equal-order pairs

In this paper the performance of various stabilized mixed finite element methods based on the lowest equal-order polynomial pairs (i.e., P (1) - P (1) or Q (1) - Q (1)) are numerically investigated for the stationary Stokes equations: penalty, regular, multiscale enrichment, and local Gauss integration methods. Comparisons between them will be carried out in terms of the critical factors: stabilization parameters, convergence rates, consistence, and mesh effects. It is numerically drawn that the local Gauss integration method is a favorite method among these methods.