As the velocity space is P2, while your mesh is P1, there is not a 1-to-1 map, such as the P1 space has with dof_to_vertex_map and vertex_to_dof_map as described here.
For a P2 space, you can use compute_vertex_values as shown in: Save the result of mpirun in a numpy array - #2 by dokken