基于DES的二维和三维空腔流动特性研究.pdfVIP

第40卷 第 1期 航 空 计 算 技 术 Vo1．40No．1 2010年 1月 AeronauticalComputingTechnique Jan．2010 基于 DES的二维和三维空腔流动特性研究 谭玉婷 ，伍贻兆 ，田书玲 (南京航空航天大学 航空宇航学院，江苏 南京210016) 摘 要 ：空腔流动广泛存在于航空航天工程中，对其流动特性的研究具有十分重要的工程意义。采 用数值方法对 比模拟了二维和三维空腔流动 ，控制方程采用N—S方程，空间离散采用有限体积方 法，对流通量计算采用Roe格式，非定常时间离散采用双时间步长方法，湍流粘性计算采用基于 sA 模型的DES方法。数值计算所得的三维空腔底面压强分布与实验结果一致 ，所得 的二维和三维空 腔 内的流动结构与相关文献中的分析相吻合。将二维和三维的计算结果进行比较 ，发现空腔底部 压强分布、空腔内的压强振荡、空腔内部流线图及声压级等均有所不同，得 出在空腔宽度有限时，横 向流动简化 了流动结构 ．削弱了振荡幅值 ，但使振荡过程复杂化 ，其三维效应不能忽略。 关键词：脱体涡模拟；Navier—Stokes方程；空腔流动；声压级；数值模拟 中图分类号：V221．3 文献标识码：A 文章编号：1671—654X(2010)01—0067—04 引言 型 3J，后来也有人采用一方程模型和两方程模型，但是 结果并不理想。近几年来，大涡模拟 LES、直接数值模 在航空航天工程 中常常存在空腔流动的问题 ，如 拟DNS也被运用到空腔流动的数值模拟中，但这两种 飞行器的内埋式弹舱，起落架舱等。高速气流流过空 方法需要 占用巨大的计算机 内存，耗费大量的计算时 腔，流动将出现 自持振荡，产生剧烈的压力、速度等脉 动 ，并产生强烈 的噪声 ，这些非定常振荡对于飞行器是 间。而DES对附体的边界层完全采用 RANS湍流模 型模拟，利用RANS可以有效处理附体边界层流动的 很不利的，因此对空腔流动特性 的研究具有深刻 的工 程意义。对于空腔流激励振荡机理的存在两种理解方 优点，不需要耗费很大的计算机资源；在远离物面的分 式，一种方式是把流动激起的振荡归因于空腔 口处的 离区域 内，对小尺度涡采用亚格子模型模拟，对大尺度 涡进行直接模拟，利用LES可以比较有效地模拟分离 剪切层与内部流动的相互作用 1j，这种作用导致 了空 流动的优点。通过结合 RANS和 LES各 自的优点， 腔后缘角处流体周期性的流入与流出，这种理解的出 DES可以比较快速而有效的模拟空腔流动。 发点是压力波的传播。另一种理解方式是由Rossit 本文以文献 [4]中的空腔为研究对象，采用基于 提出的，指出空腔 口处的剪切层 由一系列涡组成 ，这些 Spalart—Allmaras方程模型的DES方法，数值求解非定 涡从空腔前缘处脱落，并以一定的速度 向下游运动，一