The first inf-sup stable finite element spaces that yield divergence-free approximations for incompressible fluid models on simplicial triangulations were given in the classical paper by Scott and Vogelius [1]}. Their lowest order pair, defined on a two-dimensional domain, used quartic Lagrage finite elements and piecewise cubic polynomials for the discrete velocity and pressure spaces, respectively; since then, there has been several efforts to find lower order inf-sup stable and conforming finite element spaces that produce divergence-free approximations [2]}, [3]}, [4]}, [5]}, [6]}, [7]}, [8]}, [9]}, [10]}.
