Burgers方程的一类交替分组方法.pdfVIP

立章龋号：1000端杆(2∞4)02．让13’08 Burgers方程的～类交替分组方法。 王文洽 (山东大学数学与系统科学学院．济南250100) (李继彬推荐) 摘要：对于Bursers方程给出了一组新的‰11l’yev型非对称差分格式，并用这些差分格式构造了 求解非线性Burgers方程的交替分组四点方法．该算法把剖分节点分成若干组，在每组上构造能够 独立求解的整分方程．因此算法具有并行本性。能直接在并行计算机上使甩．文章还证明了所给 算法线性绝对稳定．数值试验表明，该方法使用简便，稳定性好，有很好的精度． 关蕾词：Burgers方程；Saul’yev型非对称格式；交替分组四点格式； 线性绝对稳定； 并行计算 中圈分类号：0241 文献标识码：A 引 言 我们考虑非线性Burgers方程 碧+“再3u：e生，u5L，u‘2’’8oo$L，oI?t，so(1) 万+“瓦2。一c3X2’ 初始条件和边界条件是 o*L， Ⅱ(Ⅳ，O)=“*)， (2) u(0，t)=91(‘)，Ⅱ(L，t)=92(t)，0(‘T，(3) 这里e是常数． 些性质，因此讨论(1)的数值解法还有求解方程本身以外的学术价值．求解(1)的有限差分方 法和有限元方法巴经很多，Evans和张宝琳等人针对扩散方程和对流．扩散方程设计了一类有 Explicit)AGE和(．Alternating 并行本性的(AlternatingGroup Group 文献[6，9】给出的是(1)的全隐格式，而用的是交替分组显式迭代求解．在这篇文章里，我们用 第二类Saul’yev型非对称格式，构造了方程(1)的交替分组四点方法，其中对流项的离散化不 同于文献[6一lO]．该方法线性绝对稳定，使用简便，适合在高性能多处理器的并行机上直接 使用．数值试验表明交替分组四点方法有很好的精度． ·收稿日期：2001．05—08；修订日期：2003-06—27 基金疆目： 国家教育部博士点专嘎基金资助项目；山东省自然科学基金资助项目 (Y2003A04) on)． 213 芏 文 治 1交替分组四点格式 令h和At分别表示*和t方向的网格步长，h=L／m．m为正整数．n?表示问题(1)． t¨／2)。并定义下面的记号 丸“?=(“0l一．t)／h，也“?=(“?一n卫1)／h。 也Ⅱ?=(“Il—u曼J)／2h，盈ⅡP=(u：+1一ut)／△t． 则方程(1)在节点(％，I。+l／2)处的4个第二类Saul’yev型格式是 为了对(4)一(7)进行线性化，式中u：“／2按文[9]的思想处理，如图1．取 Ⅱ{“／2=“；≈——■再—兰l_—一_口?， (8) 1+嵩(u}一ulI’ 记r=At／(2h2)，于是(4)一(7)式可以改写成下面的形式 (1+m)u1”一r(e一口批／2)U“l+i= (9) r(2e+口孤)M量l+(1—3re)“?+r(e一口弘／2)吨l， 一r(E+口孤／2)U¨i-：+(1+，￡)u1+1= (10) r(e+口胁／2)ull+(1—3m)“?+r(2￡一a}h／2)“ll， 一r(e