Efficient energy-stable schemes for the hydrodynamics coupled phase-field model

Abstract In this article, several efficient and energy-stable semi–implicit schemes are presented for the Cahn–Hilliard phase-field model of two-phase incompressible flows. A scalar auxiliary variable (SAV) approach is implemented to solve the Cahn–Hilliard equation, while a splitting method based on pressure stabilization is used to solve the Navier–Stokes equation. At each time step, the schemes involve solving only a sequence of linear elliptic equations, and computations of the phase-field variable, velocity, and pressure are totally decoupled. A finite-difference method on staggered grids is adopted to spatially discretize the proposed time-marching schemes. We rigorously prove the unconditional energy stability for the semi-implicit schemes and the fully discrete scheme. Numerical results in both two and three dimensions are obtained, which demonstrate the accuracy and effectiveness of the proposed schemes. Using our numerical schemes, we compare the SAV, invariant energy quadratization (IEQ), and stabilization approaches. Bubble rising dynamics and coarsening dynamics are also investigated in detail. The results demonstrate that the SAV approach is more accurate than the IEQ approach and that the stabilization approach is the least accurate among the three approaches. The energy stability of the SAV approach appears to be better than that of the other approaches at large time steps.

[1]  Xiaofeng Yang,et al.  Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends , 2016, J. Comput. Phys..

[2]  T. Hughes,et al.  Isogeometric analysis of the Cahn–Hilliard phase-field model , 2008 .

[3]  Lei Zhang,et al.  Pore scale simulation of liquid and gas two-phase flow based on digital core technology , 2015 .

[4]  S. M. Wise,et al.  Unconditionally Stable Finite Difference, Nonlinear Multigrid Simulation of the Cahn-Hilliard-Hele-Shaw System of Equations , 2010, J. Sci. Comput..

[5]  Jie Shen,et al.  Efficient energy stable numerical schemes for a phase field moving contact line model , 2015, J. Comput. Phys..

[6]  Xiaofeng Yang,et al.  Numerical approximations of Allen-Cahn and Cahn-Hilliard equations , 2010 .

[7]  James J. Feng,et al.  A diffuse-interface method for simulating two-phase flows of complex fluids , 2004, Journal of Fluid Mechanics.

[8]  Xiaofeng Yang,et al.  Numerical approximations for a phase-field moving contact line model with variable densities and viscosities , 2017, J. Comput. Phys..

[9]  Min Gao,et al.  A gradient stable scheme for a phase field model for the moving contact line problem , 2012, J. Comput. Phys..

[10]  Min Gao,et al.  An efficient scheme for a phase field model for the moving contact line problem with variable density and viscosity , 2014, J. Comput. Phys..

[11]  D. Jacqmin.  Regular Article: Calculation of Two-Phase Navier–Stokes Flows Using Phase-Field Modeling , 1999 .

[12]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[13]  Zhixue Sun,et al.  The numerical simulation of thermal recovery based on hydraulic fracture heating technology in shale gas reservoir , 2016 .

[14]  Franck Boyer,et al.  Numerical schemes for a three component Cahn-Hilliard model , 2011 .

[15]  Jie Shen,et al.  A Phase-Field Model and Its Numerical Approximation for Two-Phase Incompressible Flows with Different Densities and Viscosities , 2010, SIAM J. Sci. Comput..

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

[17]  Ronald Fedkiw,et al.  A review of level-set methods and some recent applications , 2018, J. Comput. Phys..

[18]  Guangpu Zhu,et al.  Pore-Scale Investigation of Carbon Dioxide-Enhanced Oil Recovery , 2017 .

[19]  Axel Voigt,et al.  Benchmark computations of diffuse interface models for two‐dimensional bubble dynamics , 2012 .

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

[21]  Shuyu Sun,et al.  Thermodynamically consistent simulation of nonisothermal diffuse-interface two-phase flow with Peng-Robinson equation of state , 2017, J. Comput. Phys..

[22]  Jia Zhao,et al.  Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method , 2017, J. Comput. Phys..

[23]  Lili Ju,et al.  Efficient linear schemes with unconditional energy stability for the phase field elastic bending energy model , 2017 .

[24]  Tao Zhang,et al.  Energy Stability Analysis of Some Fully Discrete Numerical Schemes for Incompressible Navier–Stokes Equations on Staggered Grids , 2018, J. Sci. Comput..

[25]  Jie Shen,et al.  Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[26]  Charles M. Elliott,et al.  The global dynamics of discrete semilinear parabolic equations , 1993 .

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

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

[29]  Shuyu Sun,et al.  The Superconvergence Phenomenon and Proof of the MAC Scheme for the Stokes Equations on Non-uniform Rectangular Meshes , 2015, J. Sci. Comput..

[30]  Jie Shen,et al.  Numerical simulations of jet pinching-off and drop formation using an energetic variational phase-field method , 2006, J. Comput. Phys..

[31]  Xiaofeng Yang,et al.  Numerical approximations for a three-component Cahn–Hilliard phase-field model based on the invariant energy quadratization method , 2017 .

[32]  Shuyu Sun,et al.  A Novel Energy Stable Numerical Scheme for Navier-Stokes-Cahn-Hilliard Two-Phase Flow Model with Variable Densities and Viscosities , 2018, ICCS.

[33]  Ying Chen,et al.  Efficient, adaptive energy stable schemes for the incompressible Cahn-Hilliard Navier-Stokes phase-field models , 2016, J. Comput. Phys..

[34]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[35]  Jisheng Kou,et al.  Numerical Approximation of a Phase-Field Surfactant Model with Fluid Flow , 2018, J. Sci. Comput..

[36]  D. Kuzmin,et al.  Quantitative benchmark computations of two‐dimensional bubble dynamics , 2009 .

[37]  Yi Shi,et al.  A finite element method for the numerical solution of the coupled Cahn-Hilliard and Navier-Stokes system for moving contact line problems , 2012, J. Comput. Phys..