pMultigrid和时间隐式方法对高阶精度DGM格式加速收敛的研究.pdfVIP

第十六届全国计算流体力学会议 p-Multigrid 和时间隐式方法对高阶精度 DGM 格式 加速收敛的研究 1 1,2 2 2 刘伟 ，张来平 ，李明 ，赫新 （1. 中国空气动力研究与发展中心 空气动力学国家重点实验室，四川 绵阳 621000 ） （2. 中国空气动力研究与发展中心 计算空气动力学研究所，四川 绵阳 621000 ） 摘 要：本文采用两层 V-循环的 p-Multigrid 策略，分别在顶层的 p 阶多项式和底层的 0 阶多项式进行光滑迭代， 光滑算子采用作者已经发展的基于 Newton 和 Gauss-Seidel 迭代的隐式方法。利用该方法，数值模拟了无粘的 NA CA0012 绕流、粘性的圆柱型绕流。计算结果表明，与显式的 Runge-Kutta 时间格式，以及单纯的隐式方法相比， 采用 p-Multigrid 技术和隐式相结合的策略，所需的迭代步数和 CPU 时间均在很大程度上得到减少，计算效率能 够提高明显。 关键词：间断Galerkin 有限元；时间隐式；p-Multigrid 技术 法[4]、无矩阵运算（Matrix-free ）Krylov 方法[5]、 0 引 言 LU-SGS 方法[6]等。作者此前将 Newton 迭代和 近年来，高阶精度计算方法越来越受到 CFD Gauss-Seidel 迭代相结合，发展了高阶 DG 格式的 工作者的高度重视。高阶精度格式具有较小的数值 时间隐式方法[7]，与显式时间格式相比，收敛速度 色散和耗散，适合求解多尺度流动问题，如湍流、 显著提高。 气动声学问题等等。在众多的高阶精度计算方法中， 另一方面，多重网格（Multigrid ）技术作为另 外一种计算收敛策略，也备受人们关注。多重网格 间断 Galerkin 有限元（DG ）方法，因其提高精度 技术分为基于网格的“h-Multigrid ”、基于精度的 不需要模板、计算精度高、并行算法容易实现、具 “p-Multigrid ”以及两种相结合的“hp-多重网格”。 有良好的数学特性等特点，已成为 CFD 领域卓有 近年来，国内外对多重网格研究已经比较深入，方 成效的数值方法之一，备受关注[1-2]。 法也渐趋成熟。C.M. Klai[8]等采用 h-Multigrid 技术 然而，随着空间精度的提高，DG 格式对应的 求解基于时空间断 Galerkin 有限元的 N-S 方程， 稳定性条件将越来越严格，时间推进步长受到明显 Krzysztof J. Fidkowski 等[9]采用 p-Multigrid 技术， 的限制。特别对于定常流场的计算问题，缓慢的 Cristian R. Nastase 等[10]采用 hp-Multigrid 加速 DG Runge-Kutta 时间推进需要花费大量的时间才能使 格式，都取得了不错的加速收敛效果。p-Multigrid 得问题收敛，计算效率低下。因此发展高阶格式的