FULLY DECOUPLED ARBITRARILY HIGH-ORDER CONSERVATIVE SCHEMES FOR THE SCHRÖDINGER–BOUSSINESQ EQUATION

Abstract A family of arbitrarily high-order energy-preserving methods are developed to solve the coupled Schrödinger–Boussinesq (S-B) system. The system is a nonlinear coupled system and satisfies a series of conservation laws. It is often difficult to construct a high-order decoupling numerical algorithm to solve the nonlinear system. In this paper, the original system is first reformulated into an equivalent Hamiltonian system by introducing multiple auxiliary variables. Next, the reformulated system is discretized by the Fourier pseudo-spectral method and the implicit midpoint scheme in the spatial and temporal directions, respectively, and a second-order conservative scheme is obtained. Finally, the scheme is extended to arbitrarily high-order accuracy by means of diagonally implicit symplectic Runge–Kutta methods or composition methods. Rigorous analyses show that the proposed methods are fully decoupled and can precisely conserve the discrete invariants. Numerical results show that the proposed schemes are effective and can be easily extended to other nonlinear partial differential equations.