1 ABSTRACT
A numerical two-dimensional model for river morphology is extended with bank erosion. The model is formulated in curvilinear coordinates. This allows the use of a boundary-fitted computational grid, suited for rivers with a curved centre-line and a non-uni form width. Dealing with non-homogeneous bank erodibili ty is shown to require a non-orthogonal coordinate system. This enables the use of simple algebraic grid generators, but complicates the transformed equations.
A main feature of the present study is that all coordinate derivatives in the transformed equations are replaced by grid properties that are invariant for grid rotation, such as local grid skewness and grid line divergences. This facilitates the physical interpretation of transformation terms and the estimation of their contribution to the solution.
The model is verified by comparing its results with results from the original model, and by comparing computations on non-orthogonal grids with a computation on an non-orthogonal grid. An examination of the truncation error reveals the importance of a smooth grid point distribution for accuracy.