带有阻尼项的4阶非线性薛定谔方程的显式辛格式.pdfVIP

第 37卷 第 3期 江西师范大学学报 (自然科学版) Vo1．37No．3 2013年5月 JournalofJiangxiNormalUniversity(NaturalS‘cience) Mav2013 文章编号：1000-5862(2013)03-0244-05 带有阻尼项的4阶非线性薛定谔方程的显式辛格式 徐 远 ，孔令华 ，王 兰，黄晓梅 (江西师范大学数学与信息科学学院，江西南昌 330022) 摘要 ：把带有阻尼项的4阶薛定谔方程写成标准的哈密尔顿系统，将该哈密尔顿系统分裂成2个哈密尔 顿子系统．一个子系统是可分的，可以构造显式的辛格式；而另一个子系统由点点的质量守恒可以精确求 解．这样得到的数值格式整体上是辛格式，而且避免了通常辛格式需要迭代的弊端，提高了计算效率． 关键词：4阶非线性薛定谔方程；显式辛格式；哈密尔顿系统 中图分类号：0241．8 文献标志码 ：A 于JM(，t)l∈R的有界实函数，v(x)是有界实函 0 引言 数，称之为势函数．u(x，t)是一个复数值函数， M。()是一个给定的复数值函数．经过简单的微积 自1984年在北京举行的国际 “双微 ”会议上冯 分计算可知，该定解问题具有质量守恒律 L 康院士首先系统地提出哈密尔顿系统的辛几何算法 r Q(f)：Jlu(x，t)ldx=Jlu(x，0)l。dx= 以来 J，辛几何算法得到了长足发展，构造出了多 U U 种数值方法，如生成函数法、(块)Runge—Kutta(RK) Q(0)， (4) 方法等 引．现在不仅有哈密尔顿常微分方程的辛 和能量守恒律 算法，而且已发展到哈密尔顿偏微分方程多辛算 E(f)=I(I I+G(I(，t)l)+ 法_4引．由RK方法需满足bi6—bia一66=0这一 U 条件才称为辛格式，易知显式的RK方法不可能为 ()lM(，t)l。)dx=E(0)． (5) 辛的．对于可分的哈密尔顿系统，可以构造出显式的 对于周期初边值问题 (1)～(3)的数值解，已经有一 块RK辛算法．这样对于不可分的哈密尔顿系统只 些研究，主要包含完全显式的蛙跳格式，但是此格式 能建立完全隐式的辛格式，这将大大增加计算量和 不是 自启动的，不保持辛守恒，同时也不保持守恒律 存储量，使得格式的计算效率不高．基于辛映射的复 (4)式或(5)式．完全隐式的辛格式能够保持辛守恒 合还是辛映射的理论 ，可以把不可分的哈密尔顿系 但是计算效率不高．为此，在辛算法中引入分裂步方 统分裂成一些可分的哈密尔顿子系统，然后为各个 法的基本思想：把原问题分裂成一些更简单的子问 子系统构造辛格式，最后进行复合运算．基于上述原 题，通过对子问题依次数值求解或者解析求解来达 理，本文考虑在量子力学、非线性光学中有重要应用 到求解原问题的目的． 的带有阻尼项的4阶非线性薛定谔 (FNSLET)方程 的周期初边值问题 1 Ej1『方程哈密尔顿形式 i『+ m-0 令 u(x，t)=p(x，t)+iq