二维三温辐射扩散方程组两层预条件子自适应求解.pdfVIP

计 算 物 理 Vo1．29．NO．4 第29卷 第4期 2012年 7月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS July，2012 文章编号 ：100l-246X(2012)04—0475-09 二维三温辐射扩散方程组两层预条件子的自适应求解 周志阳 ， 徐小文 ， 舒 适 ， 冯春生 ， 莫则尧 (1．湘潭大学数学与计算科学学院科学工程计算与数值仿真湖南省重点实验室，湘潭市 411105； 2．北京应用物理与计算数学研究所 高性能计算 中心 计算物理重点实验室 ，北京 100094) 摘 要 ：针对实际应用 中若干典型三温线性系统，分析求解二维三温辐射扩散方程离散线性系统的代数两层预条 件子 (PCTL)的算法效率 ．结果表明，PCTL的算法效率与三个温度之间的耦合强度 以及单温子系统对角 占优性强 弱程度有很大关系．为此 ，通过刻画三温线性系统的耦合强度和单温子系统对角 占优性特征，提 出一种 PCTL中子 系统的 自适应求解算法 ．数值结果表 明，可 以显著改善 PCTL的算法效率，对于实际数值模拟应用 中37个典型三 温线性系统 ，相对于经典 AMG算法 ，算法整体加速 2．5倍．数值实验表 明算法具有很强的鲁棒性 ． 关键词 ：二维三温方程 ；辐射扩散 ；代数多重网格(AMG)；预条件子 ；自适应 中图分类号 ：O241．6； TP338．6 文献标识码 ：A 0 引言 在惯性约束聚变等高能量密度物理领域 ，数值模拟是主要研究手段 I2．其 中，求解辐射流体力学方程 组 的三温辐射扩散方程线性解法器是该类数值模拟 中的重要组成部分 。．在很多情况下 ，三温线性解法 器的CPU时间占90％以上．近年来 ，很多作者针对二维三温 (2D3T)方程做了大量工作 ，发展 了很多求解算 法，推动了该课题的研究 ．随着该类应用大规模数值模拟的开展，线性解法器仍然是制约该类数值模拟 性能的主要瓶颈 ’ ． 代数多重网格 (AMG)方法 具有最优性质 ，是科学计算领域求解大规模稀疏线性代数方程组最有效的方法 之一．对于标量偏微分方程的离散系统 ，AMG算法设计相对成熟．然而，对于三温辐射扩散方程，由于三个温度相 互耦合 ，使得经典AMG(C—AMG)方法难以适应．在文[3]中，作者提出一种特殊的粗化策略，将电子温度视为粗变 量，光子和离子温度视为细变量．基于此，通过构造特殊的插值算子和光滑子 ，设计了一种基于物理量粗化的代数 两层预条件子 (PCTL)．数值实验表明，PCTL具有 良好的算法可扩展性．在 PCTL算法中，三温耦合方程化为三个 单温方程和一个粗网格方程共四个子系统的求解，其算法效率取决于这些子系统的求解． 在文[3]中，没有考虑子系统求解对 PCTL迭代效率的影响．本文基于辐射流体力学实际数值模拟中的 三温线性系统 ，测试并分析了PCTL的算法效率．结果表 明，子系统的求解精度与 PCTL的收敛速度之间需 要取得平衡．通过对三温系统 的矩阵性质和结构进行分析发现，各子系统的对角 占优强弱程度和各耦合项的 耦合强弱程度对 PCTL的效率有直接影响．本文提出耦合 因子和对角 占优因子的概念 ，用于刻画三温系统的 耦合强弱和对角 占优强弱的程度．基于此 ，提出一种 自适应方法 ，用于 PCTL中的子系统求解．针对辐射流 体力学数值模拟应用中的37个典型三温线性系统，数值结果表 明本文算法的效率． 1 二维三温辐射扩散方程离散线性系统及其两层预条件子 1．1 三温辐射扩散方程及其离散线性系统 辐射扩散方程组 收稿 日期 ：2011—09—20；修回 日期 ：2012一】【1—17 基金项 目：国家重点基础研究专项经费(2011CB309702)，国家 自然科学基金 l117128110935003及湖南 省科学技术厅科技计划重点项 目(2011FJ2011)资助项 目 作者简介 ：周志 阳(1983一)，男 ，博士生 ，主要从事多重网格法与区域分解法的研究，E-mail：peghoty@ 163．coln 476