非线性薛定谔方程的平均向量场方法非线性薛定谔方程的平均向量场方法.pdfVIP

2013年 2月 计 算 数 学 第 35卷第 1期 Feb．，2013 MATHEMATICA NUMERICA SINICA Vo1．35．No．1 非线性薛定谔方程的平均向量场方法 ) 李吴辰 孙建强 骆思宇 (海南大学信息科学技术学院数学系，海口570228) 摘 要 利用平均向量场方法 (AVF)对非线性薛定谔方程进行求解，在理论上得到了一个保非线性 薛定谔方程描述的系统能量守恒的AVF格式，再分别用非线性薛定谔方程的AVF格式和辛格 式数值模拟孤立波的演化行为，并比较两个格式是否保系统能量守恒特性．数值结果表明，AVF 格式也能很好地模拟孤立波的演化行为，并且比辛格式更能保持系统的能量守恒． 关键词：AVF格式；保能量算法；非线性薛定谔方程 MR (2000)主题分类：65D17 1．引 言 在数学和物理中有一大类微分方程，微分方程所描述的系统具有能量守恒的特性．在数值 模拟中，设计精确地保持系统能量守恒特性的数值格式对准确地模拟系统的行为具有重要的 意义 l【12]．早在上个世纪八十年代，冯康及其研究小组提出了Hamilton系统的辛几何算法， 能长时间精确地模拟系统的演化行为，并能近似地保持系统的能量守恒特性 3【I4J．Bridges和 Reich等人在辛几何算法的基础上提出了偏微分方程的多辛算法，能近似地保持系统的整体 和局部能量守恒特性 5【)6]．近年来，Quispel和 McLachlan等人提出了精确地保Hamilton系统 的平均向量场方法 [7-10】．唐贻发和洪佳林等人在平均向量场方法的基础上提出了精确地保 多辛结构偏微分方程整体能量守恒的数值格式 [11,12】．在这里我们直接利用平均向量场方法 和辛算法求解具有能量守恒特性的非线性薛定谔方程，并且比较两个格式的保能量守恒特性． 非线性薛定谔方程是孤立子理论中最重要的可积模型，广泛应用在物理中的许多领域，包 括非线性光学和等离子体物理 引【．设非线性薛定谔方程 初始条件 u(x，0)=uo(x)，X∈R， 0是常数．方程 (1)具有保持系统能量守恒特性 E：厂。。(1Iz一导luI)d：Eo， (2) J—o。 其中 为常数 本文首先介绍了平均向量场方法的基本思想，随后利用平均向量场方法构造了薛定谔方 程 (1)的能量守恒格式，最后分别用平均向量场方法和辛算法对薛定谔方程求解，比较两种方 法是否保薛定谔方程的整体能量守恒特性． 2012年 6月 25日收到． )基金项目：本项 目获国家 自然科学基金项 目、海南大学科研启动基金项 目(kyqd1053)和海南 大学青年基金项 目(qnjj1022)资助． 计 算 数 学 2．平均向量场方法 在这里我们简单介绍下平均向量场方法及其保系统能量守恒特性 [7--12]．给定常微分方 程 )， (0) z∈R“- (3) 假设 H：R R是一个不变量，或者是系统的能量，即 _01VH( dz (4) 系统 (3)常常可以写为 dz = S(z) ㈤ ： (矿 n ． (5) 系统 (5)又称为斜梯度系统，即S ：一S． 反之，假设存在反对称矩阵S(z)，使系统 (3)斜对称 f(z)= (z)H(z)，于是