Finite Element Modeling of the Galvanic Corrosion of Aluminum at Engineered Copper Particles

Abstract

Finite element modeling based on solving the Nernst-Planck equation was used to describe the evolution of current densities and pH distribution at the surface of a bimetallic system. This system consisted of five single copper particles fabricated on isolated, thin-film aluminum electrodes exposed to dilute aqueous chloride solutions. Excess anodic and cathodic currents flowing during exposure were used to validate the model. The corrosion of the bimetallic system exhibited a passive and an active stage. The model was used firstly to quantify the influence of pH (modified by alkalization by O 2 reduction at the cathodes and acidification by Al 3+ hydrolysis) on the passive-active transition and secondly to verify the anodic control of the active corrosion. EQCM (Electrochemical Quartz Crystal Microbalance) was used to obtain an experimental relationship between Al dissolution rate and pH. In the conditions of the model, Al 3+ hydrolysis was found to activate Al dissolution in the early instants of exposure (≈10 s) whereas OH − would have an effect only at longer times (hundreds of s). The pH-related destabilization appeared not to be sufficient enough to trigger active dissolution and it was proved that this stage was controlled by the anodic activity and not O 2 diffusion.

A Review of the Application of Finite Element Method (FEM) to Localized Corrosion Modeling

“…[90][91][92][93][94][95][96][97][98] The accepted mechanism for the location of the initiation sites is based on microgalvanic couples in which the alloy matrix can be either the anode or cathode with respect to the adjacent IMC. One can experimentally study the microgalvanic corrosion using engineered samples such as engineered copper galvanically coupled with Al at microscale 99 to systematically investigate electrochemical distributions and chemistry evolution over time at the coupling interface. Note that these previous studies nevertheless utilized the pure metal substitute (e.g., Cu) which may not truly reflect the electrochemical behavior of IMC during the corrosion.…”

Numerical Simulation of Galvanic Corrosion between Carbon Steel and Low Alloy Steel in a Bolted Joint

“…The publications of Murer et al [10][11][12] and Shi and Kelly [13] 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 [14], who applied the mathematical approach of Yan et al [15] to the modeling of deposit formation under seawater conditions, clearly introduce a possible way of a useful model built up for the mentioned purpose.…”

A method for determining the polarization state of the metal with the middle potential in order to calculate the corrosion rate of multi‐metals complicated galvanic couple

“…Firstly, it is to get the potential distribution by solving the Laplace equation; secondly, it is to obtain the current density distribution by the relationship between the potential and current density of polarization curves; thirdly, it is to calculate the galvanic corrosion rate by Faraday's law. According to the numerical method of solving Laplace equation, it develops three kinds of galvanic corrosion numerical simulation methods including finite difference method (FDM), finite element method (FEM), and boundary element method (BEM) . Doig et al evaluated the distribution of the potential within the electrolyte above a galvanic corrosion couple with semi‐infinite electrodes by FEM .…”