An improved scalar auxiliary variable (SAV) approach for the phase-field surfactant model

Abstract In this work, we develop a new linear, decoupled numerical scheme for the typical phase-field surfactant model. An improved scalar auxiliary variable (SAV) approach is used to discretize the governing equations in time. Different from the classical SAV approach, this improved form can calculate the phase field function ϕ, surfactant function ψ, and auxiliary variables in a step-by-step manner, i.e., the auxiliary variables are treated totally explicitly, thus we can directly calculate ϕ and ψ instead of computing the inner products. At each time step, the surfactant ψ can be directly obtained by an explicit way, then ϕ is updated by solving a linear system with constant coefficient. Therefore, the implementation of this improved SAV approach is easier than the classical SAV approach. The energy stability in first-order case can be analytically proved by using the our method. The numerical experiments show that our proposed method not only achieves desired first- and second-order accuracy but also satisfies the desired discrete energy dissipation law even if some larger time steps are used. Furthermore, the coarsening dynamics with different average concentrations can be well simulated by using our method. The co-continuous and drop patterns are generated in the even compositions case and uneven compositions case, respectively.

[1]  Jian-Jun Xu,et al.  A level-set method for two-phase flows with moving contact line and insoluble surfactant , 2014, J. Comput. Phys..

[2]  Charles M. Elliott,et al.  Numerical analysis of the Cahn-Hilliard equation with a logarithmic free energy , 1992 .

[3]  Lili Ju,et al.  Linear and unconditionally energy stable schemes for the binary fluid–surfactant phase field model , 2017, 1701.07446.

[4]  J. Lowengrub,et al.  A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant , 2004 .

[5]  Junxiang Yang,et al.  An unconditionally stable second-order accurate method for systems of Cahn-Hilliard equations , 2020, Commun. Nonlinear Sci. Numer. Simul..

[6]  Arshad Riaz,et al.  STUDY OF HEAT AND MASS TRANSFER IN THE EYRING–POWELL MODEL OF FLUID PROPAGATING PERISTALTICALLY THROUGH A RECTANGULAR COMPLIANT CHANNEL , 2019, Heat Transfer Research.

[7]  Francisco Guillén-González,et al.  On linear schemes for a Cahn-Hilliard diffuse interface model , 2013, J. Comput. Phys..

[8]  W. Goddard,et al.  New surfactant classes for enhanced oil recovery and their tertiary oil recovery potential , 2010 .

[9]  Harald Garcke,et al.  Diffuse interface modelling of soluble surfactants in two-phase flow , 2013, 1303.2559.

[10]  Rahmat Ellahi,et al.  Numerical study of heat transfer and Hall current impact on peristaltic propulsion of particle-fluid suspension with compliant wall properties , 2019 .

[11]  Zhifeng Weng,et al.  A Fourier spectral method for fractional-in-space Cahn–Hilliard equation , 2017 .

[12]  Hui Zhang,et al.  An energy-stable finite-difference scheme for the binary fluid-surfactant system , 2014, J. Comput. Phys..

[13]  Mohamed Laradji,et al.  The effect of surfactants on the dynamics of phase separation , 1992 .

[14]  Haihu Liu,et al.  Phase-field modeling droplet dynamics with soluble surfactants , 2010, J. Comput. Phys..

[15]  Jaemin Shin,et al.  First and second order numerical methods based on a new convex splitting for phase-field crystal equation , 2016, J. Comput. Phys..

[16]  J. Lowengrub,et al.  Conservative multigrid methods for Cahn-Hilliard fluids , 2004 .

[17]  A Numerical Method for Solving of the Boundary Value Problems for Ordinary Differential Equations , 2009 .

[18]  Junseok Kim Numerical simulations of phase separation dynamics in a water-oil-surfactant system. , 2006, Journal of colloid and interface science.

[19]  Shuyu Sun,et al.  Efficient energy-stable schemes for the hydrodynamics coupled phase-field model , 2019, Applied Mathematical Modelling.

[20]  Hiroya Kodama,et al.  TWO-ORDER-PARAMETER MODEL FOR AN OIL-WATER-SURFACTANT SYSTEM , 1997 .

[21]  M. M. Bhatti,et al.  Effects of magnetic Reynolds number on swimming of gyrotactic microorganisms between rotating circular plates filled with nanofluids , 2020 .

[22]  Xiaofeng Yang,et al.  Efficient, second oder accurate, and unconditionally energy stable numerical scheme for a new hydrodynamics coupled binary phase-field surfactant system , 2020, Comput. Phys. Commun..

[23]  J. Baret Surfactants in droplet-based microfluidics. , 2012, Lab on a chip.

[24]  I-Liang Chern,et al.  Simulating binary fluid-surfactant dynamics by a phase field model , 2012 .

[25]  Qian Yu,et al.  An unconditional stable compact fourth-order finite difference scheme for three dimensional Allen-Cahn equation , 2019, Comput. Math. Appl..

[26]  D. Myers Surfactant Science and Technology , 1988 .

[27]  R. V. D. Sman,et al.  Diffuse interface model of surfactant adsorption onto flat and droplet interfaces , 2006 .

[28]  Junseok Kim,et al.  An efficient linear second order unconditionally stable direct discretization method for the phase-field crystal equation on surfaces , 2019, Applied Mathematical Modelling.

[29]  Shuyu Sun,et al.  Linearly Decoupled Energy-Stable Numerical Methods for Multicomponent Two-Phase Compressible Flow , 2018, SIAM J. Numer. Anal..

[30]  Suchuan Dong,et al.  Numerical approximation of incompressible Navier-Stokes equations based on an auxiliary energy variable , 2018, J. Comput. Phys..

[31]  Jie Shen,et al.  Decoupled, Energy Stable Schemes for Phase-Field Models of Two-Phase Incompressible Flows , 2015, SIAM J. Numer. Anal..

[32]  Wenqiang Feng,et al.  An energy stable fourth order finite difference scheme for the Cahn-Hilliard equation , 2017, J. Comput. Appl. Math..

[33]  Ping Lin,et al.  A numerical method for the quasi-incompressible Cahn-Hilliard-Navier-Stokes equations for variable density flows with a discrete energy law , 2014, J. Comput. Phys..

[34]  Jiang Yang,et al.  The scalar auxiliary variable (SAV) approach for gradient flows , 2018, J. Comput. Phys..

[35]  Shuyu Sun,et al.  Decoupled, energy stable schemes for a phase-field surfactant model , 2018, Comput. Phys. Commun..

[36]  Ping Lin,et al.  Mass conservative and energy stable finite difference methods for the quasi-incompressible Navier–Stokes–Cahn–Hilliard system: Primitive variable and projection-type schemes , 2017, 1703.06606.

[37]  Philip Broadbridge,et al.  Nonclassical Symmetry Solutions for Fourth-Order Phase Field Reaction-Diffusion , 2018, Symmetry.

[38]  Weiqing Ren,et al.  Derivation of a continuum model and the energy law for moving contact lines with insoluble surfactants , 2014 .

[39]  Shuyu Sun,et al.  Fully discrete energy stable scheme for a phase-field moving contact line model with variable densities and viscosities , 2018, Applied Mathematical Modelling.