Skip to main content
SpringerLink
Log in

Numerical simulation of Swift–Hohenberg equation by the fourth-order compact scheme

  • Published:
Computational and Applied Mathematics Aims and scope Submit manuscript

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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Log in via an institution

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11
Fig. 12

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

    Article  MathSciNet  Google Scholar 

  • Cheng M, Warren JA (2008) An efficient algorithm for solving the phase field crystal model. J Comput Phys 227:6241–6248

    Article  MathSciNet  Google Scholar 

  • Christov CI, Pontes J (2002) Numerical scheme for Swift–Hohenberg equation with strict implementation of lyapunov functional. Math Comput Model 35:87–99

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • Cross M, Greenside H (2009) Pattern formation and dynamics in nonequilibrium systems. Cambridge University Press, Cambridge

    Book  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • Elder KR, Grant M (2004) Modeling elastic and plastic deformations in nonequilibrium processing using phase field crystals. Phys Rev E 70:051605–051623

    Article  Google Scholar 

  • Elder KR, Viñals J, Grant M (1992) Ordering dynamics in the two-dimensional stochastic Swift–Hohenberg equation. Phys Rev Lett 68:3024–3027

    Article  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • Kim JS, Kang K, Lowengrub J (2004) Conservative multigrid methods for Cahn–Hilliard fluids. J Comput Phys 193:511–543

    Article  MathSciNet  Google Scholar 

  • Lee HG (2017) A semi-analytical Fourier spectral method for the Swift–Hohenberg equation. Comput Math Appl 74:1885–1896

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • Li YB, Choi YH, Kim JS (2017) Computationally efficient adaptive time step method for the Cahn–Hilliard equation. Comput Math Appl 73:1855–1864

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • Marconi UMB, Tarazona P (1999) Dynamic density functional theory of liquids. J Chem Phys 110:8032–8044

    Article  Google Scholar 

  • 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

    MATH  Google Scholar 

  • Mohammadi V, Mohebbi A, Asgari Z (2009) Fourth order compact solution of the nonlinear Klein–Gordon equation. Numer Algorithms 52(4):523–540

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • Nikolay NA, Ryabov PN (2016) Analytical and numerical solutions of the generalized dispersive Swift–Hohenberg equation. Appl Math Comput 286:171–177

    MathSciNet  MATH  Google Scholar 

  • Staliunas K, Sánchez-Morcillo VJ (1998) Dynamics of phase domains in the Swift–Hohenberg equation. Phys Lett A 241:28–34

    Article  Google Scholar 

  • Swift J, Hohenberg PC (1977) Hydrodyamic fluctuations at the convective instability. Phys Rev A 15:319–328

    Article  Google Scholar 

  • Trottenberg U, Oosterlee C, Schüller A (2001) Multigrid. Academic Press, New York

    MATH  Google Scholar 

  • 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

    Article  Google Scholar 

  • 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

    Article  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • Xi H, Viñals J, Gunton JD (1991) Numerical solution of the Swift–Hohenberg equation in two dimensions. Phys A 177:356–365

    Article  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

  • 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

    Article  MathSciNet  Google Scholar 

Download references

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

  1. School of Mathematics and Statistics, Xi’an Jiaotong University, Xi’an, 710049, China

    Jian Su, Qian Yu & Yibao Li

  2. Department of Computational Science and Engineering, Yonsei University, Seoul, 03722, Republic of Korea

    Weiwei Fang

Authors
  1. Jian Su
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Weiwei Fang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Qian Yu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yibao Li
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Yibao Li.

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1007/s40314-019-0822-8

Keywords

Mathematics Subject Classification

Search

Navigation