Hamiltn ODEs的高效辛算法和Hamilton PDEs的多辛算法.pdfVIP

哈尔滨理工大学理学硕士论文 HamiltonODES的高效辛 算法和HamiltonPDEs的多辛算法 摘要 这篇硕士论文集中了作者在攻读硕士学位期间的主要研究成果.具体内容是 我们关于HamiltonianODEs的高效辛算法和HamiltonianPEDs的多辛积分子 的一些结果，以及它们在实际中的应用. 针对HamiltonianODE。的数值计算，我国已故科学家冯康和美国科学家 Ruth于80年代创立了辛几何算法.相关的理论分析和大量数值结果表明该方 法对HamiltonODES的计算是适合和有效的，尤其是对它们进行长期跟踪计算 时，优势更为明显.我们分析了辛算法在实际中的应用.一般说来，生成函数法可 能会出现高阶导数项，这使得它在实际应用中受到很大的限制.因此现在实际应 用中，人们很大程度上依赖于RI{(PRK)方法，生成函数法只是提供理论上的工 具.我们知道，对于不可分Hamilton系统而言，辛RK(PRK)法必然为隐式的， 而隐式RK方法在实际应用中，运算量都是非常惊人的，尤其是在高阶的水平 上，更是令人难以实现.所以辛RK(PRK)方法的在实际的应用，一直以来都局限 在低阶的方法上，高阶方法大都以理论结果给出.但是如果一个隐式RK(PRK) 方法的系数矩阵仅有实特征值，一个基于Butcher和Bickart的技巧可以大大的 减少其在实际应用中的运算量.基于上述原因，我们考虑了这种具有实特征值的 高效辛RK法和辛PRK法.我们首先证明了一个s级的上述方法的阶pGs+1, 尤其是一个s级的辛RK法，当s是偶数时，它不可能达到、+1阶，而当s是 奇数时，则没有此阶障碍的限制.同时，我们也指出一个s级的辛PRK法它的 阶完全可以达到s+1。接着，我们又证明了，在高阶的水平上，无论是这种辛 RK法还是辛PRK法，它们的运算效率都不能和相应阶的辛组合方法相媲美。 我们最近的研究工作又表明，组合方法又会带来诸如稳定性方面的问题.最后， 我们推荐了一些低阶的高效格式。 建立在Bridges意义上的多辛Hamilton系统的多辛算法被证明深刻反映 了系统的本质.在 Reich研究工作的基础上，我们把PRK方法应用到这种多 辛HamiltonianPDEs上，给出了它是多辛的一些充分条件，同时，对这种多辛 PRK法关于线性能量和动量的保持也作了相应的探讨，这些都在一定程度上发 展和拓广了现有的结果. 一I一 哈尔滨理工大学理学硕士论文 关键词 HamiltonianODES,HamiltonianPDEs,Runge-Kutta法 (RK)块 Runge-Kutta法 P〔RK)，辛积分子，多辛积分子 一11一 哈尔滨理工大学理学硕士论文 EfficientSymplecticMethodsforHamiltonianODES andMulti-SymplecticMethodsforHamiltonianPDEs Abstract Thisdissertationisacollectionoftheauthorsresearchworkduringhisstudy pursuingamasterdegree.ThecontentsareourresultsontheefficientSymplectic methodsforHamiltonianODESandmulti-SymplecticmethodsforHamiltonian PDEs,andtheirapplicationsinpractice. ForHamiltonianODES,thesymplecticmethod,foundedbythelateChinese scientistFengKangandAmericanscientistRuthinthe1980s,hasbeenproved bymanytheoriticalresultsandnumerousnumericalexperimentst