Numerical Approximations for the Cahn-Hilliard phase field model of the binary fluid-surfactant system

TL;DR

This paper develops and analyzes linear, decoupled time-stepping schemes for the binary fluid-surfactant phase field model, ensuring unconditional energy stability and accuracy through the Invariant Energy Quadratization approach.

Contribution

It introduces new linear, decoupled, and energy-stable time-stepping schemes for the coupled Cahn-Hilliard equations using the IEQ method and explicit-implicit techniques.

Findings

First order scheme is unconditionally energy stable.

Numerical simulations confirm stability and accuracy.

Linear system is well-posed with symmetric positive definite operator.

Abstract

In this paper, we consider the numerical approximations for the commonly used binary fluid-surfactant phase field model that consists two nonlinearly coupled Cahn-Hilliard equations. The main challenge in solving the system numerically is how to develop easy-to-implement time stepping schemes while preserving the unconditional energy stability. We solve this issue by developing two linear and decoupled, first order and a second order time-stepping schemes using the so-called "Invariant Energy Quadratization" approach for the double well potentials and a subtle explicit-implicit technique for the nonlinear coupling potential. Moreover, the resulting linear system is well-posed and the linear operator is symmetric positive definite. We rigorously prove the first order scheme is unconditionally energy stable. Various numerical simulations are presented to demonstrate the stability and the accuracy thereafter.

Full text

See pages 1-last of Surfactant_Jan_1.pdf