三种典型翼型边界层稳定性对比分析.pdfVIP

第29卷 第4期 空 气 动 力 学 学 报 Vo1．29。 No．4 2011年 8月 ACTA AERoDYNAMICA SDⅡCA Aug．，2011 文章编号：0258·1825(2011)o4-o48]-05 三种典型翼型边界层稳定性对 比分析 王 菲，额日其太，王 强，苏沛然 (北京航空航天大学 能源与动力工程学院，北京 100191) 摘 要：为了便于工程上翼型的选取及优化，本文针对三个典型翼型边界层稳定性进行对比分析。首先采用e方 法对翼型SD8020进行转捩预测，并采用萘升华实验法检验数值计算的准确性。再利用该数值方法对比了三种典型 翼型上表面的转捩位置、不同弦长雷诺数下的压力系数分布、扰动增长率以及最不稳定扰动波的频率。分析结果表 明，在相同弦长雷诺数下，NACA0012最先发生转捩 ，而NACA64-204和RAE2822都保持着较长的层流区；压力梯度 对扰动增长有很大影响；同一翼型随弦长雷诺数增长，转捩点雷诺数变大，而不是工程上常采用的固定值。 关键词：线性稳定性；e方法；转捩；翼型；压力梯度 中圈分类号：V211．4 文献标识码：A 0 引 言 ，1 数值模拟方法 摩擦阻力可 占亚音速运输机总阻力的50％，而层 1．1 边界层平均流求解 流的阻力比相 同雷诺数湍流的阻力要小90％ILl，因此 采用无粘流和边界层方法求解基本流，利用 F—S法 减小表面摩擦阻力可以有效地降低飞行器的总阻力 。 变换基本流方程 ，取 =( i )，，’，这里u。表示 自 为了计算机翼的摩擦阻力从而进一步更合理地预测机 由流速度。相对于流函数 (，，，)的无量纲流函数为： 翼飞行特性，转捩位置的预测成为了重要前提工作。对 于高空飞行外界扰动较小的特点，线性稳定性方法的 叼)： (1) √“。 e方法仍然是 目前工程应用中最常用的转捩预测方 由此可得连续方程和动量方程为 ： 法。L．J．Runyan和D．George-Falvy通过修正 Mack的 稳定性方法计算了翼型在不同工况下的转捩位置并与 ( )十 ∥ + 1一 )] 实验对照，得到相应的最大增长率的值 ，为之后的相 关研究提供了丰富的数据。朱万琳等人在二维不可压 = (厂善一) ㈤ 流动的稳定性数值分析基础上进行了转捩的判定 ]。 边界条件为 =0 =厂=O；= 。’厂=1。其中导数表 为了便于工程上翼型的选取及优化，针对多个典型翼型 u,dx√一=u 示对 ’7求导，m表示无量纲梯度 ：m=羔 H ， = 边界层稳定性的对比分析具有重要的参考意义。首先 本文采用萘升华实验方法验证对称翼型SD8020数值模 手二-，，88mJ=：：_，，b：I++88mJ=：。。 拟的准确性，并拟定稳定性计算的幅值放大因子的值， 然后针对常规翼型(NACA0012)，自然层流翼型 (NA— 1．2 线性稳定性方程 CA64-204)和超临界翼型(RAE2822)这三种典型翼型 扰动的发展模式可以分成时间模式和空间模式，时 采用线性理论方法计算在不同弦长雷诺数情况下转捩