A highly efficient and accurate new scalar auxiliary variable approach for gradient flows. (English) Zbl 1451.65210
Summary: We present several essential improvements to the powerful scalar auxiliary variable (SAV) approach. Firstly, by using the introduced scalar variable to control both the nonlinear and the explicit linear terms, we are able to reduce the number of linear equations with constant coefficients to be solved at each time step from two to one, so the computational cost of the new SAV approach is essentially half of the original SAV approach while keeping all its other advantages. This technique is also extended to the multiple SAV approach. Secondly, instead of discretizing the dynamical equation for the auxiliary variable, we use a first-order approximation of the energy balance equation, which allows us to construct high-order unconditionally energy-stable SAV schemes with uniform and, more importantly, variable time step sizes, enabling us to construct, for the first time, high-order unconditionally stable adaptive time-stepping backward differentiation formula schemes. Representative numerical examples are provided to demonstrate the improved efficiency and accuracy of the proposed method.
MSC:
| 65N35 | Spectral, collocation and related methods for boundary value problems involving PDEs |
| 65N22 | Numerical solution of discretized equations for boundary value problems involving PDEs |
| 65M06 | Finite difference methods for initial value and initial-boundary value problems involving PDEs |
| 65F05 | Direct numerical methods for linear systems and matrix inversion |
| 35J05 | Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation |
References:
| [1] | S.M. Allen and J.W. Cahn, A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening, Acta Metallurgica, 27 (1979), pp. 1085-1095. |
| [2] | M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman, and E.A. Cornell, Observation of Bose-Einstein condensation in a dilute atomic vapor, Science, 269 (1995), pp. 198-201. |
| [3] | W. Bao and Q. Du, Computing the ground state solution of Bose-Einstein condensates by a normalized gradient flow, SIAM J. Sci. Comput., 25 (2004), pp. 1674-1697. · Zbl 1061.82025 |
| [4] | J.W. Barrett and J.F. Blowey, Finite element approximation of a model for phase separation of a multi-component alloy with non-smooth free energy, Numer. Math., 77 (1997), pp. 1-34. · Zbl 0882.65129 |
| [5] | J.W. Barrett and J.F. Blowey, Finite element approximation of a model for phase separation of a multi-component alloy with a concentration-dependent mobility matrix, IMA J. Numer. Anal., 18 (1998), pp. 287-328. · Zbl 0899.76262 |
| [6] | J.W. Cahn and J.E. Hilliard., Free energy of a nonuniform system. I. Interfacial free energy, J. Chem. Phys., 28 (1958), pp. 258-267. · Zbl 1431.35066 |
| [7] | J.W. Cahn and J.E. Hilliard, Free energy of a nonuniform system. III. Nucleation in a two-component incompressible fluid, J. Chem. Phys., 31 (1959), pp. 688-699. |
| [8] | E. Celledoni, V. Grimm, R.I. McLachlan, D.I. McLaren, D. O’Neale, B. Owren, and G.R.W. Quispel, Preserving energy resp. dissipation in numerical PDEs using the “Average Vector Field” method, J. Comput. Phys., 231 (2012), pp. 6770-6789. · Zbl 1284.65184 |
| [9] | W. Chen, X. Wang, Y. Yan, and Z. Zhang, A second order BDF numerical scheme with variable steps for the Cahn-Hilliard equation, SIAM J. Numer. Anal., 57 (2019), pp. 495-525. · Zbl 1435.65142 |
| [10] | Q. Cheng and J. Shen, Multiple scalar auxiliary variable (MSAV) approach and its application to the phase-field vesicle membrane model, SIAM J. Sci. Comput., 40 (2018), pp. A3982-A4006. · Zbl 1403.65045 |
| [11] | Q. Cheng, J. Shen, and X.F. Yang, Highly efficient and accurate numerical schemes for the epitaxial thin film growth models by using the SAV approach, J. Sci. Comput., 78 (2019), pp. 1467-1487. · Zbl 1444.76020 |
| [12] | C.M. Elliott and A.M. Stuart, The global dynamics of discrete semilinear parabolic equations, SIAM J. Numer. Anal., 30 (1993), pp. 1622-1663. · Zbl 0792.65066 |
| [13] | D.J. Eyre, Unconditionally gradient stable time marching the Cahn-Hilliard equation, in Computational and Mathematical Models of Microstructural Evolution (San Francisco, CA, 1998), Warrendale, PA, 1998, pp. 39-46. |
| [14] | F. Guillen-Gonzalez and G. Tierra, On linear schemes for a Cahn-Hilliard diffuse interface model, J. Comput. Phys., 234 (2013), pp. 140-171. · Zbl 1284.35025 |
| [15] | Y. Gong, J. Zhao, and Q. Wang, Arbitrarily high-order unconditionally energy stable schemes for thermodynamically consistent gradient flow models, SIAM J. Sci. Comput., 42 (2020), pp. 135-156. · Zbl 1447.65099 |
| [16] | Y. Gong, J. Zhao, and Q. Wang, Arbitrarily high-order unconditionally energy stable SAV schemes for gradient flow models, Comput. Phys. Commun., 249 (2020), pp. 107033. · Zbl 1447.65099 |
| [17] | Y. He, Y. Liu, and T. Tang, On large time-stepping methods for the Cahn-Hilliard equation, Appl. Numer. Math., 57 (2006), pp. 616-628. · Zbl 1118.65109 |
| [18] | X.L. Li and J. Shen, Error Analysis of the SAV-MAC Scheme for the Navier-Stokes Equations, e-print, https://arxiv.org/abs/1909.05131, 2019. |
| [19] | X.L. Li, J. Shen, and H.X. Rui, Stability and error analysis of a second-order SAV scheme with block-centered finite differences for gradient flows, Math. Comp., 88 (2019), pp. 2047-2068. · Zbl 1422.65165 |
| [20] | L.L Lin, Z.G. Yang, and S.C. Dong, Numerical approximation of incompressible Navier-Stokes equations based on an auxiliary energy variable, J. Comput. Phys., 388 (2019), pp. 1-22. · Zbl 1452.76093 |
| [21] | F. Lou, T. Tang, and H. Xie, Parameter-free time adaptivity based on energy evolution for the Cahn-Hilliard equation, Commun. Comput. Phys., 19 (2016), pp. 1542-1563. · Zbl 1388.65060 |
| [22] | G.R.W. Quispel and D.I. McLaren, A new class of energy-preserving numerical integration methods, J. Phys. A, 41 (2008), 045206. · Zbl 1132.65065 |
| [23] | J. Shen, J. Xu, and J. Yang, The scalar auxiliary variable (SAV) approach for gradient flows, J. Comput. Phys., 353 (2018), pp. 407-416. · Zbl 1380.65181 |
| [24] | J. Shen, J. Xu, and J. Yang, A new class of efficient and robust energy stable schemes, SIAM Rev., 61 (2019), pp. 474-506. · Zbl 1422.65080 |
| [25] | J. Shen, T. Tang, and L.L. Wang, Spectral Methods: Algorithms, Analysis and Applications, Springer, Berlin, 2011. · Zbl 1227.65117 |
| [26] | J. Shen and X.F. Yang, Numerical approximations of Allen-Cahn and Cahn-Hilliard equations, Discrete Contin. Dyn. Syst., 28 (2010), pp. 1669-1691. · Zbl 1201.65184 |
| [27] | X.F. Yang, Linear, first and second-order, unconditionally energy stable numerical schemes for the phase field model of homopolymer blends, J. Comput. Phys., 327 (2016), pp. 294-316. · Zbl 1373.82106 |
| [28] | Z.G. Yang and S.C. Dong, A roadmap for discretely energy-stable schemes for dissipative systems based on a generalized auxiliary variable with guaranteed positivity, J. Comput. Phys., 404 (2020), 109121. · Zbl 1453.65276 |
| [29] | Z.G. Yang and S.C. Dong, An unconditionally energy-stable scheme based on an implicit auxiliary energy variable for incompressible two-phase flows with different densities involving only precomputable coefficient matrices, J. Comput. Phys., 393 (2019), pp. 229-257. · Zbl 1452.76244 |
| [30] | Z.G. Yang, L.L Lin, and S.C. Dong, A family of second-order energy-stable schemes for Cahn-Hilliard type equations, J. Comput. Phys., 383 (2019), pp. 24-54. · Zbl 1451.65129 |
| [31] | Z. Zhang and Z. Qiao, An adaptive time-stepping strategy for the Cahn-Hilliard equation, Commun. Comput. Phys., 11 (2012), pp. 1261-1278. · Zbl 1388.65072 |
| [32] | J. Zhao, Q. Wang, and X.F. Yang, Numerical approximations for a phase field dendritic crystal growth model based on the invariant energy quadratization approach, Internat. J. Numer. Methods Engrg., 110 (2017), pp. 279-300. · Zbl 1365.74138 |
| [33] | J.Z. Zhu, L.Q. Chen, J. Shen, and V. Tikare, Finite element approximation of a model for phase separation of a multi-component alloy with a concentration-dependent mobility matrix, Phys. Rev. E, 60 (1999), pp. 35-64. |
| [34] | Q.Q. Zhuang and J. Shen, Efficient SAV approach for imaginary time gradient flows with applications to one-and multi-component Bose-Einstein condensates, J. Comput. Phys., 396 (2019), pp. 72-88. · Zbl 1452.65116 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
