Debris flow has caused severe human casualties and economic losses in landslide-prone areas around the globe. A comprehensive understanding of the morphology and deposition mechanisms of debris flows is crucial to delineate the extent of a debris flow hazard. However, due to inherent complex field topography and varying compositions of the flowing debris, coupled with a lack of fundamental understanding about the factors controlling the geomaterial flow, interparticle interactions and its final settlement resulted in a limited understanding of the flow behaviour of the landslide debris. In this study, a physical model was set up in the laboratory to simulate and calibrate the debris flow using PFC, a distinct element modelling-based software. After calibration, a case study of the Varunavat landslide was taken to validate the developed numerical model. Following validation with an acceptable level of confidence, several models were generated to evaluate the effect of slope height, slope angle, slope profile, and grain size distribution of the dislodged geomaterial in the rheological properties of debris flow. Both qualitative and quantitative analysis of the landslide debris flow was performed. Finally, the utility of retaining wall and their effect on debris flow is also studied with different retaining wall positions along the slope surface.