二维对流扩散方程的一类交替分组方法.pdfVIP

2007年6月 山东师范大学学报(自然科学版) Jtm．狮 d Nmmal Sciel黜) V01．22No．2 第22卷第2期 Journal ShandongUniversity(Natural 二维对流扩散方程的一类交替分组方法* 冯青华1’ 王文洽2’一 (1)山东理工大学数学与信息科学学院，255049，山东淄博； 2)山东大学数学与系统科学学院，250100，济南∥第一作者29岁，男。硕士生) 摘要利用第二类Saul’yev型非对称格式给出了二维对流扩散方程的一类交替分组方法，该方法具有并行本性，易于程序实 现，并且是绝对稳定的．数值试验结果表明本方法具有较高的求解精度． ， 关键词交替分组方法；并行计算；绝对稳定 中图分类号0241 对流扩散方程是一类重要的偏微分方程，可以描述许多物理现象，研究其数值解法有重要的现实意义．数值解法中显式 方法稳定性差，隐式方法稳定性好，但不适合并行计算．因此，构造具有稳定性好且具有并行本性的差分方法就有着非常重要 的理论意义和现实意义．本文利用第二类Saul’yev型非对称格式，给出了二维对流扩散方程的一类交替分组格式，并分析了该 方法的稳定性．数值试验表明，该方法有较高的精度，适合于并行计算． 考虑二维对流扩散方程的初边值问题： l磊+如瓦+2荀2c1石≯+钇字了，，u≤菇《r2，u≤，，≤zu≤。≤1 f警+‰襄+k：考=e：寺+s：雾，，o≤菇《2，o≤，，≤2，o≤t≤r {“(菇，y,0)=，(菇，Y) (1) l H(o，，，，1)=gI(，，，t)，u(2，Y，I)=92(Y，1) L“(茹，0，￡)=^(鼻，t)，“(菇，2，t)=五(菇，I)， 为(i，．『，n)，精确解为u(x‘，竹，k)，记，=r，h2． l交替分组方法 利用第二类Saul’yev型非对称格式，构造本文的16种基本格式(图l一图16)． 令at=一号+了khr，啦=一号一警．设岛=．|}2=k，cI=e：=e． (2) (1+er)遽31+口luiga‘++^l，，+口lUI，．，I+1．1=一2n2ⅡllJ一2a2“：』一l+(1—3er)“：j一口l“：』+l一口1珥^+lJ， a2嵋二1+(1+2er)Ⅱni++11．，+2al／／‘n…+I，，+dIⅡn‘++11．』+l=一2n2珂：I，，．I—a2ujn．，+(1—2er)u11。j一口lu‘n+l，』+I， (3) 2az“a‘++11．，+(1+2er)um“+2l。J+口lunj++31．，+alMa“+2t，』+I=一20aⅡ12．卜I+(1—2￡r)n一+2．J一口l雎^+3．』一口In12．，+I， (4) az破：：．j+(1+￡r)遽：；．』+aI (5) (6) 啦筹1+(1+2er)uua‘，+』+il+口l嵋：：。』+l+2aI醒二：2=一2啦l矗l，，+l+(1—2er)畦，+I一口2畦』一口l以I，，+I，