Adaptive finite element methods for Cahn-Hilliard equations

  • Authors:
  • L'ubomír Baňas;Robert Nürnberg

  • Affiliations:
  • Department of Mathematics, Imperial College, London SW7 2AZ, UK;Department of Mathematics, Imperial College, London SW7 2AZ, UK

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2008

Quantified Score

Hi-index 7.29

Visualization

Abstract

We develop a method for adaptive mesh refinement for steady state problems that arise in the numerical solution of Cahn-Hilliard equations with an obstacle free energy. The problem is discretized in time by the backward-Euler method and in space by linear finite elements. The adaptive mesh refinement is performed using residual based a posteriori estimates; the time step is adapted using a heuristic criterion. We describe the space-time adaptive algorithm and present numerical experiments in two and three space dimensions that demonstrate the usefulness of our approach.