The conforming Crouzeix-Raviart element P_2^+-P_1^\text{disc} is a stable mixed pairing for discretising the Stokes and associated problems, see e.g.

• M. Crouzeix and P.-A. Raviart, Conforming and nonconforming finite element methods for solving the stationary Stokes equations I, Revue française d'automatique, informatique, recherche opérationnelle. Analyse numérique, 1973, 7, 33-75
• Rolf Krahl and Eberhard Bänsch, Computational Comparison between the Taylor-Hood and the Conforming Crouzeix-Raviart element, Proceedings of ALGORITMY 2005, 2005, 369-379
• Cedric Thieulot and Wolfgang Bangerth, On the choice of finite element for applications in geodynamics. Part II: A comparison of simplex and hypercube elements, 2024, preprint