On the Computational Solution of a Corrosion Model

We consider a phase field model for studying the corrosion of a 304 stainless steel metal immersed in a sodium chloride solution. Phase field models have been widely used to simulate moving interface problems in a variety of contexts. On one hand, these models have the benefit of treating implicitly the moving boundary by introducing an auxiliary variable. On the other hand, their nature is highly stiff and their computational cost is very high. In this talk we consider convenient numerical techniques for the efficient solution of the considered problem with a focus on their Matlab implementation.