Aqueous organic redox flow batteries (AORFB) have become the promising pitch in large scale energy storage facilities, which offer fast electrochemistry kinetics, and relatively low system cost [1]. Nonetheless, finding the perfect redox couples in AORFB cells among all the possible candidates remains challenging, while the electrochemistry mechanisms behind the cell operation remain generally poorly understood.Here we present a novel mesoscale kinetic Monte Carlo (kMC) algorithm combined with a dynamic electric double layer (EDL) approach allowing us to investigate in depth the electrochemical kinetics in a methyl viologen (MV) based AORFB cell. The resulting three-dimensions computational model, adapted from our previous works in different applications [2,3,4,5,6], simulates the stochastic processes happening at the electrode/anolyte interface from a molecular level viewpoint by assigning different rates to different events. In this study, three key events have been included in the model: molecular motion, electrochemical reaction, and dimerization. Both diffusion and electromigration are considered as the driving forces for molecular motion, while the electromigration effects are calculated through the EDL approach. The electrochemistry reaction rate considers the reorganization energy, which varies as a function of the electronic tunneling distance between the electrode and the anolyte. The dimerization event allows capturing MV's capacity degradation in agreement with experimental knowledge.The model has been used to simulate the galvanostatic discharging process with different input current densities and electrolyte concentrations. Along the simulation time, the model simulates the system's electrochemical response while providing insights on the dynamic EDL structure evolution and potential dynamics. The calculated observables are in good agreement with empirical knowledge advancing the understanding of these complex electrochemical interfaces' behavior. This mesoscale kMC model paves the way towards a tool able to scale up computational screening results arising from Density Functional Theory calculations into kinetics simulation of AORFB electrochemical interfaces, as being done by us in the context of the EU-funded project "SONAR" [7].