The development of slope land in Taiwan is rapidly growing, so is the potential hazard for debris flow. Therefore, identification of these potential hazard area will be an important task. The flowbility of debris flow is a control factor in estimation these area. This paper develops a numerical model utilizing nonlinear constitutive relationship and fluid conservation laws. Model can be used in field and the values of all parameters used in the model can be obtained through simple in situ tests. Results from numerical modeling fits those from flume tests very well.