积分微分方程的区域分裂并行算法

来源 :山东大学 | 被引量 : 0次 | 上传用户:zk0529
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
区域分裂方法是并行求解大型偏微分方程的有效方法,因为这种方法可以把大型计算问题分解成小型问题,从而简化了计算,上个世纪50年代,在并行机出现之前,区域分裂方法已经在串行机上得到了应用。现在,随着并行计算机和并行算法的发展,自上世纪80年代开始,区域分解算法开始蓬勃发展起来,现在高性能并行计算机已经广泛地应用于能源部门如核工业和实验、石油工业,生物,基因以及气象(天气预报及模拟)等等。 区域分裂方法通常用于以下两种情况:第一,可以通过区域分解的方法把大型问题转化为小型问题,实现问题的并行求解,缩短求解时间;第二,许多问题在不同的区域表现为不同的数学模型,那么可以在不同的区域对数学模型采取不同的方法进行求解,从而自然的引入区域分解方法,实现了并行计算。因为区域分解方法可以把大型问题分解成小型问题,复杂边值问题分解为简单边值问题,串行问题分解成并行问题,因此这种方法的研究十分活跃,具体方法也多种多样。 区域分裂算法是把计算区域分解成若干子域,于是原问题的求解转化为在域上求解,它的优越性表现在:第一,它把大问题化为若干小问题缩小计算规模;第二,允许使用局部拟一致网格,无需用整体拟一致网格,各子域可以用不同离散方法进行计算。这对于形态极不规则问题如锅炉燃烧问题有很大的灵活性;第三,允许不同子域选用不同的数学模型,以便整体模型适合于工程物理实际情况。如油、气藏模拟、气体绕飞体流动等;第四,算法是高度并行的,在各个子域内独立进行。 当用区域分裂方法来数值求解数学模型问题时,我们首先会根据问题的特性或问题求解区域的几何特点对区域进行划分,把整个求解区域划分为若干个子区域,然后在每个子区域上分别求解独立的子问题,实现并行计算。当求解抛物型偏微分方程时,一般情况下我们需要知道方程的初边值条件,但由于区域时认为划分的,那么对于子域而言,至少有一测度非零的边界条件是未知的,即相邻子区域相交内边界上的边界条件是未知的,我们的工作就是给出子区域的边界格式条件,或者虽然不能在求解子问题时明确给出子区域的内边界条件,但是在算法的求解过程中给出子区域相交内边界条件。 显隐格式区域分裂方法就是以显式的格式给出相邻子区域间相交内边界的边界条件的一种方法.我们知道,如果采用显示方法来求解问题,那么方法可以自然实现并行.但是由于显格式式条件稳定的,对于迭代步长有限制条件,增加了迭代步骤,从而使得计算机求解时间增长,因此一般在求解抛物方程时,普遍会采用隐式方法,隐式方法绝对稳定,对时间没有限制条件,但是在每一步计算时,都要求解一个大型方程组,当剖分加细时,方程组随之变大,求解这样一个方程组耗时较长,从而影响了计算时间。显隐格式区域分裂方法综合了二者的优点,借助前一层数值解的信息,给出在这一层的子问题的未知边界条件,把一个整体区域上的问题化为若干个子区域上的子问题,在每个子区域上用隐式方法求解,从而实现了并行。从计算角度而言,就是把一个整体的大型方程组分解成若干个小型方程组实现了并行。由于给出子区域间内边界条件的方法利用了前层数值解的信息,具有显性性质,导致了算法需要一个稳定性条件,但这个稳定性条件没有显式方法那么严格。 关于显隐格式区域分解方法,前人已经做了很多工作。Dawson,Q.Du和Dupont在文献[3]中采用显隐格式区域分解有限元差分方法,由显式向前差分格式给出子区域间内边界条件,得到了最优阶的ι∞模误差估计.基于这种方法,Q.Du等人在文献[9]采用一种高效显隐格式区域分解方法,用多层显格式给出内边界条件,最终得到最优的ι∞模误差估计。Dawson和Dupont在文献[1]中采用了显隐格式区域分裂方法,定义了一个函数,用这个函数给出了区域间内边界条件,得到了L2模误差估计。在此基础上,Dawson和Dupont在文献[10]中采用了基于中心显隐格式有限差分区域分裂方法,得到了ι2。模误差估计。本文共分为两部分,第一章为预备知识,第二章研究积分微分方程区域分裂方法,内边界上的函数值由上一层函数值得到,实现算法的并行,误差估计上采用了新的方法,定义一个新的包含边晃信息的算子,从而删去了原来误差估计中的日H1/2,得到了最优估计。
其他文献
集值优化理论在不动点、变分学、微分包含、最优控制、数理经济学等领域有着广泛的应用,是目前应用数学领域中备受关注的热点之一,而集值优化问题在各种解意义下的最优性条件是
种群生态学是研究种群数量动态与环境相互作用关系的科学,它起源于人口统计学,应用昆虫学和水产学.Lotka-Volterra(1925,1926)的模型理论是理论生态学的一个里程碑,生态学并由此进
证券收益率的统计规律或分布形式是金融市场的基本性质之一。大量实际的高频金融数据表明,收益率的分布远远偏离正态分布,具有尖峰、厚尾特征。在研究过程中,人们逐步发现稳定分
本文首先考虑奇异半线性椭圆问题{-△u=u-γ+g(x,u), x∈Ω,u>0, x∈Ω, (1)u=0, x∈()ΩQ,其中,Ωc RN(N≥3)是具有光滑边界()Ω的有界区域,γ>0是一个正常数,g:Ω×R→R是一个Carathéodo
变分不等式理论己有较突出的地位,其最重要也很有趣的内容是设计有效的数值计算法来寻求近似解.鉴于此,本文从以下几个方面讨论: 1.简述变分不等式理论的历史背景和研究现状。