A numerical method was developed for treating the reaction-diffusion process in multi-component systems which exhibited limited solubility in their compounds. The method combined a description, of local thermodynamic equilibrium at the moving interfaces, with a phenomenological description of the diffusion process in each phase; taking their specific solubility features into account. There were no formal restrictions on the number of components or phases, and the algorithms could be used for the systematic study of the effect of variables such as free energies of formation and mobilities upon the stable diffusion path, as well as for the prediction of the actual diffusion path.

A Model for the Prediction of Reaction Diffusion Paths in Multicomponent Systems with Limited Solubility. B.Hallemans, P.Wollants, J.R.Roos, B.Blanpain: Materials Transactions, 2001, 42[12], 2630-6