论文部分内容阅读
随着科学的快速发展,越来越多的物理、化学和生物过程可以用非线性发展方程或者电磁场方程来描述.在许多情况下这些系统是保守的,因此如何为这些系统设计高效且能保持系统守恒量的算法一直是计算科学的研究热点.本博士学位论文致力于几类非线性发展方程的局部保结构算法(多辛算法,局部能量、动量守恒算法)和三维Maxwell方程高效保结构算法的研究及其数值分析.主要研究成果包括：1)辛算法以及全局保能量、动量等传统保结构算法在处理偏微分方程时,除了要考察方程是否是保守系统外,还必须考虑边界条件是否适当,只有在适当的边界条件下才可以使用这些保结构算法.为了增加保结构算法的适用范围,我们以耦合非线性Schrodinger系统、Boussinesq系统和Klein-Gordon-Schrodinger系统为研究对象,为其构造了一系列的局部保结构算法,包括多辛算法,局部能量、动量守恒算法.这些局部保结构算法能在任何时、空区域上保持离散的局部守恒律.当边界条件适当时,这些局部保结构算法自然就是全局保结构算法,反之不然.除此之外,我们还分析了其中一些局部保结构算法的非线性稳定性和收敛性.数值实验表明局部保结构算法不但能得到较好的数值解,而且确实能保持系统的局部守恒律和全局守恒律.通过与文献中已有的数值算法比较,展现了本文算法的优点.2)三维Maxwell方程具有双哈密尔顿结构.应用谱方法离散哈密尔顿函数和哈密尔顿算子,再对所得到的有限维哈密尔顿系统用平均向量场方法求积,从而我们得到求解三维Maxwell方程的两个格式(下称“AVF(2)”和"AVF(4)"). AVF(2)和AVF(4)自动保持系统的两个哈密尔顿.我们证明了AVF(2)和AVF(4)保持离散的能量、动量和散度.数值色散分析表明它们是无条件稳定的且无耗散的.严格的误差分析表明AVF(2)和AVF(4)分别在时问方向具有二阶和四阶收敛性,在空间方向都具有谱精度.数值实验很好地证实了理论分析结果.3)AVF(2)和AVF(4)是通过直接离散Maxwell方程得到的,计算时编程实现比较复杂.为了设计更高效的保能量算法,我们利用指数算子分裂和组合的技巧分别得到逼近Maxwell方程的时间二阶和四阶分裂方法.这些分裂模型的每个子问题都是哈密尔顿系统且和原问题具有相同的哈密尔顿函数.对每个子问题,分别应用谱方法离散哈密尔顿函数及哈密尔顿算子,进而应用平均向量场方法求积所得的有限维哈密尔顿系统,从而得到时间二阶和四阶分裂格式(下称"S-AVF(2)"和"S-AVF(4)").我们证明了S-AVF(2)和S-AVF(4)能同时保持4个离散能量且是无条件稳定的.此外,利用离散的傅里叶变换,我们还可以将所得格式写成显式形式.借助于能量分析方法,我们得到了S-AVF(2)和S-AVF(4)的误差估计.数值实验证实了理论分析结果.