Decoupled, energy stable schemes for a phase-field surfactant model

被引:58
作者
Zhu, Guangpu [1 ]
Kou, Jisheng [2 ]
Sun, Shuyu [3 ]
Yao, Jun [1 ]
Li, Aifen [1 ]
机构
[1] China Univ Petr East China, Sch Petr Engn, Res Ctr Multiphase Flow Porous Media, Qingdao 266580, Peoples R China
[2] Hubei Engn Univ, Sch Math & Stat, Xiaogan 432000, Hubei, Peoples R China
[3] King Abdullah Univ Sci & Technol, Div Phys Sci & Engn, Computat Transport Phenomena Lab, Thuwal 239556900, Saudi Arabia
关键词
Phase-field; Cahn-Hilliard; Fluid-surfactant; Scalar auxiliary variable; Energy stability; DIFFUSE-INTERFACE METHOD; MOVING CONTACT LINE; 2-PHASE FLOWS; NUMERICAL APPROXIMATIONS; LATTICE-BOLTZMANN; SIMULATION; CAHN; STABILITY; DYNAMICS; TRACKING;
D O I
10.1016/j.cpc.2018.07.003
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
In this paper, we construct several linear, decoupled and energy stable schemes for a phase-field surfactant model, in which the free energy functional contains a fourth-order Ginzburg-Landau double well potential, a logarithmic Flory-Huggins potential and two nonlinear coupling terms. Several scalar auxiliary variables (SAV) are introduced to transform the governing system into an equivalent form, allowing the nonlinear potentials to be treated efficiently and semi-explicitly. At each time step, the schemes involve solving only two linear elliptic differential equations, and computations of two phase-field variables are totally decoupled. Moreover, the local concentration of surfactants can be obtained in an "explicit" way. We further establish a rigorous proof of unconditional energy stability for the semi-implicit schemes. Numerical results in both two and three dimensions are obtained, which demonstrate that the proposed schemes are accurate, efficient, easy-to-implement and unconditionally energy stable. (C) 2018 Elsevier B.V. All rights reserved.
引用
收藏
页码:67 / 77
页数:11
相关论文
共 48 条
[1]   Benchmark computations of diffuse interface models for two-dimensional bubble dynamics [J].
Aland, S. ;
Voigt, A. .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, 2012, 69 (03) :747-761
[2]  
[Anonymous], ARXIV13032559
[3]  
[Anonymous], 2017, J LIGHTWAVE TECHNOL
[4]   A hybrid numerical method for interfacial fluid flow with soluble surfactant [J].
Booty, M. R. ;
Siegel, M. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2010, 229 (10) :3864-3883
[5]   Applications of semi-implicit Fourier-spectral method to phase field equations [J].
Chen, LQ ;
Shen, J .
COMPUTER PHYSICS COMMUNICATIONS, 1998, 108 (2-3) :147-158
[6]   NUMERICAL-ANALYSIS OF THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC FREE-ENERGY [J].
COPETTI, MIM ;
ELLIOTT, CM .
NUMERISCHE MATHEMATIK, 1992, 63 (01) :39-65
[7]   On Diffuse Interface Modeling and Simulation of Surfactants in Two-Phase Fluid Flow [J].
Engblom, Stefan ;
Do-Quang, Minh ;
Amberg, Gustav ;
Tornberg, Anna-Karin .
COMMUNICATIONS IN COMPUTATIONAL PHYSICS, 2013, 14 (04) :879-915
[8]   Surfactants in foam stability: A phase-field model [J].
Fonseca, Irene ;
Morini, Massimiliano ;
Slastikov, Valeriy .
ARCHIVE FOR RATIONAL MECHANICS AND ANALYSIS, 2007, 183 (03) :411-456
[9]   An efficient scheme for a phase field model for the moving contact line problem with variable density and viscosity [J].
Gao, Min ;
Wang, Xiao-Ping .
JOURNAL OF COMPUTATIONAL PHYSICS, 2014, 272 :704-718
[10]   Transient coarsening behaviour in the Cahn-Hilliard model [J].
Garcke, H ;
Niethammer, B ;
Rumpf, M ;
Weikard, U .
ACTA MATERIALIA, 2003, 51 (10) :2823-2830