论文部分内容阅读
守恒和耗散系统广泛来源于天体力学,分子动力学,电路模拟,量子力学和电磁学等科学领域。能量是刻画守恒和耗散系统的最重要的物理量之一。从保结构算法的角度来看,好的数值方法必须尽可能地保持来源于原连续系统的离散的物理和几何结构。保能量或能量耗散方法是一类特殊的保结构算法。他们可以分别保持守恒或耗散系统的首次积分或李雅普诺夫函数。大量的理论和实验结果表明数值方法的能量守恒或能量耗散性质可以保证线性误差增长,对原系统定性性质的正确模拟和强有力的数值稳定性等优良的数值特性。针对守恒或耗散的常微分和偏微分方程,本文致力于构造和分析具有优良几何性质和极高代数精度的保能量或能量耗散方法。第一章简要地介绍了保能量方法及其相关方法的背景知识和我已被接收或发表的关于保能量算法的原创性工作。本文剩余内容可分为三部分:第二章考虑了具有线性主部的守恒或耗散的常微分系统y’=Q(My+▽U(y)), y(t0)=y0∈Rd,其中Q是一个反对称或半负定矩阵,M是一个对称矩阵.对于此系统,这一章提出并具体分析了一个二阶对称保首次积分或李雅普诺夫函数H(y)=1/2yT My+U(y)的指数型积分子。如果||QM||远大于||QHessU(y)||,该新方法比标准的保能量或能量耗散方法的容许步长更大,数值解的精度更高。针对一般的哈密顿系统y’=J-1▽H(y), y(t0)=y0∈Rd,这里J是一个常值辛矩阵,第三章构造了对称的函数适应型保能量方法。通过扩大函数适应空间,这一章在理论上证明了这类方法可以达到任意高阶。往函数适应空间中添加三角函数对{sin(ωw),cos(ωt)}可以自然地导出了任意高阶的三角适应型保能量方法。在处理以单个频率ω振荡的哈密顿常微分系统时,新提出的三角适应型保能量方法比标准的保能量方法效率更高。第四章考虑了一般的带一个时间变量t和两个空间变量x,y的哈密顿多辛偏微分方程Mzt+Kzx+Lzy=▽zS(z,x,y),z=z(x,y,t)∈Rd,这里M,K,和L都是反对称矩阵。分别用连续型Runge-Kutta方法和Gauss-Legendre/拟谱方法离散时间和空间方向,这一章提出了一个一般的框架来构造任意高阶的保某种离散局部能量律的方法。该离散守恒律是原方程局部能量守恒律(?)t(S-182z τKzx-1/2zτZzy)+(?)x(1/2zτKzt)+(?)y(182zτLzt)=的近似。除此之外,本文还展示了包括刻画三原子分子的不可分哈密顿系统,风诱导的振荡问题,α-Fermi-Pasta-Ulam系统,扰动的Kepler]问题,Duffing方程,高频振荡的Fermi-Pasta-Ulam问题,含一个或两个空间变量的(耦合)非线性薛定谔方程等一系列数值实验,来证明文中新方法的优秀的数值特性。