DES与RANS模型在压气机叶栅流动中的研究.pdfVIP

中国航空学会2007年学术年会 气动专题 71 DES 和 RANS 模型在压气机叶栅流动中的研究 陈美兰 李雪松 顾春伟 清华大学热能工程系, 北京 100084 摘要 本文对某压气机二维叶栅的内部流动及尾迹掺混进行数值模拟,讨论了 Spalart-Allmaras 模型在不同网格不 同计算格式下的模拟能力，在来流均匀，定常边界条件下，Spalart-Allmaras 模型标准 MUSCL 计算格式在 300 万 网格下定常计算得到的总压损失分布与实验符合良好；讨论了DES 模型常数对大涡模拟粘性计算的影响；采用基 C 于 Spalart-Allmaras 模型的分离涡模型（DES ）计算时， des 常数取为 0.2，得到了较好的计算结果。 C 关键词 RANS 模型 DES 模型 网格 des 常数 总压损失沿程分布 1 引言 近年来，基于传统的雷诺平均 Navier-Stokes 方程(RANS)的计算方法和湍流模型得到了极大 的发展，基本上能够成功的模拟叶轮机械内部的二维、准三维和三维真实流动。从应用范围看， 这些湍流模型为从精度较低但简单的代数模型到复杂、更为精确的、基于输运方程的一方程和二 方程模型。目前，这些模型仍然广泛应用于工程实际中。 但是，这些模型虽然大体反映了叶轮机械内部流动的基本结构，但都不能准确预测转捩、分 离、激波、尾迹、附面层和漩涡等流动结构以及他们之间流动的强相互作用[1]。 为了准确地模拟透平和压气机内部复杂流动现象，理论上只有采用直接数值模拟（DNS）和 大涡模拟方法(LES)方法。直接数值模拟方法不做任何假设，对 Navier-Stokes 方程直接求解，可 以精确模拟湍流，包括转捩，然而计算量太大，目前还不可能用于复杂的实际流动计算中。大涡 模拟方法对大尺度涡直接计算，对小尺度涡建立模型，由于小尺度涡的近似各向同性，易于建立 普适性的模型，可以准确模拟湍流与预测转捩位置，而且计算量介于雷诺平均方法与直接数值模 拟方法之间，在工程上可以接受，然而，大涡模拟方法在处理边界层内的流动时需要更细的网格， 处理能力并不比雷诺平均方法好。RANS/LES 的一种混合方法分离涡方法（DES）在主流区使用 大涡模拟方法，而在边界层区使用 RANS 方法，避免了大涡模拟方法在边界层区需要细网格的要 求，比大涡模拟方法经济，却有着与大涡模拟方法相近的精度。 基于以上选题背景和意义，本文将在压气机叶栅流动中研究目前应用广泛的雷诺平均模型和 亟待研究发展的分离涡模型。Spalart-Allmaras 模型是目前应用较为广泛的一种雷诺时均模型， 本文将比较讨论 24 万、54 万、300 万三种不同粗密的网格下以及标准 MUSCL 和 DM- MUSCL （同时适用于结构非结构网格）计算格式下Spalart-Allmaras 模型的模拟能力，为分离涡模拟做 准备，最终在分离涡模拟中选定 300 万网格设置、标准 MUSCL 格式进行计算。分离涡模型中的 模型参数直接关系着计算的精度，本文将讨论模型常数对分离涡模拟计算粘性的影响，并通过数 值实验验证模型常数 0.2 在本算例中的合理性。 2 数值方法 本文的计算采用自主开发的程序 HPPD。该程序提供了多种湍流模型，包括S-A ，k-g，k-ω， 大涡模拟 LES 模型和分离涡 DES 模型，采用有限体积法进行计算，可压不可压统一求解，非定常 计算采用双时间步长法加速收敛，采用 MUSCL、DM-MUSCL （可用于非结构网格）方法重构提 1 中国航空学会2007年学术年会 气动专题 71 高计算格式的精度。程序基于 MPI 实现并行化，并行模式采用区域分区办法。 本文在计算中选用了 S-A 一方程模型和 DES 分离涡模型，计算格式上采用 MUSCL 和 DM-MUSCL 两种重构方法。进行 3