Improving the accuracy of convexity splitting methods for gradient flow equations

Karl B Glasner, Saulo Orizaga

Research output: Contribution to journalArticle

17 Scopus citations

Abstract

This paper introduces numerical time discretization methods which significantly improve the accuracy of the convexity-splitting approach of Eyre (1998) [7], while retaining the same numerical cost and stability properties.A first order method is constructed by iteration of a semi-implicit method based upon decomposing the energy into convex and concave parts. A second order method is also presented based on backwards differentiation formulas. Several extrapolation procedures for iteration initialization are proposed. We show that, under broad circumstances, these methods have an energy decreasing property, leading to good numerical stability.The new schemes are tested using two evolution equations commonly used in materials science: the Cahn-Hilliard equation and the phase field crystal equation. We find that our methods can increase accuracy by many orders of magnitude in comparison to the original convexity-splitting algorithm. In addition, the optimal methods require little or no iteration, making their computation cost similar to the original algorithm.

Original languageEnglish (US)
Pages (from-to)52-64
Number of pages13
JournalJournal of Computational Physics
Volume315
DOIs
Publication statusPublished - Jun 15 2016

    Fingerprint

Keywords

  • Cahn-Hilliard
  • Eyre splitting
  • Gradient flows
  • Phase-field crystal
  • Spectral methods

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy (miscellaneous)

Cite this