AbstractThis paper demonstrates the effectiveness of a high-resolution numerical scheme in boundary-fitted curvilinear coordinates for simulation of flow at channel bends. Two-dimensional shallow water equations, which are transformed from physical domain to computational domain, are solved by a robust high-resolution scheme particularly developed for boundary-fitted curvilinear coordinates. Finite-volume method is adopted for discretization of governing equations in computational domain. Turbulence shear stresses are included in the model by eddy-viscosity approach. Moreover, the effect of secondary flow, which is significant in channel bends, has been included in the model and the results have been compared with the results achieved by ne...