Add support for bubble enhanced P2 element in 2D

The conforming Crouzeix-Raviart element P_2^+-P_1^\text{disc} for triangles in 2D employs for the velocity component continuous quadratic Lagrange elements extended by a cubic bubble function. This MR adds to HOG support for this P_2^+ part.