非定常求解的内迭代初值对计算效率的影响研究.pdf

２０１６年２月 西 北 工 业 大 学 学 报 Ｆｅｂ． ２０１６ 第３４卷第 １期 Ｊｏｕｒｎａｌ ｏｆ Ｎｏｒｔｈｗｅｓｔｅｒｎ Ｐｏｌｙｔｅｃｈｎｉｃａｌ Ｕｎｉｖｅｒｓｉｔｙ Ｖｏｌ．３４ Ｎｏ．１ 非定常求解的内迭代初值对计算效率的影响研究 贡伊明，张伟伟，刘溢浪 （西北工业大学 航空学院，陕西 西安　 ７１００７２） 摘　 要：基于非定常流场的双时间求解方法，提出了一种提高非定常流场求解效率的有效策略。 通过 对前几个时刻的流场信息进行外插来确定下一时刻的迭代初值，使之更接近于收敛解，降低内迭代初 始残值，进而提高了非定常流场的求解效率。 将流场中每个点的守恒量在时间方向上进行泰勒级数 展开，设计了若干种外插格式。 采用绕圆柱非定常流动的求解来验证本方法的计算效果，并研究了不 同初值外插格式、空间离散格式、时间步长和收敛标准对初值外插方法效果的影响。 研究表明，在双 时间步法基础上，采用初值外插策略可普遍提高计算效率，其中交替外插策略可以普遍将求解效率提 高１倍左右。 相比于迎风格式，该方法对中心格式的求解效率提高更显著，并且对于不同的收敛标准 和时间步数均有非常明显的效果。 关　 键　 词：初值外插；非定常流场求解；泰勒级数展开；双时间步法 中图分类号：Ｖ２１１．３　 　 　 文献标志码：Ａ　 　 　 文章编号：１０００⁃２７５８（２０１６）０１⁃００１１⁃０７ 　 　 近些年来，随着计算机性能的提高，计算流体力 （ｐｈｙｓｉｃａｌ ｔｉｍｅ ｓｕｂｉｔｅｒａｔｉｏｎ）。 双时间步法引入了伪 学（ｃｏｍｐｕｔａｔｉｏｎａｌ ｆｌｕｉｄ ｄｙｎａｍｉｃｓ，简称 ＣＦＤ）在很多 时间项，对于每一个物理时刻看做定常流动来求解。 工作中得到广泛运用。 在实际的工程应用中，很多 双时间步法的好处是原来用于定常计算的预处理、 问题涉及到非定常流动，并且其计算量远大于定常 当地时间步、对角化、多重网格都可以应用到计算 流动，比如，气动弹性动力学问题、涡轮叶片绕流问 中；真实时间步长可以取得很大；双时间法可以降低 题、分离涡脱落问题、动导数的求解等等。 由此而促 通量线性化误差和近似因子分解误差，放宽了稳定 进了各种加速收敛技术和时间推进格式的发展。 ［９］ 性限制 。 这些方法都能有效的提高收敛速度。 为了提高计算效率，前人发展了很多加速收敛 流场求解效率不仅仅与收敛速度有关，还与初 措施，如当地时间步长、隐式残值光顺方法和多重网 值与收敛解间的差距有关。 因此，在基于双时间方 ［１］ 格方法 。 隐式时间推进方法可以在保证稳定性 法的非定常流场求解中，为了提高计算效率，通常将 的同时大幅提高时间步长，进而提高计算效率，常用 上一个时刻的流场作为下一个时刻流场求解的初 ［２］ 的有近似因子分解 （ＡＦ） 、超松弛迭代方法 值。 至于从初值角度进一步提高计算效率，相关研 ［３］ ［４］ （ＳＳＯＲ） 、直接求解方法和Ｎｅｗｔｏｎ⁃Ｋｒｙｌｏｖ方法 。 究鲜有报道。 本文从该研究角度出发，将前２～３个 在近似因子分解方法中应用最广的是 ＬＵＳＧＳ 算 物理时刻的流场信息外插求解出待求时刻的流场作 ［５］ ［６］ 法 。 后来Ｃｈｅｎ和Ｗａｎｇ 对ＬＵＳＧＳ进行了改进， 为伪时间步迭代的初场，这样就进一步缩小了初值 在有限体积方法中发展了ＢＬＵＳＧＳ 算法。 Ｎｅｗｔｏｎ⁃ 与收敛解的差距，以期减小迭代步数，提高计算