Numerical solutions to a two-component Camassa–Holm equation
作者:
Highlights:
•
摘要
In the present paper, a three-step iterative algorithm for solving a two-component Camassa–Holm (2CH) equation is presented. In the first step, the time-dependent equation for the horizontal fluid velocity with nonlinear convection is solved. Then an inhomogeneous Helmholtz equation is solved. Finally, the equation for modeling the transport of density is solved in the third step. The differential order of 2CH equation has been reduced in order to facilitate numerical scheme development in a comparatively smaller grid stencil. In this study, a fifth-order spatially accurate upwinding combined compact difference scheme (UCCD5) which differs from that in Sheu et al. (2011) is developed in a four-point grid stencil for approximating the first-order derivative term. For the purpose of retaining long-time Hamiltonians in the 2CH equation, the time integrator (or time-stepping scheme) chosen is symplectic. Various numerical experiments such as the single peakon, peakon–antipeakon interaction and dam-break problems are conducted to illustrate the effectiveness of the proposed numerical method. It is shown that both the Hamiltonians and Casimir functions are conserved well for all problems.
论文关键词:Two-component Camassa–Holm equation,Inhomogeneous Helmholtz equation,Combined compact difference scheme,peakon–antipeakon,Hamiltonians,Casimir function
论文评审过程:Received 20 November 2015, Revised 27 May 2016, Available online 8 January 2018, Version of Record 5 February 2018.
论文官网地址:https://doi.org/10.1016/j.cam.2017.12.043