抛物型偏微分方程几类新并行算法及其收敛性分析.pdfVIP

原创性声明 本人郑重声明：所呈交的学位论文，是本人在导师的指导下，独立进 行研究所取得的成果。除文中已经注明引用的内容外，本论文不包含任何 其他个人或集体已经发表或撰写过的科研成果。对本文的研究作出重要贡 献的个人和集体，均已在文中以明确方式标明。本声明的法律责任由本人 承担。 论文作者签名： 尴 日 期： 型垒：』：!。 关于学位论文使用授权的声明 本人完全了解山东大学有关保留、使用学位论文的规定，同意学校保 留或向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅 和借阅：本人授权山东大学可以将本学位论文的全部或部分内容编入有关 数据库进行检索，可以采用影印、缩印或其他复制手段保存论文和汇编本 学位论文。 (保密论文在解密后应遵守此规定) 论文作者签名：谭导师签名：弛期： 山东大学硕士学位论文 抛物型偏微分方程的几类新的 并行算法及其收敛性分析 杨 建 华 (山东大学数学与系统科学学院，济南250100) 摘 要 当今实践中提出的问题，计算规模越来越大，精度要求越来越高，而单机计算 的速度已接近极限，因此不可避免地促进了并行计算和区域分解算法的发展． Schwarz算法可以把复杂区域分解为若干相互覆盖的子区域，在子区域上可以 用快速算法求解。所谓加性Schwarz算法的发展，又可克服交替方法的串行性，更 利于并行处理．本文我们给出半线性对流一扩散方程的一种加性Schwarz算法， 这种算法是基于区域分解的．由误差估计，我们知收敛率与每个时间步长上的迭 代次数无关． 加性Schwarz方法是解椭圆问题及其它静态问题的一种有效的迭代方法，在 这篇文章中我们将用此方法解决下面的抛物型对流，扩散问题： 象一耋嘉掣磐圳川) ∽眯Q×J “=0 (z．t)∈aQ×J “(z，0)=“o(z) z∈Q ，I●●●●●●●●●●，、●●●●●●●●●●L 首先，我们从纯代数的角度给出我们的第一类并行算法．由于此算法是基于区 域分解和子区域校正的，所以我们需要对区域Q进行剖分，在每个小空间W(n。) 上定义相关的投影算子： 砰：伊_÷W．许进超[5]给出了一种子域校正形式的 并行算法，可简单表述为： u“e”=札d4+∑耋。龟我们将利用许进超的并行方法 考虑我们的对流一扩散问题，从而构造我们的第一类并行算法． 其次，在前面剖分的基础上，引进粗空间VH，并且定义投影算子珊：V“_÷VH 山东大学硕士学位论文 在每个小区域上解边值问题，得霞：(?=1．2．…．。，)再解得粗空间VH上的有限 元解舀o．则我们可得一平均值面=了玎1∑a。．以此为基础我们可以构造第二类并行 12U 算法． 我们还需要给出一些条件： 条件(A)对任意的z∈而存在一个开集合职和i∈{1，2，…，Ⅳ)使得X∈D。 且D。nn∈Q： 条件(B)对任意的z∈豆至多存在‰个包含。的子区域． 在文章的第三章和第四章详细介绍了两种算法的构造，并分别证明了解的收敛性． 针对两算法得到了不同的收敛性定理，即： 解．当At=DfH2)，h=o(H)时，有 c-， “n—w“％≤c(^≈+1+△t+(会；+鲁)号)． 其中C为依赖于u但与迭代次数№空间网格剖分参数h和H与时问步长△t 无关的常数． 法的解，则有 其中C为依赖于u但与迭代次数m，空间网格剖分参数^．H以及时间步长△t 无关的常