LAUSR

We develop a second-order accurate, energy stable, and linear numerical method for a ternary Cahn–Hilliard (CH) model. The proposed scheme is an extension of typical Lagrange multiplier approach for binary CH system. The second-order backward difference formula (BDF2) is applied to construct time discretization. We theoretically prove the mass conservation, unique solvability, and energy stability of the proposed scheme. We efficiently solve the resulting discrete linear system by using a multigrid algorithm. The numerical solutions demonstrate that the proposed scheme is practically stable and second-order accurate in time and space. Moreover, we can use the proposed scheme as an effective solver to calculate the ternary CH equations in ternary phase-field fluid systems.