非线性扩散占优偏泛函微分方程时间离散化方法比较及若干相关研究

来源 :湘潭大学 | 被引量 : 0次 | 上传用户:zliang_1981
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
对于非线性扩散占优Volterra型偏泛函微分方程的时间离散化,以往最常用的方法是Volterra泛函微分方程(VFDE)的隐式Euler方法及梯形方法,有时候也使用隐式中点法及Bellen和Zennaro所建议的二级Lobatto IIIC Runge-Kutta方法。为了避免出现Bellen和Zennaro所指出的“order failure”和“stability failure”现象,本文约定恒使用VFDE的正则Runge-Kutta法(Canonical Runge-Kutta methods,abbr.CaRK),不使用连续Runge-Kutta法(abbr.CRK)。为方便计,我们将上述四种不同的时间离散化方法依次简记为CaIE、CaTr、CaIM及CaLo2。本文的第一项主要工作是通过理论分析和数值试验,对上述四种不同的时间离散化方法进行了仔细比较(参见第三章)。我们发现其中计算效率最高的是CaIM方法。这一发现为选择时间离散化方法提供了新的重要科学依据。其重要性在于对于求解上述类型的强非线性问题,有可能导致在不久的将来CaIM方法逐渐上升为使用最广的时间离散化方法之一,而CaTr方法有可能使用得越来越少,甚至有可能逐渐被淘汰。本文的第二项主要工作是证明了CaIM方法在一致时间网格上是二阶最佳B-收敛的(参见第二章)。这一新的阶结果比以往文献中已有的关于CaIM方法的阶结果提高了一阶。正是由于这一创新成果,才使人们意识到CaIM方法有可能比CaTr方法更好,因为前者是B-稳定且二阶最佳B-收敛的,而后者即使是对于求解常微分方程问题,也仅仅是A-稳定且二阶经典收敛的,而不是B-稳定和B-收敛的,因此对于求解上述类型的强非线性问题,后者有可能出现不稳定现象,甚至有可能导致整个计算失败。应当指出,CaIM方法的最佳B-收敛阶为2这一性质尽管要求时间网格均匀,但对于分段均匀的时间网格同样是成立的。而这正是实际应用中出现得最多的情形。此外,作为本项研究的延伸,我们构造了一类分数阶偏泛函微分方程的高精度算法(参见第四章)。
其他文献
铁电薄膜因其优越的性能越来越受到大家的关注。利用铁电薄膜制作的电子元器件,工作在辐射环境中时会受到辐照的影响而导致性能下降。当器件工作在航天、航空等领域时,由于在航空和航天环境中包含不同的高能射线粒子,这些射线会对工作在其中的电子系统、电子元器件等产生各种损伤效应,最终导致航天器出现功能失效、甚至会发生航天器坠毁,所以要求电子器件拥有较强的抗辐射性能。铁电抗辐射性能研究,目前主要集中在宏观性能随着
多铁材料同时具有铁弹、铁电、铁磁等多种序参量,且这些序参量通过耦合能够产生一些新的效应,比如磁电耦合效应,使得其在传感器、多态存储、自旋电子器件等领域具有广阔的应用前景。对于多铁纳米材料,它不仅能在纳米尺度上呈现出磁电耦合效应,更能促进器件实现多功能化、集成化和微型化,近年来受到广泛的关注和研究。对纳米尺度多铁材料的力电磁耦合行为进行研究,不但可以加深对多铁材料复杂的物理和力学现象的理解,而且能为
铁电材料是一种功能材料,其具有的优异电学和光学性能孕育出了它广阔的应用前景。电畴翻转是铁电材料显示宏观非线性本构行为的微观物理机制,铁电畴的取向直接决定着铁电材料的物理性质和应用方向,而电畴翻转后不同取向铁电畴体积分数的变化则直接与铁电材料及其器件的效率和稳定性挂钩。因此,确定铁电材料的畴取向及其体积分数(分布特性)对铁电器件工程至关重要。多铁性材料作为一种特殊的铁电材料,除具有铁电性外,还兼具有
软件定义网络(Software-Defined Networking,SDN)是一种新型网络体系结构,它实现了网络设备中控制逻辑与数据转发功能的分离,采用集中式的控制方式控制整个网络,并向上层应用提供开放的可编程接口。这种设计模式极大地简化了网络策略的部署,缩短了网络应用的开发周期,在数据中心和云计算中得到了广泛的应用。但是,这种网络体系结构也引入了新的安全威胁。本文将对SDN体系结构中的若干安全
热障涂层(Thermal Barrier Coatings,TBCs)因为耐高温、高隔热、抗腐蚀等优异性能,已成为航空发动机涡轮叶片等热端部件不可缺少的热防护材料。然而,在含有杂质颗粒的高温、高速燃气环境下服役时,会因为界面氧化、颗粒冲蚀、CMAS腐蚀(钙镁铝硅等金属氧化物的混合物,简称CMAS)等多种因素导致涂层剥落。其中,CMAS腐蚀指高温下钙镁铝硅混合氧化物熔融、渗透到涂层中,引起涂层结构、
本文利用谱配置方法求解一类带弱奇异核的和非线性的Volterra型积分微分方程,并且构造高精度算法,着重分析该方法的误差估计和收敛性,并进行数值实验,验证所给方法的有效性。第二章,用Legendre谱配置方法对非线性Volterra-Fredholm-Hammerstein积分方程进行求解。并且对于所给的数值格式和误差进行了分析,即当核函数充分光滑时,计算所得数值格式的L2范数和L∞范数误差呈指数
半拉格朗日(半拉氏)方法没有CFL稳定性条件限制且比欧拉方法更加稳定。在文献当中,水平集重新初始化方程的半拉氏方法是复杂的,这是没有必要的。因为重新初始化过程是辅助的,我们提出了耦合有改善紧邻交界面网格点处精度的投影技术的一阶半拉氏格式。标准的二阶半拉氏格式被用于推进水平集对流方程。实现是简单的,包括在块结构自适应网格上。半拉氏方法的有效性通过大量的数值算例被证实,其中包括在给定速度场中做被动对流
复杂网络已经吸引了来自科学技术不同领域研究者的大量关注。继复杂网络中的小世界特性与无标度性质之后,自相似性已经成为复杂网络的第三大基本特征并且在近年来得到了广泛的研究。分形分析可以有效地揭示一些分形对象的自相似性。然而,重分形分析是一个用来系统性地刻画理论和实验复杂分形对象空间异质性的更加强大而有效的工具。尽管已有一些研究复杂网络自相似性的分形和重分形分析算法,但是这些算法的效率并不高而且已有的算
本文主要研究了阶数为α∈ (0,1)的具有Caputo导数的时间分数阶Navier-Stokes方程,这类方程可以用来模拟分形介质中的反常扩散现象.我们主要讨论它的适定性.本文的工作主要分成五个部分:在第二章中我们首先通过把Helmholtz投影子作用到目标方程上消去压强项,将其转化为抽象形式的发展方程;然后根据分数阶发展方程理论给出其相应的积分方程并定义适度解;利用半群理论和不动理论得到方程的全
本论文主要研究了M-模糊化凸结构理论,(L,M)-模糊凸结构理论,及(L,M)-模糊凸结构(子)范畴理论.在M-模糊化凸结构理论中,首先,我们引入并刻画了M-模糊化(弱)JHC凸空间、M-模糊化几何区间空间、M-模糊化(弱)Peano、Pasch、沙漏区间空间,并讨论了它们之间的关系。其次,我们引入了M-模糊化基点序的概念,并用其刻画了 M-模糊化几何性质、M-模糊化Peano性质、M-模糊化Pa