In this work we use Pourbaix diagrams as a means for simulating a polarization curve at a corroding iron surface. Pourbaix diagrams show the boundaries between the changing thermodynamically stable species at a metal surface in an aqueous environment as a function of the applied potential and pH. The position of these boundaries can therefore be used to model the onset of the corresponding oxidation and reduction reactions by combining the equation of the appropriate boundary line with Butler-Volmer kinetics. At the same time, the change in pH local to the metal surface is monitored by simulating the flux of protons generated during the oxidation process, and the impact of this on the corrosion potentials and rate is taken into account. This is of great importance as the corrosion rate and the corrosion product varies according to the pH at the metal surface.
In this way, we show a simple means for the simple simulation of a polarisation curve at an iron surface, which is in excellent agreement with an experimentally recorded curve under the same conditions. This same method can then be applied to more complex metal alloys such as stainless steels, by combining the Pourbaix diagrams for the appropriate alloy components. This allows the model to be used as a standalone analytical tool for the prediction of corrosion behaviour of novel alloys before they are developed, as well as for the validation of experimental data obtained from existing samples.