Numerical simulation and analysis of the Swift-Hohenberg equation by the stabilized Lagrange multiplier approach. (English) Zbl 07453282

Summary: In this work, we develop linear, first- and second-order accurate, and energy-stable numerical scheme for the Swift-Hohenberg (SH) equation. An auxiliary variable (Lagrange multiplier) is used to control the nonlinear term so that the linear temporal scheme can be easily constructed. To further achieve the accuracy with large time steps, a proper stabilized term is adopted. For the first-order time-accurate scheme, the backward Euler approximation is adopted. The Crank-Nicolson (CN) and explicit Adams-Bashforth (AB) approximations are applied to achieve temporally second-order accuracy. We analytically perform the estimations of the semi-discrete solvability, the energy stability with respect to the original and pseudo-energy functionals, and the convergence error. To numerically solve the resulting discrete system of equations, we use an efficient linear multigrid method. We present various two- (2D) and three-dimensional (3D) computational examples to demonstrate the accuracy and energy stability of the proposed scheme.


34D20 Stability of solutions to ordinary differential equations
37M05 Simulation of dynamical systems
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
65N22 Numerical solution of discretized equations for boundary value problems involving PDEs
Full Text: DOI


[1] Abbaszadeh, M.; Khodadadian, A.; Parvizi, M.; Dehghan, M.; Heitzinger, C., A direct meshless local collocation method for solving stochastic Cahn-Hilliard-Cook and stochastic Swift-Hohenberg equations, Eng Anal Bound Elem, 98, 253-264 (2019) · Zbl 1404.65207
[2] Christov, CI; Pontes, J., Numerical scheme for Swift-Hohenberg equation with strict implementation of Lyapunov functional, Math Comput Model, 35, 87-99 (2002) · Zbl 1197.65120
[3] Church, JM; Guo, Z.; Jimack, PK; Madzvamuse, A.; Promislow, K.; Wetton, B.; Wise, SM; Yang, F., High accuracy benchmark problems for Allen-Cahn and Cahn-Hilliard dynamics, Commun Comput Phys, 26, 947-972 (2019) · Zbl 1473.65099
[4] Dehghan, M.; Abbaszadeh, M., 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 (2017) · Zbl 1403.74285
[5] Dong, L.; Wang, C.; Zhang, H.; Zhang, Z., A positivity-preserving second-order BDF scheme for the Cahn-Hilliard equation with variable interfacial parameters, Commun Comput Phys, 28, 967-998 (2020) · Zbl 07419121
[6] Elder, KR; Viñals, J.; Grant, M., Ordering dynamics in the two-dimensional stochastic Swift-Hohenberg equation, Phys Rev Lett, 68, 3024-3027 (1992)
[7] Feng, W.; Guo, Z.; Lowengrub, JS; Wise, SM, A mass-conservative adaptive FAS multigrid solver for cell-centered finite difference methods on block-structured locally-cartesian grids, J Comput Phys, 352, 463-497 (2018) · Zbl 1380.65402
[8] Gao, Y.; Li, R.; Mei, L.; Lin, Y., A second-order decoupled energy stable numerical scheme for Cahn-Hilliard-Hele-Shaw system, Appl Numer Math, 157, 338-355 (2020) · Zbl 1446.65110
[9] Gomez, H.; Nogueira, X., A new space-time discretization for the Swift-Hohenberg equation that strictly respects the Lyapunov functional, Commun Nonlinear Sci Numer Simul, 17, 4930-4946 (2012) · Zbl 1352.76109
[10] Guillén-González, F.; Tierra, G., On linear schemes for the Cahn-Hilliard diffuse interface model, J Comput Phys, 234, 140-171 (2013) · Zbl 1284.35025
[11] Haken, H., Advanced synergetics (1983), Berlin: Springer, Berlin · Zbl 0521.93002
[12] Hutt, A.; Atay, FM, Analysis nonlocal neural fields for both general and gamma-distributed connectivities, Physica D, 203, 30-54 (2005) · Zbl 1061.92020
[13] Kim, J.; Kang, K.; Lowengrub, J., Conservative multigrid methods for the Cahn-Hilliard fluids, J Comput Phys, 193, 511-543 (2004) · Zbl 1109.76348
[14] Kudryashov, NA; Sinelshchikov, DI, Exact solutions of the Swift-Hohenberg equation with dispersion, Commun Nonlinear Sci Numer Simul, 17, 26-34 (2012) · Zbl 1245.35095
[15] Lee, HG, A semi-analytical Fourier spectral method for the Swift-Hohenberg equation, Comput Math Appl, 74, 1885-1896 (2017) · Zbl 1397.65205
[16] Lee, HG, An energy stable method for the Swift-Hohenberg equation with quadratic-cubic nonlinearity, Comput Methods Appl Mech Engrg, 343, 40-51 (2019) · Zbl 1440.65150
[17] Li, H.; Ju, L.; Zhang, C.; Peng, Q., Unconditionally energy stable linear schemes for the diffuse interface method with Peng-Robinson equation of state, J Sci Comput, 75, 993-1015 (2018) · Zbl 1398.65305
[18] Li, N.; Lin, P.; Gao, G., Energy law preserving finite element scheme for the Cahn-Hilliard equation with dynamic boundary conditions, Commun Comput Phys, 26, 1490-1509 (2019) · Zbl 1474.76039
[19] Liu, H.; Yin, P., Unconditionally energy stable DG scheme for the Swift-Hohenberg equation, J Sci Comput, 81, 789-819 (2019) · Zbl 1427.65372
[20] Nikolay, NA; Ryabov, PN, Analytical and numerical solutions of the generalized dispersive Swift-Hohenberg equation, Appl Math Comput, 286, 171-177 (2016) · Zbl 1410.35174
[21] Sarmiento, AF; Espath, LFR; Vignal, P.; Dalcin, L.; Parsani, M.; Calo, VM, An energy-stable generalized-\( \alpha\) method for the Swift-Hohenberg equation, J Comput Appl Math, 344, 836-851 (2018) · Zbl 1457.65080
[22] Shivanian, E.; Abbasbandy, S., Pseudospectral meshless radial point interpolation for generalized biharmonic equation in the presence of Cahn-Hilliard conditions, Comput Appl Math, 39, 148 (2020) · Zbl 1449.65336
[23] Staliunas, K.; Sánchez-Morcillo, VJ, Dynamics of phase domain in the Swift-Hohenberg equation, Phys Lett A, 241, 28-34 (1998)
[24] Su, J.; Fang, W.; Yu, Q.; Li, Y., Numerical simulation of Swift-Hohenberg equation by the fourth-order compact scheme, Comput Appl Math, 38, 54 (2019) · Zbl 1438.65192
[25] Swift, J.; Hohenberg, PC, Hydrodynamic fluctuation at the convective instability, Phys Rev A, 15, 319-328 (1977)
[26] Viñals, J.; Hernández-Garca, E.; San, MM; Toral, R., Numerical study of the dynamical aspects of pattern selection in the stochastic Swift-Hohenberg equation in one dimension, Phys Rev A, 44, 1123-1133 (1991)
[27] Wise, S.; Kim, J.; Lowengrub, J., Solving the regularized, strongly anisotropic Cahn-Hilliard equation by an adaptive nonlinear multigrid method, J Comput Phys, 226, 414-446 (2007) · Zbl 1310.82044
[28] Xi, H.; Viñals, J.; Gunton, JD, Numerical solution of the Swift-Hohenberg equation in two dimensions, Phys A, 177, 356-365 (1991)
[29] Xu, Z.; Yang, X.; Zhang, H.; Xie, Z., Efficient and linear schemes for anisotropic Cahn-Hilliard model using the Stabilized-Invariant Energy Quadratization (S-IEQ) approach, Comput Phys Commun, 238, 36-49 (2019)
[30] Yang, X., Linear, first and second order and unconditionally energy stable numerical schemes for the phase field model of homopolymer blends, J Comput Phys, 327, 294-316 (2016) · Zbl 1373.82106
[31] Yang, J.; Kim, J., An unconditionally stable second-order accurate method for systems of Cahn-Hilliard equations, Commun Nonlinear Sci Numer Simulat, 87 (2020) · Zbl 1455.76138
[32] Yang, X.; Zhang, G., Convergence analysis for the Invariant Energy Quadratization (IEQ) schemes for solving the Cahn-Hilliard and Allen-Cahn equations with general nonlinear potential, J Comput Sci, 82, 55 (2020) · Zbl 1434.65222
[33] Zhang, J.; Yang, X., Numerical approximations for a new \(L^2\)-gradient flow based phase field crystal model with precise nonlocal mass conservation, Comput Phys Commun, 243, 51-67 (2019)
[34] Zhu, G.; Kou, J.; Sun, S.; Yao, J.; Li, A., Numerical approximation of a phase-field surfactant model with fluid flow, J Sci Comput, 80, 223-247 (2019) · Zbl 1448.76171
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.