Abstract
In this paper, we present a linearly implicit energy-preserving scheme for the Camassa–Holm equation by using the multiple scalar auxiliary variables approach, which is first developed to construct efficient and robust energy stable schemes for gradient systems. The Camassa–Holm equation is first reformulated into an equivalent system by utilizing the multiple scalar auxiliary variables approach, which inherits a modified energy. Then, the system is discretized in space aided by the standard Fourier pseudo-spectral method and a semi-discrete system is obtained, which is proven to preserve a semi-discrete modified energy. Subsequently, the linearized Crank–Nicolson method is applied for the resulting semi-discrete system to arrive at a fully discrete scheme. The main feature of the new scheme is to form a linear system with a constant coefficient matrix at each time step and produce numerical solutions along which the modified energy is precisely conserved, as is the case with the analytical solution. Several numerical results are addressed to confirm accuracy and efficiency of the proposed scheme.
Similar content being viewed by others
References
Cai, W., Jiang, C., Wang, Y., Song, Y.: Structure-preserving algorithms for the two-dimensional sine-Gordon equation with Neumann boundary conditions. J. Comput. Phys. 395, 166–185 (2019)
Camassa, R., Holm, D.: An integrable shallow water equation with peaked solitons. Phys. Rev. Lett. 71, 1661–1664 (1993)
Camassa, R., Holm, D., Hyman, J.: A new integrable shallow water equation. Adv. Appl. Mech. 31, 1–33 (1994)
Chen, J., Qin, M.: Multi-symplectic Fourier pseudospectral method for the nonlinear Schrödinger equation. Electron. Trans. Numer. Anal. 12, 193–204 (2001)
Cheng, Q., Shen, J.: Multiple scalar auxiliary variable (MSAV) approach and its application to the phase-field vesicle membrane model. SIAM J. Sci. Comput. 40, A3982–A4006 (2018)
Cohen, D., Raynaud, X.: Geometric finite difference schemes for the generalized hyperelastic-rod wave equation. J. Comput. Appl. Math. 235, 1925–1940 (2011)
Constantin, A.: On the scattering problem for the Camassa–Holm equation. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 457, 953–970 (2001)
Dahlby, M., Owren, B.: A general framework for deriving integral preserving numerical methods for PDEs. SIAM J. Sci. Comput. 33, 2318–2340 (2011)
Eidnes, S., Li, L., Sato, S.: Linearly implicit structure-preserving schemes for Hamiltonian systems (2019). arXiv preprint arXiv:1901.03573
Gong, Y., Cai, J., Wang, Y.: Multi-symplectic Fourier pseudospectral method for the Kawahara equation. Commun. Comput. Phys. 16, 35–55 (2014)
Gong, Y., Wang, Y.: An energy-preserving wavelet collocation method for general multi-symplectic formulations of Hamiltonian PDEs. Commun. Comput. Phys. 20, 1313–1339 (2016)
Gong, Y., Zhao, J., Yang, X., Wang, Q.: Fully discrete second-order linear schemes for hydrodynamic phase field models of binary viscous fluid flows with variable densities. SIAM J. Sci. Comput. 40, B138–B167 (2018)
Holden, H., Raynaud, X.: Convergence of a finite difference scheme for the Camassa–Holm equation. SIAM J. Numer. Anal. 44, 1655–1680 (2006)
Hong, Q., Gong, Y., Lv, Z.: Linear and Hamiltonian-conserving Fourier pseudo-spectral schemes for the Camassa–Holm equation. Appl. Math. Comput. 346, 86–95 (2019)
Jiang, C., Wang, Y., Gong, Y.: Arbitrarily high-order energy-preserving schemes for the Camassa–Holm equation. Appl. Numer. Math. 151, 85–97 (2020)
Kalisch, H., Lenells, J.: Numerical study of traveling-wave solutions for the Camassa–Holm equation. Chaos Solitons Fractals 25, 287–298 (2005)
Lenells, J.: Traveling wave solutions of the Camassa–Holm equation. J. Differ. Equ. 271, 393–430 (2005)
Matsuo, T., Yamaguchi, H.: An energy-conserving Galerkin scheme for a class of nonlinear dispersive equations. J. Comput. Phys. 228, 4346–4358 (2009)
Miyatake, Y., Matsuo, T.: Energy-preserving \(H^1\)-Galerkin schemes for shallow water wave equations with peakon solutions. Phys. Lett. A 376, 2633–2639 (2012)
Shen, J., Tang, T.: Spectral and High-Order Methods with Applications. Science Press, Beijing (2006)
Shen, J., Xu, J., Yang, J.: The scalar auxiliary variable (SAV) approach for gradient. J. Comput. Phys. 353, 407–416 (2018)
Shen, J., Xu, J., Yang, J.: A new class of efficient and robust energy stable schemes for gradient flows. SIAM Rev. 61, 474–506 (2019)
Xu, Y., Shu, C.-W.: A local discontinuous Galerkin method for the Camassa–Holm equation. SIAM J. Numer. Anal. 46, 1998–2021 (2008)
Yang, X., Zhao, J., Wang, Q.: Numerical approximations for the molecular beam epitaxial growth model based on the invariant energy quadratization method. J. Comput. Phys. 333, 104–127 (2017)
Zhao, J., Yang, X., Gong, Y., Wang, Q.: A novel linear second order unconditionally energy stable scheme for a hydrodynamic-tensor model of liquid crystals. Comput. Methods Appl. Mech. Eng. 318, 803–825 (2017)
Zhu, H., Song, S., Tang, Y.: Multi-symplectic wavelet collocation method for the Schrödinger equation and the Camassa–Holm equation. Comput. Phys. Commun. 182, 616–627 (2011)
Acknowledgements
The authors would like to express sincere gratitude to the referees for their insightful comments and suggestions. Chaolong Jiang’s work is partially supported by the National Natural Science Foundation of China (Grant No. 11901513), the Yunnan Provincial Department of Education Science Research Fund Project (Grant No. 2019J0956) and the Science and Technology Innovation Team on Applied Mathematics in Universities of Yunnan. Yuezheng Gong’s work is partially supported by the Natural Science Foundation of Jiangsu Province (Grant No. BK20180413), the National Natural Science Foundation of China (Grant No. 11801269) and the Foundation of Jiangsu Key Laboratory for Numerical Simulation of Large Scale Complex Systems (Grant No. 202002). Wenjun Cai’s work is partially supported by the National Natural Science Foundation of China (Grant No. 11971242) and the National Key Research and Development Project of China (Grant Nos. 2018YFC0603500, 2018YFC1504205). Yushun Wang’s work is partially supported by the National Natural Science Foundation of China (Grant No. 11771213).
Author information
Authors and Affiliations
Corresponding author
Additional information
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
Jiang, C., Gong, Y., Cai, W. et al. A Linearly Implicit Structure-Preserving Scheme for the Camassa–Holm Equation Based on Multiple Scalar Auxiliary Variables Approach. J Sci Comput 83, 20 (2020). https://doi.org/10.1007/s10915-020-01201-4
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10915-020-01201-4
Keywords
- Multiple scalar auxiliary variables approach
- Linearly implicit scheme
- Energy-preserving scheme
- Camassa–Holm equation