An extended finite-element model coupled with level set method for analysis of growth of corrosion pits in metallic structures

Mass balance between metal and electrolytic solution, separated by a moving interface, in stable pit growth results in a set of governing equations which are solved for concentration field and interface position (pit boundary evolution). The interface experiences a jump discontinuity in metal concentration. The extended finite-element model (XFEM) handles this jump discontinuity by using discontinuous-derivative enrichment formulation, eliminating the requirement of using front conforming mesh and re-meshing after each time step as in the conventional finite-element method. However, prior interface location is required so as to solve the governing equations for concentration field for which a numerical technique, the level set method, is used for tracking the interface explicitly and updating it over time. The level set method is chosen as it is independent of shape and location of the interface. Thus, a combined XFEM and level set method is developed in this paper. Numerical analysis for pitting corrosion of stainless steel 304 is presented. The above proposed model is validated by comparing the numerical results with experimental results, exact solutions and some other approximate solutions. An empirical model for pitting potential is also derived based on the finite-element results. Studies show that pitting profile depends on factors such as ion concentration, solution pH and temperature to a large extent. Studying the individual and combined effects of these factors on pitting potential is worth knowing, as pitting potential directly influences corrosion rate.

[1]  Gad Frankel,et al.  Pitting Corrosion of Metals. A Review of the Critical Factors , 1998 .

[2]  G. Frankel Pitting Corrosion of Metals A Review of the Critical Factors , 1998 .

[3]  Frédéric Gibou,et al.  Second-Order Accurate Computation of Curvatures in a Level Set Framework Using Novel High-Order Reinitialization Schemes , 2008, J. Sci. Comput..

[4]  A. D. Bisgard,et al.  Determining the potential independent critical pitting temperature (CPT) by a potentiostatic method using the Avesta Cell , 1996 .

[5]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[6]  Petrus Christiaan Pistorius,et al.  Growth of corrosion pits on stainless steel in chloride solution containing dilute sulphate , 1992 .

[7]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[8]  A. Atrens Environmental Conditions Leading to Pitting/Crevice Corrosion of a Typical 12% Chromium Stainless Steel at 80 C , 1983 .

[9]  S. Osher,et al.  High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations , 1990 .

[10]  J. Sethian,et al.  Numerical Schemes for the Hamilton-Jacobi and Level Set Equations on Triangulated Domains , 1998 .

[11]  Olivier Coussy,et al.  Propagation Fronts during Calcium Leaching and Chloride Penetration , 2000 .

[12]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[13]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[14]  Chi-Wang Shu Total-variation-diminishing time discretizations , 1988 .

[15]  R. Newman,et al.  Perforated Covers for Propagating Pits , 1998 .

[16]  Stefan Scheiner,et al.  Stable pitting corrosion of stainless steel as diffusion-controlled dissolution process with a sharp moving electrode boundary , 2007 .

[17]  S. Osher,et al.  A Level Set Formulation of Eulerian Interface Capturing Methods for Incompressible Fluid Flows , 1996 .

[18]  D. Landolt,et al.  ROTATING DISC ELECTRODE STUDY OF ANODIC DISSOLUTION OR IRON IN CONCENTRATED CHLORIDE MEDIA , 1975 .

[19]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

[20]  Ted Belytschko,et al.  Discontinuous enrichment in finite elements with a partition of unity method , 2000 .

[21]  T. Belytschko,et al.  On the construction of blending elements for local partition of unity enriched finite elements , 2003 .

[22]  Wang Hai-bing,et al.  High-order essentially non-oscillatory schemes for Hamilton-Jacobi equations , 2006 .

[23]  J. Galvele Transport processes in passivity breakdown—II. Full hydrolysis of the metal ions , 1981 .

[24]  T. Belytschko,et al.  The extended finite element method (XFEM) for solidification problems , 2002 .

[25]  J. Sethian Evolution, implementation, and application of level set and fast marching methods for advancing fronts , 2001 .

[26]  Tsuguo Suzuki,et al.  Critical Potential for Growth of Localized Corrosion of Stainless Steel in Chloride Media , 1972 .

[27]  R. Newman,et al.  Pit growth studies in stainless steel foils. II. Effect of temperature, chloride concentration and sulphate addition , 2002 .

[28]  D. Chopp,et al.  Modelling crack growth by level sets , 2013 .

[29]  S. Osher,et al.  A PDE-Based Fast Local Level Set Method 1 , 1998 .

[30]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[31]  S. Osher,et al.  Regular Article: A PDE-Based Fast Local Level Set Method , 1999 .

[32]  P. Ernst,et al.  Pit growth studies in stainless steel foils. I. Introduction and pit growth kinetics , 2002 .

[33]  Ted Belytschko,et al.  Arbitrary discontinuities in finite elements , 2001 .

[34]  John E. Dolbow,et al.  On strategies for enforcing interfacial constraints and evaluating jump conditions with the extended finite element method , 2004 .

[35]  C. Sukenik,et al.  Siloxane-based thin films for corrosion protection of stainless steel in chloride media , 2011 .

[36]  David L. Chopp,et al.  A hybrid extended finite element/level set method for modeling phase transformations , 2002 .

[37]  Petrus Christiaan Pistorius,et al.  Metastable pitting corrosion of stainless steel and the transition to stability , 1992, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

[38]  N. Sato The stability of localized corrosion , 1995 .

[39]  H. S. Isaacs,et al.  In Situ X‐Ray Microprobe Study of Salt Layers during Anodic Dissolution of Stainless Steel in Chloride Solution , 1995 .

[40]  R. Newman,et al.  Temperature dependence of pitting potentials for austenitic stainless steels above their critical pitting temperature , 1998 .

[41]  T. P. Hoar,et al.  The production and breakdown of the passivity of metals , 1967 .

[42]  S. Osher,et al.  A Simple Level Set Method for Solving Stefan Problems , 1997, Journal of Computational Physics.

[43]  Stéphane Bordas,et al.  Enriched finite elements and level sets for damage tolerance assessment of complex structures , 2006 .

[44]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[45]  B. Šarler Stefan's work on solid-liquid phase changes , 1995 .