Abstract
Formation process of gully-headcut has been simulated through numerical methods. 2-D depth-averaged flow model with non-orthogonal boundary fitted grid system is employed to describe discharge distribution around gully-headcut. The heterogeneity of unit discharge along headcut promotes gully development process. Simulated result shows initial shape of the gully-headcut is very important for its further development and advancement. Headcut convex shape from up-stream side is a necessary condition to develop gully, while the concave shape is disadvantageous in gully development process. Higher convexity, especially around gully-head, makes the gully development process faster.
