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

来源 :山东大学 | 被引量 : 0次 | 上传用户:zshihao
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
诸多物理现象都可用泊松方程和扩散方程来描述,泊松方程是数学中一个常见于静电学,机械工程和理论物理的偏微分方程;扩散方程常见于化学扩散,热传导,医学,生化方面和一定的生物反应过程.由许多偏微分方程的标准离散化得到的用于求解这两类方程的若干快速计算方法,受到了人们的密切关注.目前,大规模的科学工程计算需要高速度,大容量的并行机和有效的数值并行方法和并行算法.考虑到成本与速度,并行机的使用起到了非常重要的作用.同时,区域分解方法是一个很有力的工具,可将整个问题化为各个子域问题,然后并行求解,由于其构造简单,速度快,得到了广泛应用.在现代数值方法中,有限差分方法是最早且很完美的求解方法,所以对于研究抛物型和椭圆型问题的有限差分方法,备受人们的关注和重视.有限差分方法对于椭圆型问题的逼近往往需要求解较大的稀疏矩阵问题[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阶,以及也适用于二维变系数扩散问题.
其他文献
本文通过分析高职院校科研成果应用推广的现状,剖析科研成果应用推广存在的问题与原因。提出加强科研成果推广应用应转变意识、营造氛围,成立专门机构,深化校内制度改革等关措施
如果说儒为我们指向了不无高远的“天下理想”的话,那么,禅则以极其直切的“活在当下”为自身的旨归。然而,为中国文化所共有的“根身性”却决定了二者看似分歧,实际上却难掩
本文对砂轮修整技术领域的技术发展进行了梳理,描述了砂轮修整技术由手工到设备、由通用到专用的技术演进路线,对我国国内砂轮修整技术的发展进行了思考和展望。
采用每秒100幅(LOCAM型)高速扫描摄影方法,现场拍摄了北京亚运会男子三级跳远银牌得主邹四新和1991年全国田径冠军赛金牌获得者赵本宏在同一次比赛中的三跳起跳技术影片。经
我国经济和科学技术的快速发展,促进了我国建筑行业的快速发展,与此同时,建筑工程施工效率、施工进度、施工质量都得到了有效的提高。各种各样的机电设备被广泛的应用在建筑
随着我国教育事业的快速发展,学生的认知水平与学习能力逐渐成为影响教学质量的关键因素。当代教学活动是面向全体学生的,重视学生的知识基础,提高教学活动的针对性,才能促进
【正】融合创新理论,是当下国际经济界普遍认同的经济创新理论。这种理论称将各种创新要素通过创造性地融合,使各创新要素之间相互匹配,从而使创新系统的整体功能发生质的飞
本文通过估算2011-2020年我国成品油需求,进而估计原油总需求及由此可能产生的能源缺口。研究结果表明,我国汽车和交通运输业对成品油的需求量较大、增长较快。如果基于能源
<正>2018年冲压、封头、精冲行业技术、经济数据调研工作于2019年8月结束。本次调研针对51家冲压企业、13家封头企业、13家精冲企业进行了技术、经济数据的采集和分析。调研
经济的发展离不开道路的支撑,而边坡是最为常见的一种道路安全事故类型,会给来往行人与车辆带来极大的安全隐患,所以必须要通过各种手段对边坡进行监测,并对监测结果进行分析