动态预条件技术及其在ICF 数值模拟中的应用 - 中物院高性能数值模拟 .PDFVIP

技术报告 IAPCM-SCNS-HPC201508 动态预条件技术及其在ICF 数值模拟中 的应用 徐小文 2015/10/7 北京应用物理与计算数学研究所 中物院高性能数值模拟软件中心 技术报告 声明 本技术报告引用方式： “徐小文，动态预条件技术及其在ICF 数值模拟中的应用，中物院高性能数值模拟软件中心技 术报告，NO. IAPCM-SCNS-HPC201508，2015 年10 月.” 在正式发表之前，本报告内容可能修改。 技术报告 动态预条件技术及其在ICF数值模拟中的应用1 徐小文2 北京应用物理与计算数学研究所 摘要 在基于偏微分方程求解的数值模拟应用中，隐式离散后，最终通常归结为求解一个由 大量线性代数方程组构成的方程组序列。在实际模拟中，求解该方程组序列的线性解法器通 常采用静态预条件策略，即：在模拟全过程采用一个固定预条件子。考虑到迭代收敛的健壮 性，静态预条件策略通常选择Setup-based 型预条件，这类预条件需要一个额外的Setup 阶段， 其并行可扩展性较差，是主要的性能瓶颈。本文针对实际数值模拟应用中线性方程组序列的 求解，设计了动态预条件策略。对于方程组序列中每个线性方程，该策略按如下步骤执行： 首先，采用Setup-free 型预条件进行求解，并自动监控预条件收敛行为，通过设计切换条件， 在Setup-free 型预条件无法有效收敛时，切换到更健壮的Setup-based 型预条件子。其次，如 果需要切换到Setup-based 预条件，则重用方程组序列中上一个方程的Setup 过程，通过设计 刷新条件，在重用无法有效收敛时，基于当前矩阵重新Setup，构造新的预条件子。相对于实 际数值模拟应用中普遍采用静态策略，本文的动态策略可在保持健壮性的同时，尽可能减低 Setup-based 型预条件中的Setup 开销，提高并行可扩展性。针对ICF 数值模拟二维总体程序 和三维辐射流体界面不稳定性程序的典型模型，测试了动态预条件技术的效果，测试结果表 明，动态策略可有效减低Setup 开销，可明显改善解法器性能和数值模拟效率。 关键词 线性解法器，预条件技术，代数多重网格，ICF 数值模拟 1 引言 实际应用数值模拟通常按时间积分向前推进。对于非线性多物理耦合应用问题，时空离散 后，每个时间步包含多个物理过程的求解，在单个物理过程，又涉及非线性迭代，每个非线 性迭代需要求解一个线性代数方程组。通常，典型多物理数值模拟应用的一次完整模拟需要 求解的方程组数目可能达数万、数十万、甚至数百万，这些方程组形成一个待求解的方程组 序列： {A x = b, i = 1, …, n } (1.1) i i i eq 其中neq 为该序列中方程组数目。在很多大规模数值模拟中，求解该方程组序列的CPU 时间 1 本文得到国家自然科学基金 、国家973 项目( 2011CB309702)资助. 2 通信作者：xwxu@iapcm.ac.cn 1 技术报告 比例高达80% 以上，成为数值模拟性能的主要瓶颈。 受计算和存储复杂度的限制，迭代方法是当前大规模计算中求解该类方程的主要方法，例 如，Krylov 子空间迭代方法[1] （包括：CG、GMRES、BiCGS