Abstract
A high-order accurate compact scheme for the Swift–Hohenberg equation is presented in this paper. We discretize the Swift–Hohenberg equation by a fourth-order compact finite difference formula in space and a backward differentiation with second-order accurate in time, respectively. A stabilized splitting scheme is presented and a Newton-type iterative method is introduced to deal with the nonlinear term. Therefore, a large time step can be used. The resulting discrete systems are solved by a fast and efficient nonlinear multigrid solver. Adaptive time step method is implemented to reduce the computational cost. Various numerical simulations including a convergence test of the proposed scheme, comparison with second-order scheme, a test of the stability of the proposed scheme, extension of the adaptive time step method, comparison with the phase field crystal equation, a study of the effect of computational domain and boundary condition, and an evolution of Swift–Hohenberg equation in three dimensions, are performed to demonstrate the efficiency of our proposed method.
Similar content being viewed by others
References
Abbaszadeh M, Khodadadian A, Parvizi M, Dehghan M, Heitzinger C (2019) A direct meshless local collocation method for solving stochastic Cahn–Hilliard Cook and stochastic Swift–Hohenberg equations. Eng Anal Bound Elem 98:253–264
Cheng M, Warren JA (2008) An efficient algorithm for solving the phase field crystal model. J Comput Phys 227:6241–6248
Christov CI, Pontes J (2002) Numerical scheme for Swift–Hohenberg equation with strict implementation of lyapunov functional. Math Comput Model 35:87–99
Christov CI, Pontes J, Walgraef D, Velarde MG (1997) Implicit time splitting for fourth-order parabolic equations. Comput Methods Appl Mech Eng 148:209–224
Clayton JD, Knap J (2016) Phase field modeling and simulation of coupled fracture and twinning in single crystals and polycrystals. Comput Methods Appl Mech Eng 312:447–467
Cross M, Greenside H (2009) Pattern formation and dynamics in nonequilibrium systems. Cambridge University Press, Cambridge
Dehghan M, Abbaszadeh M (2017) The meshless local collocation method for solving multi-dimensional Cahn–Hilliard, Swift–Hohenberg and phase field crystal equations. Eng Anal Bound Elem 78:49–64
Dehghan M, Mohammadi V (2015) The numerical solution of Cahn–Hilliard (CH) equation in one, two and three-dimensions via globally radial basis functions (GRBFs) and RBF differential quadrature (RBFDQ). Eng Anal Bound Elem 51:74–100
Elder KR, Grant M (2004) Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals. Phys Rev E 70:051605–051623
Elder KR, Viñals J, Grant M (1992) Ordering dynamics in the two-dimensional stochastic Swift–Hohenberg equation. Phys Rev Lett 68:3024–3027
Gomez H, Nogueira X (2012) A new space-time discretization for the Swift–Hohenberg equation that strictly respects the Lyapunov functional. Commun Nonlinear Sci Numer Simul 17:4930–4946
Kim JS, Kang K, Lowengrub J (2004) Conservative multigrid methods for Cahn–Hilliard fluids. J Comput Phys 193:511–543
Lee HG (2017) A semi-analytical Fourier spectral method for the Swift–Hohenberg equation. Comput Math Appl 74:1885–1896
Lee C, Jeong D, Shin J, Li YB, Kim JS (2014) A fourth-order spatial accurate and practically stable compact scheme for the Cahn–Hilliard equation. Phys A 409:17–28
Lee HG, Shin J, Lee JY (2015) First and second order operator splitting methods for the phase field crystal equation. J Comput Phys 299:82–91
Li YB, Kim JS (2017) An efficient and stable compact fourth-order finite difference scheme for the phase field crystal equation. Comput Methods Appl Mech Eng 319:194–216
Li YB, Lee HG, Xia B, Kim JS (2016) A compact fourth-order finite difference scheme for the three-dimensional Cahn–Hilliard equation. Comput Phys Commun 200:108–116
Li YB, Kim J, Wang N (2017) An unconditionally energy-stable second-order time-accurate scheme for the Cahn–Hilliard equation on surfaces. Commun Nonlinear Sci Numer Simul 53:213–227
Li YB, Choi YH, Kim JS (2017) Computationally efficient adaptive time step method for the Cahn–Hilliard equation. Comput Math Appl 73:1855–1864
Li YB, Luo C, Xia B, Kim J (2019) An efficient linear second order unconditionally stable direct discretization method for the phase-field crystal equation on surfaces. Appl Math Model 67:477–490
Lloyd DJB, Sandstede B, Avitabile D, Champneys AR (2008) Localized hexagon patterns of the planar Swift–Hohenberg equation. SIAM J Appl Dyn Syst 7:1049–1100
Marconi UMB, Tarazona P (1999) Dynamic density functional theory of liquids. J Chem Phys 110:8032–8044
Mohammadi V, Dehghan M (2010) High-corder solution of one dimensional sine Gordon equation using compact finite difference and DIRKN methods. Math Comput Model 51(5):537–549
Mohammadi V, Mohebbi A, Asgari Z (2009) Fourth order compact solution of the nonlinear Klein–Gordon equation. Numer Algorithms 52(4):523–540
Mohebbi A, Dehghan M (2010) High-order compact solution of the one dimensional heat and advection diffusion equations. Appl Math Model 34(10):3071–3084
Nikolay NA, Ryabov PN (2016) Analytical and numerical solutions of the generalized dispersive Swift–Hohenberg equation. Appl Math Comput 286:171–177
Staliunas K, Sánchez-Morcillo VJ (1998) Dynamics of phase domains in the Swift–Hohenberg equation. Phys Lett A 241:28–34
Swift J, Hohenberg PC (1977) Hydrodyamic fluctuations at the convective instability. Phys Rev A 15:319–328
Trottenberg U, Oosterlee C, Schüller A (2001) Multigrid. Academic Press, New York
Vignal P, Dalcin L, Brown DL, Collier N, Calo VM (2015) An energy-stable convex splitting for the phase-field crystal equation. Comput Struct 158:355–368
Viñals J, Hernández-Garca E, San MM, Toral R (1991) Numerical study of the dynamical aspects of pattern selection in the stochastic Swift–Hohenberg equation in one dimension. Phys Rev A 44:1123–1133
Wise SM, Wang C, Lowengrub JS (2009) An energy stable and convergent finite-difference scheme for the phase field crystal equation. SIAM J Numer Anal 47(3):2269–2288
Xi H, Viñals J, Gunton JD (1991) Numerical solution of the Swift–Hohenberg equation in two dimensions. Phys A 177:356–365
Yang X, Han D (2017) Linearly first-and second-order, unconditionally energy stable schemes for the phase field crystal model. J Comput Phys 330:1116–1134
Zhao J, Yang X, Li J, Wang Q (2016) Energy stable numerical schemes for a hydrodynamic model of nematic liquid crystals. SIAM J Sci Comput 38(5):3264–3290
Zhao J, Yang X, Shen J, Wang Q (2016) A decoupled energy stable scheme for a hydrogynamic phase-field model of mixtures of nematic liquid crystals and viscous fluids. J Comput Phys 305:539–556
Zouraris GE (2018) An IMEX finite element method for a linearized Cahn–Hilliard Cook equation driven by the space derivative of a space time white noise. Comput Appl Math 37(5):5555–5575
Acknowledgements
J. Su is supported by National Natural Science Foundation of China (no. 91630206, no. 11771348). Y. B. Li is supported by National Natural Science Foundation of China (no. 11601416, no. 11631012) and by the China Postdoctoral Science Foundation (no. 2018M640968).
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Jorge X. Velasco.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Su, J., Fang, W., Yu, Q. et al. Numerical simulation of Swift–Hohenberg equation by the fourth-order compact scheme. Comp. Appl. Math. 38, 54 (2019). https://doi.org/10.1007/s40314-019-0822-8
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s40314-019-0822-8
Keywords
- Swift–Hohenberg equation
- Fourth-order compact scheme
- Nonlinear stabilized splitting scheme
- Adaptive time step method