二维双曲型方程的分组并行格式及其数值实验.pdfVIP

第27卷第4期 贵州大学学报(自然科学版) V01．27No．4 ofGuizhou Journal University(Natural 2010年8月 Sciences) Aug．2010 文章编号1000-5269(2010)04-0009-06 二维双曲型方程的分组并行格式及其数值实验 刘轶中’，张大凯 (贵州大学理学院，贵州贵阳550025) 摘 要：本文构造了求解二维双曲型方程u。+口“。+bu，=0的初边值问题的一组分组并行算法 验证了理论结果． 关键词：二雏双曲型方程；分组显式格式：稳定性：截断误差 中图分类号：0241．82文献标识码：A 设数学模型为： ——+口——+D——2塑+口塑+6塑：0 0≤舅≤L，0≤y≤L；t≥0，a0，60 挑 缸 耖 u(o，Y，t)=八0，，，，f) 0≤Y≤L；t≥0 M(z，0，t)=g(x，0，f) 0≤茗≤三：t≥0 u(x，Y，0)=h(x，Y，0)0≤茗≤L；0≤y≤工 由于对此方程的计算具有极强的方向性且仅 下面针对．，奇偶性，设计如下几种分组并行的 具有单边边界条件，故对二维双曲型串行差分格式 计算方法： 的并行化是一件很不容易的事情．从已有文献看 1．1 GE格式 当．，为偶数时，为了设计分组并行差分格式， (见文[1]、[2]、[3]、[4])，尚未发现二维双曲型 方程的并行化格式．本文利用一个二维显格式、两 用如下的四个格式构造逼近式(1)的并行差分方 个二维显隐格式和一个二维隐格式构造了一组分 程组： 组显式格式，格式的局部截断误差阶一般为D(r+ fu譬1=rz‘0一。+兀。乞。J+(1—2r)耻0 h)，稳定性条件为0r≤1． 1构造并行差分格式 设问题(1)的解u(x，Y，t)充分光滑，以缸、 【一九‘／“t+ll-，一九‘2：l+(1+2r)Ⅱ甜√+I=Ⅱ11J+I △，，和r分别为x、)痢t方向的网络步长，其中：缸 (2) r r =三，Ay=三，m和n为正整数，在网点(毛，乃， m 乃 其中r=簪于是在第k+1时间层上的4个函数 ，‘ t。)处的网格函数l‘(毛，乃，t。)记为[1‘]乞，其近似 值可由第k时间层上的8个函数值显式地进行计 值记为u：，．其中≈=ih(i=0，1，…，m)；乃=正(， 算： =0，1，…，，1)；气=打(后=0，1，…)．为简单起见， 本文取口=b且m=n，则Ax=Ay=h，其中h= 量 m‘ 收穑日期：2010一07—30 基金项目：贵州省省长优秀科技人才项目(黔省专(2008)19)；贵州大学引进人才科研项目(XOS5024) 作者简介：刘轶中(1973一)，男．湖南衡东人，硕士，实验师，研究方向：微分方程数值解，Emsil：liuyiz2004@126．toni． ·通讯作者：刘轶中，Entail：liuyiz2004@126．corn． ．10· 贵州大学学报(自然科学版) 第27卷 u21=ru。kJ．1+rz‘I．1-，+(1—2r)M乞 (3) u搿．，=南uo南u‰+南