Arbitrarily High-Order Unconditionally Energy Stable Schemes for Thermodynamically Consistent Gradient Flow Models
Abstract
We present a systematic approach to developing arbitrarily high-order, unconditionally energy stable numerical schemes for thermodynamically consistent gradient flow models that satisfy energy dissipation laws. Utilizing the energy quadratization method, we formulate the gradient flow model into an equivalent form with a corresponding quadratic free energy functional. Based on the equivalent form with a quadratic energy, we propose two classes of energy stable numerical approximations. In the first approach, we use a prediction-correction strategy to improve the accuracy of linear numerical schemes. In the second approach, we adopt the Gaussian collocation method to discretize the equivalent form with a quadratic energy, arriving at an arbitrarily high-order scheme for gradient flow models. Schemes derived using both approaches are proved rigorously to be unconditionally energy stable. The proposed schemes are then implemented in four gradient flow models numerically to demonstrate their accuracy and effectiveness. Detailed numerical comparisons among these schemes are carried out as well. These numerical strategies are rather general so that they can be readily generalized to solve any thermodynamically consistent PDE models.
1. . Gutierrez-Santacreu, Finite element approximation of nematic liquid crystal flows using a saddle-point structure , J. Comput. Phys. , 230 ( 2011 ), pp. 1686 -- 1706 .
2. R. Bird, C. F. Curtiss, R. C. Armstrong, and O. Hassager, Dynamics of Polymeric Liquids: Kinetic Theory, Vol. 2, John Wiley, New York, 1987.
3. , Regularized linear schemes for the molecular beam epitaxy model with slope selection , Appl. Numer. Math. , 128 ( 2018 ), pp. 138 -- 156 .
4. , Applications of semi-implicit Fourier-spectral method to phase-field equations , Comput. Phys. Commun. , 108 ( 1998 ), pp. 147 -- 158 .
5. , Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential , J. Comput. Phys. , 3 ( 2019 ), 100031 .
6. , A linear iteration algorithm for a second-order energy stable scheme for a thin film model without slope selection , J. Sci. Comput. , 59 ( 2014 ), pp. 574 -- 601 .
7. , An energy stable fourth order finite difference scheme for the Cahn--Hilliard equation , J. Comput. Appl. Math. , 362 ( 2019 ), pp. 574 -- 595 .
8. , A general framework for deriving integral preserving numerical methods for PDEs , SIAM J. Sci. Comput. , 33 ( 2011 ), pp. 2318 -- 2340 .
9. , Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals , Phys. Rev. E , 70 ( 2004 ).
10. , Modeling elasticity in crystal growth , Phys. Rev. Lett. , 88 ( 2002 ).
11. , A simple and efficient scheme for phase field crystal simulation , ESAIM Math. Model Numer. Anal. , 47 ( 2013 ), pp. 1413 -- 1432 .
12. , Unconditionally gradient stable time marching the Cahn-Hilliard equation , Comput. Math. Models Microstruct. Evol. , 529 ( 1998 ), pp. 39 -- 46 .
13. , A uniquely solvable, energy stable numerical scheme for the functionalized Cahn--Hilliard equation and its convergence analysis , J. Sci. Comput. , 76 ( 2018 ), pp. 1938 -- 1967 .
14. , A second-order energy stable backward differentiation formula method for the epitaxial thin film equation with slope selection , Numer. Methods Partial Differential Equations , 34 ( 2018 ), pp. 1975 -- 2007 .
15. , Stabilized Crank-Nicolson/Adams-Bashforth schemes for phase field models , East Asian J. Appl. Math. , 3 ( 2013 ), pp. 59 -- 80 .
16. , Improving the accuracy of convexity splitting methods for gradient flow equations , J. Comput. Phys. , 315 ( 2016 ), pp. 52 -- 64 .
17. , Provably unconditionally stable, second-order time-accurate, mixed variational methods for phase-field models , J. Comput. Phys. , 230 ( 2011 ), pp. 5310 -- 5327 .
18. H. Gomez and X. Nogueira, An unconditionally energy-stable method for the phase field crystal equation, Comput. Methods Appl. Mech. Engrg., 249--252 (2012), pp. 52--61.
19. , Fully discretized energy stable schemes for hydrodynamic equations governing two-phase viscous fluid flows , J. Sci. Comput. , 69 ( 2016 ), pp. 921 -- 945 .
20. , An energy stable algorithm for a quasi-incompressible hydrodynamic phase-field model of viscous fluid mixtures with variable densities and viscosities , Comput. Phys. Commun. , 219 ( 2017 ), pp. 20 -- 34 .
21. , Linear second order in time energy stable schemes for hydrodynamic models of binary mixtures based on a spatially pseudospectral approximation , Adv. Comput. Math. , 44 ( 2018 ), pp. 1573 -- 1600 .
22. , Second order fully discrete energy stable methods on staggered grids for hydrodynamic phase field models of binary viscous fluids , SIAM J. Sci. Comput. , 40 ( 2018 ), pp. B528 -- B553 .
23. , Second-order linear schemes for hydrodynamic phase field models of viscous fluid flows with variable densities , SIAM J. Sci. Comput. , 4 ( 2018 ), pp. B138 -- B167 .
24. , A thermodynamically consistent numerical method for a phase field model of solidification , Commun. Nonlinear Sci. Numer. Simul. , 19 ( 2014 ), pp. 2309 -- 2323 .
25. , Second order schemes and time-step adaptivity for Allen-Cahn and Cahn-Hilliard models , Comput. Math. Appl. , 68 ( 2014 ), pp. 821 -- 846 .
26. , An H$2$ convergence of a second-order convex-splitting, finite difference scheme for the three-dimensional Cahn-Hilliard equation , Commun. Math. Sci. , 14 ( 2016 ), pp. 489 -- 515 .
27. , Mass conservative and energy stable finite difference methods for the quasi-incompressible Navier-Stokes-Cahn-Hilliard system: Primitive variable and projection-type schemes , Comput. Methods Appl. Mech. Engrg. , 326 ( 2017 ), pp. 144 -- 174 .
28. E. Hairer, C. Lubich, and G. Wanner, Geometric Numerical Integration: Structure-preserving Algorithms for Ordinary Differential Equations, Springer Ser. Comput. Math. 31, Springer, New York, 2006.
29. , A second order in time uniquely solvable unconditionally stable numerical schemes for Cahn-Hilliard-Navier-Stokes equation , J. Comput. Phys. , 290 ( 2015 ), pp. 139 -- 156 .
30. , Energy stability and error estimates of exponential time differencing schemes for the epitaxial growth model without slope selection , Math. Comp. , 87 ( 2018 ), pp. 1859 -- 1885 .
31. , Epitaxial growth without slope selection: Energetics, coarsening and dynamic scaling , J. Nonlinear Sci. , 14 ( 2004 ), pp. 429 -- 451 .
32. , Thin film epitaxy with or without slope selection , European J. Appl. Math. , 14 ( 2003 ), pp. 713 -- 743 .
33. , Characterizing the stabilization size for semi-implicit Fourier-spectral method to phase field equations , SIAM J. Numer. Anal. , 54 ( 2016 ), pp. 1653 -- 1681 .
34. , A class of conservative phase field models for multiphase fluid flows , J. Appl. Mech. , 81 ( 2014 ), 021004 .
35. , Stabilized semi-implicit spectral deferred correction methods for Allen-Cahn and Cahn-Hilliard equations , Math. Methods Appl. Sci. , 38 ( 2013 ), pp. 4564 -- 4575 .
36. , Reciprocal relations in irreversible processes. I , Phys. Rev. , 37 ( 1931 ), pp. 405 -- 426 , https://doi.org/10.1103/PhysRev.37.405.
37. , Second-order convex splitting schemes for gradient flows with Ehrlich--Schwoebel type energy: Application to thin film epitaxy , SIAM J. Numer. Anal. , 50 ( 2012 ), pp. 105 -- 125 .
38. , Stabilized predictor-corrector schemes for gradient flows with strong anisotropic free energy , Commun. Comput. Phys. , 24 ( 2018 ), pp. 635 -- 654 .
39. J. Shen, J. Xu, and J. Yang, A New Class of Efficient and Robust Energy Stable Schemes for Gradient Flows, arXiv:1710.01331, 2017.
40. , The scalar auxiliary variable (SAV) approach for gradient flows , J. Comput. Phys. , 353 ( 2018 ), pp. 407 -- 416 .
41. , Numerical approximations of Allen-Cahn and Cahn-Hilliard equations , Discrete Contin. Dyn. Syst. Ser. A , 28 ( 2010 ), pp. 1669 -- 1691 .
42. , A phase-field model and its numerical approximation for two-phase incompressible flows with different densities and viscosities , SIAM J. Sci. Comput. , 32 ( 2010 ), pp. 1159 -- 1179 .
43. , Efficient energy stable numerical schemes for a phase field moving contact line model , J. Comput. Phys. , 284 ( 2015 ), pp. 617 -- 630 .
44. , Unconditionally stable methods for gradient flow using convex splitting runge-kutta scheme , J. Comput. Phys. , 347 ( 2017 ), pp. 367 -- 381 .
45. , An energy stable convex splitting for the phase-field crystal equation , Comput. Struct. , 158 ( 2015 ), pp. 355 -- 368 .
46. , Unconditionally stable schemes for equations of thin film epitaxy , Discrete Contin. Dyn. Syst. , 28 ( 2010 ), pp. 405 -- 423 .
47. , An energy stable and convergent finite-difference scheme for the modified phase field crystal equation , SIAM J. Numer. Anal. , 49 ( 2011 ), pp. 945 -- 969 .
48. , An energy-stable and convergent finite-difference scheme for the phase field crystal equation , SIAM J. Numer. Anal. , 47 ( 2009 ), pp. 2269 -- 2288 .
49. , A second-order energy stable BDF numerical scheme for the Cahn-Hilliard equation , Commun. Comput. Phys , 32 ( 2018 ), pp. 572 -- 602 .
50. , Linearly first- and second-order, unconditionally energy stable schemes for the phase field crystal equation , J. Comput. Phys. , 333 ( 2017 ), pp. 1116 -- 1134 .
51. , Hydrodynamic theories for flows of active liquid crystals and the generalized onsager principle , Entropy , 18 ( 2016 ), 202 .
52. , On linear and unconditionally energy stable algorithms for variable mobility Cahn-Hilliard type equation with logarithmic Flory-Huggins potential , Commun. Comput. Phys , 25 ( 2019 ), pp. 703 -- 728 .
53. , Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method , J. Comput. Phys. , 333 ( 2017 ), pp. 102 -- 127 .
54. , Numerical approximations for a three component Cahn-Hilliard phase-field model based on the invariant energy quadratization method , Math. Models Methods Appl. Sci. , 27 ( 2017 ), pp. 1993 -- 2023 .
55. , A linearly decoupled energy stable scheme for phase-field models of three-phase incompressible flows , J. Sci. Comput. , 70 ( 2017 ), pp. 1367 -- 1389 .
56. , Numerical approximations to a new phase field model for two phase flows of complex fluids , Comput. Method Appl. Mech. Engrg. , 310 ( 2016 ), pp. 77 -- 97 .
57. , A novel linear second order unconditionally energy-stable scheme for a hydrodynamic Q tensor model for liquid crystals , Comput. Methods Appl. Mech. Engrg. , 318 ( 2017 ), pp. 803 -- 825 .
58. , A general strategy for numerical approximations of non-equilibrium models-part I: Thermodynamical systems , Int. J. Numer. Anal. Model. , 15 ( 2018 ), pp. 884 -- 918 .
