Numerical modeling of gravitational erosion in rill systems