In order to optimise the electrochemical behaviour of CLs, it is important to quantify the effect of fabrication parameters on the microstructure. In the process of proposing a novel DEM model to probe the co-aggregation of Pt/C and ionomer dispersed in a solvent phase, we perform two crucial approximations. We represent the Pt/C catalyst as an effective colloidal particle. Similarly the ionomer is approximated as a collection of colloidal beads. We use LAMMPS as the Newtonian dynamics engine to track the evolution of particle positions with time. The DLVO theory proposes calculating interparticle conservative interactions as a combination of attractive van der Waals forces and repulsive electrostatic forces.