In this article, we construct an efficient numerical scheme to solve the coupled nonlinear Klein–Gordon system in a multi-dimensional unbounded domain Rd(d=1,2,3). We first use the SAV method to equivalently deform the original system, then we treat the nonlinear part of the system explicitly in temporal discretization. Thanks to the SAV method, we ensure the unconditional energy conservation of the discrete system, as a cost, we only need to solve two completely linearized decoupled systems. To avoid the errors caused by domain truncation and the construction of artificial boundary conditions, we apply the Hermite–Galerkin spectral method with scaling factor for spatial discretization. In order to save computational costs, we utilize two time adaptive strategies to adjust the time step size, through which we can improve computational efficiency without sacrificing accuracy. Finally, we presented numerical experiments to demonstrate the accuracy, efficiency, and energy conservation properties of our algorithm, and applied it to simulate the interactions of two-dimensional and three-dimensional vector solitons.
Read full abstract