This report details the development and implementation of a Discrete Element Method (DEM) simulation designed to model the behavior of deformable two dimensional (2D) elliptical particles. The simulation implements particle-particle and particle-wall interactions, with the foundational contact model based on Hertzian theory, extended with modifications to account for geometric nonlinearities. A central focus of this project is on the deformation of elliptical particles under various loading conditions and the subsequent evolution of the deformations post-collision.