Stable and unstable cross-grid PkQlPkQl mixed finite elements for the Stokes problem

In this paper we develop and analyze a family of mixed finite element methods for the numerical solution of the Stokes problem in two space dimensions. In these schemes, the pressure is interpolated on a mesh of rectangular elements, while the velocity is approximated on a triangular mesh obtained by dividing each rectangle into four triangles by its diagonals. Continuous interpolations of degrees kk for the velocity and ll for the pressure are considered, so the new finite elements are called cross-grid PkQlPkQl. A stability analysis of these approximations is provided, based on the macroelement technique of Stenberg. The lowest order P1Q1P1Q1 and P2Q1P2Q1 cases are analyzed in detail; in the first case, a global spurious pressure mode is shown to exist, so this element is unstable. In the second case, however, stability is rigorously proved. Numerical results obtained in these two cases are also presented, which confirm the existence of the spurious pressure mode for the P1Q1P1Q1 element and the stability of the P2Q1P2Q1 element.