跨音速绕流的数值模拟和验证.pdfVIP

LHD 2011 年度夏季学术研讨会 157 跨音速绕流的数值模拟和验证 郑洪伟，杨国伟 （中国科学院力学研究所高温气体动力学国家重点实验室 （筹），北京海淀区 100190） 摘要 高雷诺数粘性流的的气动升阻力的精确计算对于工程问题具有很重要的意义，因此有必要对 CFD 结果进行验证。本文将工程上广泛使用的几个跨音速流场算例（RAE2822 翼型和 DLR-F6 翼身组 合体）作为标准算例， 对我们自己开发的“基于任意形状的多面体单元的混合网格”的高性能气 动仿真软件平台 (PolySim)进行了验证。具体对压力分布、激波位置及升阻力与实验结果进行了细 致的比较，结果表明，用PolySim 计算所得结果均与实验结果较吻合，差异较小。 关键词 Spalart-Allmaras ，RAE2822 翼型，DLR-F6 翼身组合体，跨音速绕流 模型，因此不具普遍性，对复杂工程问题计算 引 言 结果不是很好，而且不适合非定场问题的计 计算流体力学 (CFD) 在过去的半个多世纪 算。与 RANS 不一样，大涡模拟(LES)求解的是 已经取得了很大的进展。它已经成为人们在很 空间平均量。它对大尺度涡进行直接计算，而 多工程应用的设计和优化过程中的一种不可或 认为小尺度涡对大涡的作用具有普遍均匀性和 缺的工具。然而对于高雷诺数粘性流的的气动 各向同向性，一般用亚格子 (sub-grid)模型来 升阻力的精确计算仍是一个艰巨的任务。国际 模拟。小尺度涡从大尺度涡那里获得的能量， 上对于这方面的研究, 主要有欧洲的计算空气 最后以湍流粘性耗散掉。但是 LES 对网格的分 动 力 学 研 究 项 目 ECARP[1] （European 辨率的要求很高，特别是壁面层的网格基本与 Computational Aerodymamics Research DNS 差别不是太大，因此计算成本也很高。混 Project）和 AIAA 应用气动技术委员会（AIAA 合法是结合了 RANS 和 LES 两种方法的优点， Applied Aerodynamics Technical 也就是：RANS 对于附着流的计算成本低和 LES Committee ）召开的四次会议 DPW[2] (Drag 对于分离涡的优势。混合方法一般是对于靠近 Prediction Workshop)。近年来，国内也先后 壁面地方采用非定常雷诺平均 (URANS)来计算, 组织了几次全国范围内的 CFD 软件可信度确认 而远离壁面采用 LES 来计算。 因此, 混合法 工作，如 2004 年 8 月的“CFD 统一算例研讨” 很适合流场中有大分离的壁面湍流问题。对于 活动和 2005 年 6 月的 CT-1 大攻角气动特性研 附着流的工程问题, 目前一般还是采用RANS。 讨会。 根据我们多年的积累，我们已经初步开发 影响气动升阻力的计算精度的因素大概有 了一套“基于任意形状的多面体单元的混合网 三个方面, 一个方面是湍流模型, 另外一个方 格”的高性能气动仿真软件平台 PolySim。本 面是网格, 最后一个是数值方法。关于湍流的 文应用两个跨音速绕流问题 (RAE2822 翼型和 模型和计算方法，大致可以分成这四种：直接 DLR-F6 翼身组合体)的数值模拟来考察我们自 数值模拟 (DNS)、雷诺平均 (RANS)、大涡