Communications in Theoretical Physics ›› 2021, Vol. 73 ›› Issue (2): 025005. doi: 10.1088/1572-9494/abd0e5

• Mathematical Physics • Previous Articles     Next Articles

Numerical simulation of the soliton solutions for a complex modified Korteweg-de Vries equation by a finite difference method

Tao Xu1,2(),Guowei Zhang2,Liqun Wang2,Xiangmin Xu2,Min Li3,   

  1. 1State Key Laboratory of Heavy Oil Processing, China University of Petroleum, Beijing 102249, China
    2College of Science, China University of Petroleum, Beijing 102249, China
    3Department of Mathematics and Physics, North China Electric Power University, Beijing 102206, China
  • Received: 2020-08-22 Revised: 2020-10-13 Accepted: 2020-12-06 Published: 2021-02-01
  • Contact: Tao Xu E-mail:xutao@cup.edu.cn

Abstract:

In this paper, a Crank-Nicolson-type finite difference method is proposed for computing the soliton solutions of a complex modified Korteweg-de Vries (MKdV) equation (which is equivalent to the Sasa-Satsuma equation) with the vanishing boundary condition. It is proved that such a numerical scheme has the second-order accuracy both in space and time, and conserves the mass in the discrete level. Meanwhile, the resulting scheme is shown to be unconditionally stable via the von Nuemann analysis. In addition, an iterative method and the Thomas algorithm are used together to enhance the computational efficiency. In numerical experiments, this method is used to simulate the single-soliton propagation and two-soliton collisions in the complex MKdV equation. The numerical accuracy, mass conservation and linear stability are tested to assess the scheme’s performance.

Key words: complex modified Korteweg-de Vries equation, finite difference method, soliton solutions