A programed numerical tool that allows two-dimensional simulation of debris flows is presented in this dissertation. Cundall´s Distinct Element Method (DEM) is used to this purpose, which was developed in 1979. Following this methodology, motion equations are simultaneilly solved by the numerical method of Dynamic Relaxation (MRD) for each distinct particle. This method is used in order to the transient behavior of this particular problem. The numerical modeling is based on the SAND program, developed at PUC (2002) for a numerical application on sand production for petroleum extraction process considering biphasic flow motion. Two new features incorporated in this kind of analyses are the graphical representation of walls with cubic spline curves segmentation and the implementation of the Munjiza´s method for contact detection. They attemp to implement irregular curves that are closely to represent real sliding surface and to decrease the total program executation time respectively. The idoneousness of the programed numerical tool for the representation of the most caracteristic phisical mechanisms of these kind of flows is tested using several curves configurations as variation on the inicial parametrics values of the numerical model. Moreover, main features and associated terminologies for identification and classification of debris flows and similar events are largely described here because of the usual confusion in the use of them within scientific literatures and professional communications between technician of Geology, Geography and Geotechnical Engineering.