热传导方程三层并行差分格式初始条件的计算.pdfVIP

计 算 物 理 Vo1．28．No．4 第 28卷 第 4期 2011年 7月 CHINESE JOURNAL OF COMPUTATIONALPHYSICS Ju|．，20l1 文章编号：1001—246X(2011)04-0488—05 热传导方程三层并行差分格式初始条件的计算 左风丽 ， 崔 霞 ， 袁光伟 (1 北京应用物理与计算数学研究所高性能计算中心 ，北京 100088； 2．北京应用物理与计算数学研究所计算物理重点实验室，北京 100088) 摘 要 ：给出二维热传导问题的三层差分格式初始条件 的一种显式计算方法 ，对于 由此形成 的内边界预估校正三 层并行差分算法 ，证明稳定性和收敛性定理．并行数值试验表 明，方法稳定 ，且与通常采用隐式格式计算初始 条件的方法相 比，易于程序实现 ；与 已有 的扰动算法相 比，能大幅度减小误差． 关键词 ：三层差分格式；初始条件 ；稳定显式计算方法 中图分类号 ：0241．82；0246 文献标识码 ：A 0 引言 三层时间差分格式是数值求解抛物型微分方程 的离散方法之一 ¨ ，如 Richardson格式 、DuFort-Frankel (DF)和蛙跳格式等．一些并行差分算法也利用了三层差分格式 ，如袁光伟等利用 DF格式构造 了解一维热 传导方程的分段隐式并行差分格式 ]．在此基础上 ，Gtinter等给出极坐标系下求解 电磁等离子体中强各 向 异性二维热传导方程 的显隐混合差分格式 J．近年来，袁光伟等又提出了基于内边界预估校正的并行差分 格式 卜 ，其 内边界预估利用了前两个时间步值的线性组合．对于三层差分格式，必须预知前两个时间步的 值．记最初两个时间步为t。和t，其上的值称为三层差分格式的初始条件．文 [5，8]对初值的处理是由已知 t。时刻的值利用隐式格式求解 t时刻的值，此时需要在整个空间区域上求解代数方程组，而其后续各时间步 的计算仅需要独立并行求解各子区域上的代数方程组；这样虽然求解精度有保证 ，但对于实际问题 的求解 ， 程序结构及其在并行机上实现均比较复杂．文[6]对极坐标下二维热传导方程 的DF三层差分格式，给出了 初始条件的扰算动法 ，其中引入 了一个辅助时间层 t卜”，即将t。之前的t卜”(即一△ )和t。作为计算的最初 两个时间步，数值结果表明，采用不适当的初始条件处理将会导致较大的数值误差． 本文对文[6]中t卜”时刻初始条件 的取法进行改进 ，即利用已知的t。时刻的值和该时刻值满足的隐式 差分格式 ，求出t卜”时刻的值．该方法能以显式形式计算 t卜”时刻的值 ，无需在整个区域上解代数方程组 ， 在初始化程序中即可完成初始条件的计算．与通常采用隐式格式计算t时刻值作为初始条件 的方式相 比，在 程序实现上要容易得多．本文对二维热传导问题的三层并行差分格式 ，分别采用显式计算方法 、扰动算法和 精确解赋值方法计算 t卜”时刻的初始条件 ，然后利用文[7—8]提出的内边界预估校正并行差分算法进行后 续时间步的计算 ，并且基于块数据结构在 JASMIN框架上进行 了并行数值试验．结果表明，该数值方法是稳 定的，其误差明显 比扰动算法低 ，几乎能达到由精确解作为初始条件所得到的数值解的精度． 1 初始条件的计算方法 1．1 问题和记号 以二维热传导方程为例，考虑如下初边值问题 ： 亡Iu at ()+()+ ， c 。，∈(0,TI， 收稿 日期 ：2010—08—12；修回 日期 ：2011—02—18 基金项 目：国防基础科研项 目(B1520110011)、中国工程物理研究院科学技术基金 (2010A0202010)、计算物理实验室基金、国家 自然科学 基金61033009)、863课题 (2009AA01A134，2010AA012303)及973课题 (2011CB309702)资助项 目 作者简介 ：左风丽(1971一)，女，河南，副研究员，硕士，主要从事大规模科学与工程并行计算研究． {通讯作者 E—mail：cuixia09@