资源描述：
《多段翼型的一种分析方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、多段翼型的一种分析方法2000年8月第18卷第3期西北工业大学JournalofNorthwesternPolytechniea[UniversityAug.2000Vo1.18No.3一【多段翼型的一种分析方法詹浩,华俊张仲寅一————一(西北工业大学飞机工程系.陕西西安710072)V2ii,V2f1.摘要:采用定常二堆败拉方程和积分形式的可压缩附面屡方程作为基本方程,井将附面屡方程耦合到求解酞拉方程的奎局性牛横解法中,形成了一种求解多段翼型高升力系统的数值计算方法计算蛄果与风洞试验蛄果进行了
2、时比,两者基本吻合.计算实践表明,谊方法具有适用性强,收敛快的优点,是一种工程实用的计算方法美键词:!竺垫垦:附面层,兰篮墨型中田分类号:V211.3文献标识码:A为了提高运输机的起飞着陆性能,现代运输机普遍采用了多段翼型技术,对多段翼型的高升力性能进行分析便成为气动设计的一个重要课题.由于目前国内求解多段翼型的已有工程方法是面元法加附面层惨正,精度不高且计算结果不稳定因此,设计部门对多段翼型的气动设计还主要依靠风洞试验.由于求解多段翼型的有粘流动比单段翼型复杂得多,一般用于求解单段翼型粘性流动的
3、有粘/无粘迭代方法不能很好地应用到多段翼型上i又由于多段翼型结构罔格生成的复杂性限制了用NS方程求解多段翼型i因此,急需新的高精度气动分析方法来代替目前国内外求解多段翼型的全速势方程加附面层惨正的方法.本文根据文献[1]提供的方法,采用欧拉方程与可压缩附面层方程耦台求解,形成了一种求解多段翼型的工程实用的数值计算方法.1方法描述采用定常二维欧拉方程作为基本方程.对于封闭的控制体,积分形式的定常二维欧拉方程为质量方程f阳?山一.(1)动量方程P(q?;)q+户一o(2)可压编附方程参影j]c文章编号
4、:1000—2758(2000)03—0438-04能量方程中阳?^,ds一0(3)式中,为单位外法向矢量,为驻点焙.在流线网格中使用有限体积法对欧拉方程求解.由于单元的两侧边界为流线,因此没有质量通量穿过该两边,只有式(2)中的压力项有贡献这表明,只在另外两侧边界上需定义密度和速度.速度的方向必须与当地的几何形状相一致,即和单元的两个侧边是流线的规定相一致采用全局性牛顿方法[解算欧拉方程.用守恒型差分格式,对超音速区使用了人工粘性,以便正确捕捉激渡.考虑粘性影响的办法是,认为附面层和尾流把无粘流
5、动从物体表面推开,推开的量等于"位移厚度".本法采用积分形式的可压缩附面层方程,采用着名的卡门动量积分关系式百d8一吾字一[等+一]妄丽du,2c)百一丁一l可十一;j丽J式中的动量厚度口,位移厚度'和表面摩擦系数c,定义如下一肌一]老却(5a占Cff[一收藕日期:1998—12—03作者筒介借l惦(1972一),男,上海市人.西北工业大学博士生,主要从事计算及设计空气动力学的研究.(sb)(5c)第3期詹特等:多殷翼型的一种分析方法动能积分关系式为d面O"=zc.一[等+s一]繁(6)式中的能量
6、厚度',密度厚度d～和耗散系数C.定义如下.一[一]却一肛1一P-l--ud~(7b)d一一Jl一J(7bc.一士嵩d(7c)再定义三个形状参数H=等H'=鲁H～一等㈣式(4)和式(6)可改写成吾鬈=吾譬一(H+2一)du,(9)∈dH'2C0Cf矿F一百万一百一[等+l—Hi~du,(10)式(10)称为形状参数方程.在此基础上构造层流和紊流封闭关系式,附面层转捩准则是基于臭尔一桑耿弗尔德方程计算的e准则嘲.把附面层求解耦合到求解欧拉方程的全局性牛顿解决中.在选代求解时,按附面层位移厚度逐敬惨改
7、流线形状,即惨改流线网格.最终得出翼顶压力分布,升,阻力系数,以及与之相应的流场流线型网格.2算例及结果分析首先采用RAE一2822翼型作为算例,验证了本文方法的可行性和精度.计算状态为来流马赫敷为0.73雷诺敷为6.5×10'图l蛤出了本文方法计算结果,NS方程计算结果和风洞试验结果的压力分布比较.表1比较了三种结果的升力系数和阻力系数.由图1和表l可见,本文方法的计算结果与风洞试验结果吻合较好}与采用NS方程的计算结果也十分接近,且对阻力系数的计算更准确,表明本文的方法能较好的模拟真实的牯性流
8、动,是一种高精度的数值分析方法.在本算倒中,采用本文方法的计算所用时间是NS方程计算所用时间的近十分之一,这表明本文方法具有收敛快的特点,是一种工程实用的气动分析方法.X/C图1压力分布寰1培果比较本文采用文献[7]提供的GAW—l翼型与30弦长襟翼的两段翼型构型作为算例.计算状态为来流迎角为12.8.襟翼偏角为20.0.来流马赫数为0.13雷诺数为2+2×10'图2给出了压力分布的计算结果与风洞试验结果的比较,两者基本吻合.X/C图2压力分布?440?西北工业大学第18卷由于本文
