This paper is concerned with three-dimensional numerical analysis of turbulent flows around submerged spur-dikes. A finite volume method is employed to discretize the governing equations. A nonlinear k-ε model proposed by Shih et al. and a linear k-ε model modified by Zhu-Shih are applied to two cases of flow with different spur dikes. The comparison with measured data shows that the nonlinear k-ε model reproduces mean velocity field better than the modified linear k-ε model in the two cases. The reason of the superiority of the non-linear model is discussed. It is shown, however, that further improvement in the numerical model is required to reproduce the complex flow phenomena accurately.