By fixing the anodic current to the protonator (WE1), a flux of protons was applied at the surface of the protonators, and the flux was assumed to be proportional to the current applied to the electrodes. The initial pH within the experimental domain was set to pH 7. The proton diffusion coefficient used for the simulation was 9.31x10 -5 cm 2 s -1 . Once the initial static diffusion simulations were undertaken, the addition of a single-phase