Stability and Error Estimate of the Operator Splitting Method for the Phase Field Crystal Equation
Abstract
In this paper, we propose a second-order fast explicit operator splitting method for the phase field crystal equation. The basic idea lied in our method is to split the original problem into linear and nonlinear parts. The linear subproblem is numerically solved using the Fourier spectral method, which is based on the exact solution and thus has no stability restriction on the time-step size. The nonlinear one is solved via second-order strong stability preserving Runge–Kutta method. The stability and convergence are discussed in -norm. Numerical experiments are performed to validate the accuracy and efficiency of the proposed method. Moreover, energy degradation and mass conservation are also verified.
References
[1]
Elder KR, Katakowski M, Haataja M, and Grant M Modeling elasticity in crystal growth Phys. Rev. Lett. 2002 88 24 245701
[2]
Elder KR and Grant M Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals Phys. Rev. E 2004 70 5 051605
[3]
Berry J, Provatas N, Rottler J, and Sinclair CW Defect stability in phase-field crystal models: stacking faults and partial dislocations Phys. Rev. B 2012 86 224112
[4]
Provatas N, Dantzig J, Athreya B, Chan P, Stefanovic P, Goldenfeld N, and Elder K Using the phase-field crystal method in the multi-scale modeling of microstructure evolution JOM 2007 59 83-90
[5]
Stolle J and Provatas N Characterizing solute segregation and grain boundary energy in binary alloy phase field crystal models Comput. Mater. Sci. 2014 81 493-502
[6]
Swift J and Hohenberg PC Hydrodynamic fluctuations at the convective instability Phys. Rev. A 1977 15 319
[7]
Wang C, Wise SM, and Lowengrub JS An energy-stable and convergent finite-difference scheme for the phase field crystal equation SIAM J. Numer. Anal. 2009 47 2269-2288
[8]
Wang C and Wise SM An energy stable and convergent finite-difference scheme for the modified phase field crystal equation SIAM J. Numer. Anal. 2011 49 945-969
[9]
Baskaran A, Hu ZZ, Lowengrub JS, Wang C, Wise SM, and Zhou P Energy stable and efficient finite-difference nonlinear multigrid schemes for the modified phase field crystal equation J. Comput. Phys. 2013 250 270-292
[10]
Baskaran A, Lowengrub JS, Wang C, and Wise SM Convergence analysis of a second order convex splitting scheme for the modified phase field crystal equation SIAM J. Numer. Anal. 2013 51 2851-2873
[11]
Li Q, Mei LQ, and You B A second-order, uniquely solvable, energy stable BDF numerical scheme for the phase field crystal model Appl. Numer. Math. 2018 134 46-65
[12]
Xia BH, Mei CL, Yu Q, and Li YB A second order unconditionally stable scheme for the modified phase field crystal model with elastic interaction and stochastic noise effect Comput. Methods Appl. Mech. Eng. 2020 363 112795
[13]
Cheng, K.L., Wang, C., Wise, S.M.: An energy stable BDF2 Fourier pseudo-spectral numerical scheme for the square phase field crystal equation. arXiv:1906.12255 [math.NA]
[14]
Yang XF and Han DZ Linearly first- and second-order, unconditionally energy stable schemes for the phase field crystal model J. Comput. Phys. 2017 330 1116-1134
[15]
Zhao J, Wang Q, and Yang XF Numerical approximations for a phase field dendritic crystal growth model based on the invariant energy quadratization approach Int. J. Numer. Methods Eng. 2017 110 279-300
[16]
Li Q, Mei LQ, Yang XF, and Li YB Efficient numerical schemes with unconditional energy stabilities for the modified phase field crystal equation Adv. Comput. Math. 2019 45 1551-1580
[17]
Liu ZG and Li XL Two fast and efficient linear semi-implicit approaches with unconditional energy stability for nonlocal phase field crystal equation Appl. Numer. Math. 2020 150 491-506
[18]
Liu ZG and Li XL Efficient modified stabilized invariant energy quadratization approaches for phase-field crystal equation Numer. Algorithms 2020 85 107-132
[19]
Li XL and Shen J Stability and error estimates of the SAV Fourier-spectral method for the phase field crystal equation Adv. Comput. Math. 2020 46 48
[20]
Li, X.L., Shen, J.: Efficient linear and unconditionally energy stable schemes for the modified phase field crystal equation. arXiv:2004.04319 [math.NA]
[21]
Goldman D and Kaper T th-order operator splitting schemes and nonreversible systems SIAM J. Numer. Anal. 1996 33 349-367
[22]
Strang G On the construction and comparison of difference schemes SIAM J. Numer. Anal. 1968 5 506-517
[23]
Yanenko NN The Method of Fractional Steps: The Solution of Problems of Mathematical Physics in Several Variables 1971 New York Springer
[24]
Holden H, Karlsen KH, and Risebro NH Operator splitting methods for generalized Korteweg-de Vries equations J. Comput. Phys. 1999 153 203-222
[25]
Cheng YZ, Kurganov A, Qu ZL, and Tang T Fast and stable explicit operator splitting methods for phase-field models J. Comput. Phys. 2015 303 45-65
[26]
Zhai SY, Weng ZF, and Feng XL Fast explicit operator splitting method and time-step adaptivity for fractional non-local Allen–Cahn model Appl. Math. Model. 2016 40 1315-1324
[27]
Zhai SY, Wu LY, Wang JY, and Weng ZF Numerical approximation of the fractional Cahn–Hilliard equation by operator splitting method Numer. Algorithms 2020 84 1155-1178
[28]
Li X, Qiao ZH, and Zhang H Convergence of a fast explicit operator splitting method for the epitaxial growth model with slope selection SIAM J. Numer. Anal. 2017 55 265-285
[29]
Bao WZ, Li HL, and Shen J A generalized-Laguerre–Fourier–Hermite pseudospectral method for computing the dynamics of rotating Bose–Einstein condensates SIAM J. Sci. Comput. 2009 31 3685-3711
[30]
Shen J and Wang ZQ Error analysis of the Strang time-splitting Laguerre–Hermite/Hermite collocation methods for the Gross–Pitaevskii equation Found. Comput. Math. 2013 13 99-137
[31]
Zhang C, Huang JF, Wang C, and Yue XY On the operator splitting and integral equation preconditioned deferred correction methods for the “Good” Boussinesq equation J. Sci. Comput. 2018 75 687-712
[32]
Zhang C, Wang H, Huang JF, Wang C, and Yue XY A second order operator splitting numerical scheme for the “good” Boussinesq equation Appl. Numer. Math. 2017 119 179-193
[33]
Lubich C On splitting methods for Schrödinger–Poisson and cubic nonlinear Schrödinger equations Math. Comput. 2008 77 2141-2153
[34]
Thalhammer M Convergence analysis of high-order time-splitting pseudospectral methods for nonlinear Schrödinger equations SIAM J. Numer. Anal. 2012 50 3231-3258
[35]
Zhai SY, Wang DL, Weng ZF, and Zhao X Error analysis and numerical simulations of Strang splitting method for space Fractional nonlinear Schrödinger equation J. Sci. Comput. 2019 81 965-989
[36]
Lee HG, Shin J, and Lee JY First and second order operator splitting methods for the phase field crystal equation J. Comput. Phys. 2015 299 82-91
[37]
Shen J, Tang T, and Wang LL Spectral Methods Algorithms: Analyses and Applications 2010 1 Berlin Springer
[38]
Gottlieb S and Shu CW Total variation diminishing Runge–Kutta schemes Math. Comput. 1998 67 73-85
[39]
Mishra S and Svärd M On stability of numerical schemes via frozen coefficients and the magnetic induction equations BIT Numer. Math. 2010 50 85-108
[40]
Tadmor E The exponential accuracy of Fourier and Chebyshev differencing methods SIAM J. Numer. Anal. 1986 23 1-10
[41]
Canuto C, Quarteroni A, Hussaini MY, and Zang TA Spectral Methods: Fundamentals in Single Domains 2006 Berlin Springer
[42]
Gottlieb S and Wang C Stability and convergence analysis of fully discrete Fourier collocation spectral method for 3-D viscous Burgers’ Equation J. Sci. Comput. 2012 53 102-128
[43]
Gottlieb S, Tone F, Wang C, Wang X, and Wirosoetisno D Long time stability of a classical efficient scheme for two-dimensional Navier–Stokes equations SIAM J. Numer. Anal. 2012 50 126-150
Recommendations
Numerical approximation of the fractional Cahn-Hilliard equation by operator splitting method
AbstractIn this paper, we consider a fast explicit operator splitting method for a fractional Cahn-Hilliard equation with spatial derivative (α ∈ (1,2]), where the choice α = 2 corresponds to the classical Cahn-Hilliard equation. The original ...
A Second-Order Operator Splitting Fourier Spectral Method for Models of Epitaxial Thin Film Growth
In this paper, we develop an operator splitting Fourier spectral method for models of epitaxial thin film growth with and without slope selection. A main idea of the method is to split the original equation into linear and nonlinear parts, and then to ...
First and second order operator splitting methods for the phase field crystal equation
In this paper, we present operator splitting methods for solving the phase field crystal equation which is a model for the microstructural evolution of two-phase systems on atomic length and diffusive time scales. A core idea of the methods is to ...
Comments
Please enable JavaScript to view thecomments powered by Disqus.Information & Contributors
Information
Published In
© The Author(s), under exclusive licence to Springer Science+Business Media, LLC part of Springer Nature 2021.
Publisher
Plenum Press
United States
Publication History
Published: 01 January 2021
Accepted: 10 December 2020
Revision received: 02 November 2020
Received: 14 August 2020
Author Tags
Author Tags
Qualifiers
- Research-article
Funding Sources
Contributors
Other Metrics
Bibliometrics & Citations
Bibliometrics
Article Metrics
- 0Total Citations
- 0Total Downloads
- Downloads (Last 12 months)0
- Downloads (Last 6 weeks)0
Reflects downloads up to 22 Nov 2024
Other Metrics
Citations
View Options
View options
Login options
Check if you have access through your login credentials or your institution to get full access on this article.
Sign in