Grand-potential based multiphase-field model is extended to include surface diffusion. Diffusion is elevated in the interface through a scalar degenerate term. In contrast to the classical Cahn-Hilliard-based formulations, the present model circumvents the related difficulties in restricting diffusion solely to the interface by combining two second-order equations, an Allen-Cahn-type equation for the phase field supplemented with an obstacle-type potential and a conservative diffusion equation for the chemical potential or composition evolution. The sharp interface limiting behavior of the model is deduced by means of asymptotic analysis. A combination of surface diffusion and finite attachment kinetics is retrieved as the governing law. Infinite attachment kinetics can be achieved through a minor modification of the model, and with a slight change in the interpretation, the same model handles the cases of pure substances and alloys. Relations between model parameters and physical properties are obtained which allow one to quantitatively interpret simulation results. An extensive study of thermal grooving is conducted to validate the model based on existing theories. The results show good agreement with the theoretical sharp-interface solutions. The obviation of fourth-order derivatives and the usage of the obstacle potential make the model computationally cost-effective.

Owing to its unique thermodynamical description, phase separation has largely been modeled in the Cahn-Hilliard framework. In the present work, as a computationally efficient alternative, a multicomponent, multiphase-field model operating in Allen-Cahn framework is presented and subsequently employed to simulate spinodal decomposition. Stability analysis shows that the formulation can cover the effect of phase separation while simplifying the extension to multiple phases and incorporation of additional driving forces. Computational efficiency of the proposed approach is compared with the conventional technique by modeling intercalation in a representative one-dimensional domain. Moreover, intercalation within a multigrain system involving multiparticle interaction is studied. Our results suggest initiation of phase transformation at higher order junctions as well as a grain-by-grain intercalation behavior in a two-phase cathode material such as the well-studied ${\text{LiFePO}}_4$.

Existing grand-potential based multicomponent phase-field model is extended to handle systems with interstitial sublattice. This is achieved by treating the concentration of alloying elements in site-fraction. Correspondingly, the chemical species are distinguished based on their lattice positions, and their mode of diffusion, interstitial or substitutional, is appropriately realised. An approach to incorporate quantitative driving-force, through parabolic approximation of CALPHAD data, is introduced. By modelling austenite decomposition in ternary Fe–C–Mn, albeit in a representative microstructure, the ability of the current formalism to handle phases with interstitial components, and to distinguish interstitial diffusion from substitutional in grand-potential framework is elucidated. Furthermore, phase transformation under paraequilibrium is modelled to demonstrate the limitation of adopting mole-fraction based formulation to treat multicomponent systems.