Special Issues

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

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

    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:

  • 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.
  • Reader Comments
  • © 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)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1010) PDF downloads(159) Cited by(1)

Article outline

Figures and Tables

Figures(6)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog