基于DES的并列空腔流动特性数值模拟研究.pdfVIP

3O 2010年第 3期 基于 DES的并列空腔流动特性数值模拟研究 谭玉婷 李 青 (中航 工业一飞院，西安) NumericalSimulationofTwinCavityFlowsUsingDES TAN Yuting LIQing 摘 要 ：采用数值方法模拟了单个空腔及并列空腔的跨音速流动 ，控制方程采用 N—S方程 ，空间离散采用有限体积 方法，对流通量计算采用Roe格式，非定常时间离散采用双时间步长方法，湍流粘性计算采用基于 A模型的DES 方法。数值计算所得 的单个空腔底面压强分布与试验结果一致 ，所得 的并列空腔流动特性与相关文献 中的分析相 吻合 。 关键词 ：脱体涡模拟 (DES)；N—S方程 ；空腔流动 ；数值模拟 空腔流动广泛存在于航空工程中，如美国的F／ Stokes)方程作为控制方程 ，其积分形式为 A一22、JSF，俄罗斯的米格 1．44、苏一47等的内埋式弹 w d力+ (Fc-- ds— o (1) 舱。当武器被投放时，气流流过空腔会产生强烈的 加 自持性振荡 ，对控制武器 的电子系统及结构疲劳造 式中：W、F 和 F 分别为守恒矢量、对流通量和粘 成很大的危险性 。 性通量，它们的具体含义见文献[3]。 目前 ，关于空腔流动特性 的数值模拟研究主要 流动控制方程空问离散采用有 限体积方法 ，其 是针对单个空腔的情况 ，对并列空腔流动特性的数 中对流通量计算采用 Roe格式 ，时间离散采用双时 值模拟研究 比较少 。但是，在工程实际中，往往存在 间步长方法 。 多个空腔并列 的情况L1]，当两个这样 的并列空腔相 层流粘性系数由Sutherland公式计算 ，紊流粘性 对距离 比较近时，它们的流动不稳定性会相互影响， 系数 由基于 A(Spalart-Allmaras)紊流模 型 DES 自持性振荡会更加剧烈 。因此 ，在对单个空腔流动 (DetachedEddySimulation)方法计算 。 特性的研究逐步成熟的同时，进一步探讨并列空腔 1．2 DES方法 的流动特性具有重要的意义 。文献E2]对并列空腔 A方程模型的积分形式是 ： 自持性振荡的模态进行了试验研究，指 出当两个并 列空腔的位置 比较靠近时，会 出现耦合的 自持性振 dn+ 。 (．一．)dS=f。Qd (2) 荡 ，两个并列空腔中会产生耦合频率 (耦合频率定义 式中： 为并列的两个空腔中同时出现的频率)，与单个空腔 F = (3日) — 、 的自持性振荡有明显 的不同。因此，这里进行 了一 一 些数值模拟研究，并对单个空腔及两个空腔并列的 ， (+)(m + +也 ) 情况分别进行了数值模拟研究，对两种情况下流动