A positivity-preserving and convergent numerical scheme for the binary fluid-surfactant system

Abstract

In this paper, we develop a first order (in time) numerical scheme for the binary fluid surfactant phase field model. The free energy contains a double-well potential, a nonlinear coupling entropy and a Flory-Huggins potential. The resulting coupled system consists of two Cahn-Hilliard type equations. This system is solved numerically by finite difference spatial approximation, in combination with convex splitting temporal discretization. We prove the proposed scheme is unique solvable, positivity-preserving and unconditionally energy stable. In addition, an optimal rate convergence analysis is provided for the proposed numerical scheme, which will be the first such result for the binary fluid-surfactant system. Newton iteration is used to solve the discrete system. Some numerical experiments are performed to validate the accuracy and energy stability of the proposed scheme.

References in corpus (6)

• Numerical Approximations for a three components Cahn-Hilliard phase-field Model based on the Invariant Energy Quadratization method
• Linear and Unconditionally Energy Stable Schemes for the binary Fluid-Surfactant Phase Field Model
• Efficient and accurate numerical schemes for a hydrodynamically coupled phase field diblock copolymer model
• A positivity-preserving, energy stable scheme for a Ternary Cahn-Hilliard system with the singular interfacial parameters
• A positivity-preserving second-order BDF scheme for the Cahn-Hilliard equation with variable interfacial parameters
• A Positive and Energy Stable Numerical Scheme for the Poisson-Nernst-Planck-Cahn-Hilliard Equations with Steric Interactions