The objective of this study is to develop an accurate and robust two-dimensional finite element method for turbulence simulation in open channels. The model is based on Streamline Upwind/Petrov-Galerkin finite element method and Boussinesq's eddy viscosity theory. The method developed in the study is depth-averaged mixing length model which assumes anisotropic and local equilibrium state of turbulence. The model calibration and validation were performed by comparing with analytical solutions and observed data. Several numerical simulations were carried out, which examined the performance of the turbulence model for the purpose of sensitivity analysis. The uniform channels that appear horizontal flow and vertical flow were carried out. The model was also applied to the Han river was in for the applicability test. The results were compared with the observed data. The suggested model displayed reasonable flow distribution compare to the observed data in natural river flow. As a result of this study, the two-dimensional finite element model provides a reliable results for flow distribution based on the turbulence simulation in open channels.