A Numerical Method for 2-D Bed Morphology Calculations

A fully coupled two-dimensional sub-critical and/or supercritical, free-surface flow numerical model is developed to calculate bed variations in alluvial channels. Vertically averaged free-surface flow equations in conjunction with sediment transport equation are numerically solved using an explicit finite-volume scheme in integral form. The capabilities of the proposed method are first demonstrated by analyzing supercritical flow in an expansion channel. Thereafter, one and two-dimensional applications referring to aggradation and scouring are reported. For each of these test cases, computed results compare satisfactorily with measurements as well as with other numerical solutions. The method is stable, reliable and accurate, although time consuming, handling a variety of sediment transport equations with rapid changes of sediment transport at the boundaries.