HOME JOURNALS CONTACT

Journal of Applied Sciences

Year: 2009 | Volume: 9 | Issue: 13 | Page No.: 2457-2462
DOI: 10.3923/jas.2009.2457.2462
Mathematical Modeling of Corrosion Inhibition Behavior of Low Carbon Steel in HCl Acid
A.A. Khadom, Aprael S. Yaro, Abdul Amir H. Kadum and Ahmed S. Altaie

Abstract: In this study, the nonlinear region of polarization curve near the corrosion potential can be discussed depending on data of over potential as a function of current densities. These data can be analyzed by suggestion a mathematical model take into account the effect of mass transfer on activation process. This model was used to evaluate the values of polarization resistance PR in presence and absence of phenylthiourea (PTU) as a corrosion inhibitor of low carbon steel in hydrochloric acid. The values of polarization resistance obtained by mathematical model compared with these obtained by Tafel parameters. The suggested β-model was very useful in evaluating the values of PR. The maximum value of polarization resistance was 775.5 mΩ cm2 at 30°C and 1M HCl in presence of 1 g L-l PTU, while the minimum value was 3.1 mΩ cm2 at 60°C and 5 M HCl in absence of PTU. These results indicate that the polarization resistance values increased with increase in PTU concentration and decreased with increase in both temperature and acid concentration.

Fulltext PDF Fulltext HTML

How to cite this article
A.A. Khadom, Aprael S. Yaro, Abdul Amir H. Kadum and Ahmed S. Altaie, 2009. Mathematical Modeling of Corrosion Inhibition Behavior of Low Carbon Steel in HCl Acid. Journal of Applied Sciences, 9: 2457-2462.

Keywords: polarization resistance, Current density, phenylthiourea, over potential and Stern and Geary equation

INTRODUCTION

Many researches and literature concentrated on corrosion inhibition, but few of them deal with the effect of inhibitor concentration on polarization resistance. The polarization resistance is a measure of a metal to undergo corrosion in a certain environment, it is a slope of the potential-current (E-i) curve at the corrosion potential, i.e.,

The shape of this curve around the corrosion potential and hence PR, is determined by the sum of the partial E-i curve both anodic and cathodic half reaction (Bardal, 2004) The shape of both curves, in turn, depends on the respective kinetic parameters of their reaction. The current-potential curve of corroding metal is rather complex non-linear equation; hence a general analytical solution for PR is equally complex. Under some conditions, the E-i relation reduces to a simple form and hence a simple analytical solution for PR can be obtained. The Stern and Geary equation (Bardal, 2004) which can be written in equation as:

(1)

where, η is over potential (i.e., E-Ecorr), icorr is corrosion current density (A cm-2), ba and bc are anodic and cathodic Tafel slopes (mV). Equation 1 assumed that the applied current is a linear function of electrode potential (about 10 mV around corrosion potential) and slope of this linear relationship is PR. Equation 1 will have high percent of error because of the non-linearity of (η-i) relationship near the corrosion potential. According to Stern, the quotient (∂η/∂i) or (ΔE/Δiapp) is called the polarization resistance by analogy to resistance = (E/l) from conventional electricity. Other has suggested that the inverse of (Δiapp/ΔE), is called the corrodance by obvious analogy to the inverse of electrical resistance, which is called the conductance (Allen and Faulkner, 2001).

Recently, Mansfeld (2005) stated that, as new experimental technique mature and are more commonly used, various misconception and misnomers are sometimes perpetuated, which, when used often enough, are considered to be truth. Some of these myths, misconception and misnomers are the linear polarization technique. There exists the polarization resistance (PR) technique and linear polarization technique, but there is no Linear Polarization Resistance (LPR) technique (Mansfeld, 2005). As pointed out many years ago (El-Feki and Walter, 2000) polarization curves recorded in the vicinity of the corrosion potential (Ecorr) are linear only in the exceptional case that the anodic (ba) and cathodic (bc) Tafel slopes are equal, i.e., b c= ba. In the general case there can be appreciable curvature at Ecorr. It is surprising that the little use is made of experimental data recorded for determining PR values to calculate Tafel slopes and icorr values from data in pre-Tafel region (Sastri et al., 2007).

The term polarization resistance (PR) must be not confused with resistance polarization (R). PR is polarization resistance at electrode/solution interface, knowing PR one can calculate electrochemical reactions rates (Sastri, 2001). PR represents the charge transfer resistance resulting from the separation of charges across the solid interface to the outside edge of the double layer. While R represents the resistance through the bulk solution (Rs) and the resistance of any films present on the working electrode surface (Rf) (i.e., R = Rs+Rf) (Trethewey and Chamberlain, 1996). The goal of this research is to study the effect of temperature, acid concentration and inhibitor concentration on polarization resistance. Two models suggested to evaluate the values of PR, one of them take into account the effect of mass transfer on the activation control corrosion reaction and the deviation of PR from linearity. While, the other model is Stern and Geary equation which neglect the effect of mass transfer and assuming that the corrosion reaction is completely activation control.

MATERIALS AND METHODS

Specimens of cylindrical shape with outside diameter of 1.96 cm and a length of 0.5 cm, of low carbon steel were used. The specimens were first decreased with analar benzene and acetone and then annealed in a vacuum at 600°C for 1 h and cooled to room temperature. Samples were abraded in sequence under running tap water using emery paper of grade number 220, 320, 400 and 600 then washed with running tap water followed by distilled water, dried with clean tissue, immersed in acetone and benzene, kept in desiccators over silica gel bed until use. Corrosion cell has four necks, one was fitted with working electrode, one for immersing a thermometer in order to observe the test temperature and the other one had a spherical joint for manipulating the lugging capillary probe. The probe was adjusted to be in a distance not more than 2 mm from the working electrode. The fourth necks input graphite electrode. All potential values were measured in reference to Saturated Calomel Electrode (SCE). DC power supply which is equipped with current and voltage limiters offers better stability and control, its type is 6236B, triple output supply Hewlet Packard, USA. In present study, the corrosion reaction of low carbon steel in 1, 3 and 5 M HCl, at 30, 40, 50 and 60°C, in presence of 0.1, 0.75 and 1 g L-l phenylthiourea (PTU) as a corrosion inhibitor were studied. The low carbon steel working electrode specimens have the following chemical composition C 0.041 wt.%, Mn 0.311%, P 0.05%, S 0.007% and the remainder is iron.

Mathematical model derivation: The generalized equation for the polarization curve for metal corroding is given in Eq. 2 (Bockris et al., 1993):

(2)

where, λ is mass-transport correction factor,

and the dimensionless factor (icorr/il) can be substituted by β. The cathodic part of the equation is mathematically equivalent to those of Albaya et al. (1973) and Flitt and Schweinsberg (2005), but it is expressed in a more general and easily applicable form. Massiani et al. (1984) use a correction factor of (1-i/il) in place of λ. The expression of Britz and Hougaard (1983) can be shown to be an approximation of Eq. 2 when (icorr/il) <<1.

If the exponential terms of Eq. 2 are expanding in series, we come to the equation of polarization curve by the way of degree polynomial:

(3)

where, Co, C1, C2, ...Cn are constants. Equation 3 is a non-linear equation of (η-i) relationship. The coefficient of η in Eq. 3 can be obtained by a non-linear estimation method using a statistical computer program from the set of experimental points of current densities vs. over potentials for the corrosion of low carbon steel in 1, 3 and 5 M HCl acid in absence and presence of 0.75 and 1 g L-l phenylthiourea (PTU) as a corrosion inhibitor, at different temperatures of 30, 40, 50 and 60°C. Equation 3 can be written in the form of the Maclaurin formula as (Andrei and Alexander, 2007):

(4)

By comparison of Eq. 3 with Eq. 4, the polynomial coefficient can be written as:

(5)

The dimension less factor (icorr/il) or (β) can be used to determine the corrosion mechanism. The values of β have two limiting cases; β →ÿ 0 the corrosion kinetic is subjected by slow charge transportation; at β = 1 the limiting stage of total process is the depolarizer transport to electrode surface. It was suggested that in the range of β<0.05 and β>0.95, it is possible to use the simplified form of Eq. 2 and 5, (for activation and diffusion mechanism) with out precision of calculations evaluation. The approach that suggested here was to use the values of the second coefficient of β-model to evaluate the value of PR, i.e., by using Eq. 5. It is clear that when β = 0, C1 = (1/PR). This analysis may take into account the nonlinearly of polarization curve near the corrosion potential.

RESULTS AND DISCUSSIONS

Polarization experiments were carried out in acid solutions at different inhibitor concentration for each temperature and acid molarities as shown in Table 1 which shows the values of PR at different conditions. Table 1 shows 48 experimental tests.

The effect of inhibitor concentration on values of PR was calculated as mentioned previously from the values of C1 obtained from β-model. These values compared with the values obtained using Tafel extrapolation method depending on the values of ba, bc and icorr (i.e., PR = babc/(2.303 icorr (ba +bc))). As shown in Fig. 1-6, the values of PR in presence of inhibitor were higher than the values in absence of inhibitor. The values of polarization resistance decrease with both acid concentration and temperature increasing.

Figure 1 and 2 shows the effect of inhibitor concentration on polarization resistance. The values of PR in presence of inhibitor were higher than the values in absence of inhibitor, which indicates that the resistance to the corrosion was increased. The values of polarization resistance PR represent the charge resistance resulting from the separation of the charges across the solid interface to the outside edge of the double layer, these values were increased by increasing of inhibitor’s concentration.

Fig. 1: Variation of polarization resistance with inhibition concentration at 30°C at different acid concentration

Table 1: Polarization resistance of steel in HCl at different inhibitor concentration, temperature and acid concentration

Fig. 2: Variation of polarization resistance with inhibition concentration at 40°C at different acid concentration

Fig. 3: Variation of polarization resistance with acid concentration at different PUT concentration at 50°C

This suggests that as inhibitor concentration increase the fraction of surface of metal which covered by inhibitor will increas and this may prevent or reduce the charge transfer from the metal to the solution and hence, increase the value of PR. Mahmoud (2005) studied the inhibition effect of Pyrazol on the corrosion of metal in 1 M hydrochloric acid. The values of PR are deduced for the working electrode from the slopes of the straight lines obtained from potential-current plots (ΔE/Δi), he found that the values of PR for the working electrode increase by increasing the concentration of inhibitors, which is the same as in present study.

Figure 3 and 4 show the selective effect of acid molarities on polarization resistance. In absence of inhibitor the values of PR reduced as acid concentration increased. Normally, as the concentration of corrosive acid media is increased, the corrosion rate is likewise increased.

Fig. 4: Variation of polarization resistance with acid concentration at different PUT concentration at 40°C

Fig. 5: Variation of polarization resistance with temperature at different PUT concentration at 5 M HCl

This is primarily due to the fact that the amounts of hydrogen ions, which are active species, are increased. It is known that corrosion process consists of two or more partial reactions. For the corrosion of steel in acids, on of the reactions being the anodic reaction consists of an oxidation of the metal so that it changes from metallic state to an ionic state, Fe →Fe++ + 2e, while, the cathodic reaction occurs at different site from the anodic one, as follows 2H+ + 2e →H2. Both cathodic and anodic reactions proceed at the same rate during electrochemical corrosion. Therefore, increasing in hydrogen ions concentration will encourage the hydrogen evolution or cathodic reaction and at the same time, this will leads to increase in anodic reaction, i.e., increasing in separation of metallic ions. This behavior will reduce the resistance of metal to the corrosion as acid concentration increased. The same behavior observed in presence of inhibitor, the values of PR also reduced as acid concentration increased.

Fig. 6: Variation of polarization resistance with temperature at different PUT concentration at 1 M HCl

The values of PR were higher in presence of inhibitor due to the formation of protective layer which reduces the corrosion process.

Temperature increases the rate of almost all chemical reaction (Shreir et al., 2000). Like most chemical reactions, the rates of corrosion of metal in aqueous acid solutions increase with temperature increasing and this leads to reduction in metal resistance. Therefore, as shown in Fig. 5 and 6, the polarization resistances decrease by increasing in temperature at any acid concentration. In presence of inhibitor, the effect of increasing in temperature on inhibited acid-metal reaction can be highly complex. The rate of the uninhibited, acid-metal, heterogeneous reaction increased and the fraction of the metal surface covered by adsorbed inhibitor changes, usually decrease (Philippe, 2002). These two reasons increase the corrosion rates and decrease the resistance of metal. Therefore, the same behavior of PR observed in the presence of different concentration of PTU at different acid concentration. Generally, the values of PR obtained by β-model were lower than that obtained using Tafel parameters. This β-model behavior explain the mixed influence of mass transfer and activation effect on corrosion process which leads to decrease the metal resistance, while the second model neglecting the mixed effect and hence a higher values of polarization resistance were obtained.

CONCLUSION

The values of PR were obtained from the set of data of over potential as a function of current densities. The values of PR increased by increasing in inhibitor concentration and this may be, due to the formation of protective film which increases the resistance of metal to corrosion. The values of PR decreased with both temperature and acid concentration increase. This attributed to the increasing in corrosion rates with hydrogen ion concentration and with temperature increase, which lead to reduction in metal resistance. The highest value of polarization resistance was 775.5 mΩ cm2 at 30°C and 1 M Hcl in presence of 1 g L-l PTU, while the lower value was 3.1 mΩ cm2 at 60°C and 5 M HCl in absence of PTU. These results indicating the improvement of polarization resistance by increasing inhibitor’s concentration and decreasing of temperature.

ACKNOWLEDGMENT

The authors are gratefully acknowledge Universiti Kebangsaan Malaysia for assist and support of this study (Grant No. UKM-GUP-BTT-07-25-170).

REFERENCES

  • Albaya, H.C., O.A. Cobo and J.B. Bessone, 1973. Some consideration in determining corrosion rates from linear polarization measurements. Corrosion Sci., 13: 287-287.


  • Bardal, E., 2004. Corrosion and Protection. 11th Edn., Springer-Verlag, London, pp: 199-201


  • Bockris, J., N.G. Smart and M. Gamboa-Aldeco, 1993. Corrosion mechanisms of iron in concentrated acidic zinc chloride media. Corrosion Sci., 34: 759-777.


  • Britz, D. and G. Hougaard, 1983. Effects of mixed activation and transport control of the cathodic process on corrosion current measurements. Corrosion Sci., 23: 987-987.


  • Shreir, L.L., R.A. Jarman and G.T. Burstein, 2000. Corrosion, Metal/Environment Reactions. 3rd Edn., Butterworth-Heinemann, Britain


  • Philippe, M., 2002. Corrosion Mechanisms in Theory and Practice. 2nd Edn., Marcel Dekker, Inc., USA., ISBN: 0 8247 0666-8


  • Flitt, H.J. and D.P. Schweinsberg, 2005. Evaluation of corrosion rate from polarization curves not exhibiting a Tafel region. Corrosion Sci., 47: 3034-3052.
    CrossRef    


  • Andrei, D. and V.P. Alexander Manzhirov, 2007. Handbook of Mathematics for Engineers and Scientists. Taylor and Francis Group, LLC., ISBN: 1584885025


  • Mahmoud, S.S., 2005. Corrosion behavior and its inhibition of a copper-iron alloy in hydrochloric acid solution. Corrosion, 12: 1127-1135.
    Direct Link    


  • Mansfeld, F., 2005. Editorial: Misconceptions and misnomers. Corrosion, 61: 1019-1019.


  • Sastri, V.S., E. Ghali and M. Elboujdaini, 2007. Corrosion Prevention and Protection. Johan Wiley and Sons Publisher, USA


  • Massiani, Y., J.P. Crousier, J. Galea and R. Romanetti, 1984. Determination des parametres des courbes de polarisation a l'aide d'un programme de calcul par regression non lineaire. Application a la corrosion d'un metal en milieu acide en presence d'oxygene dissous. Electrochim Acta, 29: 1679-1679.


  • El-Feki, A.A. and G.W. Walter, 2000. Corrosion rate measurements under conditions of mixed charge transfer plus diffusion control including the cathodic metal ion deposition partial reaction. Corrosion Sci., 42: 1055-1070.
    CrossRef    


  • Allen, J.B. and L.R. Faulkner, 2001. Electrochemical Methods Fundamentals and Applications. 2nd Edn., John Wiley and Sons, Inc., London


  • Sastri, V.S., 2001. Corrosion Inhibitors Principles and Applications. Johan Wiley and Sons Publisher, New York, ISBN: 0-471-97608-3


  • Trethewey, K.R. and J. Chamberlain, 1996. Corrosion for Science and Engineering. 2nd Edn., Longman, USA., ISBN: 0582238692

  • © Science Alert. All Rights Reserved