涡轮叶栅流动分离涡模拟验证研究的论文.pdfVIP

中国工程热物理学会 热机气动热力学 学术会议论文 编号：082034 涡轮叶栅流动分离涡模拟验证研究 顾春伟 李 艳 李雪松 清华大学热能工程系 北京：100084 TEL：010 E-mail：gcw@mails.tsinghua.edu.cn 摘要：本文采用 DDES 模型分析透平叶栅内部亚音速和跨音速流动。DDES 模型能有效且较精确地 计算亚音速流动，但计算跨音速流动时，出现修正函数f d “保护”边界层功能失效的问题，极大地降 低了计算精度；人为设定 RANS 和 LES 计算域，保证边界层完全被 RANS 计算域覆盖，计算结果得到很 大改善。 关键词：S-A模型 DDES分离涡模型 修正函数f d 激波/尾迹干扰 １．引言 一般认为，湍流的基本物理性质可由 Navier-Stokes 方程来描述，原则上可以通过 数值方法直接求解该方程而解决所有湍流问题，但受计算机计算能力、湍流流动复杂性 和工程应用的时效性要求等限制，在目前乃至可预见的将来，直接数值模拟( DNS) 实 际湍流仍不易实现，采用湍流模型仍然是研究工程流动最便利可行的技术途径。当前广 泛应用的是求解 Reynolds 平均的Navier-Stokes方程组即 RANS 方法。RANS湍流模型能 够较为准确预测无分离或小分离的边界层状态，但是在研究具有明显非定常气动特性的 模拟中，并不能保证有足够的计算精度和普遍性。大涡模拟（LES）是一种介于雷诺时 均（RANS）模式与直接数值模拟之间的高级湍流模拟方法，它可以精确模拟分离流动以 及几何相关的大尺度非定常运动，所花费的计算资源仅需 DNS 方法的很小部分。但是， [1] LES 方法模拟边界层时，由于网格要求高，需要的计算资源和直接数值模拟相当 。 近年来，由 Spalart[2]提出的分离涡（DES）模型得到发展，基本思想是在边界层使 用 RANS 模型，分离区与主流区使用 LES 模型，它既能很好地模拟主流湍流特性，获得 和 LES 近似的计算精度，但又将计算量减少到有可能接受的水平。因此，对 DES模型的 研究非常有意义，所以本文以涡轮叶栅内部流动为研究对象，进行了二维分离涡模拟的 验证研究。 ２． 数值计算方法和边界条件 本文的计算采用自主开发的程序 HPPD。该程序提供了多种湍流模型，包括 S-A 模型、 k-g 模型、k-ω模型等 RANS 模型、大涡模拟 LES 和分离涡 DES模型，采用有限体积法， 在控制体中心存储参数，数值格式采用可压不可压统一求解方法[3]，并用 MUSCL 重构方 法提高格式精度，非定常计算采用双时间步长法加速收敛，采用 MUSCL 方法对一阶基础 格式重构提高格式精度。程序基于 MPI 实现并行化，并行模式采用区域分区办法。 本文在计算中选用了改进的分离涡模型 DDES（Delayed Detached-Eddy Simulation），该模型基于 S-A 模型实现，计算格式上采用 MUSCL 重构方法得到了二阶 精度。DES 模型非定常计算采用双时间步长法加速收敛，物理时间步长为 2e-6，在每个 物理时间步长中进行 50 次内迭代，得到的收敛解即是当前物理时间步的瞬态解。本文 基于 MPI，采用区域分解的办法，在计算机上实现 HPPD 程序的并行化，完成并行计算。 2.1 DES 模型 DES 模型是由 Spalart 等人[2]于 1997 年提出的一种 RANS 模型和大涡模拟 LES 的混 合模型。最初 DES 模型基于 S-A 模型，它引入长度尺度d min (d , CDES Δ)，替代 S-A 涡