This article focuses on the development of high-order energy stable schemes for the multi-length-scale incommensurate phase-field crystal model which is able to study the phase behavior of aperiodic structures. These high-order schemes based on the scalar auxiliary variable (SAV) and spectral deferred correction (SDC) approaches are suitable for the $ L^2 $ gradient flow equation, i.e., the Allen-Cahn dynamic equation. Concretely, we propose a second-order Crank-Nicolson (CN) scheme of the SAV system, prove the energy dissipation law, and give the error estimate in the almost periodic function sense. Moreover, we use the SDC method to improve the computational accuracy of the SAV/CN scheme. Numerical results demonstrate the advantages of high-order numerical methods in numerical computations and show the influence of length-scales on the formation of ordered structures.
Citation: Kai Jiang, Wei Si. High-order energy stable schemes of incommensurate phase-field crystal model[J]. Electronic Research Archive, 2020, 28(2): 1077-1093. doi: 10.3934/era.2020059
This article focuses on the development of high-order energy stable schemes for the multi-length-scale incommensurate phase-field crystal model which is able to study the phase behavior of aperiodic structures. These high-order schemes based on the scalar auxiliary variable (SAV) and spectral deferred correction (SDC) approaches are suitable for the $ L^2 $ gradient flow equation, i.e., the Allen-Cahn dynamic equation. Concretely, we propose a second-order Crank-Nicolson (CN) scheme of the SAV system, prove the energy dissipation law, and give the error estimate in the almost periodic function sense. Moreover, we use the SDC method to improve the computational accuracy of the SAV/CN scheme. Numerical results demonstrate the advantages of high-order numerical methods in numerical computations and show the influence of length-scales on the formation of ordered structures.
[1] | Phenomenological theory of icosahedral incommensurate ("quasiperiodic") order in Mn-Al alloys. Phys. Rev. Lett. (1985) 54: 1517-1519. |
[2] | C. Corduneanu, Almost Periodic Functions, 2$^{nd}$ edition, Chelsea Publishing Company, New York, 1989. |
[3] | Simultaneous Diophantine approximation. Duke Math. J. (1946) 13: 105-111. |
[4] | Spectral deferred correction methods for ordinary differential equations. BIT Numer. Math. (2000) 40: 241-266. |
[5] | Semi-implicit spectral deferred correction method based on the invariant energy quadratization approach for phase field problems. Commun. Comput. Phys. (2019) 26: 87-113. |
[6] | Long-range icosahedral orientational order and quasicrystals.. Phys. Rev. Lett. (1985) 55: 607-610. |
[7] | Stability of three-dimensional icosahedral quasicrystals in multi-component systems. Philos. Mag. (2020) 100: 84-109. |
[8] | K. Jiang, J. Tong, P. Zhang and A.-C. Shi, Stability of two-dimensional soft quasicrystals in systems with two length scales, Phys. Rev. E, 92 (2015), 042159. |
[9] | Numerical methods for quasicrystals. J. Comput. Phys. (2014) 256: 428-440. |
[10] | K. Jiang, P. Zhang and A.-C. Shi, Stability of icosahedral quasicrystals in a simple model with two-length scales, J. Phys.: Condens. Matter, 29 (2017), 124003. doi: 10.1088/1361-648X/aa586b |
[11] | X. Li and J. Shen, Stability and error estimates of the SAV Fourier-spectral method for the phase field crystal equation, Advances in Computational Mathematics, arXiv: 1907.07462, 46 (2020), Article number: 48. doi: 10.1007/s10444-020-09789-9 |
[12] | Soft quasicrystals–Why are they stable?. Philos. Mag. (2007) 87: 3021-3030. |
[13] | Theoretical model for Faraday waves with multiple-frequency forcing. Phys. Rev. Lett. (1997) 79: 1261-1264. |
[14] | Multiple-scale structures: From Faraday waves to soft-matter quasicrystals. IUCrJ (2018) 5: 247-268. |
[15] | The scalar auxiliary variable (SAV) approach for gradient flows. J. Comput. Phys. (2018) 353: 407-416. |
[16] | Numerical approximations of Allen-Cahn and Cahn-Hilliard equations. Discrete Contin. Dyn. Syst. (2010) 28: 1669-1691. |
[17] | An energy-stable and convergent finite-difference scheme for the phase field crystal equation. SIAM J. Numer. Anal. (2009) 47: 2269-2288. |
[18] | 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. |