UNIFORM ERROR ESTIMATES OF FINITE DIFFERENCE METHODS FOR THE NONLINEAR SCHRODINGER EQUATION WITH WAVE OPERATOR

被引:174
|
作者
Bao, Weizhu [1 ,2 ]
Cai, Yongyong [1 ]
机构
[1] Natl Univ Singapore, Dept Math, Singapore 119076, Singapore
[2] Natl Univ Singapore, Ctr Computat Sci & Engn, Singapore 119076, Singapore
关键词
nonlinear Schrodinger equation with wave operator; energy method; error estimates; conservative Crank-Nicolson finite difference method; semi-implicit finite difference method; MODELING LIGHT BULLETS; KLEIN-GORDON; SPLITTING METHODS; NONRELATIVISTIC LIMIT; NUMERICAL-SOLUTION; PRESERVING SCHEME; SINE-GORDON; SPACE; ENERGY; CONVERGENCE;
D O I
10.1137/110830800
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We establish uniform error estimates of finite difference methods for the nonlinear Schrodinger equation (NLS) perturbed by the wave operator (NLSW) with a perturbation strength described by a dimensionless parameter epsilon (epsilon is an element of (0, 1]). When epsilon -> 0(+), NLSW collapses to the standard NLS. In the small perturbation parameter regime, i.e., 0 < epsilon << 1, the solution of NLSW is perturbed from that of NLS with a function oscillating in time with O(epsilon(2))-wavelength at O(epsilon(4)) and O(epsilon(2)) amplitudes for well-prepared and ill-prepared initial data, respectively. This high oscillation of the solution in time brings significant difficulties in establishing error estimates uniformly in epsilon of the standard finite difference methods for NLSW, such as the conservative Crank-Nicolson finite difference (CNFD) method, and the semi-implicit finite difference (SIFD) method. We obtain error bounds uniformly in epsilon, at the order of O(h(2) + tau) and O(h(2) + tau(2/3)) with time step tau and mesh size h for well-prepared and ill-prepared initial data, respectively, for both CNFD and SIFD in the l(2)-norm and discrete semi-H-1 norm. Our error bounds are valid for general nonlinearity in NLSW and for one, two, and three dimensions. To derive these uniform error bounds, we combine epsilon-dependent error estimates of NLSW, epsilon-dependent error bounds between the numerical approximate solutions of NLSW and the solution of NLS, together with error bounds between the solutions of NLSW and NLS. Other key techniques in the analysis include the energy method, cut-off of the nonlinearity, and a posterior bound of the numerical solutions by using the inverse inequality and discrete semi-H-1 norm estimate. Finally, numerical results are reported to confirm our error estimates of the numerical methods and show that the convergence rates are sharp in the respective parameter regimes.
引用
收藏
页码:492 / 521
页数:30
相关论文
共 50 条