Cathodic Protection of Steel in Reinforced Concrete
Cathodic protection (CP) is a common strategy for retarding the corrosion of reinforcing steel in concrete structures, such as bridges and parking garages. By the use of CP the potential of the corroding surface is lowered, thereby decreasing the rate of undesired anodic corrosion reactions.
This example models the cathodic protection of a steel reinforcing bar (rebar) in concrete. The corrosion cell consists of a zinc anode, the concrete, acting as electrolyte, and the steel surface. Iron oxidation, water reduction (hydrogen evolution), and oxygen reduction are considered on the steel surface, whereas oxygen and charge transport are accounted for in the concrete electrolyte.
The anode and the steel surface are connected electrically via a potentiostat that controls the cell voltage.
Concrete is a porous material, and hence its transport properties for ions and gases vary with the moisture content. Therefore the electrolyte conductivity and oxygen diffusion coefficient are modeled to vary with the concrete pore saturation level using empirical data.
The corrosion rate for various moisture contents is investigated in this model.
The model example is based on a paper by Muehlenkamp and others. (Ref. 1).
Model Definition
Geometry
Figure 4 shows the model geometry. The geometry represents a two dimensional cross section of a repeating unit cell in a larger structure where three symmetry planes (top, bottom, and right) have been used to reduce the model geometry. The zinc anode has been coated onto the concrete by thermal spraying and is assumed to be permeable to air.
Figure 4: Model geometry. One electrolyte domain and two electrode surfaces.
Concrete domain equations
A Tertiary Current Distribution, Nernst–Planck interface is used to model the current distribution in the cell. The electrolyte conductivity depends on the pore saturation level according to Figure 5.
Figure 5: Electrolyte conductivity (S/m) as function of the concrete pore saturation level.
The oxygen diffusivity depends on the pore saturation level according to Figure 6.
Figure 6: Oxygen diffusivity (m2/s) in the concrete as function of pore saturation.
Boundary Conditions
Choose the electric potential of the Zn anode as ground for the system. By assuming the kinetics of the Zn anode to be very fast, the polarization is neglected in the model and sets the electrolyte potential to
where Eeq,Zn is an empirically measured equilibrium potential for the thermally sprayed zinc/concrete electrode. In this model Eeq,Zn is set to 0.68 V.
The concentration of oxygen at the Zn anode is set to atmospheric conditions according to
Three different electrode reactions are considered on the steel rebar boundary: iron oxidation, oxygen reduction, and hydrogen evolution:
The reaction kinetics for these reactions is modeled with an Electrode Surface node in the Tertiary Current Distribution, Nernst–Planck interface, on which the external electric potential of the steel bar, ϕs,steel, is set to the applied cell potential of 1 V.
The electrode kinetics of the steel bar reactions are described by Tafel expressions according to
using the parameters shown in the table, where the overpotential for each reaction is calculated as
O2
H2
The oxygen reduction reaction causes a flux of oxygen at the steel surface according to Faraday’s law. The flux is set up using the stoichiometric coefficients entered in the Electrode Reaction node for the oxygen reduction reaction.
Symmetry conditions of the Insulation type apply on all other boundaries.
Atmospheric concentration is used as the initial value for the oxygen concentration variable.
Study
The model is solved using a parametric sweep over a stationary study step and solving for a range of pore saturation values from 0.2 to 0.8.
Results and Discussion
Figure 7 shows the electrolyte potential for a pore saturation level of 0.8. The electrolyte potential is lower toward the back (the right side) of the rebar.
Figure 7: Electrolyte potential for a pore saturation (moisture) level of 0.8.
Figure 8 on the next page shows the oxygen concentration in the electrolyte for a pore saturation level of 0.8. The concentration is very low close to the rebar, indicating that the oxygen reduction kinetics should be mass transport limited for this pore saturation level. The concentration is lower toward the back of the rebar.
Figure 8: Oxygen concentration for a pore saturation (moisture) level of 0.8.
An important factor for the corrosion rate of the rebar is the operating electrode potential, which is the difference between the electric potential (here the potential applied by the potentiostat) and the electrolyte potential. Figure 9 on the next page shows the operating electrode potential for various pore saturation levels for three different points (front, middle, and back) of the rebar surface. The potential drops considerably at a pore saturation level of 0.65.
Figure 9: Operating electrode potential for three points at the rebar-concrete interface.
Figure 10 shows the local oxygen concentration at the rebar for various pore saturation levels. The concentration drops significantly toward higher saturation levels. This is an effect of the decreasing diffusivity of oxygen in the concrete for higher saturation levels.
Figure 10: Local oxygen concentration at the rebar-concrete interface.
The local oxygen reduction current densities at the rebar are shown in Figure 11 on the next page. The magnitude of the reduction current is highest around a pore saturation level of 0.6–0.65. Up to this point the current densities are increasing due to increased electrolyte conductivity, but for higher pore saturation levels the current densities decrease due to decreased oxygen diffusivity.
Figure 11: Local oxygen reduction current densities at the rebar-concrete interface.
The hydrogen evolution currents are shown in Figure 12 on the next page. Hydrogen evolution is very limited below a PS level of 0.65, which is the saturation level at which the electrode potential gets below the equilibrium potential (1.03 V) for the hydrogen evolution reaction, see Figure 9.
Figure 12: Local hydrogen evolution current densities at the rebar-concrete interface.
Finally, the iron oxidation currents are shown in Figure 13. Corrosion current densities are higher for low PS levels, which is in line with the higher electrode potential for low PS levels (Figure 9). The magnitude of the iron oxidation current density is considerably smaller than the oxygen reduction and hydrogen evolution current densities at the steel rebar. This low magnitude indicates that the zinc coating applied onto the concrete is effective in protecting the steel rebar from being corroded.
Figure 13: Iron corrosion current densities at the rebar-concrete interface.
Reference
1. E.B. Muehlenkamp, M.D. Koretsky, and J.C. Westall, “Effect of Moisture on the Spatial Uniformity of Cathodic Protection of Steel in Reinforced Concrete,” Corrosion, vol. 61, no. 6, pp. 519–533, 2005.