Uniform Error Estimates of a Finite Difference Method for the Klein-Gordon-Schrödinger System in the Nonrelativistic and Massless Limit Regimes
Weizhu Bao,Chunmei Su
DOI: https://doi.org/10.3934/krm.2018040
2018-01-01
Kinetic and Related Models
Abstract:We establish a uniform error estimate of a finite difference method for the Klein-Gordon-Schrodinger (KGS) equations with two dimensionless parameters 0 < gamma <= 1 and 0 < epsilon <= 1, which are the mass ratio and inversely proportional to the speed of light, respectively. in the simultaneously nonrelativistic and massless limit regimes, i.e., gamma similar to epsilon and epsilon -> 0(+), the KGS equations converge singularly to the Schrodinger-Yukawa (SY) equations. When 0 < epsilon << 1, due to the perturbation of the wave operator and/or the incompatibility of the initial data, which is described by two parameters alpha >= 0 and beta >= -1, the solution of the KGS equations oscillates in time with O(epsilon)-wavelength, which requires harsh meshing strategy for classical numerical methods. We propose a uniformly accurate method based on two key points: (i) reformulating KGS system into an asymptotic consistent formulation, and (ii) applying an integral approximation of the oscillatory term. Using the energy method and the limiting equation via the SY equations with an oscillatory potential, we establish two independent error bounds at O(h(2) + tau(2)/epsilon) and O(h(2) + tau(2) + tau epsilon(alpha)* + epsilon(1+alpha)*) with h mesh size, tau time step and alpha* = Min{1,alpha, 1+beta}. This implies that the method converges uniformly and optimally with quadratic convergence rate in space and uniformly in time at O(tau(1/3)) and O(tau(1+alpha*/2+alpha*)) for well-prepared alpha* = 1) and ill-prepared (0 <= alpha* < 1) initial data, respectively. Thus the e-scalability of the method is tau = O(1) and h = O(1) for 0 < epsilon <= 1, which is significantly better than classical methods. Numerical results are reported to confirm our error bounds. Finally, the method is applied to study the convergence rates of KGS equations to its limiting models in the simultaneously nonrelativistic and massless limit regimes.