微分方程保能量算法的研究

来源 :南京大学 | 被引量 : 0次 | 上传用户:hongyu203311
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
守恒和耗散系统广泛来源于天体力学,分子动力学,电路模拟,量子力学和电磁学等科学领域。能量是刻画守恒和耗散系统的最重要的物理量之一。从保结构算法的角度来看,好的数值方法必须尽可能地保持来源于原连续系统的离散的物理和几何结构。保能量或能量耗散方法是一类特殊的保结构算法。他们可以分别保持守恒或耗散系统的首次积分或李雅普诺夫函数。大量的理论和实验结果表明数值方法的能量守恒或能量耗散性质可以保证线性误差增长,对原系统定性性质的正确模拟和强有力的数值稳定性等优良的数值特性。针对守恒或耗散的常微分和偏微分方程,本文致力于构造和分析具有优良几何性质和极高代数精度的保能量或能量耗散方法。第一章简要地介绍了保能量方法及其相关方法的背景知识和我已被接收或发表的关于保能量算法的原创性工作。本文剩余内容可分为三部分:第二章考虑了具有线性主部的守恒或耗散的常微分系统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问题,含一个或两个空间变量的(耦合)非线性薛定谔方程等一系列数值实验,来证明文中新方法的优秀的数值特性。
其他文献
本文将现存文献中有关具有泊松到达过程交换网中的最优尺寸标度问题发展到具有一般连续时间更新到达过程的交换网中。具体地:我们考虑一具有N个队列的交换排队网络,数据包在网
随着全球化的加强,国家之间的交流合作逐渐增多。尤其在面对全球性的问题时,不同国家和地区间的学者会进行互访调研,以寻求最优解决方案。在调研的过程中,采访是一种最直接、效率最高的调研方法。在欧盟环境署提出的“以自然为本的创新”项目中,作者担任采访口译。在释意学派翻译程序理论的指导下,作者结合在采访中的口译录音转录完成了本篇实践报告。释意学派口译理论中提出了“理解——脱离语言外壳——表达”的翻译程序。该
随着Higgs玻色子的发现,人们已经发现了粒子物理标准模型预言的所有基本粒子,这是物理学史上的一次伟大壮举。但是目前对于强子结构的认识仍然存在很多问题。在传统的夸克模
目的通过本临床研究,观察天地解郁汤治疗阴虚痰瘀阻络型缺血性脑卒中后抑郁的临床疗效,为天地解郁汤更好的应用于临床提供理论依据。方法将阴虚痰瘀阻络型缺血性脑卒中后抑郁患者62例,随机分为观察组和对照组,每组31例。两组在常规治疗基础上,均予以舍曲林治疗,观察组联合天地解郁汤治疗。观察周期为2个疗程,共12周。记录治疗前后HAMD、NIHSS、ADL和中医证候积分变化。结果1.两组临床疗效总有效率比较,
在运筹学中排序作为一个非常重要的部分,有十分重要的研究意义。传统的经典排序中加工时间一般是一定的常数。但是在我们的生活中,任务的实际加工时间经常会受基本加工时间、
本文研究可压缩气体在一个缓慢扩张球中的解的整体存在性和长时间行为。气体由Boltzmann方程描述。假设扩张的球是Ωt={x∈R3:|x|
本文主要研究了离散Lomax和离散Rayleigh分布。Lomax分布被视为指数和Gamma分布的混合分布。由于Lomax分布包含了单调递减和单调递增的失效率,使其在研究医学、生物和工程等
磁性多层膜所具有的巨磁电阻效应自问世以来对磁电子学、凝聚态物理学和材料学等学科的发展进步做出了巨大的贡献。磁性多层膜广泛地应用于存储、家电、工业器械等方面,极大
本文主要整理了C*-代数与有限维逼近[1]中的一些基本定义和定理,重点整理了Stinespring定理,Arveson扩张定理以及Voiculescu定理,主要参考了Brown和Ozawa的书[1].但对其中一
本文主要介绍了1.5keV电子穿越玻璃直管/锥管的时间演化的二维角分布,主要包括2个部分1对低能带负电荷的电子束穿透绝缘高绷硅玻璃管的过程进行动力学研究。观察到了低能电子