A two-dimensional vertical (2DV), Eulerian two-phase model or complete two-fluid model of the free surface flow was developed to simulate water-sediment flow in a local scour hole. In the model, the complete forms of the vertical, two-dimensional, two-fluid Navier-Stokes equations were discretized using a finite volume scheme. This discretization was done based on a standard staggered grid system using a curvilinear network system in compliance with the bed boundaries and water level. At the beginning of the computational cycle, the equations goving the fluid phase were solved based on the two-step projection method with a pressure-correction technique. In the first step, the intermediate fluid velocities were obtained by solving different phases of the momentum equations of the fluid phase using the time-splitting technique. In the second step, pressure was obtained and fluid velocities were updated. In this step a simple discretization method was applied for decreasing the computational complexity. After obtaining all the fluid phase variables at a new time step, the sediment phase momentum equa-tions were solved using the time-splitting technique and sediment velocities were obtained. Then, at the end of the computational cycle, the sediment phase mass equation was solved and the concentrations of both phases were updated. At last, the capacity of the model for simulating of the longitudinal fluid velocity and sediment concentration in a local scour hole was evaluated. Numerical results were found to be in good agreement with experimental data.