积分微分方程的区域分裂并行算法-计算数学专业论文.docxVIP

原创性声明本人郑重声明：所呈交的学位论文，是本人在导师指导下，独立 原创性声明 本人郑重声明：所呈交的学位论文，是本人在导师指导下，独立 进行研究所取得的成果。除文中已经注明引用的内容外，本论文不 包含任何其他个人或集体已经发表或撰写过的科研成果。对本论文 的研究作出重要贡献的个人和集体，均已在文中以明确方式标明。 本声明的法律责任由本人承担。 论文作者签名： 懿 日 期： 砂多、亏，≯z，蒯、叮，≯z， 关于学位论文使用授权的声明 本人完全了解山东大学有关保留、使用学位论文的规定，同意学 校保留或向国家有关部门或机构送交论文的复印件和电子版，允许 论文被查阅和借阅；本人授权山东大学可以将本学位论文全部或部 分内容编入有关数据库进行检索，可以采用影印、缩印或其他复制 手段保存论文和汇编本学位论文。 (保密的论文在解密后应遵守此规定) 论文作者签名：趟导师签名： 积分微分方程的区域分裂并行算法王庆超 积分微分方程的区域分裂并行算法 王庆超 (山东大学数学与系统科学学院，济南，山东，250loo) 中文摘要 区域分裂方法是并行求解大型偏微分方程的有效方法，因为这 种方法可以把大型计算问题分解成小型问题，从而简化了计笄，上 个世纪50年代，在并行机出现之前，区域分裂方法已经在串行机上 得到了应用．现在，随着并行计算机和并行算法的发展，自上世纪80 年代开始，区域分解算法开始蓬勃发展起来，现在高性能并行计算 机已经广泛地应用于能源部门如核工业和实验、石油工业，生物， 基因以及气象(天气预报及模拟)等等。 区域分裂方法通常用于以下两种情况：第一，可以通过区域分解 的方法把大型问题转化为小型问题，实现问题的并行求解，缩短求 解时间；第二，许多问题在不同的区域表现为不同的数学模型，那 么可以在不同的区域对数学模型采取不同的方法进行求解，从而自 然的引入区域分解方法，实现了并行计算。因为区域分解方法可以 把大型问题分解成小型问题，复杂边值问题分解为简单边值问题， 串行问题分解成并行问题，因此这种方法的研究十分活跃，具体方 法也多种多样。 区域分裂算法是把计算区域分解成若干子域，于是原问题的求 解转化为在域上求解，它的优越性表现在：第一，它把大问题化为 若干小问题缩小计算规模；第二，允许使用局部拟一致网格，无需 用整体拟一致网格，各子域可以用不同离散方法进行计算。这对于 形态极不规则问题如锅炉燃烧问题有很大的灵活性；第三，允许不 同子域选用不同的数学模型，以便整体模型适合于工程物理实际情 山东大学硕士学位论文况。如油、气藏模拟、气体绕飞体流动等；第四，算法是高度并行 山东大学硕士学位论文 况。如油、气藏模拟、气体绕飞体流动等；第四，算法是高度并行 的，在各个子域内独立进行。 当用区域分裂方法来数值求解数学模型问题时，我们首先会根 据闻题的特性或问题求解区域的几何特点对区域进行划分，把整个 求解区域划分为若干个子区域，然后在每个子区域上分别求解独立 的子问题，实现并行计算。当求解抛物型偏微分方程时，一般情况 下我们需要知道方程的初边值条件，但由于区域时认为划分的，那 么对于子域而言，至少有一测度非零的边界条件是未知的，即相邻 子区域相交内边界上的边界条件是未知的，我们的工作就是给出子 区域的边界格式条件，或者虽然不能在求解子问题时明确给出子区 域的内边界条件，但是在箅法的求解过程中给出子区域相交内边界 条件。 显隐格式区域分裂方法就是以显式的格式给出相邻子区域间相 交内边界的边界条件的一种方法．我们知道，如果采用显示方法来求 解问题，那么方法可以自然实现并行．但是由于显格式式条件稳定 的，对于迭代步长有限制条件，增加了迭代步骤，从而使得计算机求 解时间增长，因此一般在求解抛物方程时，普遍会采用隐式方法， 隐式方法绝对稳定，对时间没有限制条件，但是在每一步计算时， 都要求解一个大型方程组，当剖分加细时，方程组随之变大，求解 这样一个方程组耗时较长，从而影响了计算时间。显隐格式区域分 裂方法综合了二者的优点，借助前一层数值解的信息，给出在这一 层的子问题的未知边界条件，把一个整体区域上的问题化为若干个 子区域上的子问题，在每个子区域上用隐式方法求解，从而实现了 并行。从计算角度而言，就是把一个整体的大型方程组分解成若干 个小型方程组实现了并行。由于给出子区域闻内边界条件的方法利 用了前层数值解的信息，具有显性性质，导致了算法需要一个稳定 性条件，但这个稳定性条件没有显式方法那么严格。 关于显隐格式区域分解方法，前人已经做了很多工作。Dawson，Q．Du 2 山东大学硕士学位论文和Dupont在文献【31中采用显隐格式区域分解有限元差分方法，由显 山东大学硕士学位论文 和Dupont在文献【31中采用显隐格式区域分解有限元差分方法，由显 式向前差分格