非线性Schrodinger方程的显式辛算法研究.pdfVIP

摘 要 在这篇论文里，我们用显式辛格式数值求解非线性Schrodinger方程，并 验证其对系统不变量的保持情况。数值结果表明L-L-N分裂方法是模拟非线性 Schrodinger方程的一个又快又好的方法. 首先，我们把非线性Schrodinger方程通过空间离散后，写成一个Hamilton 系统，我们证明这个离散系统的解，当步长充分小时，与原初始方程的解是充 分接近的。初始的非线性Schrodinger方程有一系列的守恒量，我们用中心差分 给出前六个相应的逼近量。 然后，我们发现这个Hamilton函数可以分成三部分，每一部分对应的Hamil- ton系统是可以精确求解的.我们通过这三部分的精确解，给出原Hanmilton 系统的两个显式辛格式S1和S2，它们分别为一阶精度和两阶精度的，而且S2 是时间可逆的.我们分别给出其形式能量展开式的前几项。 最后，对一个孤立子、两个孤立子和三个孤立子情形，我们用这两个格式 分别进行数值模拟，结果表明两阶精度的格式S2可以很好地保持原方程的孤 立子解、各守恒量及辛格式的形式能量. q.t引言 1 引言 我们考虑带初值的非线性Schrodinger方程 (NLSE): iW+W..+川W12W二0 (1) W(x,0)=WO(x) xER,a为常数，W(X,t)是复函数二这个方程出现在许多物理系统的描 述中，例如非线性光学和等离子物理学间[171.对于方程(1)，两种流行的空 间离散是: i州1)+Wp+i〕一2W(1)+W(1-1)+a1W(1)12W(1)=0 (2) h2 W(i+1)一2W(1)+W(I-1) i司1)+ h2 +鲁IWL()I2W(一“，，+W(1+1))一0 3() 这里，h是空间步长，w(l)(t)=W(lh,t),l=一，一1，0，1，… 不同的初值条件WO(x)决定不同的运动.例如，Wf-=。时，为明孤立 子运动[[41;IWo(x)I=P0o，为暗孤立子运动[51;周期的WO(x)导致周期 运动 1【1.在这篇文章里，我们只研究明孤立子的情况. 许多种NLSE方程(1)的数值格式被提出(例如[31，[61，[71，[8]， 191，[101，[ill，[121，[131，[141，[151，[16])，包括一些离散能量和载荷 的守恒格式([3]，[61，[71，[11]，[151，[161)和辛算法([81，[91，[131， [14])。他们或者用空间离散和时间离散同时，或者先空间离散，后时间离散。 在!3}中，Delfour，Fortin和Payre(1981)提出一种有限微分法来近似 比(1)更一般一些的Schrodinger方程 ‘这种方法是非线性隐式的，在空间和 时间上同时离散，这个方法的主要前景是它满足一些重要守恒量的离散分析， 例如能量，载荷等.他们给出了一个孤立子解和两个孤立子解的实验. Sanz-Serna[10](1984)和Guo[6](1986)分别分析了NLSE(1)的一些数值 方法的收敛性. 在这篇文章里我们描述怎么用辛的分裂格式来积分NLSE(1)，这个格式 是完全显式的。 2空间离散 2 空间离散 让我们考虑 (2)，我们有下面的定理: 定理 1:设初值条件we是平方可积的，即 J_-Ivvotx)l-ax十CC 则空间离散方程(2)的解收敛于原始NLSE(a0)的解，当h--30. if明设，(x,t)是方程(1)的解，w(l)(t)=二(th,t)，则有 :、{)+‘,w(t+1)一2w(0+w(1-1)+a}二(L),2w(L)=M