轴对称变几何进气道初步研究.pdfVIP

1.1 数值仿真模型及网格 为了验证本文所采用的数值仿真方法对模拟圆锥凹腔流动及其外部流场的准确性，首先对该数值仿真 [5,6,7] 方法进行校核。选用日本东京大学航空航天系Yusuke Maru等人所做的圆锥凹腔实验数据 ，图1给出 了该圆锥凹腔实验模型示意图。实验模型总长为 594mm，底部圆柱直径为 120mm，圆锥的半锥角为8 度， 距锥顶下游295mm处开前后壁面均垂直于圆锥轴线的凹腔，其底面与锥面平行。凹腔的几何特征由凹腔前 壁面深度D和前后壁面间的轴向长度L描述，实验模型凹腔前壁面深度D为15mm，不同尺寸的凹腔实验模 型，凹腔距前缘距离固定不变，通过改变凹腔长度 L，使 L/D 分别为 0.5，1.0，3.7。实验时，在距凹腔 后壁面下游404mm处布置总压探针，该总压探针可沿垂直于锥面高度0~7mm 内移动以测量该截面气流总压 分布。 图1 单凹腔圆锥实验模型示意图 图2 数值仿真网格及计算边界条件 图2给出了单凹腔圆锥实验模型的结构化网格及计算边界条件，为准确模拟凹腔内部和附面层内部流 动，对近壁面及凹腔内网格进行了加密处理，保证附面层内部区域至少有3个网格节点，整个计算域网格 总数约为3.3万。采用基于有限体积法的二维雷诺平均N-S方程求解器，标准k 两方程模型，标准壁 面函数处理。设置压力远场，压力出口，固壁和轴边界等边界条件，考虑流动的对称性，选取了整个流动 区域的一半作为计算域。采用与风洞实验条件一致的来流参数 （见表1）进行数值模拟。 表1 风洞实验条件 来流马赫数 来流总压 来流总温 来流雷诺数 测试持续时间 3.0 0.441MPa 288K 3.48 10 7 90 seconds 1.2 圆锥凹腔数值仿真结果与实验数据对比 按上述数值仿真方法，对单凹腔圆锥流动进行了数值模拟。图3给出了圆锥凹腔实验的纹影图和数值 仿真马赫数等值图对比。从图可以看出，凹腔的前壁面和后壁面肩点处分别产生了一道扰动波并且在凹腔 内气流与主流相接触的界面上形成了气动型面，气动型面与前后两固壁平滑连接。数值仿真显示的流场结 构与风洞试验获取流场的纹影图相似。 2 0.08 0.2 0.4 0.6 1 1.6 2 2 .4 2.6 2.8 0.06 m m / Y 0.04 0.26 0.28 0.3 0.32 0.34