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

二维泊松方程和扩散方程的一类显式并行算法

许秋燕  
【摘要】: 诸多物理现象都可用泊松方程和扩散方程来描述,泊松方程是数学中一个常见于静电学,机械工程和理论物理的偏微分方程;扩散方程常见于化学扩散,热传导,医学,生化方面和一定的生物反应过程.由许多偏微分方程的标准离散化得到的用于求解这两类方程的若干快速计算方法,受到了人们的密切关注.目前,大规模的科学工程计算需要高速度,大容量的并行机和有效的数值并行方法和并行算法.考虑到成本与速度,并行机的使用起到了非常重要的作用.同时,区域分解方法是一个很有力的工具,可将整个问题化为各个子域问题,然后并行求解,由于其构造简单,速度快,得到了广泛应用. 在现代数值方法中,有限差分方法是最早且很完美的求解方法,所以对于研究抛物型和椭圆型问题的有限差分方法,备受人们的关注和重视. 有限差分方法对于椭圆型问题的逼近往往需要求解较大的稀疏矩阵问题[1],随着现代计算机的发展,迭代方法如Jacobi, Gauss-Seidel, SOR方法[2,3]常被用在此大型计算中,也是解此大型方程组的一类很好的方法.众所周知,泊松方程的并行单元是利用Jacobi并行迭代算法[4]来解决的,因其有明显的并行性。[5]中冯慧等通过不同点的隐式差分格式之间的相互约化来建立新型迭代(Stencil)方法,此方法和Jacobi方法同样具有并行性,却比Jacobi收敛快.超松弛(SOR)迭代方法是很有效的方法,却没有本质的并行性.目前,仅有较少的方法可以实现并行迭代,如着色法[6]和点并行SOR方法[7]等,但很难拓展. 而对于求解时间依赖的扩散问题的数值方法,主要有两种类型,显式和隐式差分方法.显式差分方法很容易在并行机上实现,有好的并行性,但由于稳定性限制[8],对时间步长的限制相当苛刻。隐式差分方法是绝对稳定的,但在每一时间步都需要去求解较大的线性或非线性代数方程组,计算效率不高。因此需要构造具有良好稳定性,并行性和计算精度的新的差分方法。八十年代初,Evans和Abdullah[9-11]设计了分组显式方法保证了数值计算的稳定性,同时由于显式求解而使该方法具有很好的并行性质,它是不同类型Saul'yev非对称格式[12]的恰当组合.在此基础上,张宝琳等在[13-15]中提出利用Saul'yev非对称格式构造分段隐式格式的思想,并恰当的使用交替技术建立了多种显-隐式和纯隐式交替并行方法,取得了稳定性和并行兼顾的研究结果.另外,Peaceman与Rachford在[16]中提出的交替方向隐式方法(ADI)具有一些相当好的性质,但是对于多指令流多数据流(MIMD)并行机来说,由于每个并行处理器在奇数时间层计算一列或若干列的数值,而在偶数时间层计算一行或若干行的数值,或者是相反,因而在每个时间层大部分数据需要从一个处理器传送到另一个处理器,这是很不经济的. 在导师王文洽教授的悉心指导下,本文作者基于并行计算,区域分解和交替分组显式思想,对二维泊松方程和扩散方程构造了一类有效的显式并行算法,对方法进行了理论分析并给出了数值算例,均证实了我们在本文中所提的算法是有效的,实用的.全文共分四章,下面几段将简要介绍一下各章的主要内容. 第一章,基于区域分解思想,对二维泊松方程提出了一种新的有限差分并行迭代(FDPI)算法.首先将求解区域划分为四个子区域,并从经典的五点差分格式构造出四个新的迭代格式.然后,在迭代次数为奇数或偶数时,分别给出新算法的实现过程,且证明了其收敛性.最后给出了数值算例以验证此迭代算法的有效性和精确性,而且也适用于二维变系数椭圆型问题.此章内容已发表在《Numerical Methods for Partial Differential Equations》. 本章内容的创新之处在于虽然迭代格式是半隐式的,但可结合边界条件和分组显式思想,显式地,并行地计算。特别地,格式中添加了松弛因子ω并理论分析得到最佳松弛因子ωopt,大大的提高了收敛速率,减少了迭代次数.在与Jacobi, Gauss-Seidel迭代算法和数学Stencil[5]方法的数值结果比较中发现,我们所提出的FDPI算法具有更少的迭代次数,更短的计算时间和更快的收敛速率. 第二章,给出了求解二维泊松方程的另一种新的并行超松弛(P-SOR)迭代算法.利用四个超松弛(sOR)迭代格式完成算法,实现过程类似于第一章,最后数值试验说明此算法的优越性.本章内容已投到《International Journal of Computer Mathemat-ics》. 此章内容的创新之处在于,众所周知的SOR算法并没有明显的并行性,而我们给出的新算法(P-SOR)具有本质并行性,并且通过与Jacobi, Gauss-Seidel,数学Stencil[5]及第一章中的FDPI算法关于时间的消耗,速度及空间复杂性的比较,表明本节提出的算法更为有效。 第三章,基于区域分解和分组显式思想,给出了求解二维扩散方程的一种新的有限差分并行算法.首先构造了逼近扩散方程的四个新的差分格式,然后将求解域划分为四个子域,在奇数或偶数时间层上算法的实现过程类似于第一章,并证明了格式的无条件稳定性,且误差阶为O((?)+h2).最后数值试验说明算法的适用性,并且时间步长充分小的前提下误差关于空间达到2阶.本章内容已投到《Applied Mathematics and Computation》. 第四章,提出了求解二维扩散方程一种新的交替分组(N-AGE)算法.采用显式方法[10],使用第三章中构造的四个差分格式在奇数或偶数时间层上交替使用,来完成新的显式算法.最后数值试验通过与第三章中的并行算法和分组显式方法[10]在同时间和不同空间参数下误差的比较,表明了此算法的有效性和精确性.本章内容已投到《计算数学》. 第三,四章的创新之处在于,构造的新的差分格式虽是隐式的,但利用分组显式思想可显式求解.并且,差分格式的交替使用导致截断误差中的主要项符号相反,从而抵消.数值试验也验证了此类算法的有效性,且关于空间的收敛速率达到2阶,以及也适用于二维变系数扩散问题.


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 陈泮勤;非定常两维平流扩散方程的数值实验[J];气象学报;1985年01期
2 张忠诚;;Zn-MnO_2干电池扩散方程的求解[J];山东大学学报(工学版);1989年02期
3 陈军;李秋姝;;一类非线性扩散方程解的估计及非负解的存在性[J];贵州大学学报(自然科学版);1990年01期
4 张庆红;郁国瑞;;扩散方程的具有振动边值的Dirichlet问题[J];河北师范大学学报(自然科学版);2006年03期
5 李静然;高志明;马逸尘;;广义拉格朗日乘子的虚拟区域法在泊松问题中的应用[J];西安交通大学学报;2006年08期
6 陈昊;钟尔杰;;扩散方程的高精度恒稳定差分格式[J];电子科技大学学报;2007年S1期
7 梁志辉;李之杰;;柱体三维电磁场的泊松方程的数值计算[J];赤峰学院学报(自然科学版);2009年03期
8 杨传富;;扩散方程系数的半逆问题[J];数学年刊A辑(中文版);2011年01期
9 杨志峰;周同明;;含源扩散方程的一种高精度差分方法[J];北京师范大学学报(自然科学版);1992年01期
10 范学刚;一类扩散方程逆问题的稳定性[J];焦作师范高等专科学校学报;2003年01期
11 邓滨,施保昌,王广超;求解含源项扩散方程的两种格子BGK模型比较[J];华中科技大学学报(自然科学版);2004年12期
12 张新东;胡月宏;;同伦分析方法求具有边界条件的扩散方程的精确解[J];山东大学学报(理学版);2008年12期
13 王盘贞;应用细胞神经网络求解泊松方程的新途径[J];南京邮电学院学报;1995年01期
14 孙亮;马东军;秦丰华;孙德军;;二维泊松方程的两步预估校正格式[J];力学与实践;2010年01期
15 刘建新;;三维泊松方程的一种解析解法[J];华北电力大学学报;1987年02期
16 马昌前;扩散方程的岩浆动力学意义[J];地质科技情报;1988年02期
17 苏秀娟;用K理论研究南京大厂地区污染物扩散[J];南京气象学院学报;1989年03期
18 田振夫;泊松方程的高精度三次样条差分方法[J];西北师范大学学报(自然科学版);1996年02期
19 张海亮,于鸣歧;一类非线性反应扩散方程解的Blow-up问题[J];数学杂志;1997年04期
20 梅家斌;分形在电镀中的应用[J];中南民族大学学报(自然科学版);2002年01期
中国重要会议论文全文数据库 前10条
1 苏娟;沙宪政;;复杂毛细管分布组织中氧及葡萄糖传输数学模型及计算机模拟的研究[A];中国生物医学工程学会成立30周年纪念大会暨2010中国生物医学工程学会学术大会报告论文[C];2010年
2 蒋丽丽;孙建国;;基于泊松方程的正交曲网格生成技术[A];中国地球物理学会第二十四届年会论文集[C];2008年
3 王青平;白武明;王洪亮;;多重网格在二维泊松方程有限元分析中的应用[A];中国科学院地质与地球物理研究所第11届(2011年度)学术年会论文集(上)[C];2012年
4 桑苏玲;郑璐石;任可亮;;用差分线法求解泊松方程[A];“力学2000”学术大会论文集[C];2000年
5 吴振翔;缪柏其;;MCMC方法在金融中的一点应用[A];2003中国现场统计研究会第十一届学术年会论文集(下)[C];2003年
6 张庆红;谭亦丽;韩伟;;基于扩散方程的上解与解的研究[A];数学·力学·物理学·高新技术研究进展——2006(11)卷——中国数学力学物理学高新技术交叉研究会第11届学术研讨会论文集[C];2006年
7 金蔚青;;材料制备的界面质量输运现象[A];全国第三届溶胶—凝胶科学技术学术会议论文摘要集[C];2004年
8 杨水龙;;一个化学反应扩散方程奇异行波摄动解[A];数学·力学·物理学·高新技术研究进展——2002(9)卷——中国数学力学物理学高新技术交叉研究会第9届学术研讨会论文集[C];2002年
9 袁光伟;盛志强;;大变形网格上扩散方程计算格式的设计[A];全国计算物理学会第六届年会和学术交流会论文摘要集[C];2007年
10 袁光伟;;扩散方程计算方法研究进展[A];中国工程物理研究院科技年报(2009年版)[C];2010年
中国博士学位论文全文数据库 前10条
1 许秋燕;二维泊松方程和扩散方程的一类显式并行算法[D];山东大学;2010年
2 张勇;薛定谔—泊松方程组的数值计算和分析及其应用[D];清华大学;2012年
3 张瀛;一维分数阶扩散方程的若干反问题研究[D];复旦大学;2012年
4 苏丽娟;分数阶扩散方程的几种数值解法[D];山东大学;2010年
5 李灿;几类分数阶反常扩散方程的数值分析[D];兰州大学;2012年
6 靳朝辉;离子交换动力学的研究[D];天津大学;2004年
7 赵强;扩散方程有限体积格式及守恒型并行差分格式研究[D];中国工程物理研究院;2011年
8 李颖花;具周期位势或周期源的高阶扩散方程[D];吉林大学;2008年
9 金春花;具非线性源的非散度型扩散方程[D];吉林大学;2008年
10 吕龙进;分数阶奇异扩散方程的几种解法及其应用[D];复旦大学;2012年
中国硕士学位论文全文数据库 前10条
1 徐旭光;基于有限差分法的泊松方程算法研究与软件实现[D];电子科技大学;2011年
2 张方方;柱坐标系下时间分数阶波扩散方程及分析[D];山东大学;2011年
3 陈寅;基于泊松方程的浮雕处理技术研究[D];国防科学技术大学;2010年
4 徐倩倩;扩散方程高阶格式的分组迭代法[D];山东大学;2011年
5 史蓉;扩散方程的势对称和Degasperis-Procesi方程的精确解[D];西北大学;2010年
6 郭宾;二维分数阶反应子扩散方程的有限差方法[D];湘潭大学;2010年
7 张庆红;扩散方程的Dirichlet问题的推广[D];河北工业大学;2003年
8 贾云涛;对流扩散方程并行算法的若干研究[D];浙江大学;2006年
9 肖玲玲;非线性分数阶反应子扩散方程数值解法[D];湘潭大学;2010年
10 谷文娟;时间分数阶对流—扩散方程反问题研究[D];山东理工大学;2011年
中国重要报纸全文数据库 前1条
1 顾恺;GPRS的计费[N];网络世界;2001年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978