时空二阶精度隐式迭代NND算法的改进应用.pdfVIP

时空二阶精度隐式迭代NND 算法的改进与应用 *,+,1) *,+ * * 谢昱飞 ，袁先旭 ，陈琦，司芳芳 * + （中国空气动力研究与发展中心，绵阳 621000） （空气动力学国家重点实验室，绵阳 621000) [1] 摘要：在文献 所发展时空二阶精度的隐式迭代NND算法的基础上，吸收隐式Beam-Warming格式、隐式单/ 双时间步法的优点，参考Hsu 和 Jameson 发展的ADI-BDF 显隐混合格式，提出用双时间步法增强隐式迭代 NND算法的思路，以兼顾非定常计算的时间精度和效率。典型算例考核表明，改进后的隐式迭代NND算法有 较高的精度和效率，适合于非定常流动的计算研究。并运用该算法，针对细长体非对称涡流诱导的自激锥转 运动、细长方形截面飞行器非对称涡流诱导的自激滚转运动等有工程应用背景的问题，开展了非定常数值模 拟研究，取得了若干有价值的结论，表明改进后的算法是可行的，有一定的工程应用价值。 主题词：隐式迭代NND 算法；锥转运动；方形截面飞行器；数值模拟 前 言 随着现代飞行器机动性的增强，采用非定常CFD技术研究复杂非定常流动及动态特性的需求日益增长。 AntonyJameson提出的dual-time-steppingscheme极大地促进了非定常流动CFD模拟技术的发展，并得 [1] 到了广泛应用。文献 等也提出了一种时空二阶精度的隐式迭代NND算法，该方法采用GAUSS-SEIDEL迭代， 方程的右端项RHS只计算一次（n时刻），效率比较高，但不能消除线化、边界条件显式处理等带来的误差， 从而限制了大物理时间步长的选取。 [1] 本文在文献 所发展的时空二阶精度的隐式迭代NND算法的基础上，吸收隐式Beam-Warming格式、隐 式单/双时间步法的优点，参考Hsu和Jameson发展的ADI-BDF 显隐混合格式，提出用双时间步法增强隐 式迭代NND算法的思路，以兼顾非定常计算的时间精度和效率。首先，用隐式迭代NND算法获得一个名义 二阶的流场增量；然后用 dual-time-stepping 迭代方法消除线化、边界条件显式处理等带来的误差。典 型算例考核表明，改进后的时空二阶精度的隐式迭代NND算法有较高的精度和效率，适于非定常流动的计 算研究。针对细长飞行器的动态特性实际工程问题，开展应用计算研究，取得了较好的结果。 1 计算方法 在曲线坐标系下，基于完全气体假设且不考虑质量力的三维非定常Navier-Stokes方程组为： U E F G E F G v v v      (1) t       控制方程 （1）经过空间离散后，可写成下面的半离散形式： 1 J U R U 0 (2)   t