Research Article | Open Access
Rachid Radouani, Younes Echcharqy, Mohamed Essahli, "Numerical Simulation of Galvanic Corrosion between Carbon Steel and Low Alloy Steel in a Bolted Joint", International Journal of Corrosion, vol. 2017, Article ID 6174904, 10 pages, 2017. https://doi.org/10.1155/2017/6174904
Numerical Simulation of Galvanic Corrosion between Carbon Steel and Low Alloy Steel in a Bolted Joint
The galvanic corrosion of a bolt joint combining carbon steel end plate and low alloy steel bolt was investigated electrochemically in a 1 M HCl solution. The corrosion parameters of the joint components were used for numerical simulation using Comsol Multiphysics software to analyze the galvanic corrosion behavior at the contact zone between the head bolt and the end plate. In this research work we evaluate the variation of the corrosion rate in the steel end plate considered as the anode, in order to determine the lifetime of the bolted assembly used in steel structures. Three materials (20MnCr5, 42CrMo4, and 32CrMoV13) and three bolts (M12, M16, and M20) were tested in two thicknesses of electrolyte 1 M HCl ( = 1 mm, = 20 mm). It is found that the corrosion rate of the anode part (end plate) is higher for 32CrMoV13 materials and it increases if both diameter of the bolt and thickness of the electrolyte increase (Cr(M20) > Cr(M16) > Cr(M12) and Cr( = 20 mm) > Cr( = 1 mm)). This corrosion rate is higher in the contact area between the bolt head and the end plate, and it decreases if we move away from this contact area.
Galvanic corrosion can simply be defined as the corrosion that occurs as a result of one metal being in contact with another in a conducting, corrosive environment. The corrosion is stimulated by the potential difference that exists between the two metals: the more noble material acting as a cathode where some oxidizing species is reduced and the more active metal, which corrodes, acting as the anode. The anodic reaction is, by definition, some form of metal dissolution; the cathodic reaction is, in the majority of practical situations, either oxygen reduction or hydrogen evolution, or a combination of both. Many factors affecting galvanic corrosion are already discussed to determining whether or not galvanic corrosion will occur in a particular instance and if so at what rate; it is important when considering the theory of galvanic corrosion to be aware of these factors including electrode potential, reaction kinetics, alloy composition, protection film characteristics, bulk solution environment, total geometry, and type of joint .
There is a high incidence of past scientists taking an interest in corrosion to understand what causes it and what limits or accelerates the process. Numerous studies have been conducted; some take a more global outlook , whereas some take a more focused approach . The study conducted in  looked at many different galvanic couples commonly used in seawater applications. The study focused on developing reasonable models for systems experiencing varying periods of exposure to the corrosive environment.
Simulation of galvanic corrosion between magnesium and aluminum has been performed by Lacroix et al. , Deshpande [5–7], Jia et al. , and Trinh et al.  who have studied the corrosion of magnesium alloys in contact to mild steel under static conditions. The publications of Murer et al. [10–12] and Shi and Kelly  in this context also gave an extended insight into the topic, especially in the very important choice of boundary conditions. New studies of Sun et al. , who applied the mathematical approach of Yan et al.  to the modeling of deposit formation under seawater conditions, clearly introduce a possible way of a useful model built up for the mentioned purpose. The following studies and results are based on the progress achieved by them. Basic galvanic current density computations were modified by layer growth aspects leading to time dependent variations in the electrochemical response of the electrodes.
In the work of Höchez , simulation of galvanic corrosion between Mg and Al on a geometry like a self-pierce punch rivet joint has been performed. Starting from the initial model setup and related to the scientific challenge, the mathematical model requires fundamental assumptions. As always in modeling action they define the quality and accuracy of simulation results. The analysis via computer simulations showed the ability to study parameters influencing corrosion performance like the affected surface fraction by a reaction product.
In the works of Johnson and Abbott  and Xu et al. , the impact of the nature of the materials on the corrosion rate of mild steel in a galvanic coupling was studied. These articles show that the corrosion rate of mild steel changes if the material is changed in the galvanic corrosion. The addition of Cr in steel alloys improves the corrosion resistance of the steel. It is found that as the Cr content increases, the corrosion rate of the mild steel decreases.
Finite element analysis (FEA) of potential and current distributions in galvanic systems has long been studied in the literature [21–27]. Such studies are often carried out to investigate fundamental effects of electrolyte geometry [26, 28, 29], electrode kinetics , and unique part geometries . Early modeling attempts have provided an analytical solution [26, 30] for the galvanic current and potential distribution around a galvanic junction involving two dissimilar metals kept side by side and covered by a corroding solution. These theoretical predictions were also later validated with experimental measurement of the potential and current density across the junction. Recently, semianalytical and numerical methods [29, 31] were developed to solve for current density distribution on anode and cathode in a galvanic couple and similar electrochemical phenomena such as electrodeposition. Galvanic corrosion is found to be better modeled with a continuum approach than a lattice model since the curvature of the interface is an important parameter in determining the electrochemical dissolution rate.
In our work the finite element method was used to perform a parametric study of the galvanic corrosion of a column-beam bolted assembly used in steel structures. This study takes into consideration the material of the cathode “bolt,” the thickness of the electrolyte “,” the dimensions of the bolt head, and the distance from the contact area between the head bolt and the end plate in order to determine the parameters which can cause the maximum corrosion rate of the anode “end plate.”
2. Materials and Methods
2.1. Bolted Joint and Materials
Three types of materials used for the bolt were tested to predict the corrosion rate of the end plate. Thus three types of galvanic connections were studied in a bolted joint. The material composition of the elements of the assembly is stated in Table 1. The bolt was made by three types of low alloy steel which have the higher corrosion potential, and the metal of the end plate is the carbon steel (S235JR), Table 2. Galvanic corrosion therefore took place in the contact zone between end plate and head bolt in each bolted joint. In other words, the end plate acted as anode, and the bolt becomes cathode.
2.2. Geometry and Modeling
The bolted assembly is column-beam type and it is used in steel structures, Figure 1. The developed model is Cartesian 2D type and it is used for a representation of the plane modeling of the elements entering into the galvanic corrosion, Figure 2. The modeled part is the area of contact between the bolt head, the end plate, and the electrolyte. The width of the anode (end plate) is = 26.5 mm and the dimensions of the bolt head depend on the size of the bolt (Table 3).
2.3. Study Parameters
The analysis of galvanic corrosion is parametric, taking into account the effect of the bolt material, the bolt size, the electrolyte thickness, and the immersion time on the corrosion rate of the end plate, Table 4 and Figure 3.
Three types of low alloys were used for the bolt (20MnCr5, 42CrMo4, and 32CrMoV13) and will represent the cathode. For the bolt size three diameters were taken in this study (M12, M16, and M20); the dimensions that will be taken into consideration during the galvanic analysis are the dimensions of the bolt head as it is the part that will be in contact with the end plate. Galvanic corrosion will be studied in two thicknesses of the 1 M HCl electrolyte ( = 1 mm, = 20 mm) and for the following period range: 1 month, 6 months, and 12 months.
3. Method and Boundaries Conditions
3.1. Analysis of the Electrode Potential Distribution
The present analysis is based on a model galvanic corrosion couple consisting of two elements shown in Figure 2. the couple consists of a cathodic element (bolt) of a width “” in the + direction and a height “” in + direction; the width of the anode (plate) is “” in the − direction. The couple is covered by an electrolyte of depth “” in the + direction. The electrolyte is bounded by perfect insulators at , , and = +
The various current versus electrode potential relationships for the anodic and cathodic reactions are assumed to be subject to activation control with logarithmic (Tafel) polarization behavior. Thus the net cathodic current density per unit length of corrosion couple on the bolt at an electrode potential iswhere is the free corrosion potential of element , is the free corrosion current density, and and are the Tafel parameters for the anodic and cathodic reactions, respectively. Similarly for the end plate, the net anodic current density per unit length of corrosion couple on end plate at an electrode potential iswhere is the free corrosion potential of end plate, is the free corrosion current density, and and are the Tafel parameters for the anodic and cathodic reactions, respectively.
The equation governing the potential distribution in the electrolyte becomes Laplace’s equation (for constant conductivity): The governing equation for the potential must be solved using software based on finite element method (FEM); this equation must be subjected to appropriate boundary conditions. At any electrode-electrolyte interface, both anodic and cathodic processes take place simultaneously and the current density through the interface is the result of the electronic exchanges in both processes. In a galvanic couple, the cathodic process dominates in the more noble metal (the bolt) while the anodic process is dominant in the less noble metal (the end plate).
3.2. Boundary Conditions
For = 0 the normal derivative of electrode potential is given bywhere is the conductivity of the corrosion electrolyte.
Since there is no current flow normal to the insulating boundaries,The numerical solution of the Laplace equation requires the definition of the electrochemical parameters of the anode, the cathode, and the electrolyte in a 1 M HCl electrolyte whose characteristics are cited in the Table 5 . These parameters are determined from experimental tests already carried out in the literature  and they are cited in Table 6.
4. Results and Discussion
4.1. Electrolyte Potential
Once the electrochemical parameters were defined, the geometry was applied, the boundary conditions and governing equation were applied, and the appropriate mesh was found for the geometries, the model for galvanic corrosion was solved. The resolution of Laplace’s equation enables us to present electrolyte potential on the surfaces of electrodes bolt/end plate. Figures 5, 6, and 7 show us the surface plot of electrolyte potential with deformed geometry at different bolt sizes and different materials of bolt. It can be seen that the potential of the electrolyte increases if the bolt size increases: (M20) > (M16) > (M12). These potential values also depend on the bolt material: (32CrMoV13) > (42CrMo4) > (20MnCr5).
The potential of the electrolyte also depends on the thickness of the electrolyte; it increases if the thickness of the electrolyte increases, and this applies to the two parameters: bolt size and bolt material, ( = 20 mm) > ( = 1 mm), Figures 8, 9, and 10.
4.2. Corrosion Rate
The corrosion rate was calculated using the following relation : where is the corrosion rate (cm/yr), is the corrosion current density (A/cm2), is the molar mass (Fe) = 55.85 g/mol, is the valence of iron , is the faraday constant = 96500 A·s/mol, and is the density of steel .
Table 7 and Figures 11 and 12 show that the corrosion rate depends on the bolt size, the bolt material, and the electrolyte thickness. It can be seen that the corrosion rate of the end plate increases if the bolt size increases: (M20) > (M16) > (M12). These corrosion rate values depend also on the bolt material: (32CrMoV13) > (42CrMo4) > (20MnCr5).
The corrosion rate of the end plate depends also on the thickness of the electrolyte; it increases if the thickness of the electrolyte increases, and this applies to the two parameters: bolt size and bolt material, 20 mm) > 1 mm).
The choice of bolts in galvanic bolted joints of steel structures must take into account the effect of the bolt size and bolt material. It is shown in this study that as the size of the bolt increases, the corrosion rate increases. Also if the steel alloys of the bolt change, the corrosion rate changes. The increasing of the corrosion rate in the process of galvanic corrosion of the end plate leads to a reduction in lifetime of the bolted assembly. Hence a great interest must be shown by engineers/designers in the choice of the bolts when designing the steel structures to ensure long lifetime.
Conflicts of Interest
The authors declare no conflicts of interest.
- J. W. Oldfield, “Electrochemical theory of galvanic corrosion,” in Galvanic Corrosion, H. P. Hack, Ed., pp. 5–22, American Society for Testing and Materials, Philadelphia, Pa, USA, 1988, ASTM STP 978.
- H. P. Hack and J. R. Scully, “Galvanic corrosion prediction using long- and short-term polarization curves,” Corrosion, vol. 42, no. 2, pp. 79–90, 1986.
- J. F. Yan, S. N. Pakalapati, T. V. Nguyen, and R. E. White, “Mathematical modeling of cathodic protection using the boundary element method with a nonlinear polarization curve,” Journal of the Electrochemical Society, vol. 139, no. 7, pp. 1932–1936, 1992.
- L. Lacroix, C. Blanc, N. Pébère, B. Tribollet, and V. Vivier, “Localized approach to galvanic coupling in an aluminum-magnesium system,” Journal of The Electrochemical Society, vol. 156, no. 8, pp. C259–C265, 2009.
- K. B. Deshpande, “Validated numerical modelling of galvanic corrosion for couples: Magnesium alloy (AE44)-mild steel and AE44-aluminium alloy (AA6063) in brine solution,” Corrosion Science, vol. 52, no. 10, pp. 3514–3522, 2010.
- K. B. Deshpande, “Experimental investigation of galvanic corrosion: Comparison between SVET and immersion techniques,” Corrosion Science, vol. 52, no. 9, pp. 2819–2826, 2010.
- K. B. Deshpande, “Numerical modeling of micro-galvanic corrosion,” Electrochimica Acta, vol. 56, no. 4, pp. 1737–1745, 2011.
- J. X. Jia, G. Song, and A. Atrens, “Experimental measurement and computer simulation of galvanic corrosion of magnesium coupled to steel,” Advanced Engineering Materials, vol. 9, no. 1-2, pp. 65–74, 2007.
- D. Trinh, P. Dauphin Ducharme, U. Mengesha Tefashe, J. R. Kish, and J. Mauzeroll, “Influence of edge effects on local corrosion rate of magnesium alloy/mild steel galvanic couple,” Analytical Chemistry, vol. 84, no. 22, pp. 9899–9906, 2012.
- N. Murer, N. Missert, and R. Buchheit, Towards the Modeling of Microgalvanic Corrosion in Aluminum Alloys: the Choice of Boundary Conditions, COMSOL Users Conference, 2008.
- N. Murer, R. Oltra, B. Vuillemin, and O. Néel, “Numerical modelling of the galvanic coupling in aluminium alloys: A discussion on the application of local probe techniques,” Corrosion Science, vol. 52, no. 1, pp. 130–139, 2010.
- N. Murer, N. A. Missert, and R. G. Buchheit, “Finite element modeling of the galvanic corrosion of aluminum at engineered copper particles,” Journal of The Electrochemical Society, vol. 159, no. 6, pp. C265–C276, 2012.
- Y. Shi and R. G. Kelly, ECS Transactions, 41, 155, 2012.
- W. Sun, G. Liu, L. Wang, T. Wu, and Y. Liu, “An arbitrary Lagrangian-Eulerian model for studying the influences of corrosion product deposition on bimetallic corrosion,” Journal of Solid State Electrochemistry, vol. 17, no. 3, pp. 829–840, 2013.
- J. F. Yan, T. V. Nguyen, R. E. White, and R. B. Griffin, “Mathematical modeling of the formation of calcareous deposits on cathodically protected steel in seawater,” Journal of The Electrochemical Society, vol. 140, 733 pages, 1993.
- V. G. Artemov, A. A. Volkov, N. N. Sysoev, and A. A. Volkov, “Conductivity of aqueous HCl, NaOH and NaCl solutions: Is water just a substrate?” EPL (Europhysics Letters), vol. 109, no. 2, 2015.
- Okba Belahssen. Thesis: Etude comparative du comportement chimique et tribologique des aciers nitrurés.
- D. Höchez, “Simulation of corrosion product deposit layer growth on bare magnesium galvanically coupled to aluminum,” Journal of The Electrochemical Society, vol. 162, no. 1, pp. C1–C11, 2015.
- K. E. Johnson and J. S. Abbott, “Bimetallic corrosion effects on mild steel in natural environments,” British Corrosion Journal, vol. 9, no. 3, pp. 171–176, 1974.
- Q. Xu, K. Gao, W. Lv, and X. Pang, “Effects of alloyed Cr and Cu on the corrosion behavior of low-alloy steel in a simulated groundwater solution,” Corrosion Science, vol. 102, pp. 114–124, 2016.
- J. Scully and H. P. Hack, ASTM STP 978 Galvanic Corrosion, 136, 1988.
- C. Wagner, “Theoretical analysis of the current density distribution in electrolytic cells,” Journal of The Electrochemical Society, vol. 98, no. 3, pp. 116–128, 1951.
- F. J. Presuel-Moreno, H. Wang, M. A. Jakab, R. G. Kelly, and J. R. Scully, “Computational modeling of active corrosion inhibitor release from an Al-Co-Ce metallic coating,” Journal of The Electrochemical Society, vol. 153, no. 11, Article ID 002611JES, pp. B486–B498, 2006.
- J. S. Lee, M. L. Reed, and R. G. Kelly, “Combining rigorously controlled crevice geometry and computational modeling for study of crevice corrosion scaling factors,” Journal of The Electrochemical Society, vol. 151, no. 7, pp. B423–B433, 2004.
- E. Kennard and J. T. Waber, “Mathematical Study of Galvanic Corrosion: Equal Coplanar Anode and Cathode with Unequal Polarization Parameters,” Journal of The Electrochemical Society, vol. 117, no. 7, pp. 880–885, 1970.
- J. T. Waber and B. Fagan, “Mathematical Studies on Galvanic Corrosion IV. Influence of Electrolyte Thickness on the Potential and Current Distributions over Coplanar Electrodes Using Polarization Parameters,” Journal of The Electrochemical Society, vol. 103, no. 1, pp. 64–72, 1956.
- J. W. Fu, “Technical note: a finite element analysis of corrosion cells,” Corrosion, vol. 38, no. 5, pp. 295-296, 1982.
- J. T. Waber, “Mathematical Studies of Galvanic Corrosion VI. Limiting Case of Very Thin Films,” Journal of The Electrochemical Society, vol. 103, no. 10, pp. 567–570, 1956.
- P. Doig and P. E. J. Flewitt, “A finite difference numerical analysis of galvanic corrosion for semi-infinite linear coplanar electrodes,” Journal of The Electrochemical Society, vol. 126, no. 12, pp. 2057–2063, 1979.
- E. McCafferty, “Distribution of potential and current in circular corrosion cells having unequal polarization parameters,” Journal of The Electrochemical Society, vol. 124, no. 12, pp. 1869–1878, 1977.
- A. C. West, S. Mayer, and J. Reid, “A superfilling model that predicts bump formation,” Electrochemical and Solid-State Letters, vol. 4, no. 7, pp. C50–C53, 2001.
- N. Perez, Electrochemistry and Corrosion Science, Kluwer Academic Publishers, New York, NY, USA, 2004.
Copyright © 2017 Rachid Radouani et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.