Comparison of different 2 nd order formulations for the solution of the 2 D groundwater flow problem over irregular triangular meshes