Premium
Development of a numerical phase optimized upwinding combined compact difference scheme for solving the C amassa– H olm equation with different initial solitary waves
Author(s) -
Yu C. H.,
Sheu Tony W. H.,
Chang C. H.,
Liao S. J.
Publication year - 2015
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/num.21965
Subject(s) - mathematics , upwind scheme , stencil , camassa–holm equation , helmholtz equation , partial differential equation , mathematical analysis , order of accuracy , finite difference , peakon , ordinary differential equation , finite difference method , compact finite difference , differential equation , method of characteristics , boundary value problem , integrable system , computational science , discretization
In this article, the solution of Camassa–Holm (CH) equation is solved by the proposed two‐step method. In the first step, the sixth‐order spatially accurate upwinding combined compact difference scheme with minimized phase error is developed in a stencil of four points to approximate the first‐order derivative term. For the purpose of retaining both of the long‐term accurate Hamiltonian property and the geometric structure inherited in the CH equation, the time integrator used in this study should be able to conserve symplecticity. In the second step, the Helmholtz equation governing the pressure‐like variable is approximated by the sixth‐order accurate three‐point centered compact difference scheme. Through the fundamental and numerical verification studies, the integrity of the proposed high‐order scheme is demonstrated. Another aim of this study is to reveal the wave propagation nature for the investigated shallow water equation subject to different initial wave profiles, whose peaks take the smooth, peakon, and cuspon forms. The transport phenomena for the cases with/without inclusion of the linear first‐order advection term κu x in the CH equation will be addressed. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1645–1664, 2015