收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

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

张守慧  
【摘要】: 抛物型方程是偏微分方程中基本方程之一.在自然科学的众多领域中,许多现象是用抛物型方程或者方程组描述的[1],例如热传导以及其它扩散现象、化学反应、粒子的运输等等.另外在一些问题的数值处理中也经常出现抛物型偏微分方程.在现代数值方法中,最早为人们所注意且理论分析完善的是有限差分法,因此抛物型偏微分方程的有限差分方法一直是人们关心的焦点.随着并行机的问世和发展,传统的有限差分方法在不同方面暴露出各自的弱点.例如,古典显式虽然适合于并行计算,但它是条件稳定的,特别是多维问题中计算步长受到严格的限制;古典隐式和Crank-Nicolson格式是绝对稳定的,但需要求解联立方程组,不便于直接在并行机上应用.因此需要构造具有良好稳定性、并行性和计算精度的新的差分方法. 七十年代初,Miranker指出用有限差分逼近偏微分方程时,主要是组织传统差分方法的并行实现,至于设计新算法,推动力是很小的;之后的十几年高阶差分格式方面的研究得到了发展[3-7].八十年代初,上述情况由于Evans和Abdullah的工作[8-12]而发生了变化,他们设计的分组显式方法保证了数值计算的稳定性,同时由于显式求解而使该方法具有很好的并行性质.它是不同类型Saul'yev非对称格式[13]的恰当组合.由于不同的Saul'yev格式的截断误差中某些项绝对值相等,符号相反,在同一时间层和不同时间层上连续交替使用不同的非对称格式,可带来截断误差的部分抵消,从而提高方法的计算精度.这些非对称格式都是隐格式,但由于它们之间的巧妙结合,可以显式求解,这就是Evans-Abdullah分组显式(GE).这项工作说明了建立满足上述要求的新的差分格式是可能的.但是在将分组显式思想应用于变系数问题时,稳定性的证明遇到了困难. 在此基础上,张宝琳等在[14-16]中提出利用Saul'yev非对称格式构造分段隐式的思想,并恰当的使用交替技术建立了多种显-隐式和纯隐式交替并行方法,取得了稳定性和并行兼顾的研究成果.之后又将方法推广到变系数问题,并用能量法证明了方法的绝对稳定性.在数值试验中发现,分段或分块并行计算的结果一般都比原来相应的未加分裂时的结果精确.所以通过分而治之的策略来建立新算法,不但可以用于并行,还可以提高精度.之后涌现出大量的并行差分算法的研究成果,韩臻在[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 ofComputer Mathematics》. 第二章的创新之处在于对于变系数的抛物问题给出和Crank-Nicolson格式相匹配的交替分组显格式,并用能量方法证明了稳定性. 接下来的三章内容中,主要借鉴了[14,23-26,46-48]中交替分组(段)格式的思想,将它与高阶差分格式[27-32]相结合,建立了高阶的交替格式;经证明方法都是绝对稳定的,且具有并行性质;在时间步长足够小时,空间的局部截断误差可达到O(h~4) 第三章中首先给出高阶的显、隐差分格式,在隐格式的基础上构造了四种非对称格式,通过它们之间的巧妙组合建立了交替分段显隐格式;由[33-34]中的Kellogg引理证明了方法的无条件稳定性;得到了方法的局部截断误差;数值算例证实了方法的实用性,并且可以达到O(h~4)的误差精度.本章内容已被《计算物理》接受. 第三章的创新之处在于将高阶差分格式与交替分段显隐格式的思想相结合,构造出高阶的交替分段显隐格式.方法具有良好的数值稳定性,空间误差阶可以达到O(h~4)阶. 第四、五章是在第三章的基础上,引入高阶Crank-Nieolson差分格式并适当变形,构造了八个非对称逼近方程;通过交替使用这些差分格式建立了两种不同的数值计算方法. 第四章是单独应用八个非对称差分格式构造了交替分组显格式,通过Kellogg引理证明方法的稳定性,通过误差分析得到两层抵消部分误差后的误差可以达到O(Τh)阶;数值试验说明了方法的实用性,并且时间步长充分小的前提下误差对于空间来说可以达到四阶.本章内容发表在《山东大学学报》(理学版). 第五章利用交替分段格式思想,引入了在非对称格式之间插入对称格式的思想,从而建立交替分段Crank-Nicolson格式.本格式同样具有数值稳定性、可并行性质,误差分析时由于插入中心对称格式,所以在这些点处两层之间部分抵消后误差较小,这在数值试验中得到了证实.交替分段Crank-Nicolson格式的结果比相应的交替分组显格式结果要好.本章内容已投稿到《应用数学与力学》. 后两章的创新之处在于将高阶差分格式与交替分组、交替分段Crank-Nicolson思想充分的结合,建立与第三章不同的非对称Saul'yev格式,在此基础上构造相应的交替方法.这些方法都具有绝对稳定性、可并行性质;并且在空间上截断误差可以达到O(h~4)阶;最后给出数值算例说明方法的适用性.


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 刘轶中;李物兰;方春华;曾诚;;双曲型方程的一类高阶算法及其数值实验[J];湖南师范大学自然科学学报;2011年03期
2 胡劲松;王玉兰;;广义非线性Sine-Gordon方程的一个隐式差分格式及其Richardson外推[J];西华大学学报(自然科学版);2011年04期
3 安静;罗振东;;三维抛物方程基于POD基的差分格式及后验误差估计[J];数学物理学报;2011年03期
4 张天德;左进明;段伶计;;广义improved KdV方程的守恒差分格式[J];山东大学学报(理学版);2011年08期
5 胡劲松;王玉兰;;广义正则长波方程的拟紧致守恒差分格式[J];东北师大学报(自然科学版);2011年03期
6 段国东;蒋建;牛成虎;;半差分格式在支付交易费用的欧式看涨期权定价模型中的应用[J];广东石油化工学院学报;2011年04期
7 郭瑞;阿布都热西提·阿布都外力;;基于不同分裂的修正局部Crank-Nicolson方法对一维热传导方程的应用[J];兵团教育学院学报;2011年02期
8 李军锋;李国明;李如山;;地震波传播数值模拟中的频散校正补偿[J];价值工程;2011年19期
9 李英杰;张顶立;刘保国;;考虑变形模量劣化的应变软化模型在FLAC~(3D)中的开发与验证[J];岩土力学;2011年S2期
10 刘睿;曾现洋;;抛物型方程的人为解问题[J];聊城大学学报(自然科学版);2011年01期
11 童慧峰;唐志平;;靶结构对烧蚀模式激光推进效应影响的数值模拟[J];高压物理学报;2011年03期
12 苏克勤;程俊芳;;结合矢通量分裂的NND格式及其应用[J];河南大学学报(自然科学版);2011年05期
13 刘相国;郝江锋;;一类双曲型方程的有限差分法[J];宜宾学院学报;2011年06期
14 王巧云;;考虑二次梯度影响的非线性渗流模型的预估校正法[J];数学理论与应用;2011年02期
15 孙佳慧;;有限体积法和有限元方法之间的比较[J];长春师范学院学报;2011年08期
16 ;[J];;年期
17 ;[J];;年期
18 ;[J];;年期
19 ;[J];;年期
20 ;[J];;年期
中国重要会议论文全文数据库 前10条
1 高智;;双曲和对流扩散方程的时空高精度摄动差分格式[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
2 张新军;梁昌洪;金谋平;;MEI结合有限差分法二维应用的几个技术[A];1997年全国微波会议论文集(下册)[C];1997年
3 王双虎;;Euler-Lagrange方程的Jet保结构算法[A];中国工程物理研究院科技年报(1998)[C];1998年
4 陈光南;李德元;万正苏;王一博;;求解三维扩散方程的两种差分格式[A];中国工程物理研究院科技年报(2002)[C];2002年
5 秦孟兆;;多辛几何差分格式[A];Structure Preserving Algorithm and Its Applications--Proceedings of CCAST (World Laboratory) Workshop[C];1999年
6 冉政;胡文沛;;差分格式的群不变特性与激波捕捉问题的研究[A];计算流体力学研究进展——第十一届全国计算流体力学会议论文集[C];2002年
7 高树椿;韩隆恒;赫新;;关于《差分格式的群不变特性与激波捕捉问题的研究》一文的错误[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
8 李晓勇;;基础沉降的差分格式近似解分析[A];2004年度上海市土力学与岩土工程学术年会论文集[C];2004年
9 曲凯;邢耀国;张旭东;李兵尚;;基于NND差分格式某固体火箭发动机启动过程一维内流场计算[A];中国宇航学会2005年固体火箭推进第22届年会论文集(发动机分册)[C];2005年
10 刘玉侠;;适应于复杂流场的二阶Godunov型差分格式[A];中国工程物理研究院科技年报(2005)[C];2005年
中国博士学位论文全文数据库 前10条
1 张守慧;抛物型方程的几种可并行的有限差分方法[D];山东大学;2009年
2 邱建贤;流体动力学方程的差分格式及其收敛性研究[D];南京航空航天大学;2001年
3 王姗姗;若干非线性Schr(o|¨)dinger方程及其相关问题的数值研究[D];南京航空航天大学;2011年
4 刘伟;抛物问题的时空局部网格加密方法[D];山东大学;2006年
5 杭旭登;偏微分方程迭代并行解法与网格优化方法[D];中国工程物理研究院;2004年
6 顾培英;基于应变模态技术的结构损伤诊断直接指标法研究[D];河海大学;2006年
7 李志涛;多孔及变形介质中混溶及非混溶驱动问题的数值方法[D];山东大学;2007年
8 崔艳芬;线性传输方程和KdV方程满足两个守恒律的差分格式[D];上海大学;2008年
9 刘全;Z-Pinch磁流体力学方程组的数值模拟方法[D];中国工程物理研究院北京研究生部;2002年
10 黎芳;几何分析方法在图像处理中的应用[D];华东师范大学;2007年
中国硕士学位论文全文数据库 前10条
1 申霞;基于POM的近海三维水质模型研究及其应用[D];河海大学;2006年
2 王爽;求解KdV的方程的两种差分格式[D];新疆大学;2005年
3 盛志明;对流扩散方程的并行差分格式[D];中国工程物理研究院;2010年
4 朱蓓佳;亚式期权及其编程计算[D];华东师范大学;2008年
5 马菊香;泥沙扩散问题数值解法的研究[D];新疆大学;2006年
6 费荔枝;二维四阶抛物型方程差分格式[D];郑州大学;2009年
7 郭瑞;修正局部Crank-Nicolson方法对KDV方程的应用[D];新疆大学;2010年
8 陈丽文;城市交通噪声分析及其系统研究[D];内蒙古工业大学;2006年
9 张波;基于流体动力学的交通流模型及数值模拟[D];吉林大学;2007年
10 刘晓蕾;局部相平衡模型(LPEM)方程数值解法的实现与分析[D];中国地质大学(北京);2007年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978