基于线性抛物化稳定性方程的后掠翼边界层内横流稳定性研究.pdfVIP

第27卷 第5期 计 算 物 理 CHINESE JOURNALOFCOMPUTATIONALPHYSICS V01．27．NO．5 2010年 9月 Sep．，2010 文章编号 ：1001—246X(2010)050665—06 基于线性抛物化稳定性方程 的后掠翼 边界层 内横流稳定性研究 左岁寒， 杨 永， 李 栋 (西北工业大学翼型叶栅国家重点实验室 ，陕西 西安 710072) 摘 要 ：采用有限差分法求解考虑模 型 曲率 的线性抛物化稳定性方程 (LPSE)，分析无 限展长后掠翼边界层 内 的横流驻波不稳定 ，并与实验结果进行对 比，研究 LPSE方法 的模拟效果及其适用范围．研究表 明，在横流驻波 扰动增长的初期 ，LPSE能够准确 的预测扰动 的 e 曲线 ，较好地描述边界层 内的流动结构和扰动形态 ；当扰动 增长到足够大时，扰动 的高阶项不能再被忽略，LPSE的线性假设不再成立 ，需要采用非线性的方法 (NPSE) 来分析该状态．计算分析发现 ，模型 曲率和边界层非平行性对后掠翼边界层 内横流驻波 的稳定性分析影响很大， 影响程度与雷诺数无关．对于本文研究的模型，曲率对边界层 内的横流扰动起着稳定的作用 ，而非平行性对扰 动起不稳定 的影响． 关键词 ：线性抛物化稳定性方程 ；无 限展长后掠翼 ；横流驻波不稳定 中图分类号 ：V211．3 文献标 识码 ：A 0 引言 后掠翼边界层 的转捩主要 由横流不稳定触发 ¨，了解该不稳定进而找到有效的层流控制方法对于提高 飞行器的经济性有着重要的意义．横流不稳定包括驻波和行波两种扰动波 ，德 国DLR的Bippes等人 研 究表明，当来流湍流度较低时 ，横流不稳定 由横流驻波主导 ；而湍流度较高时，横流不稳定 由横流行波主导． 考虑到高空飞行时来流湍流度很低 ，后掠翼边界层 内的横流不稳定主要表现为横流驻波不稳定．美 国ASU 的 Saric及其同事 一¨通过实验和计算发现 ，对于后掠翼边界层内的横流驻波问题 ，传统的线性稳定性理论 (LST)得到的结果与实验结果有着很大的差异 ． 抛物化稳定性方程 (PSE)¨ u作为稳定性分析 的一种新方法 ，近年来愈发受到人们的重视．这是 由于 相对于传统的 LST方法，该方法能够考虑边界层 内非平行流和非线性的影 响，而计算难度又远小于直接数 值模拟 (DNS)方法．本文通过求解考虑模型曲率的线性抛物化稳定性方程 (LPSE)，分析无限展长后掠翼边 界层 内的横流驻波不稳定 ，并与实验结果进行对 比，研究 LPSE方法的模拟效果及其适用范围．此外 ，文中还 计算讨论 了模型曲率和边界层非平行性对稳定性分析的影响． 1 考虑 曲率 的线性抛物化方程 首先将流场中的变量都表示成平均量和脉动量两部分 M= +_五， = + ， =w+ ， P=p十 ， (1) 然后把 (u，，W，P)和 (五，，，P一)分别代人曲线坐标系的不可压 N—S方程并相减 ，略去脉动量的高次项 (非线 性项)，定义 = (，云，，)，可 以得到线性的扰动方程 厂 + +B Ot dx Oy+c警dz+D一E【ax+‘一，a+EOz】=0． (2) 本文研究的是无限展长后掠翼边界层 ，因此流场 中的平均量 (五，， ，)关于展向 的梯度都为零 ，最终简化 得到式 (2)的各个系数 ： 收稿 日期：2009—07—17；修 回日期 ：2009—12—23 作者简介：左岁寒 (1985一)，男，江苏淮安，博士生，主要从事计算流体力学和流动稳定性研究，西北工业大学 111信箱 710072 计