基于脱体涡模拟方法数值研究武器舱的流场.pdfVIP

值模拟研究一直受到研究者的广泛关注。 从20世纪50年代起，国外就开始对空腔流动进行研究。六七十年代随着高速战斗机计划的引入，炸弹 等悬挂物、起落架及武器舱问题对飞行器飞行的影响开始逐渐显现，武器舱流动问题开始成为一个热门研 究话题。研究者们对此进行了许多的实验及数值模拟研究，来流马赫数由低亚音速到超音速，空腔的几何 形状由开式空腔到闭式空腔，还出现了少量的过渡型空腔，对压力脉动、声辐射等测量的实验研究相继展 开，并出现了对流动机理、振荡及噪声预估等理论探讨和对抑制振荡技术的研究。对于来流马赫数为亚音 速、跨音速的情况， Rossiter[4]根据实验在1964年提出了预估空腔流振荡频率的模型，这个半经验模型认 为空腔内部存在一个由脱落涡与声辐射组成的反馈循环。基于以上假设，提出了半经验方程，并被Heller [5] [6] 等人 进行修改。实验方面，Maureen B. Tracy ，E. B. Plentovich 等学者利用风洞实验对各种情况下的模 型进行了研究分析，给出了很多非常有价值的实验数据和结论。 [7] 国内对空腔流动的研究发展则相对较晚。在数值模拟方面，侯中喜等 对超音速空腔流动进行了研究， 空间采用WENO格式，时间采用R-K法，湍流模型为B-L模型，把研究结果与Rossiter公式进行了一阶模态频 [8] 率对比，罗柏华 采用TVD格式及MacCormack格式对亚音速空腔流动进行了研究，湍流模型为C-S模型，定 [9] 性模拟出了空腔流程的性质。在实验方面，罗柏华 等经过研究得出了特点条件下空腔振荡以声学模态共 振为主, 并发现了可预估这类 “不深不浅”空腔振荡频率的预估公式。 尽管国内外对武器舱流动展开了许多的研究，但仍然存在一些问题值得我们去深入探讨。本文在已开 发的非定常雷诺平均N-S方程代码的基础之上，采用基于S-A模型的DES方法，通过对武器舱流场的数值模 拟分析武器舱流场流动，比较计算所得到的数值结果，与已有结果相吻合。 1.数值计算方法 1.1控制方程 在笛卡尔坐标系下，三维雷诺平均Navier-Stokes 方程的无量纲形式的表达式为： Q (F F ) (G G ) (H H )  i v  i v  i v 0 （1） t x y z 其中，Q (, u, v, w, E )T 是守恒变量向量。对流项通量向量的表达式为： u  v  w  2 u p  uv  uw     2    Fi uv , Gi v p , H i wv .       uw  vw  w2 p 