Abstract
We consider several seconder order in time stabilized semi-implicit Fourier spectral schemes for 2D Cahn–Hilliard equations. We introduce new stabilization techniques and prove unconditional energy stability for modified energy functionals. We also carry out a comparative study of several classical stabilization schemes and identify the corresponding stability regions. In several cases the energy stability is proved under relaxed constraints on the size of the time steps. We do not impose any Lipschitz assumption on the nonlinearity. The error analysis is obtained under almost optimal regularity assumptions.
Similar content being viewed by others
Explore related subjects
Discover the latest articles and news from researchers in related subjects, suggested using machine learning.References
Bourgain, J., Li, D.: Strong ill-posedness of the incompressible Euler equation in borderline Sobolev spaces. Invent. Math. 201, 97–157 (2015)
Bourgain, J., Li, D.: Strong illposedness of the incompressible Euler equation in integer \(C^m\) spaces. Geom. Funct. Anal. 25, 1–86 (2015)
Bertozzi, A., Ju, N., Lu, H.: A biharmonic-modified forward time stepping method for fourth order nonlinear diffusion equations. Discrete Contin. Dyn. Syst. 29, 1367–1391 (2011)
Cahn, J.W., Hilliard, J.E.: Free energy of a nonuniform system. I. Interfacial energy free energy. J. Chem. Phys. 28, 258–267 (1958)
Chen, L.Q., Shen, J.: Applications of semi-implicit Fourier-spectral method to phase field equations. Comput. Phys. Commun. 108, 147–158 (1998)
Chen, F., Shen, J.: Efficient energy stable schemes with spectral discretization in space for anisotropic Cahn-Hilliard systems. Commun. Comput. Phys. 13, 1189–1208 (2013)
Cheng, K., Wang, C., Wise, S., Yue, X.: A second-order, weakly energy-stable pseudo-spectral scheme for the Cahn-Hilliard equation and its solution by the homogeneous linear iteration method. J. Sci. Comput. (2016). doi:10.1007/s10915-016-0228-3
Christlieb, A., Jones, J., Promislow, K., Wetton, B., Willoughby, M.: High accuracy solutions to energy gradient flows from material science models. J. Comput. Phys. 257, 193–215 (2014)
Feng, W.M., Yu, P., Hu, S.Y., Liu, Z.K., Du, Q., Chen, L.Q.: A Fourier spectral moving mesh method for the Cahn–Hilliard equation with elasticity. Commun. Comput. Phys. 5, 582–599 (2009)
Feng, X.B., Prohl, A.: Error analysis of a mixed finite element method for the Cahn-Hilliard equation. Numer. Math. 99, 47–84 (2004)
Feng, X., Tang, T., Yang, J.: Stabilized Crank-Nicolson/Adams-Bashforth schemes for phase field models. East Asian J. Appl. Math. 3, 59–80 (2013)
Guillen-Gonzalez, F., Tierra, G.: Second order schemes and time-step adaptivity for Allen–Cahn and Cahn–Hilliard models. Comput. Math. Appl. 68, 821–846 (2014)
Gavish, N., Jones, J., Xu, Z., Christlieb, A., Promislow, K.: Variational models of network formation and ion transport: applications to perfluorosulfonate ionomer membranes. Polymers 4, 630–655 (2012)
Gomez, H., Hughes, T.J.R.: Provably unconditionally stable, second-order time-accurate, mixed variational methods for phase-field models. J. Comput. Phys. 230, 5310–5327 (2011)
Gottlieb, S., Tone, F., Wang, C., Wang, X., Wirosoetisno, D.: Long time stability of a classical efficient scheme for two dimensional Navier–Stokes equations. SIAM J. Numer. Anal. 50, 126–150 (2012)
Gottlieb, S., Wang, C.: Stability and convergence analysis of fully discrete Fourier collocation spectral method for 3-D viscous Burgers’ equation. J. Sci. Comput. 53, 102–128 (2012)
Guo, J., Wang, C., Wise, S., Yue, X.: An \(H^2\) convergence of a second-order convex-splitting, finite difference scheme for the three-dimensional Cahn–Hilliard equation. Commun. Math. Sci. 14, 489–515 (2016)
He, Y., Liu, Y., Tang, T.: On large time-stepping methods for the Cahn–Hilliard equation. Appl. Numer. Math. 57, 616–628 (2007)
Ju, L., Zhang, J., Du, Q.: Fast and accurate algorithms for simulating coarsening dynamics of Cahn–Hilliard equations. Comput. Mater. Sci. 108, 272–282 (2015)
Li, B., Liu, J.G.: Thin film epitaxy with or without slope selection. Eur. J. Appl. Math. 14, 713–743 (2003)
Li, D.: On a frequency localized Bernstein inequality and some generalized Poincaré-type inequalities. Math. Res. Lett. 20, 933–945 (2013)
Li, D., Qiao, Z., Tang, T.: Characterizing the stabilization size for semi-implicit Fourier-spectral method to phase field equations. SIAM J. Numer. Anal. 54, 1653–1681 (2016)
Li, D., Qiao, Z., Tang, T.: Gradient bounds for a thin film epitaxy equation. J. Differ. Equ. (2016, submitted)
Li, D., Qiao, Z.: On the stabilization size of semi-implicit Fourier-spectral methods for 3D Cahn–Hilliard equations. Comm. Math. Sci. (2016, submitted)
Qiao, Z., Zhang, Z., Tang, T.: An adaptive time-stepping strategy for the molecular beam epitaxy models. SIAM J. Sci. Comput. 33, 1395–1414 (2011)
Sun, Z.Z.: A second-order accurate linearized difference scheme for the two-dimensional Cahn–Hilliard equation. Math. Comput. 64, 1463–1471 (1995)
Shen, J., Yang, X.: Numerical approximations of Allen–Cahn and Cahn–Hilliard equations. Discrete Contin. Dyn. Syst. A 28, 1669–1691 (2010)
Shen, J., Wang, C., Wang, X., Wise, S.M.: Second-order convex splitting schemes for gradient flows with Ehrlich–Schwoebel type energy: application to thin film epitaxy. SIAM J. Numer. Anal. 50, 105–125 (2012)
Shen, J., Yang, X.: Decoupled energy stable schemes for phase-field models of two-phase complex fluids. SIAM J. Sci. Comput. 36, B122–B145 (2014)
Shen, J., Yang, X.: Decoupled, energy stable schemes for phase-field models of two-phase incompressible flows. SIAM J. Numer. Anal. 53, 279–296 (2015)
Wang, C., Wang, S., Wise, S.M.: Unconditionally stable schemes for equations of thin film epitaxy. Discrete Contin. Dyn. Sys. Ser. A 28, 405–423 (2010)
Xu, C., Tang, T.: Stability analysis of large time-stepping methods for epitaxial growth models. SIAM J. Numer. Anal. 44, 1759–1779 (2006)
Zhu, J., Chen, L.-Q., Shen, J., Tikare, V.: Coarsening kinetics from a variable-mobility Cahn–Hilliard equation: application of a semi-implicit Fourier spectral method. Phys. Rev. E 60(3), 3564–3572 (1999)
Acknowledgments
D. Li was supported by an Nserc discovery grant. The research of Z. Qiao is partially supported by the Hong Kong Research Council GRF Grant 15302214, NSFC/RGC Joint Research Scheme N_HKBU204/12 and the Hong Kong Polytechnic University internal Grant 1-ZE33.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Li, D., Qiao, Z. On Second Order Semi-implicit Fourier Spectral Methods for 2D Cahn–Hilliard Equations. J Sci Comput 70, 301–341 (2017). https://doi.org/10.1007/s10915-016-0251-4
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-016-0251-4
Keywords
Mathematics Subject Classification
Profiles
- Zhonghua Qiao View author profile