Numerical modelling of backward front propagation in piping erosion by DEM-LBM coupling