A finite element method for piping erosion in levees