An SDG Galerkin structure-preserving scheme for the Klein-Gordon-Schrodinger equation

摘要

In this paper, we use the Galerkin weak form to construct a structure-preserving scheme for Klein-Gordon-Schrodinger equation and analyze its conservative and convergent properties. We first discretize the underlying equation in space direction via a selected finite element method, and the Hamiltonian partial differential equation can be casted into Hamiltonian ordinary differential equations based on the weak form of the system afterwards. Then, the resulted ordinary differential equations are solved by the symmetric discrete gradient method, which yields a charge-preserving and energy-preserving scheme. Moreover, the numerical solution of the proposed scheme is proved to be bounded in the discrete L-infinity norm and convergentwith the convergence order of O(h(2)+ tau(2)) in the discrete L-2 norm without any grid ratio restrictions, where h and tau are space and time step, respectively. Numerical experiments conducted last to verify the theoretical analysis.