Numerical schemes for a three component Cahn-Hilliard model

Franck Boyer (1) , Sebastian Minjeaud (2)

1 LATP - Laboratoire d'Analyse, Topologie, Probabilités
2 IRSN - Institut de Radioprotection et de Sûreté Nucléaire

In this article, we investigate numerical schemes for solving a three component Cahn-Hilliard model. The space discretization is performed by using a Galerkin formulation and the finite element method. Concerning the time discretization,the main difficulty is to write a scheme ensuring, at the discrete level, the decrease of the free energy and thus the stability of the method. We study three different schemes and prove existence and convergence theorems. Theoretical results are illustrated by various numerical examples showing that the new semi-implicit discretization that we propose seems to be a good compromise between robustness and accuracy.

Domaines: Analyse numérique [math.NA]

Franck Boyer, Sebastian Minjeaud. Numerical schemes for a three component Cahn-Hilliard model. ESAIM: Mathematical Modelling and Numerical Analysis, 2011, 45 (4), pp 697-738. ⟨10.1051/m2an/2010072⟩. ⟨hal-00390065v2⟩