Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
4963991 | Computer Methods in Applied Mechanics and Engineering | 2017 | 25 Pages |
Abstract
In this paper, we consider the numerical solution of a binary fluid-surfactant phase field model, in which the free energy contains a nonlinear coupling entropy, a Ginzburg-Landau double well potential, and a logarithmic Flory-Huggins potential. The resulting system consists of two nonlinearly coupled Cahn-Hilliard type equations. We develop a first and a second order time stepping schemes for this system using the “Invariant Energy Quadratization” approach; in particular, the system is transformed into an equivalent one by introducing appropriate auxiliary variables and all nonlinear terms are then treated semi-explicitly. Both schemes are linear and lead to symmetric positive definite systems in space at each time step, thus they can be efficiently solved. We further prove that these schemes are unconditionally energy stable in the discrete sense. Various 2D and 3D numerical experiments are performed to validate the accuracy and energy stability of the proposed schemes.
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science Applications
Authors
Xiaofeng Yang, Lili Ju,