Flow near riverbank erosion was examined to investigate mechanism of bank erosion through experimental study of flow field of model bank erosion shape in the authors' previous study. This paper presents prediction of flow field near bank erosion. A numerical model is developed to compute two dimensional flow velocity and water depth for different model bank erosion shapes. The purpose of this study is to reproduce measured flow field of the model bank erosion shape. Curvilinear coordinate system is used in this analysis to obtain the bank erosion shape appropriately. The model could predict the flow field inside bank erosion for small erosion surface angle (equal or less then 4 degrees) with higher accuracy. It is also able to reproduce flow field inside erosion part for larger erosion surface angle (8degrees), but absolute value of the computed velocity is larger than the measured velocity.