Recent researches on debris flow is focused on understanding its movement mechanism and building a numerical simulator to predict its behavior. However, previous simulators emulating fluid-like debris flow have limitations in numerical stability, geometric modeling and application of various boundary conditions. In this study, depth integration is applied to continuity equation and force equilibrium for debris flow. Thickness of sediment, and average velocities in x and y flow direction are chosen for main variables in the analysis, which improve numerical stability in the area with zero thickness. Petrov-Galerkin formulation uses a discontinuous test function of the weighted matrix from DG scheme. Presented mechanical constitutive model combines fluid and granular behaviors for debris flow. Effects on slope angle, inducing debris height, and bottom friction resistance are investigated for a simple slope. Numerical results also show the effect of embankment at the bottom of the slope. Developed numerical simulator can assess various risk factors for the expected area of debris flow, and facilitate embankment design in order to minimize damage.