Numerical simulations of continental collision provide insights into how deformation and surface erosion shape mountain ranges. Given the high strain levels involved in orogeny and the discontinuous velocity boundary conditions characteristic of such systems, a discontinuum numerical approach can be useful for such simulations. In this study, the distinct element method has been applied to investigate collisional orogeny. The key controls on the kinematics and surface morphology of mountain belts two-dimensional experimental are incorporated in the model. These include (1) asymmetric flow of the continental crust material from the upper and lower plate sides, (2) flexural isostasy controlling large scale vertical motions, and (3) denudation at the orogen’s surface. The evolution of the model captures fundamental characteristics of collisional orogens, such as strain distribution, particle movement, and exhumation patterns.