Jump to ContentJump to Main Navigation
Show Summary Details
More options …

International Journal of Nonlinear Sciences and Numerical Simulation

Editor-in-Chief: Birnir, Björn

Editorial Board: Armbruster, Dieter / Chen, Xi / Bessaih, Hakima / Chou, Tom / Grauer, Rainer / Marzocchella, Antonio / Rangarajan, Govindan / Trivisa, Konstantina / Weikard, Rudi

8 Issues per year

IMPACT FACTOR 2017: 1.162

CiteScore 2017: 1.41

SCImago Journal Rank (SJR) 2017: 0.382
Source Normalized Impact per Paper (SNIP) 2017: 0.636

Mathematical Citation Quotient (MCQ) 2017: 0.12

See all formats and pricing
More options …
Volume 18, Issue 5


Energy Stable Interior Penalty Discontinuous Galerkin Finite Element Method for Cahn–Hilliard Equation

Ayşe Sarıaydın-Filibelioğlu
  • Corresponding author
  • Econometry Department, 100. Yıl University, Van, Turkey
  • Institute of Applied Mathematics, Middle East Technical University, Ankara, Turkey
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Bülent Karasözen
  • Institute of Applied Mathematics, Middle East Technical University, Ankara, Turkey
  • Department of Mathematics, Middle East Technical University, Ankara, Turkey
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ Murat Uzunca
Published Online: 2017-07-14 | DOI: https://doi.org/10.1515/ijnsns-2016-0024


An energy stable conservative method is developed for the Cahn–Hilliard (CH) equation with the degenerate mobility. The CH equation is discretized in space with the mass conserving symmetric interior penalty discontinuous Galerkin (SIPG) method. The resulting semi-discrete nonlinear system of ordinary differential equations are solved in time by the unconditionally energy stable average vector field (AVF) method. We prove that the AVF method preserves the energy decreasing property of the fully discretized CH equation. Numerical results for the quartic double-well and the logarithmic potential functions with constant and degenerate mobility confirm the theoretical convergence rates, accuracy and the performance of the proposed approach.

Keywords: Cahn–Hilliard equation; gradient systems; discontinuous Galerkin method; average vector field method

JEL Classification: 02.60.Cb; 02.60.Lj; 02.70.Dh


  • [1]

    Chen L. and Xu C., A time splitting space spectral element method for the Cahn–Hilliard equation, East Asian J. Appl. Math. 3 (4) (2013), 333–351.CrossrefGoogle Scholar

  • [2]

    Guo R. and Xu Y., Efficient solvers of discontinuous Galerkin discretization for the Cahn-Hilliard equations, J. Sci. Comput. 58 (2) (2014), 380–408.Web of ScienceCrossrefGoogle Scholar

  • [3]

    Cahn J. W. and Hilliard J. E., Free energy of a nonuniform system. I. Interfacial free energy, J. Chem. Phys. 28 (1958), 258–267.CrossrefGoogle Scholar

  • [4]

    Wu X., van Zwieten G. J. and van der Zee K. G., Stabilized second-order convex splitting schemes for Cahn–Hilliard models with application to diffuse-interface tumor-growth models, Int. J. Numer. Methods Biomed. Eng. 30 (2) (2014), 180–203.CrossrefGoogle Scholar

  • [5]

    Wells G.N., Kuhl E. and Garikipati K., A discontinuous Galerkin method for the Cahn–Hilliard equation, J. Comput. Phys. 218 (2006), 860–877.Google Scholar

  • [6]

    Barrett J. W., J. Blowey F. and Garcke H., Finite element approximation of the Cahn–Hilliard equation with degenerate mobility, SIAM J. Numer. Anal. 37 (1) (2000), 286–318.CrossrefGoogle Scholar

  • [7]

    Bartels S. and Müller R., Error control for the approximation of Allen–Cahn and Cahn–Hilliard equations with a logarithmic potential, Numer. Math. 119 (2011), 409–435.Web of ScienceGoogle Scholar

  • [8]

    Xia Yinhua, Yan Xu and Shu Chi-Wang, Local discontinuous Galerkin methods for the Cahn–Hilliard type equations, J. Comput. Phys. 227 (2007), 472–491.Google Scholar

  • [9]

    van der Zee K.G., Oden T.J., Prudhomme S. and Hawkins-Daarud A., Goal-oriented error estimation for Cahn–Hilliard models of binary phase transition, Numer. Methods Partial Differ. Equ. 27 (2011), 160–196.CrossrefGoogle Scholar

  • [10]

    Elliott C. and Garcke H., On the Cahn–Hilliard equation with degenerate mobility, SIAM J. Math. Anal. 27 (2) (1996), 404–423.Google Scholar

  • [11]

    Christlieb A., J. Jones, K. Promislow, B. Wetton and Willoughby M., High accuracy solutions to energy gradient flows from material science models, J. Comput. Phys. 257 (2014), 193–215.Google Scholar

  • [12]

    Kim J., A numerical method for the Cahn–Hilliard equation with a variable mobility, Commun. Nonlinear Sci. Numer. Simul. 12 (2007), 1560–1571.CrossrefWeb of ScienceGoogle Scholar

  • [13]

    Zhu Jingzhi, Chen Long-Qing, Shen Jie and Tikare Veena, Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method, Phys. Rev. E 60 (1999), 3564–3572.Google Scholar

  • [14]

    Arnold D. N., An interior penalty finite element method with discontinuous elements, SIAM J. Numer. Anal. 19 (1982), 724–760.CrossrefGoogle Scholar

  • [15]

    Rivière B., Discontinuous Galerkin methods for solving elliptic and parabolic equations, Society for Industrial and Applied Mathematics, 2008.Google Scholar

  • [16]

    Feng X. and Li Y., Analysis of interior penalty discontinuous Galerkin methods for the Allen–Cahn equation and the mean curvature flow, arXiv:1310.7504v2 [math.NA (2014).Google Scholar

  • [17]

    Tierra G. and Guillén-González F., Numerical methods for solving the Cahn–Hilliard equation and its applicability to related energy-based models, Arch. Comput. Methods Eng. 22 (2015), 1–21.CrossrefWeb of ScienceGoogle Scholar

  • [18]

    E.i Hairer and C. Lubich, Energy-diminishing integration of gradient systems, IMA J. Numer. Anal. 34 (2014), 452–461.Crossref

  • [19]

    Celledoni E., Grimm V., McLachlan R.I., McLaren D.I., O’Neale D., Owren B. and Quispel G. R. W., Preserving energy resp. dissipation in numerical {PDEs} using the “Average Vector Field” method, J. Comput. Phys. 231 (2012), 6770–6789.Web of ScienceGoogle Scholar

  • [20]

    Vemaganti K., Discontinuous Galerkin methods for periodic boundary value problems, Numer. Methods Partial Differ. Equ. 23 (2007), 587–596.CrossrefGoogle Scholar

  • [21]

    Hairer E., Energy-preserving variant of collocation methods, J. Numer. Anal. Ind. Appl. Math. 5 (2010), 73–84.Google Scholar

  • [22]

    Estep D. J. and Freund R. W., Using Krylov-subspace iterations in discontinuous Galerkin methods for nonlinear reaction-diffusion systems, Discontinuous Galerkin methods (Newport, RI, 1999), Lect. Notes Comput. Sci. Eng. 11, Springer, 2000, pp. 327–335.CrossrefGoogle Scholar

  • [23]

    Barrett J. W. and Blowey J. F., Finite element approximation of the Cahn–Hilliard equation with concentration dependent mobility, Math. Comput. 68 (1999), 487–517.CrossrefGoogle Scholar

  • [24]

    Hairer E. and Wanner G., Solving ordinary differential equations II: Stiff and differential–algebraic problems, Springer Series in Computational Mathematics: Berlin, 1996.Google Scholar

  • [25]

    Gomez H., Reali A. and Sangalli G., Accurate, efficient, and (iso)geometrically flexible collocation methods for phase-field models, J. Comput. Phys. 262 (2014), 153–171.Web of ScienceGoogle Scholar

About the article

Received: 2016-02-11

Accepted: 2017-05-04

Published Online: 2017-07-14

Published in Print: 2017-07-26

Citation Information: International Journal of Nonlinear Sciences and Numerical Simulation, Volume 18, Issue 5, Pages 303–314, ISSN (Online) 2191-0294, ISSN (Print) 1565-1339, DOI: https://doi.org/10.1515/ijnsns-2016-0024.

Export Citation

© 2017 Walter de Gruyter GmbH, Berlin/Boston.Get Permission

Citing Articles

Here you can find all Crossref-listed publications in which this article is cited. If you would like to receive automatic email messages as soon as this article is cited in other publications, simply activate the “Citation Alert” on the top of this page.

Tuğba Akman Yıldız, Murat Uzunca, and Bülent Karasözen
Applied Mathematics and Computation, 2019, Volume 347, Page 194

Comments (0)

Please log in or register to comment.
Log in