数值分析带二次流跨音速尾喷管流场.pdfVIP

中国工程热铀理学会 编号：20∞4112 燃烧学学术会议 数值分析带二次流跨音速尾喷管流场 赵坚行周琳 (南京航空航天大学动力工程系，南京·210016) ■要，车文用有限元法数值分析带二次流可调收—扩嚷管辟音遮藏动，素巍桔性采用k_一c双方程纂 丑束描述．并利用a-1．rklf．法建立有限元方霉．求鼻可压■桔眭蠢件藏动N—s方程，奠蕾■酉采用■醯 散．薯过不同曩警、不同工况下尾曩警内帮蠢栩计算．并用实t敷霸盎行对比衰嘲·廿算螬摹音曩，鲁枣可 簟，可利用有限元j缶对曩管超音燕荫瓒盎行鼓值研究，并为_凳砖机尾曩蕾优化骷P和瞢啊摄口pl用敷据． 1引言 航空发动机尾喷管可分为收敛喷譬，收敛一扩散喷管以及引射喷臂三种，其中翦者在麓 警内气流可为亚音速流动，后两种在喷譬内燃气流为跨音逮流动，本文所研究的带二次巍 可调收一扩喷瞥是一种适用于渭扇发动机的可调喷臂。它可利用谭节二次流量和喷管几何 形状来满足发动机工作状态需要．因此燃气在喷臂内的流动状况直接影响发动机特性，由于 超音遮气流流动复杂，故敦值模拟跨膏建气流是囊臂流场计算的关鲁． 有限差分法和有限元法都是目前被广瑟采用两种近似计算方法，前者应用历史较长-后 者是一种新的数值法．用它来解流动同曩时，其优点t一是对流场形状没有什么限倒I二是边 界条件容易处理，随着有限元数值分析法迅速发晨，近年来应用有限元法处理流悻力学问题 逐渐增多，为此本文采甩有限元法来教值分析带二欢流的收一扩喷管流动． 2擂值函数和四边形等参元 2．1捶擅函数 在有限元分析中，一个极为重要内容是进舞具体的单元和确定单元的位咎函散·然后便 可按照标准步翼进行运算．在计算流动问囊时，主要有三角形单元和四边形单元，由于后者 适用于结构较为规则几何形状，而且精度又较高．因此本文采用任意形状和大小的四边形单 元来蚓分计算区域． 在有限元法中，位移函数的描述，可以有两种方式：其一是采用古有若干特定系散的简 单多项式，这些系数继面变换成为相应的节点位移参散．其二可采用插值函数直接来描述， 337 本文采用后者，对于二雏流场，任—变量口．其插值函数可以写戚t ^ (1) 中“，y)膏Nl(x，y)ml+Nl(x，y)tb+…N．(x，y)o、I备Nl“，y)也 式中N“，y)由具体单元和节点来定，变量中用上式表达后，播值函数棘描述了单元中任一 点的变量．对稍黼基单元来说，所选舞的位移函数除了擒足完备和协调条件外，其特定 系数的数目还应和单元的节点自由度数目一致，以便由节赢位移值唯一地确定全部秦靛· ● 2．2四边搭善●元 ● 由于在计算区壤中四边形单元形状和大小可任意，直按对其进行积分较目难，故可对其 进行坐标转换，把任意瞬边形单元变换为正方形单元，坐标变换不奄针对整个区域丽定，而 是分别对每—单元．设(x，y)为整体坐标它适用于求解整个甘算区域而及(}，q)为局部坐标t 它只适用于每一个单元．总体坐标“，y)对局部坐标(E，1)之筒的变换具有如下的关系 ，x—N—(1+N而+NtK+N-)c- (2) bIN*+N觑+Nth+N^ 式中 Nl暑(1+皇)(1+q)／中 N．=(1一毫)(1一棚)／o， N。-(1·专)(1+q)，母 (3) NI高(1+￡)(1一n)／巾， 称为插值西羲琥丑螽敷 如把正方形单元扮造成驭筑性播值函数； 夺置b