Arbitrary High-Order Fully-Decoupled Numerical Schemes for Phase-Field Models of Two-Phase Incompressible Flows

Abstract

Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows, it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations. One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them. The resulting schemes are only first-order accurate in time, and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version. In this paper, we employ the spectral deferred correction method to improve the temporal accuracy, based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea. The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency. Within the framework of the spatially discretized local discontinuous Galerkin method, the resulting numerical schemes are fully decoupled, efficient, and high-order accurate in both time and space. Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.

Keywords

Two-phase incompressible flows / Fully-decoupled / High-order accurate / Linear implicit / Spectral deferred correction method / Local discontinuous Galerkin method

References

[1.] Caffarelli LA, Muler NE. An L ∞ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^{\infty }$$\end{document} bound for solutions of the Cahn-Hilliard equation. Arch. Ration. Mech. Anal., 1995, 133: 129-144,

[2.] Chen C, Yang X. Efficient Numerical Scheme for a dendritic Solidification Phase Field model with melt convection. J. Comput. Phys., 2019, 388: 41-62,

[3.] Cockburn B, Shu C-W. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws II: general framework. Math. Comp., 1989, 52: 411-435

[4.] Dong B, Shu C-W. Analysis of a local discontinuous Galerkin method for fourth-order time-depenedent problems. SIAM J. Numer. Anal., 2009, 47: 3240-3268,

[5.] Dutt A, Greengard L, Rokhlin V. Spectral deferred correction methods for ordinary differential equations. BIT, 2000, 40: 241-266,

[6.] Eyre, D.J.: Unconditionally gradient stable time marching the Cahn-Hilliard equation. MRS Online Proceedings Library 529, 39–46 (1998)

[7.] Gong Y, Zhao J, Wang Q. Arbitrarily high-order unconditionally energy stable SAV schemes for gradient flow models. Comput. Phys. Commun., 2020, 249,

[8.] Guermond JL, Minev P, Shen J. An overview of projection methods for incompressible flows. Comput. Methods Appl. Mech. Eng., 2006, 195: 6011-6045,

[9.] Guermond JL, Salgado A. A splitting method for incompressible flows with variable density based on a pressure Poisson equation. J. Comput. Phys., 2009, 228: 2834-2846,

[10.] Guo R, Xia Y, Xu Y. An efficient fully-discrete local discontinuous Galerkin method for the Cahn-Hilliard-Hele-Shaw system. J. Comput. Phys., 2014, 264: 23-40,

[11.] Guo R, Xia Y, Xu Y. Semi-implicit spectral deferred correction methods for highly nonlinear partial differential equations. J. Comput. Phys., 2017, 338: 269-284,

[12.] Han D, Wang X. A second order in time, uniquely solvable, unconditionally stable numerical scheme for Cahn-Hilliard-Navier-Stokes equation. J. Comput. Phys., 2015, 290: 139-156,

[13.] Han D, Wang X. A second order in time, decoupled, unconditionally stable numerical scheme for the Cahn-Hilliard-Darcy system. J. Sci. Comput., 2018, 77: 1210-1233,

[14.] Ju L, Zhang J, Du Q. Fast and accurate algorithms for simulating coarsening dynamics of Cahn-Hilliard equations. Comput. Mater. Sci., 2015, 108: 272-282,

[15.] Li X, Ju L, Meng X. Convergence analysis of exponential time differencing schemes for the Cahn-Hilliard equation. Commun. Comput. Phys., 2019, 26: 1510-1529,

[16.] Liu C, Shen J, Yang X. Decoupled energy stable schemes for a phase-field model of two-phase incompressible flows with variable density. J. Sci. Comput., 2015, 62: 601-622,

[17.] Minion ML. Semi-implicit spectral deferred correction methods for ordinary differential equations. Commun. Math. Sci., 2003, 1: 471-500,

[18.] Minjeaud S. An unconditionally stable uncoupled scheme for a triphasic Cahn-Hilliard/Navier-Stokes model. Numer. Methods Partial Differ. Equ., 2013, 29: 584-618,

[19.] Reed W, Hill T. . Triangular mesh methods for the neutron transport equation, Technical report LA-UR-73-479, 1973 Los Alamos, NM Los Alamos Scientific Laboratory

[20.] Shen J, Xu J, Yang J. The scalar auxiliary variable (SAV) approach for gradient flows. J. Comput. Phys., 2017, 353: 407-416,

[21.] Shen J, Xu J, Yang J. A new class of efficient and robust energy stable schemes for gradient flows. SIAM Rev., 2019, 61: 474-506,

[22.] Shen J, Yang X. Numerical approximations of Allen-Cahn and Cahn-Hilliard equations. Discrete Contin. Dyn. Syst., 2010, 28: 1169-1691,

[23.] Shen J, Yang X. A phase-field model for two-phase flows with large density ratio and its numerical approximation. SIAM J. Sci. Comput., 2010, 32: 1159-1179,

[24.] Shen J, Yang X. Decoupled, energy stable schemes for phase-field models of two-phase incompressible flows. SIAM J. Numer. Anal., 2015, 53: 279-296,

[25.] Shu, C.-W.: Discontinuous Galerkin methods: general approach and stability, numerical solutions of partial differential equations. In: Bertoluzza, S., Falletta, S., Russo, G., Shu, C.-W. (eds.) Advanced Courses in Mathematics CRM Barcelona, pp. 149–201. Birkhauser, Basel (2009)

[26.] Tang T, Qiao Z. Efficient numerical methods for phase-field equations. Scientia Sinica Mathematica, 2020, 48: 1-20

[27.] Wang C, Wise SM. An energy stable and convergent finite-difference scheme for the modified phase field crystal equation. SIAM J. Numer. Anal., 2011, 49: 945-969,

[28.] Xu C, Tang T. Stability analysis of large time-stepping methods for epitaxial growth models. SIAM J. Numer. Anal., 2006, 44: 1759-1779,

[29.] Xu Y, Shu C-W. Local discontinuous Galerkin methods for high-order time-dependent partial differential equations. Commun. Comput. Phys., 2010, 7: 1-46

[30.] Yan F, Xu Y. Stability analysis and error estimates of local discontinuous Galerkin method with semi-implicit spectral deferred correction time-marching for the Allen-Cahn equation. J. Comput. Appl. Math., 2020, 376,

[31.] Yang X. Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends. J. Comput. Phys., 2016, 327: 294-316,

[32.] Yang X. A new efficient fully-decoupled and second-order time-accurate scheme for Cahn-Hilliard phase-field model of three-phase incompressible flow. Comput. Methods Appl. Mech. Eng., 2021, 376,

[33.] Yang X. On a novel fully-decoupled, linear and second-order accurate numerical scheme for the Cahn-Hilliard-Darcy system of two-phase Hele-Shaw flow. Comput. Phys. Commun., 2021, 263,

[34.] Yang X, He X. A fully-discrete decoupled finite element method for the conserved Allen-Cahn type phase-field model of three-phase fluid flow system. Comput. Methods Appl. Mech. Eng., 2022, 389,

[35.] Zhao J, Han D. Second-order decoupled energy-stable schemes for Cahn-Hilliard-Navier-Stokes equations. J. Comput. Phys., 2021, 443,

[36.] Zhou L, Xu Y. Stability analysis and error estimates of semi-implicit spectral deferred correction coupled with local discontinuous Galerkin method for linear convection-diffusion equations. J. Sci. Comput., 2018, 77: 1001-1029,

[37.] Zhu J, Chen LQ, Shen J, Tikare V. Morphological evolution during phase separation and coarsening with strong inhomogeneous elasticity. Model. Simul. Mater. Sci. Eng., 2001, 9: 499-511,

Funding

National Natural Science Foundation of China(12271492); Natural Science Foundation of Henan Province(12271498)