High-order energy stable schemes of incommensurate phase-field crystal model

Kai Jiang; Wei Si; Kai Jiang; Wei Si

doi:10.3934/era.2020059

Electronic Research Archive

2020, Volume 28, Issue 2: 1077-1093. doi: 10.3934/era.2020059

Previous Article Next Article

Special Issues

High-order energy stable schemes of incommensurate phase-field crystal model

Kai Jiang ^,,
Wei Si

School of Mathematics and Computational Science, Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Xiangtan University, Hunan, Xiangtan 411105, China

Received: 01 April 2020 Revised: 01 May 2020
Primary: 35K55, 35K35; Secondary: 65M70

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.
- Incommensurate phase-field crystal model,
- high-order energy stable schemes,
- spectral deferred correction method,
- projection method,
- quasicrystals
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

Related Papers:

Abstract

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.

References

[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.

Reader Comments

Your name:*

Email:*
© 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)