论文部分内容阅读
一切耗散效应可以忽略不计的物理过程都可表示成能够保持辛几何结构不变的哈密尔顿系统的形式,它在自然界中具有普适性,也就是说大多数孤子方程都可以表示成哈密尔顿形式.现代数值计算的基本原则是尽可能保持原问题的本质特征.因此,研究保持哈密尔顿系统的辛几何结构特征的数值方法是必然的.本论文主要讨论无穷维Hamilton系统的多辛几何算法.与辛几何算法的主要差别是,多辛算法不仅在一定的边界条件下保持系统的离散空间上的辛形式之和,而且能够保持局部的辛形式,从而多辛几何算法更多的是体现在系统局部的守恒性质,更能体现系统的本质特征.
本文研究的无穷维Hamilton系统中Klein-Gordon-Zakharov(简称KGZ)方程组是一个重要的模型,它是由一个Klein-Gordon方程和一个Zakharov方程耦合而成.我们通过对KGZ方程组作正则变换后,得到了它的一个多辛方程组及其几个相关守恒律.然后用Gauss-Legendre Runge-Kutta方法对此多辛方程组离散,得到了KGZ方程组的多辛格式,证明了该格式具有离散形式的多辛守恒律.对中点格式,通过消去中间变量得到了与多辛格式等价的多辛Preissman格式.我们通过大量数值实验验证了所构造的多辛格式的有效性和长时间的数值稳定性,同时多辛格式还能很好地模拟原孤立波的波形,说明我们的理论分析是正确的.
另外,本文还通过对空间和时间方向分别用Fourier拟谱方法和中点方法离散KGZ方程组的多辛方程组,得到了非线性KGZ方程组的多辛Fourier拟谱格式.我们通过数值实验证明了该格式的有效性.数值结果表明多辛Fourier拟谱格式是正确可行的.