Numerical solutions of Schrödinger–Boussinesq system by orthogonal spline collocation method

Abstract

This paper is concerned with the numerical solutions of Schrödinger–Boussinesq (SBq) system by an orthogonal spline collocation (OSC) discretization in space and Crank–Nicolson (CN) type approximation in time. By using the mathematical induction argument and standard energy method, the proposed CN+OSC scheme is proved to be unconditionally convergent at the order O(τ2+h4) with mesh-size h and time step τ in the discrete L2-norm. We devise a new computation method based on the orthogonal diagonalization techniques (ODT) to realize the proposed CN+OSC scheme. In order to compare the performance of ODT, we devise an alternating direction implicit (ADI) method to compute the CN+OSC scheme for high spatial dimension SBq system. As an alternative implementation, the new method ODT not only exhibits more accurate numerical results, but also demonstrates stronger invariance preserving ability. Numerical results are reported to verify the error estimates and the discrete conservation laws.