A discrete numerical model of the front region in piping erosion