A discontinuous Galerkin method for the Cahn-Hilliard equation

  • Authors:

  • Garth N. Wells;Ellen Kuhl;Krishna Garikipati

  • Affiliations:

  • Faculty of Civil Engineering and Geoseiences, Delft University of Technology, Delft, Netherlands;Chair of Applied Mechanics, Technische Universität Kaiserslautern, Postfach, Kaiserslautern, Germany;Department of Mechnical Engineering and Program in Applied Physics, University of Michigan, Ann Arbor, Michigan

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2006

Quantified Score

Hi-index 31.50

Visualization

Abstract

A discontinuous Galerkin finite element method has been developed to treat the high-order spatial derivatives appearing in the Cahn-Hilliard equation. The Cahn-Hilliard equation is a fourth-order nonlinear parabolic partial differential equation, originally proposed to model phase segregation of binary alloys. The developed discontinuous Galerkin approach avoids the need for mixed finite element methods, coupled equations or interpolation functions with a high degree of continuity that have been employed in the literature to treat the fourth-order spatial derivatives. The variational formulation of the discontinuous Galerkin method, its implementation and numerical examples are presented. In this communication, it is also shown under what conditions the method is stable, and an error estimate in an energy-type norm is presented. The method is evaluated by comparison with a standard finite element treatment in which the Cahn-Hilliard equation is decomposed into two coupled partial differential equations.