带扩散项四阶抛物方程一类新的交替分组显格式.pdfVIP

第42卷第lo期 数学的实践与认识 V01．42，No．10 2012年5月 MATHEMATICSINPRACTICEANDTHEORY May，2012 带扩散项四阶抛物方程一类新的交替分组显格式 郭阁阳-，刘播z (1．天津职业技术师范大学理学院，天津300222) (2．吉林大学数学学院，吉林长春130012) 摘要：首先给出逼近带扩散项四阶抛物方程初边值问题一类非对称差分格式，利 用该组非对称格式构造了一类新的交替分组显格式算法，并给出了截断误差分析和 绝对稳定性结论，最后给出数值实验． 关键词： 带扩散项四阶抛物方程；非对称差分格式；并行计算；绝对稳定性 l引言 近年来，抛物方程初边值问题的并行数值求解逐渐活跃，相关研究工作可见文献『1-9]． 早在80年代，Evans等人针对扩散方程提出了交替分组显格式算法[1l，该算法使得通过构造 新的差分格式来实现并行计算成为可能．此后，很多学者致力于这方面的研究．90年代初，张 法都具有高度并行性、绝对稳定性和较高的计算精度．随后，袁光伟等人将上述算法推广到 具有周期边界条件的三阶色散方程[4-5】，并用能苣法给出了绝对稳定性证日月-王文恰等人分 别针对带扩散项色散方程和二三阶非线性KdV方程分别构造了交替分组显格式算法和交替分 段显隐格式算法16-7】，并利用矩阵法证明了所给算法的绝对稳定性．文献【8—9】分别针对四阶 意的结果． 本文研究带扩散项的四阶抛物方程的并行化问题．利用逼近该方程的8个非对称差分格 式和恰当的交替技术，构造了交替分组显格式算法，进而给出了该算法的截断误差分析和绝 对稳定性结论，分析表明，所给算法具有良好的性质． 考虑如下带扩散项的四阶抛物方程 。∈ 7 I z∈【u，LJ， T】J (1)【1J 窑+a筹：E筹，z∈舭】，t∈[0 瓦+a否万2E万≯， 初值条件与边值条件分别为 u(x，0)=uo(x)，z∈[0，L】 (2) u(O，t)=札。。(o，t)=“(Lt)=ux：(L，t)=0，t∈[0，T】(3) 其中uo(x)是给定函数，口、E为给定常数． (K—s)方程的高阶线性部分．K—S方程在许多实际问题的研究中都有应用，例如等离子物理问 收稿日期：2011—12—20 2011) 186 数学的实践与认识 42卷 题、热传导及氧化反应的扩散、动力学、自由膜表面的流动等问题．关于K—S方程的各种解 析解和数值解法已有很多研究．然而，由于K—S方程的复杂以及方程阶数较高，关于该方程 的并行数值解法并不多见．本文以该方程的高阶线性部分为例，针对该部分建立了并行数值 求解方法，尤其是对高阶项百04万It的处理是至关重要的，希望本文的工作对并行数值求解K—S 方程做出重要贡献． 2交替分组显格式的构造 2．1非对称差分格式的构造 近方程(1)的8个非对称差分格式(4)一(11)，示意图如图1所示． —亡=L—亡丁 乳格式(4) b．格式(5) a’．格式(6) b’．格式(7) c．格式(8) d．格式(9) c’．格式(10)d’．格式(11) 图1非对称格式示意图 三嘣一(r+三)蹦+(·+；+；)哼“ (·+．；+妒1一(r+；)嘣+互r啪nq-1