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

四阶抛物方程的一个有限差分并行格式

刘胜利  
【摘要】:在自然科学的许多领域中,很多现象是用抛物方程或方程组描述的,因此,用有限差分方法数值求解抛物偏微分方程问题具有重要的理论意义和应用价值。 随着向量机与并行机的问世与发展,并行数值算法已成为研制高性能并行计算程序的关键技术之一,到目前为止,有限差分并行算法主要是考虑二阶抛物方程问题。实际上,四阶抛物方程的并行化更有必要,相对二阶抛物方程而言,四阶抛物方程的显式格式的稳定性要求更高,达到r≤1/8,其中r=Δt/h~4;而其隐式格式的方程组是5对角的,这给大规模求解也带来了困难。 本文主要研究四阶抛物方程的并行化问题,受已有二阶抛物方程结果的启发,我们主要采用三层交替技术来实现四阶抛物方程的有限差分并行计算,得到一个可以将区域分块并行求解的差分格式,该格式绝对稳定而且局部截断误差达到O(△t+h~2)。 本文的主要工作分为以下几部分. 1.四种非对称格式 在第二章第一节中,我们考虑如下一维四阶抛物方程的初边值 计算可得,显式格式的稳定性条件为 利用中值定理,对方程(1)离散可得四种非对称格式如下, 对上面的四种非对称格式,我们可以得出如下的结论: 在第二章第二节中,我们主要利用非对称格式构造四阶抛物方 程的分组显式算法,所得算法进一步改善了误差阶, 可以计算得在四个点上的显式解,利用展开,可以计算得 显式解的局部截断误差都 定理3四种非对称格式组成的四点组的截断误差 由此可以看出,单独使用以上几种格式时,计算精度一般较差, 因为截断误差中包含了书项,即使网比r固定日寸,这一项随。趋于 零而趋于零,但是该项误差随:增大而增大,计算误差不好控制. 但是把这四种非对称格式组合起来形成如下节所给出的分组显式使 用会使计算精度大大改善,因而更具有实际应用价值. 已知。层的值,当求第层时,沿x增加或减少方向在顺序四 点上使用四点组格式(6)求解问题(1)的方法称为分组显式(Group Exphcit)方法,简记为GE方法. 现要求剖分节点是(N为自然数),此时有 要计贷.的内点.因此有川个4点组,剩余两点放在同一边,要么左 边,要么右边.剩余两点放在右边的情况,我们称为GER方法,剩 余两点放在左边的情况称为GEL方法.今以GER方法为代表进行 较仔细的讨论,GEI一方法的讨论类似,这里从略. 在第二章第三节中,我们已经看到,当单独使用GER或CEL方 法时,两种方法都是条件稳定的.这里我们把这种在 不同时间层交替使用GER和GEL方法,称为交替分组显式方法, 简称AGE方法.我们可以给出入GE方法的数学描述如下 定理5 AGE方法是绝对稳定的. AGE方法的一个特点是两个差分格式GER.和GEL的交替使用, 以连续交替计算为一个周期,不断重复进行,每一个计算周期内实 际上可以看作是使用了一个三层差分格式,以上稳定性证明就是基 于这种周期性情况做的. 在抛物方程或方程组的大型科学计算中,由于隐式格式一般具 有一好的稳定性,但由于隐式差分方程的并行求解比较困难,所以其 并行化研究可以从构造差分格式入手.受到构造分组显式格式方法 的启示,有人提出恰当地使用交替技术建立了多种显一隐式和纯隐 式交替并行方法,得到了稳定性和并行性兼顾的研究成果.我们对 厂一维问题(l)提出r有限并行解法. 泞先构造求解问题(l)的隐式段,对某个,‘1.考虑(,(、+,.,,十l)(,= 1.2.·…乙)诸点上的计算.在四个“端点”(,(。+1.l‘+l).(l()十2·11十 l).(,「,+乙一2.,;+l)和(,〔,+乙一1.,,+1)处分别用非对称格式(2).(5), (:l)和川.而在“内点”怀、+,.。十1)(,二1.2.·…L)处使用古典隐 式,即得到分段隐式. 对于问题(l)利用段隐式,一般的AsE一I方法是这样设计的:设 I一1=万LN、L为正整数,L全:.将同一奇数时间层要计算的点 分为\段,并自左而右依次按“古典显式一分段隐式一古典显式” 的规则作出安排,在下一时间层(偶数层)仍为N段,每段计算格式 交替地进行,即古典显式变为古典隐式,古典隐式变为古典显式, 非对称格式(2)变为非对称格式(3),非对称格式(5)变为非对称格式 (4),非对称格式(3)变为非对称格式(2),非对称格式(4)变为非对称 格式(匀.这样在偶数时间层上N段的计算规则就变为“分段隐式- 古典显式一分段隐式”,使得分段显一隐在不同的时间层之间交替进 行. 48 入S仁一I矩阵描述如一「: (I+,“.)U门 (l+z‘扩:)U“ 二(I一,f扩:)U 11. 二(了一,‘了.)U“- 11二0 .2.4二 明 Jr. 二 其中U’‘ 定理6由式了别所描述得交替分段显一隐方法了」乡E一l)是绝对 稳定的. 这样,我们构造了一种用于四阶抛物方程的有限差分并行格式 ASE一1.这种格式不仅绝对稳定,而且有很好的截断误差. 在L面所述抛物方程有限差分并行算法这个研究方向上,除了 我们已经提到的方法之外,目前还有许多工作正在研究之中,例如 多维变系数问题方法的稳定性和收敛性等,非线性问题在这方面的 研究更是一个十分广阔的领域.


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 赵才地;高维抛物方程Cauchy问题的构造性解法[J];温州师范学院学报;2003年05期
2 张壮志;;拟线性抛物型方程解的爆破[J];应用数学;1990年04期
3 潮兰萍;一类拟线性抛物方程的古典解[J];安庆师范学院学报(自然科学版);2002年03期
4 赵才地;抛物方程柯西问题的构造解法及机械化求解[J];华中科技大学学报(自然科学版);2004年03期
5 陈宏伟;半线性抛物方程及方程组的爆破问题[J];数学物理学报;1986年01期
6 袁洪君;一类非线性扩散方程广义源型解[J];吉林大学自然科学学报;1994年02期
7 张健,雷开泉;一类含非局部非线性项的抛物方程[J];四川师范大学学报(自然科学版);1994年04期
8 朱爱玲;;抛物方程的扩展混合体积元方法[J];山东师范大学学报(自然科学版);2006年04期
9 镇方雄;彭冬英;;一类具非线性对流项抛物方程解的渐近分析[J];佳木斯大学学报(自然科学版);2009年03期
10 何文平;周杰荣;;二阶非线性抛物方程解的定性性质[J];浙江万里学院学报;2010年04期
11 赵立纯;确定一维数理抛物方程的未知源问题[J];鞍山师范学院学报;1989年03期
12 刘振海;;关于抛物方程系数的最优控制问题[J];长沙电力学院学报(自然科学版);1990年01期
13 吴启光,李继春;奇异摄动半线性抛物方程的数值解法[J];应用数学和力学;1993年09期
14 孙志忠;解拟线性抛物方程的一类二阶差分格式[J];计算数学;1994年04期
15 姚正安,周铁强;一类拟线性抛物方程解的性质[J];数学杂志;1998年03期
16 王世祥,敬石心;一类具对流项非线性退化抛物方程古典解的存在性[J];长春大学学报;1999年01期
17 赵军生,辛洪学;具有圆盘形隔层尖灭的试井模型解析解[J];哈尔滨师范大学自然科学学报;1999年02期
18 郭金勇;;一个退化四阶抛物方程弱解的惟一性[J];广西科学;2007年02期
19 易青;邹新鹏;;R~N中有界域上的一类抛物方程非负径向解的爆破问题[J];南昌航空大学学报(自然科学版);2008年02期
20 杨金林;;二阶线性抛物方程在特征边界点C~∝解的存在性[J];包头钢铁学院学报;1991年01期
中国重要会议论文全文数据库 前10条
1 彭朝晖;;三维耦合简正波-抛物方程理论[A];中国声学学会2003年青年学术会议[CYCA'03]论文集[C];2003年
2 张林;笪良龙;卢晓亭;;基于抛物方程模型的南海中尺度涡声场分析[A];2009年全国水声学学术交流暨水声学分会换届改选会议论文集[C];2009年
3 李蔚清;吴慧中;;基于GPU的雷达盲区快速估算方法[A];中国图学新进展2007——第一届中国图学大会暨第十届华东六省一市工程图学学术年会论文集[C];2007年
4 宋俊;;抛物方程近场传播算子[A];2004年全国水声学学术会议论文集[C];2004年
5 付晓玉;;二阶偏微分方程能控能观性问题的统一处理[A];第25届中国控制会议论文集(中册)[C];2006年
6 林巨;孙军平;王欢;袁振林;;浅海内波环境下三维声传播规律仿真研究[A];2009年度全国物理声学会议论文集[C];2009年
7 陈任昭;张丹松;;带迁移因素和依赖年龄的种群生长的非线性扩散动力系统[A];1997年中国控制会议论文集[C];1997年
8 马啸;杨顶辉;;波动方程的加权近似解析离散化方法[A];中国地球物理学会第二十四届年会论文集[C];2008年
9 徐五岗;马力;;声场预报的实验研究[A];中国声学学会1999年青年学术会议[CYCA'99]论文集[C];1999年
10 彭朝晖;李风华;张仁和;;CMPE理论在水平变化浅海中的应用[A];中国声学学会2001年青年学术会议[CYCA'01]论文集[C];2001年
中国博士学位论文全文数据库 前10条
1 杜润梅;退化奇异抛物方程(组)的可控性理论[D];吉林大学;2013年
2 王昭;一类六阶抛物方程的若干问题[D];吉林大学;2013年
3 潘佳庆;一类拟线性抛物方程的定性研究[D];复旦大学;2004年
4 杨雄锋;双曲方程组与双曲—抛物方程组解的大时间状态估计[D];武汉大学;2004年
5 丁俊堂;若干非线性偏微分方程的爆破理论与最大值原理[D];山西大学;2006年
6 张志娟;具有间断系数扩散方程的局部间断Galerkin方法[D];中国工程物理研究院;2009年
7 雷玲;具备k-近似周期解的抛物系统的最优控制问题[D];浙江大学;2007年
8 吕桂霞;抛物方程有限差分并行算法理论[D];吉林大学;2004年
9 冯晓莉;椭圆和抛物方程的几类不适定问题[D];兰州大学;2010年
10 胡绘斌;预测复杂环境下电波传播特性的算法研究[D];国防科学技术大学;2006年
中国硕士学位论文全文数据库 前10条
1 贾慧美;一类抛物方程的逆源问题[D];东北师范大学;2011年
2 周莹;一类非自治的抛物方程的quench现象研究[D];湖南师范大学;2013年
3 林雪清;一类抛物方程(组)适定性的若干性质之研究[D];集美大学;2011年
4 郝芳;一类亚抛物方程的Cauchy问题[D];郑州大学;2011年
5 李振邦;六阶抛物方程行波解的不稳定性[D];吉林大学;2011年
6 张洪萍;一类高阶抛物方程的时间周期解[D];吉林大学;2013年
7 袁海峰;退化的奇异抛物方程初边值问题的解的分析[D];中南民族大学;2013年
8 刘胜利;四阶抛物方程的一个有限差分并行格式[D];吉林大学;2004年
9 赵艳冰;一类具变指数反应项的半线性抛物方程的Fujita指标[D];吉林大学;2014年
10 赵建斌;正交小波基的构造及其在抛物方程数值解中的应用[D];湘潭大学;2007年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978