一种新的重排方法在二维三温问题求解中的应用.pdfVIP

计 算 物 理 第26卷第1期 V01．26．No．1 CHINESE OFCOMPUTATIONALPHYSICS 2009年1月 JOURNAL Jan．。2009 [文章编号】1001—246X(2009)Ol一0035—07 一种新的重排方法在二维三温问题求解中的应用 张 勇1， 黄廷祝1， 刘兴平2， 谷同祥2， 李厚彪1 (1．电子科技大学应用数学学院，四川成都610054；2．北京应用物理与计算数学研究所计算物理实验室，北京lo0088) [摘要】 针对二维三温问题离散化得到的稀疏线性方程组，提出一种新的重排技术(交替hyperplane重排)， 并结合1LU分解预条件技术，在Krylov子空间迭代法下进行测试．数值实验表明，在一定的填充模式及预处理消 耗大致相同的前提下，使用交替hyperplane重排技术的迭代收敛效果明显优于红黑排序、hyperplane排序等方法． [关键词】 红黑排序；Schur补；Krylov子空间；ILU分解；填充技术 [中图分类号】TP301．6 [文献标识码] A O 引言 在二维三温辐射流体力学方程组的数值计算中，求解能量方程是一个十分重要的环节，对其求解所占 的机时比重相当大(约占整个系统计算的80％以上)．因此，寻求一个收敛快、稳定性好的二维三温能量方 程数值解法是一个值得探讨的问题．二维三温能量方程可表示为非线性抛物型方程组，在二维柱对称 Lagrange坐标系下，其标量式经简化后可得如下方程组…： 钆等+(Pc+等)警+ⅢV耻虬(卜¨+“n¨， c¨等+(肼¨O—E,)O (1) V州村。一。(to_)， ％等+(Pr+鲁)警州iv耻“卜¨． 边值条件如下： ①固壁和轴对称条件：(F。)。=0，口=e，i，r，／7,为法向； ②自由平面边界条件：(F。)。=(F；)。=0，T，=T，(石，Y，t)． P。，P；，P，分别表示电子、离子、光子的压强，q．为人为粘性压强，c”c¨，c¨分别是电子、离子和光子的等容 1 比热，P是介质密度，y=土为比重，加。∞。．分别为电子和离子、电子和光子的能量交换系数，只= P —K(1D，T。)grad 对于二维三温问题离散得到的稀疏线性方程组，已有多篇文章讨论．这些方法主要采用Krylov子空间 迭代法求解，结合ILU，稀疏近似逆等预条件方法求解．其中一类方法是利用红黑排序，采用Schur补的方 法，在保持矩阵的稀疏结构的同时，降低方程组的一半规模【2’31．在求解Schur补系统的预处理方法中，往 往在应用ILU等预处理方法之前，先采用尺度化、重排序等方法进行前预处理，形成混合预处理方法，以求 效果随着k的增大更好，因此，本文尝试研究不同的重排方法与ILU(k)结合使用的有效性，并试图寻找和 建立适合这类混合预处理方法的重排方法．我们主要关注基于结构的重排技术(如二维区域)，这样可以在 整个二维三温问题求解中采用相同的重排方法，在一定程度上有助于节省重排所花时间．本文通过一些数 【收稿日期]2007一∞一17；【修回13期]2008—02—08 【基金项目]北京应用物理与计算数学研究所计算物理国家级重点实验室合作基金(M1701100104WS01)资助项目 【作者简介】张勇(1974一)，男，四川资阳．讲师．博士，从事数值代数方面的研究． 万方数据万方数据 36 计 算 物 理 第26卷 Red-black)重排、 BRB(Block 常有效，在实现方面具有良好的并行特性． 1基于红黑排序的系统 由于二维三温问题的求解随着网格划分规模的