An approximately divergence-free 9-node velocity element (with variations) for incompressible flows

Explicit basis functions are constructed for 9-node biquadratic velocity fields which guarantee that a weak form of the continuity equation is satisfied. The corresponding pressure approximations are either piecewise constant, piecewise linear or piecewise bilinear. These results are extended to give bases for bilinear velocity/piecewise constant pressure elements and also to some three-dimensional brick elements.