To simulate 3D flow of a non-uniform and highly concentrated sediment, a numerical model using a multi-phase Eulerian method for air, water, and particles of various class size is developed. This model accounts for turbulence in pore water, particle-particle collisions, and enduring friction. To test its performance, simulations were performed for large scale debris-flow experiments. The numerical model reproduces successfully the flow property and sediment size segregation of the debris flow.