A new Phase-Field Model for Anisotropic Surface Diffusion: Anisotropic Cahn-Hilliard Equation with Improved Conservation (ACH-IC)

TL;DR

This paper introduces the ACH-IC model, a novel anisotropic Cahn-Hilliard equation that significantly improves volume conservation and convergence in simulating anisotropic surface diffusion, validated through asymptotic analysis and numerical experiments.

Contribution

The paper develops a new variational anisotropic Cahn-Hilliard model with enhanced volume conservation and second-order accuracy, advancing phase-field modeling of surface diffusion.

Findings

Achieves second-order volume conservation.

Demonstrates improved convergence to sharp-interface limit.

Shows superior performance in numerical simulations.

Abstract

As popular approximations to sharp-interface models, the Cahn-Hilliard type phase-field models are usually used to simulate interface dynamics with volume conservation. However, the convergence rate of the volume enclosed by the interface to its sharp-interface limit is usually at first order of the interface thickness in the classical Cahn-Hilliard model with constant or degenerate mobilities. In this work, we propose a variational framework for developing new Cahn-Hilliard dynamics with enhanced volume conservation by introducing a more general conserved quantity. In particular, based on Onsager's variational principle (OVP) and a modified conservation law, we develop an anisotropic Cahn-Hilliard equation with improved conservation (ACH-IC) for approximating anisotropic surface diffusion. The ACH-IC model employs a new conserved quantity that approximates a step function more effectively, and yields second-order volume conservation while preserving energy dissipation for the classical anisotropic surface energy. The second-order volume conservation as well as the convergence to the sharp-interface surface diffusion dynamics is derived through comprehensive asymptotic analysis. Numerical evidence not only reveals the underlying physics of the proposed model in comparison with the classical one, but also demonstrates its exceptional performance in simulating anisotropic surface diffusion dynamics.