Typically I have seen P2/P1 elements used for Stokes equation, but I want to use P2/(P1+P0) and P2/P0 elements because I want to ensure local mass conservation. When I say P2/(P1+P0), I simply mean the concatenation of the basis for P1 and P0 separately.
I ran simple problems in FEniCS with both element pairs and while the velocity works great in both cases (i.e., element-wise conservative), the pressure does get a little screwed up.
I have been told mixed opinions about using anything with piecewise constants for pressure. I have heard that P2/P0 is not LBB stable for 3D elements (however, I fail to see this happening because I do not get node-to-node spurious oscilliations) although this could potentially explain why my pressure looks screwed up. I have also heard that P2/P0 should be doable in 3D (hence P2/(P1+P0) should work as well).
So I guess my question is, theoretically are P2/P0 and/or P2(P1+P0) elements stable and usable for 3D Stokes and Darcy problems?