Finite element calculation of the polarisation behaviour of a metal in an aqueous solution using the dilute solution model

This article describes the calculation of the polarisation behaviour of a m etal in an aqueous environment using a computational electrochemistry code based on finite elements. The corrosion system consisted of three parallel reactions (oxygen reduction, hydrogen evolution and metal dissolution), tak ing place in the centre at the surface of the rotating disc electrode (RDE) . The selected problem is simple enough to be solved analytically, which al lows verification of the finite element calculation. Simultaneously, the ex ample is rich enough to demonstrate the potential and the advantages of fin ite element simulation for corrosion science and engineering. The analytica l and finite element calculations were in excellent agreement. The influenc e of the speed of the RDE, of the bulk concentration of dissolved oxygen an d of the bulk hydrogen ion concentration on the polarisation behaviour have been determined. The use of finite elements had the following specific adv antages. Firstly, no prior assumption on the diffusion layer thickness had to be made. Secondly, transport by diffusion, convection and electro-migrat ion were taken into account. Thirdly, the finite element approach is genera l and flexible enough to tackle more complex systems. The calculations pres ented here are a first step towards the finite element simulation of stress corrosion cracking of stainless steel in an aqueous nuclear environment. ( C) 2000 Elsevier Science Ltd. All rights reserved.