Hamiltonian boundary value methods applied to KdV-KdV systems

Abstract

In this paper, we propose a highly accurate scheme for two KdV systems of the Boussinesq type under periodic boundary conditions. The proposed scheme combines the Fourier-Galerkin method for spatial discretization with Hamiltonian boundary value methods for time integration, ensuring the conservation of discrete mass and energy. By expanding the system in Fourier series, the equations are firstly transformed into Hamiltonian form, preserving the original Hamiltonian structure. Applying the Fourier-Galerkin method for semi-discretization in space, we obtain a large-scale system of Hamiltonian ordinary differential equations, which is then solved using a class of energy-conserving Runge-Kutta methods, known as Hamiltonian boundary value methods. The efficiency of this approach is assessed, and several numerical examples are provided to demonstrate the effectiveness of the method.