During the past decades, many numerical models have been used to predict responses of asphalt mixes under different types of loading. Some of these models were simple due to practicality but overestimated the response of asphalt mixes. On the other hand, sophisticated but effective numerical models have been developed to address the shortcomings of the simpler models, and were used mostly in finite element analysis (FEA). However, these models were complicated and not user friendly. Recently, the approach of the discrete element method (DEM) was adopted. Unlike traditional FEA, DEM can simulate crack propagation by allowing the separation of elements in the simulated models. Understanding these challenges, this study was initiated to investigate the utilization of a simple visco-elasto-plastic model that had been used successfully in predicting deformation in asphalt mixes using the DEM embedded in Particle Flow Code in Two Dimensions (PFC2D) software simulations. Simulation results, when compared to flow time (FT) and number (FN) test results, showed that this model could simulate actual tests, thus predicting deformation of asphalt mixes using the DEM on a larger scale.