Parallel simulation of the Poisson–Nernst–Planck corrosion model with an algebraic flux correction method

Abstract

In this paper, a chain of scientific tools is applied to provide a flexible, efficient, and extensible toolchain that is capable of addressing limitations in corrosion simulation. A previously developed computational model of galvanic corrosion for Magnesium–Aluminium couples is considered and the necessities on the choice of implementation strategy is presented. Within this context the suitability of a novel approach, based on a recent developed algebraic flux correction method, is discussed. This suitability is argued in two directions. First, the agreement of the finite element theory with the characteristics of the problem. Second, the advantage of the method and the benefits gained by utilizing the implementation strategy. It is shown that the composition of the tools, and the finite element method successfully captures the expected model results and pose a good trade-off between simplicity, flexibility, and efficiency for a wider range of models. • Simulating complex corrosion models yields many challenges in the implementation process. • A resource consuming part of corrosion simulation is to solve the Poisson–Nernst–Planck (PNP) model. • The choice of the finite element method and an implementation tool facilitates solving the PNP model. • The developed toolchain in this work tackles the simulation challenges and speeds it up by the employed implementation strategy.