Energy Stable Numerical Schemes for Ternary Cahn-Hilliard System

Abstract

We present and analyze a uniquely solvable and unconditionally energy stable numerical scheme for the ternary Cahn-Hilliard system, with a polynomial pattern nonlinear free energy expansion. One key difficulty is associated with presence of the three mass components, though a total mass constraint reduces this to two components. Another numerical challenge is to ensure the energy stability for the nonlinear energy functional in the mixed product form, which turns out to be non-convex, non-concave in the three-phase space. To overcome this subtle difficulty, we add a few auxiliary terms to make the combined energy functional convex in the three-phase space, and this, in turn, yields a convex-concave decomposition of the physical energy in the ternary system. Consequently, both the unique solvability and the unconditional energy stability of the proposed numerical scheme are established at a theoretical level. In addition, an optimal rate convergence analysis in the $$\ell ^\infty (0,T; H_N^{-1}) \cap \ell ^2 (0,T; H_N^1)$$ norm is provided, with Fourier pseudo-spectral discretization in space, which is the first such result in this field. To deal with the nonlinear implicit equations at each time step, we apply an efficient preconditioned steepest descent (PSD) algorithm. A second order accurate, modified BDF scheme is also discussed. A few numerical results are presented, which confirm the stability and accuracy of the proposed numerical scheme.