This work set out to devise an end-to-end modeling framework to create a digital model of the PEMFC electrode fabrication process, while taking into account various common fabrication parameters. In the pursuit of resolving explicit carbonionomer distribution in CL microstructures, we propose a novel discrete element method-based aggregation model with a colloidal representation of ionomer and carbon particles. The process involves a brownian dynamics simulation of the particles, with the DLVO theory forming a basis for inter-particle conservative forces and the solvent approximated through an implicit representation using pair-wise hydrodynamic interactions.