计算物理 ›› 2005, Vol. 22 ›› Issue (4): 283-291.

• 研究论文 •    下一篇

二维三温能量方程组离散求解的两个新预处理技术

吴建平1, 刘兴平2, 王正华1, 戴自换2, 李晓梅3   

  1. 1. 国防科技大学并行与分布处理重点实验室, 湖南 长沙 410073;
    2. 北京应用物理与计算数学研究所计算物理实验室, 北京 100088;
    3. 装备指挥技术学院, 北京 101416
  • 收稿日期:2004-11-18 修回日期:2005-01-07 出版日期:2005-07-25 发布日期:2005-07-25
  • 作者简介:吴建平(1974-),男,湖南新化,博士,助理研究员,从事科学计算与并行算法方面的研究.
  • 基金资助:
    国家自然科学基金重点项目(69933030);北京应用物理与计算数学研究所计算物理实验室基金(51479040103KG0201)资助项目

Two Preconditioning Techniques for Two-dimensional Three-temperature Energy Equations

WU Jian-ping1, LIU Xing-ping2, WANG Zheng-hua1, DAI Zi-huan2, LI Xiao-mei3   

  1. 1. School of Computer, National University of Defense Technology, Changsha 410073, China;
    2. Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China;
    3. Institute of Command and Technology of Equipment, Beijing 101416, China
  • Received:2004-11-18 Revised:2005-01-07 Online:2005-07-25 Published:2005-07-25

摘要: 二维三温能量方程离散后得到的稀疏线性代数方程组中,系数矩阵各行的对角占优性相差十分悬殊,矩阵元素相差也十分大.针对前一问题,提出了改善对角占优性的一个新比例化方法.针对后一问题,利用每次舍弃前计算多个行的技术提出了多行ILUT预条件方法.最后,将对角占优性改善技术、多行ILUT与对角元比例化技术、RCM排序联合使用于实际的能量方程离散求解中,取得了较好的加速效果.

关键词: 二犀三温能量方程, 预处理, ILUT, Krylov子空间迭代

Abstract: In a sparse linear system derived from two-dimensional three-temperature energy equations, the diagonal dominan varies greatly from row to row and so is the magnitude of the elements. We provide a new scaling method to improve the diagonal dominance. As ILUT is used to the derived linear system, it reserves the number of elements in each row and several relatively large elements related to the photon are dropped due to the large difference among elements. To improve the equality of the ILUT, we provide a new method named multiple row ILUT (MRILUT), in which multiple rows are computed before dropping. The provided methods are embedded into a preconditioned Krylov subspace method to solve the actual two-dimensional energy equations with three temperatures. The number of iteration at each time step and the total computation time are both greatly reduced.

Key words: two-dimensional energy equations with three temperatures, preconditioning, ILUT, Krylov subspace iteration

中图分类号: