A finite difference method for a moving-interface diffusion-reaction problem

Abstract

A software is described which solves a diffusion-reaction problem with N diffusing elements which can react to form M non-diffusing compounds. The numerical solution of such a problem is complicated by the fact that different systems of partial differential equations hold in regions in which different compounds are present, and by the fact that the interfaces between regions are mobile. This software is applied to a problem involving the diffusion of carbon and chromium in a NiCr alloy, in which multiple chromium carbides may form and dissociate.