For single-phase and multiphase discrete element method model development, I am using following codes