A debris flow model based on the constitutive equation derived from energy conservation and preferable in debris flow and sediment laden flow was developed. The model also considered the erosion and deposition effects contributed by the concentration variation of debris flow. The model was applied to simulate the debris flow disaster that occurred in Jiaopu Creek , Hsiaolin Village during typhoon Morakot. The scenario of the Jiaopu Creek debris flow event was considered in two stages. The primary purpose in the first stage was to see if the elevation of debris flow was enough to bury bridge no. 8. The simulation result and 5m DEM made after the Hsiaolin landslide dam failure were combined as the topographic data in the second stage. The simulation result was comparable to the survey data.