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

被引:179
作者
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
相关论文
共 28 条
[1]   ON FULLY DISCRETE GALERKIN METHODS OF 2ND-ORDER TEMPORAL ACCURACY FOR THE NONLINEAR SCHRODINGER-EQUATION [J].
AKRIVIS, GD ;
DOUGALIS, VA ;
KARAKASHIAN, OA .
NUMERISCHE MATHEMATIK, 1991, 59 (01) :31-53
[2]   FINITE-DIFFERENCE DISCRETIZATION OF THE CUBIC SCHRODINGER-EQUATION [J].
AKRIVIS, GD .
IMA JOURNAL OF NUMERICAL ANALYSIS, 1993, 13 (01) :115-124
[3]  
Bao W., MATH COMP IN PRESS
[4]   Analysis and comparison of numerical methods for the Klein-Gordon equation in the nonrelativistic limit regime [J].
Bao, Weizhu ;
Dong, Xuanchun .
NUMERISCHE MATHEMATIK, 2012, 120 (02) :189-229
[5]   Comparisons between sine-Gordon and perturbed nonlinear Schrodinger equations for modeling light bullets beyond critical collapse [J].
Bao, Weizhu ;
Dong, Xuanchun ;
Xin, Jack .
PHYSICA D-NONLINEAR PHENOMENA, 2010, 239 (13) :1120-1134
[6]   Numerical solution of the Gross-Pitaevskii equation for Bose-Einstein condensation [J].
Bao, WZ ;
Jaksch, D ;
Markowich, PA .
JOURNAL OF COMPUTATIONAL PHYSICS, 2003, 187 (01) :318-342
[7]   A SINGULAR PERTURBATION PROBLEM FOR AN ENVELOPE EQUATION IN PLASMA PHYSICS [J].
BERGE, L ;
COLIN, T .
PHYSICA D-NONLINEAR PHENOMENA, 1995, 84 (3-4) :437-459
[8]   Order estimates in time of splitting methods for the nonlinear Schrodinger equation [J].
Besse, C ;
Bidégaray, B ;
Descombes, E .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2002, 40 (01) :26-40
[9]  
CHANG QS, 1995, MATH COMPUT, V64, P537, DOI 10.1090/S0025-5718-1995-1284664-5
[10]  
Colin T., 1998, Discret. Contin. Dyn. Syst, V4, P671