We use LAMMPS, an open source molecular dynamics engine to perform our simulations . In particular, we use the meso-scale particle integration methods implemented in LAMMPS to perform discrete element method simulations. This allows us to perform arbitrarily large simulations by utilising the efficient domain-decomposition to handle thousands of particles using many cores.