抛物型方程的几种可并行的有限差分方法

来源 :山东大学 | 被引量 : 0次 | 上传用户:chenchen0513
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
抛物型方程是偏微分方程中基本方程之一.在自然科学的众多领域中,许多现象是用抛物型方程或者方程组描述的,例如热传导以及其它扩散现象、化学反应、粒子的运输等等.另外在一些问题的数值处理中也经常出现抛物型偏微分方程.在现代数值方法中,最早为人们所注意且理论分析完善的是有限差分法,因此抛物型偏微分方程的有限差分方法一直是人们关心的焦点.随着并行机的问世和发展,传统的有限差分方法在不同方面暴露出各自的弱点.例如,古典显式虽然适合于并行计算,但它是条件稳定的,特别是多维问题中计算步长受到严格的限制;古典隐式和Crank-Nicolson格式是绝对稳定的,但需要求解联立方程组,不便于直接在并行机上应用.因此需要构造具有良好稳定性、并行性和计算精度的新的差分方法. 七十年代初,Miranker指出用有限差分逼近偏微分方程时,主要是组织传统差分方法的并行实现,至于设计新算法,推动力是很小的;之后的十几年高阶差分格式方面的研究得到了发展.八十年代初,上述情况由于Evans和Abdullah的工作]而发生了变化,他们设计的分组显式方法保证了数值计算的稳定性,同时由于显式求解而使该方法具有很好的并行性质.它是不同类型Saulyev非对称格式的恰当组合.由于不同的Saulyev格式的截断误差中某些项绝对值相等,符号相反,在同一时间层和不同时间层上连续交替使用不同的非对称格式,可带来截断误差的部分抵消,从而提高方法的计算精度.这些非对称格式都是隐格式,但由于它们之间的巧妙结合,可以显式求解,这就是Evans-Abdullah分组显式(GE)这项工作说明了建立满足上述要求的新的差分格式是可能的。但是在将分组显式思想应用于变系数问题时,稳定性的证明遇到了困难. 在此基础上,张宝琳等在[14-16]中提出利用saulyev非对称格式构造分段隐式的思想,并恰当的使用交替技术建立了多种显--隐式和纯隐式交替并行方法,取得了稳定性和并行兼顾的研究成果.之后又将方法推广到变系数问题,并用能量法证明了方法的绝对稳定性.在数值试验中发现,分段或分块并行计算的结果一般都比原来相应的未加分裂时的结果精确.所以通过分而治之的策略来建立新算法,不但可以用于并行,还可以提高精度.之后涌现出大量的并行差分算法的研究成果,韩臻在[17,18]中详细研究了一类纯显-隐分段和分块交替方法;[19-21]中冯慧等通过不同点的隐式差分格式之间的相互约化来建立新型迭代方法,此方法和Jacobi方法同样具有并行性,却比Jacobi收敛快.[22]中张志跃等给出变系数抛物型问题的分组显式方法,并用能量方法给出稳定性证明.[23-26]中王文洽等针对不同的方程建立了分段的显隐格式,证明方法的稳定性并给出数值算例.上述的方法在并行性和稳定性方面都有其优良的表现,但是他们都存在一个共同的问题,那就是它们都是基于二阶差分格式建立的,这直接影响数值计算中空间的误差精度.近年来,研究人员开始致力于研究高阶紧致差分格式.[27]中Sanjiva K.Lele提出了高阶的紧致差分格式,文中对格式的误差做了Fourier分析,并将它与经典差分格式做了比较.[28-32]中Mark H.Carpenter等针对不同问题提出高阶紧致格式并给出理论分析.将这些高阶差分格式与交替分组思想结合起来,是否可以得到稳定性好,可并行且精度高的数值算法呢?近年来,涌现出大量的高阶交替分组格式的研究工作[44-52]. 本文作者在王文洽教授的精心指导下,就抛物型问题的几类数学模型利用有限差分方法的技巧,构造了具有良好数值性质和计算效果的迭代方法、分组显式方法和交替分段方法,对方法做了理论分析并给出算例说明方法的适用性.本人拓广了前人的工作,不具有重复性.本文共分为五章. 第一章中主要利用[19]中冯慧提出的数值Stencil的概念,将其应用于二维对流扩散方程,建立了比Jacobi迭代收敛快的新型迭代算法.本章首先给出针对对流扩散方程的数值Stencil的概念,经过三次消元过程得到最终的数值Stencil,在此基础上建立了新型迭代算法;通过分析迭代误差证明了方法的收敛性,并与Jacobi迭代比较收敛阶;最后数值试验说明方法的适用性,证实了理论分析的结论。本章内容已被《International Journal of Computer Mathematics》接受. 第一章的创新之处在于将[19]中的方法应用到含有时间项的高维抛物型问题中,建立了收敛速度快、具有并行性质的新型迭代格式,通过分析迭代误差证明了方法的收敛性以及与古典迭代法之间收敛阶的比较;最后用实际例子说明了算法的有效性. 第二章主要运用[22,24]中的构造思想,将中心差分格式与分组显式思想相结合,针对含有变系数的对流扩散方程建立了分组显式方法,并用能量方法证明了该格式的稳定性.本章首先给出基于Crank-Nicolson差分格式的四种非对称的逼近方程,通过它们的巧妙组合建立交替分组显式方法;由于扩散项为变系数,所以采用能量法证明稳定性;最后数值试验说明方法的适用性.本章内容已投到《International Journal ofCornputer Mathematics》. 第二章的创新之处在于对于变系数的抛物问题给出和Crank-Nicolson格式相匹配的交替分组显格式,并用能量方法证明了稳定性. 接下来的三章内容中,主要借鉴了[14,23-26,46-48]中交替分组(段)格式的思想,将它与高阶差分格式[27-32]相结合,建立了高阶的交替格式;经证明方法都是绝对稳定的,且具有并行性质;在时间步长足够小时,空间的局部截断误差可达到O(h4) 第三章中首先给出高阶的显、隐差分格式,在隐格式的基础上构造了四种非对称格式,通过它们之间的巧妙组合建立了交替分段显隐格式;由[33-34]中的Kellogg引理证明了方法的无条件稳定性;得到了方法的局部截断误差;数值算例证实了方法的实用性,并且可以达到O(h4)的误差精度.本章内容已被《计算物理》接受. 第三章的创新之处在于将高阶差分格式与交替分段显隐格式的思想相结合,构造出高阶的交替分段显隐格式.方法具有良好的数值稳定性,空间误差阶可以达到O(h4)阶. 第四、五章是在第三章的基础上,引入高阶Crank-Nicolson差分格式并适当变形,构造了八个非对称逼近方程;通过交替使用这些差分格式建立了两种不同的数值计算方法. 第四章是单独应用八个非对称差分格式构造了交替分组显格式,通过Kellogg引理证明方法的稳定性,通过误差分析得到两层抵消部分误差后的误差可以达到O(Th)阶;数值试验说明了方法的实用性,并且时间步长充分小的前提下误差对于空间来说可以达到四阶.本章内容发表在《山东大学学报》(理学版) 第五章利用交替分段格式思想,引入了在非对称格式之间插入对称格式的思想,从而建立交替分段Crank-Nicolson格式.本格式同样具有数值稳定性、可并行性质,误差分析时由于插入中心对称格式,所以在这些点处两层之间部分抵消后误差较小,这在数值试验中得到了证实.交替分段Crank-Nicolson格式的结果比相应的交替分组显格式结果要好.本章内容已投稿到《应用数学与力学》. 后两章的创新之处在于将高阶差分格式与交替分组、交替分段Crank-Nicolson思想充分的结合,建立与第三章不同的非对称Saulyev格式,在此基础上构造相应的交替方法.这些方法都具有绝对稳定性、可并行性质;并且在空间上截断误差可以达到O(h4)阶;最后给出数值算例说明方法的适用性.
其他文献
武山铜矿1996年四季度三井延伸,建成投产的记破系统工艺流程见图1。溜破系统于12月投入使用起就暴露出许多问题,严重制约着生产的发展。1存在的问题武山矿石性质较为复杂,含硫高
通常的度量空间都是概率度量空间的一种特殊情况,所以研究概率度量空间中的非线性问题具有非常重要的意义。本文主要研究概率度量空间中非线性算子的理论与应用。 第一章,介
聚类分析作为多元统计分析的一个分支,用于图像分割已有相当一段历史。本文基于Gap Statistic(GS)方法从估计最佳聚类数的角度对图像存在的最佳分割进行研究。 2000年,斯坦
竖井提升中容器的过卷问题一直是困扰生产影响安全的事故隐患。虽然有多种后备保护,但过卷事故仍时有发生,轻者影响生产,重者造成人员伤亡,设备损坏。为解决这一难题.1992年版《煤
近年来,有监督字典学习由于其较好的分类效果获得了研究者的广泛关注。其主旨思想是通过对有标签训练样本的学习提取出具有类间判别性的特征,为分类提供更好的判断依据。但是对
面对国际形势越来越动荡,建立我国自主控制的卫星导航系统,开展相关的组合导航定位研究对我国的防御战略有很重要的意义。本文研究的主要内容是用线性最小二乘法对北斗/罗兰C
在这篇论文中考虑了两类中立型随机时滞微分方程的稳定性,一类是带有双时标马尔科夫切换的中立型随机泛函微分方程的指数稳定性,一类是带有中立项和变量时滞的非线性随机系统的
本文探讨了拓扑动力系统中有关熵和拓扑压中的一些问题,定义了新的熵和拓扑压,证明了它们的变分原理.本文的核心内容也就是对两个变分原理的证明. 在第一章中,主要介绍熵和拓
TheMainpointsinstallingthesingle-ropewindinghoistandtheessentialelementseffectingtheinstallingaccuracyZHANGJiping单绳缠绕式提升机主要由主轴装置、减速器、电机、
本文主要研究了Wirtinger型不等式和Lyapunov型不等式,共分为四章.  第一章为绪论,介绍了不等式的重要性.  第二章,在文献 [4],[5] 的基础上,定义新的微分算子,利用Picone恒