For the double-well potential case \(f(u)\) , we set \(\kappa =2\)  and
the phase structures captured at some moments
are presented in Figure REF , and the constant steady state \(u\equiv -1\)  is reached at around \(t=604\) .
The left picture given in Figure REF  implies
the preservation of the MBP during the whole phase transition process.
The energy evolution is plotted in the right graph of Figure REF ,
which states the energy dissipation of the process. For the Flory–Huggins potential case, we set \(\kappa =8.02\)  and
the simulated results are shown in Figures REF  and REF . We observe that
the steady state is reached at around \(t=602\)  and
the whole process of phase separation is similar to that of the double-well potential case. Those results are almost identical to
those produced using the IFRK4 scheme in [1]}.
<FIGURE><FIGURE><FIGURE><FIGURE>