Simulation of corrosion reactions has become an important tool in the analysis of novel materials for corrosion resistant materials and applications. A number of analytical, boundary element and finite element methods currently exist in the literature, which use experimental corrosion parameters as model inputs for the extraction of simulated polarization curves, oxide growth rates, and other corrosion behaviors. Here, we propose the combination of finite element simulation with thermodynamic data from Pourbaix diagrams to allow for the simulation of metal corrosion, where individual corrosion reactions can be monitored at the metal surface. (C) The Author(s) 2019. Published by ECS.