A class of conservative orthogonal spline collocation schemes for solving coupled Klein-Gordon-Schrödinger equations (Q2378991)

The paper considers one-dimensional coupled Klein-Gordon-Schrödinger equations with homogeneous Dirichlet boundary conditions and smooth initial conditions. The scheme is based on a weighted implicit three-level finite difference approximation in time and orthogonal cubic spline collocation in space. Two main results are established in the paper: (1) the proposed scheme has two invariants (in time) that are discrete analogs of the corresponding quantities of the original system of equations; (2) the scheme converges with a rate \(O(\tau^2+ h^4)\) (\(\tau\) the time step-size and \(h\) the spatial mesh-size) in \(L^2\)-norm in space and maximum norm in time. Numerical experiments showing propagation of solitary waves and collisions of two head-on solitary waves are presented and discussed.