Abstract
This article presents a two-step iterative method that uses u – P formulation to study a Degasperis-Procesi (DP) equation, with which the DP equation is decomposed into the nonlinear advection equation and the Helmholtz equation. The first-order derivative terms in the advection equation are approximated by an Optimized Compact Reconstruction Weighted Essentially Non-Oscillatory (OCRWENO) scheme that reduces dispersion and dissipation errors and suppresses oscillation near discontinuities. Besides, high-order symplectic Runge-Kutta scheme is used for time marching. A rigorous analysis of the dispersion and dissipation errors are provided for the OCRWENO scheme. Single smooth soliton solution is investigated to check the accuracy and the order of the proposed method. Peakon, peakon-peakon, peakon-antipeakon and shockpeakon solutions of the DP equation are then predicted. Finally, wave breaking phenomena of smoothed initial condition from the DP equation are addressed. In addition, two conservative quantities associated with the bi-Hamiltonian form of the DP equation are calculated to demonstrate capabilities of the present method.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call