一些非线性发展方程(组)的辛和多辛算法

来源 :中国科学技术大学 | 被引量 : 18次 | 上传用户:C_k_b
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文主要研究了一些非线性演化方程和方程组的辛和多辛算法.用辛和多辛算法研究了对称正则长波(简称SRLW)方程和Klein-Gordon-SchrSdinger (简称KGS)方程组的孤立波随时间的演化情况,及其有关的守恒律. 孤立波及其理论是现代非线性科学研究的重要组成部分,它所对应的数学模型往往是线性或者非线性的偏微分方程或者方程组,它遍布现代科学研究的各个领域和角落,在流体力学,原子分子物理,等离子体物理,光纤通讯,化学化工,生物医药等诸多科学领域有广泛应用.近几十年来,孤立波方程定解问题的精确求解及其理论研究一直是科学研究的一个热门课题,取得了丰硕的研究成果.求解孤立波方程精确解的方法除了传统的反散射方法、Hirota双线性方法、Backlüind变换方法外,近年来又涌现了很多新方法,如齐次平衡法、双曲正切法,级数展开法等.虽然已经有了这么多求解孤立波方程精确解的方法,但是由于事物的复杂性,大多数定解问题还是不能精确求解的,只能通过数值方法近似求解,特别是对于非线性的情形. 一切耗散效应可以忽略不计的物理过程都可表示成能够保持辛几何结构不变的哈密尔顿系统的形式,它在自然界中具有普适性,也就是说大多数孤子方程都可以表示成哈密尔顿形式.现代数值计算的基本原则是尽可能保持原问题的本质特征.因此,研究保持哈密尔顿系统的辛几何结构特征的数值方法是必然的.为了适应这一需要,我国计算数学的奠基人冯康院士首次在1984年系统地提出了能够保持哈密尔顿系统辛结构不变的辛几何算法。随后,辛算法成为国内外计算科学讨论的一个热门课题,在这一领域涌现了一大批的研究成果.自从冯康院士提出辛算法以后,这一算法有了两次大的飞跃和发展,一次是从有限维向无限维的推广,一次是向多辛算法的深入发展.二十世纪九十年代后期,Marsden等从变分原理的角度提出了多辛积分的概念,而Bridges,Reichs从辛几何的角度提出了多辛算法,近年来这一算法得到了迅猛的发展,它已经成功地用来解决了很多现实问题,模拟了各种物理现象. 本文第一章对辛和多辛算法的相关背景作了简单介绍,对辛空间的基本知识作了简要的回顾. 第二章简单总结了构造哈密尔顿系统辛算法的常用方法,主要有生成函数法,Runge-Kutta方法,块Runge-Kutta,可分哈密尔顿系统的显式方法和构造高阶精度格式的复合方法等.本章最后还介绍了把无限维哈密尔顿系统降低为有限维系统的Fourier拟谱方法. 第三章构造了KGS方程组的一族辛格式,分析了此族格式的守恒律和收敛速度.证明了我们构造的辛格式保持电荷守恒,而且分析了它的能量误差,证明了数值解的整体截断误差和解的收敛速度均为0(T<2>+h<2m>).数值实验表明我们所构造的辛格式具有长时间的数值模拟能力,我们的理论分析是正确的. 第四章简单介绍了多辛哈密尔顿系统及其相关的守恒律.以KGS方程组为例介绍了一些构造多辛算法的常用方法,主要有Fortrier拟谱方法,GaUSS-Legendre Runge-Kutta方法等.对KGS方程组研究了它的多辛格式的守恒律.通过分析我们发现对KGS方程组而言,Preissman格式在加权意义下电荷守恒,而多辛Fourier拟谱方法具有经典意义下的电荷守恒律,而能量表达式由于是三次多项式,Preissman格式和多辛Fourier拟谱格式均不具有能量守恒的特征,尽管如此,用多辛格式模拟所产生的误差相对较小.我们还精确分析了能量和动量的残量.数值例子表明我们所构造的格式能够模拟各种孤立波,得到了很多有趣的物理现象.同时也说明了我们的理论分析的正确性. 第五章研究了SRLW方程的多辛格式,对它构造了多辛Preissman格式和Fourier拟谱格式.数值例子说明了格式的有效性和优越性.这一章我们以SRLW方程为例简单讨论如何把多辛格式修正成局部能量守恒格式或者局部动量守恒格式. 第六章就本文的主要内容作了简单的总结.就辛和多辛算法未来的发展方向做了简单的展望,并罗列了一些公开问题及其面临的一些挑战。
其他文献
偏微分方程在工程实际和科学技术中有广泛的应用背景,研究其数值解对处理在电磁学、声学等领域中的很多物理问题具有很重要的意义。传统求解偏微分方程的数值方法都基于网格
Internet为我们提供了大量的信息资源,日渐成为我们生活中不可缺少的一部分。个性化服务是未来信息服务的发展趋势,获取用户个性化的兴趣及信息需求是开展个性化服务的基础,
基于偏微分方程(PDE)图像分析和处理方法已经越来越成为图像处理学中的研究热点,PDE方法是一种比较精细的图像分析和处理方法,可以用于图像降噪、增强和分割等领域。图像膨胀
本文主要用变分法研究了(1)和(2)两类带Hardy-Sobolev临界指数的椭圆方程组解存在性:({-div(|Du|p-2Du/|x|ap)-μ|u|p-2u/|x|(a+1)p=λ1|u|p-2u/|x|(a+1)p-c+|u|p*-2u/|x|bp*
数据包络分析(data envelopment analysis)简称DEA,是数学、运筹学、数理经济学和管理科学的一个新的交叉领域。由A.Charnes和W. W.Cooper等人于1978年开始创建,并被命名为DEA
1997年,Zheng首次提出了新的密码原语——签密。签密能够在一个合理的逻辑步骤内同时完成数字签名和公钥加密两项功能,而其计算量和通信成本都要低于传统的“先签名后加密”
基于主元分析的过程监控技术是目前过程自动化和控制领域的研究热点问题之一。然而,传统主元分析直接应用于工业过程时,由于其本身性能的限制,会出现很多问题。国内外学者针
学位
刚性延迟微分方程广泛地出现于自动控制、电力工程、生态学等领域。目前,国内外文献大都集中对求解这类问题的Runge-Kutta方法、Rosenbrock方法、线性多步法及更一般的线性方