Studies on sediment transport in rivers with steep bed slopes are facing several obstacles compared to delta rivers due to sudden changes in riverbed topography and rapidly changing flow regimes, etc. Sediment transport in rivers with steep bed slope is a complicated problem because the riverbed material is often heterogeneous and contains a variety of particles of different sizes such as soil, gravel, pebbles and boulders, etc. This study, therefore, initially develops a computer program to simulate the bed level variation processes in the steep channels based on the grain size fraction approach together with considering the horizontal movement of sediment particles to the bed channel maintenance process. The two–dimensional (2D) hydrodynamic equations and bed load sediment transport are solved by finite difference method (FDM) and written on the software language Fortran 90. The developed program is, then, applied to simulate the bed level variation in the Teabeak River section, South Korea after a severe flood event. The performance of the developed model is evaluated through the NASH and RMSE statistical indexes. The results carried out with NASH = 0.79–0.83 and RMSE = 13%–19% confirmed the appropriateness of the developed model, thereby results indicating the prospect of the model deployment in the future