The reaction environment of the CO 2 electrochemical reduction process is highly complex with several coupled phenomena operating at vastly varying length and time scales that incorporate important degrees of freedom for the overall system performance. A multi-scale approach is therefore essential in order for computational simulations to effectively develop the understanding of such a reaction system to compliment the experimental efforts in the area.