论文部分内容阅读
在数学物理的研究中,许多偏微分方程可以表示成哈密尔顿系统的辛结构形式或多辛结构形式,如复修正KdV方程,KGS方程,耦合薛定谔Boussinesq方程,BBM方程等.这些偏微分方程具有能量守恒特性,在数值计算中,尽量构造其保持能量守恒特性的数值算法在数值模拟微分方程的行为中具有重要的意义.哈密尔顿系统的辛几何算法在1984年首次由冯康院士及其研究小组提出,在辛几何算法的基础之上Bridges和Reich等人在1997年构造出了多辛算法.以上两种算法具有长时间精确计算的优点,但不足之处在于这两种算法只能近似保持方程的能量守恒.近年来,许多学者构造了保持这些偏微分方程的保能量算法.在1999年保持哈密尔顿系统能量守恒的平均向量场方法被Quispel和McLachlan等人提出了,王雨顺利用平均向量场方法构造了多辛结构偏微分方程的保能量方法.本文利用平均向量场方法和拟谱方法构造耦合和复偏微分方程的高阶保能量格式和多辛整体保能量格式,利用这些新格式对这些偏微分方程进行数值模拟,并对数值结果进行分析.在第一章,时间上利用四阶平均向量场方法,空间上利用傅里叶拟谱方法对复修正KdV方程进行离散,构造了复修正KdV方程的高阶保能量格式,利用构造的高阶保能量格式数值模拟孤立波的演化行为.数值结果表明复修正KdV方程的高阶保能量格式可以很好地模拟孤立波的演化行为,并且可以精确地保持了方程的离散能量.在第二章,对于求解耦合偏微分方程,我们在时间上利用四阶平均向量场方法,空间上利用傅里叶拟谱方法对KGS方程和CSBE方程构造高阶保能量格式,并模拟孤立波的演化行为.数值结果表明KGS方程和CSBE方程的高阶保能量格式可以很好地达到预期效果.在第三章,利用二阶平均向量场方法,拟谱方法对具有多辛结构的一维偏微分方程:BBM方程和复修正KdV方程构造多辛整体保能量格式,利用构造的多辛整体保能量格式数值模拟孤立波的演化行为,并证明了新格式能保方程离散的整体能量守恒特性.数值结果表明BBM方程和复修正KdV的多辛整体保能量格式可以很好地模拟孤立波的演化行为,并且可以精确地保持BBM方程和复修正KdV的离散整体能量守恒特性.在整体保能量守恒特性方面,BBM方程和复修正KdV方程新构造的格式比已有的经典的多辛格式更加精确,计算时间也比高阶平均向量场方法大大缩短.在第四章,对具有多辛结构的二维偏微分方程:ZK方程构造多辛整体保能量格式,利用构造的多辛整体保能量格式数值模拟孤立波的演化行为.数值结果表明整体保能量方法不仅可以长时间的模拟孤立波的演化行为,而且也可以精确地保持ZK方程的离散能量守恒.