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

色散方程的一类高精度并行算法

张青洁  
【摘要】: 色散方程作为重要的数学物理方程之一一直受到业内人士的普遍关注,在非线性波及孤立子理论的物理问题中,也一直占有相当重要的位置。鉴于色散方程在物理领域重要的应用价值,人们已经开始广泛的关注其数值解法的研究,许多专家、学者在这方面也已经做了不少工作,例如文献([1—13])。其中,文献([4])详细讨论了各种差分格式及其相应的稳定性情况。这些差分格式可分为显式格式和隐式格式两大类。我们知道显格式形式简单并适于并行计算,但其稳定性条件通常比较苛刻、不易实现。尽管人们在稳定性条件方面做过一些改进([7,8]),但是这些改进还是极其有限的。比如文献([5])讨论了一类两参数的恒稳显格式,但是参数的选择还是要满足比较复杂的条件的。相对于显格式而言,隐式格式虽然具有稳定性好的优点,然而它却不能直接应用于并行计算。随着信息时代的到来和计算机的蓬勃发展,并行计算以其快速解决大型且复杂的计算问题的特点迅速吸引了很多业内人士的普遍关注([10—37]),因此,怎样找到一个稳定的、适合并行计算的数值求解方法,便成了相关研究人员亟待解决的重要问题。 关于交替分组算法的研究是随着并行数值计算在计算机上的广泛应用而逐步深入的,目前,两类主要的并行算法就是:交替分组方法([2,4,9—24])和区域分裂算法([38—49,63,64,66])。前者是无条件稳定的,所以我们通常可以采用比较大的时间步长,而后者是条件稳定的,因此,在使用过程中,我们通常需要选取比较小的时间步长来进行计算。交替分组方法已成为行之有效的并行数值算法之一,它不但是绝对稳定的,而且还具有本质并行的特性。 例如,抛物方程的并行差分解法已经在很多文献中被广泛研究([15,23,25,63—67]),关于扩散方程和对流-扩散方程的并行差分算法的研究也已经有了不少成果([15—17,19—21]),近些年来,交替分组方法的研究和应用又逐渐扩展到了三阶色散方程、KdV方程等领域([1—13,22,50])。 不过,对于三阶色散偏微分方程而言,这样的并行差分算法并不多见。 早在1983年,Evans和Abdullah首先提出了交替分组显式算法(AGE)([15,16]),后来,张宝林提出了交替分段显隐算法(ASEI)([19])。2000年以来,朱少红又将交替分组显式算法(AGE)推广到了三阶色散方程的求解过程中来([11,12])。我们谈到的这些算法都是无条件稳定的,并且可以并行计算(随着计算机的蓬勃发展,并行计算也越来越多的被人们关注),不过,这些算法在空间上的收敛阶都只能接近2阶。 众所周知,提高数值解的精度也一直是数值解法研究人员的一个重要的目标和努力方向([50—62]),这也是我们在求解理论问题和实际应用问题过程中都不会改变的追求。 综上所述,本文作者在导师的悉心指导和精心培育下,提出了一类求解具有周期边界条件的色散方程的高精度、可并行、绝对稳定的算法。 在论文的第一章,作者介绍了色散方程的高精度并行迭代法。 在论文的第二、三章中,我们将给出四类Saul'yev型非对称差分格式来求解色散方程。基于这些Saul'yev型格式,我们又分别给出了求解带周期边界条件的色散方程的新的交替六点分组算法、新的高精度的交替显隐算法、高精度交替十二点分组算法以及4阶交替分段Crank-Nicolson算法。这四个新算法不仅具有无条件稳定和能在计算机上实现并行计算的特点,而且它们在空间上都具有4阶精度。通过数值算例,我们也容易看到,数值结果和理论分析是一致的。数值算例说明,新算法们在精度和稳定性上都优于算法AGE([11])和ASEI([12])。 论文的部分内容已在国际国内刊物上公开发表([68—72])。 全文共分为三章: 第一章介绍色散方程的高精度并行迭代法。 本章导出了一种数值求解色散方程的高精度交替分组迭代格式,此格式收敛速度快并可以在并行计算机上直接应用。 本章内容公开发表在([71])。 第二章介绍色散方程基于6点差分格式的高精度并行算法。 在第一节中,我们介绍了色散方程的高精度交替6点分组算法。 本节,我们将给出一类Saul'yev型非对称差分格式来求解色散方程。基于这些Saul'yev型格式,我们给出了求解带周期边界条件的色散方程的新的交替六点分组算法。这个新算法不仅具有无条件稳定和能在计算机上实现并行计算的特点,而且它在空间上具有4阶精度。通过数值算例,我们也容易看到,数值结果和理论分析是一致的。数值算例说明,新算法在精度和稳定性上都优于算法AGE([11])。 在第二节中,介绍了色散方程的一类新的高精度交替分组显隐算法。 本节针对色散方程提出的nAGEI新方法不但绝对稳定、本质并行,而且误差分析和数值试验表明,其数值解关于空间步长的收敛速度几乎是4阶的。通过与AGE([11])和ASEI([12])等方法的数值比较,我们容易看到本文方法确实具有更高的精度。本节内容已在《应用数学和力学》发表,请见[69]。 第三章介绍了色散方程基于12点差分格式的高精度并行算法。 在第一节中,我们介绍了色散方程的高精度交替12点分组算法。 近年来,随着并行计算机的发展,并行数值计算也越来越多的受到人们的关注和重视。像区域分裂算法一样([38—49,63,64,66]),交替分组方法也因其绝对稳定、本质并行的特点而日渐成为行之有效的并行数值方法之一。1983年,Evans首先提出了交替分组显方法(AGE)([15—16]),历经近20年的发展,交替分组算法的思想已经被成功运用到求解扩散方程([15—17,19—21])、色散方程([1—13])以及Kdv([22,50])方程等方程中去。但是,在已有交替分组方法材料中,它们的数值解在空间上都是有接近2阶的收敛速度。我们在本节给出的新算法不仅仍然具有绝对稳定、本质并行的优良特性,而且我们随后的截断误差分析和数值算例将表明新算法的数值解在空间上具有接近4阶的收敛速度。我们在数值算例中给出了本节算法与已有算法AGE([11])的数值比较。 在第二节中,我们介绍了色散方程的一类4阶交替分段Crank-Nicolson算法。 在本节,我们将给出一个新的4阶nASCN算法来求解色散方程,这个算法不仅绝对稳定,而且可以直接应用到并行计算中去。事实上,交替分组方法是随着并行计算机的发展而蓬勃发展起来的。目前,两类主要的并行算法就是:交替分组方法([2,4,9—24])和区域分裂算法([38—49,63,64,66])。前者是无条件稳定的,所以我们通常可以采用比较大的时间步长,而后者是条件稳定的,因此,在使用过程中,我们通常需要选取比较小的时间步长来进行计算。1983年,Evans率先提出了交替分组显式计算方法(AGE),此后,又有人提出了交替分组显隐算法(ASEI)以及交替分段Crank-Nicolson([19,20])(ASCN)算法。近年来,我们也开始看到交替分组方法被应用到求解色散方程、Kdv方程等方程中。不过,在已看到的交替分组算法文献中,几乎所有算法的数值解在空间上都只能接近2阶。新算法nASCN不仅格式无条件稳定,而且还具有本质并行的特点。此外,我们随后的截断误差分析和数值试验表明新算法可在空间上达到4阶收敛,这比已知的AGE([11])和ASEI([12])都精确。本节内容已发表在《Computers and Mathematics with Applications》请见[68]。 考虑到一维算法为高维算法的基础和依托,我们也正在将本文的高精度并行算法推广应用到Burgers方程、Kdv方程以及相关的二维问题的求解过程中去。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 张大凯;求解色散方程的两类阶梯形半显格式[J];数值计算与计算机应用;1991年01期
2 黎益,李北杰;关于色散方程u_t=au_(xxx)的两个显式差分格式[J];计算数学;1986年03期
3 戴伟忠;色散方程的显式与半显式差分格式[J];厦门大学学报(自然科学版);1988年01期
4 秦孟兆;色散方程u_t=au_(xxx)的差分格式[J];计算数学;1984年01期
5 王殿辉;色散方程u_t=au_(xxx)的一类具高稳定性的三层显式格式D_3[J];应用数学;1994年01期
6 王文洽,付树军;带扩散项色散方程的交替分组差分方法(英文)[J];计算物理;2005年01期
7 袁忠信;;推广的非线性色散方程的弱解[J];郑州大学学报(理学版);1984年02期
8 武蔚文;解色散方程u_t=αu_(xxx)的高精度绝对稳定的隐式差分格式[J];四川大学学报(自然科学版);1990年03期
9 明祖芬;;求解色散方程的两种高稳定指数型显格式[J];贵州大学学报(自然科学版);1991年02期
10 黎益,王红川;色散方程u_t=au_(xxx)的八点显格式[J];四川大学学报(自然科学版);1993年01期
11 黎益 ,李北杰;逼近色散方程u_t=au_(xxx)的高精度差分格式[J];四川大学学报(自然科学版);1985年04期
12 黎益,廖晓峰;色散方程的任意阶精度的显式差分格式[J];四川大学学报(自然科学版);1993年04期
13 王昌标 ,張世昌 ,刘盛纲;电子回旋脉塞色散方程[J];电子科技大学学报;1986年02期
14 曾文平;解色散方程u_t=au_(xxx)的一族绝对稳定的高精度差分格式[J];计算数学;1987年04期
15 林鹏程;色散方程u_t=au_(xxx)两类绝对稳定的两层半显式格式[J];福州大学学报(自然科学版);1992年02期
16 韩凌燕;曲小钢;;色散方程的两个隐式差分格式[J];山东科学;2008年01期
17 胡文忠;任意多个磁层的层状结构中MSFVW和MSBVW的通用色散方程[J];物理学报;1989年03期
18 王殿辉;解色散方程的一族三层显格式[J];东北大学学报(自然科学版);1992年05期
19 黎益;色散方程的四点显式差分格式[J];应用数学和力学;1993年03期
20 黎益;色散方程u_t=au_(xxx)的两类显式差分格式[J];四川大学学报(自然科学版);1994年03期
中国重要会议论文全文数据库 前10条
1 姜遵富;周文表;;采用Webber-Schafheitlin积分分析毫米波介质插入波导[A];1995年全国微波会议论文集(下册)[C];1995年
2 黄凌雄;赵斌;张戈;陈建中;黄呈辉;庄乃峰;魏勇;位民;;Yb:YP_(0.05)V_(0.95)O_4晶体折射率测量[A];大珩先生九十华诞文集暨中国光学学会2004年学术大会论文集[C];2004年
3 赵爽;吴福全;;石英晶体的色散方程及折射率温度系数[A];江苏、山东、河南、江西、黑龙江五省光学(激光)联合学术'05年会论文集[C];2005年
4 李伟;魏彦玉;刘盛纲;;有限磁场作用下填充等离子体波导中的高功率微波辐射的特性研究[A];2001年全国微波毫米波会议论文集[C];2001年
5 麻来宣;张厚;张晨新;水现辉;;等离子体干扰机理分析[A];2009年全国微波毫米波会议论文集(下册)[C];2009年
6 李伟;谢鸿全;杨福丽;魏彦玉;刘盛纲;;磁化等离子体填充的切伦可夫脉塞的特性分析[A];中国电子学会真空电子学分会第十三届学术年会论文集(上)[C];2001年
7 于新华;蒙林;鄢杨;张红卫;;波纹内导体同轴慢波结构特性的研究[A];四川省电子学会高能电子学专业委员会第四届学术交流会论文集[C];2005年
8 王龙;张玲芬;;左手介质平板波导模式色散特性分析[A];2007年全国微波毫米波会议论文集(下册)[C];2007年
9 张祖存;刘文鑫;梁正;杨梓强;;毫米波段Smith-Purcell辐射特性研究[A];第九届真空技术应用学术年会论文集[C];2006年
10 覃志东;宫玉彬;王文祥;;大功率毫米波环板行波管的小讯号增益研究[A];中国电子学会真空电子学分会第十三届学术年会论文集(下)[C];2001年
中国博士学位论文全文数据库 前10条
1 刘晓风;色散方程和退化松驰Dirichlet问题的若干问题[D];浙江大学;2001年
2 张青洁;色散方程的一类高精度并行算法[D];山东大学;2009年
3 王晓英;荷电液体雾化的研究与应用[D];江苏大学;2008年
4 殷海荣;光子晶体行波管研究[D];电子科技大学;2007年
5 徐桂香;色散方程(组)的低正则性[D];中国工程物理研究院;2006年
6 翦知渐;用右手、左手材料制成的光子晶体的光学特性的研究[D];湖南大学;2007年
7 吴奕飞;几类非线性色散方程和波动方程的低正则性理论[D];华南理工大学;2010年
8 信江波;a轴向单晶光纤模式传输特性的理论研究[D];哈尔滨工业大学;2009年
9 路志刚;矩形波导栅行波放大器的研究[D];电子科技大学;2008年
10 王彬;等离子体耦合腔行波管高频特性及电子枪的研究[D];电子科技大学;2008年
中国硕士学位论文全文数据库 前10条
1 付吉美;五阶色散方程的并行计算方法[D];山东大学;2010年
2 韩凌燕;色散方程的两个八点差分格式[D];西安建筑科技大学;2008年
3 吴蕾蕾;等离子体加载周期慢波结构[D];电子科技大学;2007年
4 张祖存;毫米波段平板矩形光栅辐射特性研究[D];电子科技大学;2007年
5 胡莉;同轴腔回旋自谐振脉塞(CARM)的线性研究[D];西南交通大学;2004年
6 李悦宝;周期性介质切伦柯夫脉塞的研究[D];电子科技大学;2007年
7 滕军;倒同轴膜片加载圆波导高频特性研究[D];电子科技大学;2008年
8 杨娜;圆柱波导TM模式回旋自动谐振脉塞(CARM)的线性特性[D];西南交通大学;2009年
9 刘洪华;三阶mKdV方程的两类并行差分算法[D];山东大学;2007年
10 王辉;时间分数阶色散方程的高阶紧致差分方法[D];中国海洋大学;2013年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978