摘要:在使用非结构网格求解器WiseCFD-UG软件对NASA第三届高升力会议JSM模型进行数值模拟研究的过程中,发现前缘缝翼与主机翼交接处的几何修型方式对失速特性的模拟结果有较大影响。针对切缝与不切缝这两种不同的几何修型方式,分别生成了同等规模的非结构网格。两套网格得到的计算结果表明:大攻角下不切缝的几何修型方式得到的压力分布及升力系数与实验值吻合的更好,能够更好的模拟低速构型的失速特性。
加入收藏
CFD技术发展到今天,已经被大量计算和实验证明是足够可靠的,可以较好地模拟典型运输机高速巡航构型的气动特性。但是对于复杂的高升力构型,CFD计算仍然存在着很多的不确定性。高升力构型的流动复杂,包括尾迹区和边界层流动的混合、分离流、翼尖涡、层流到湍流的转捩等等,这些因素共同作用无疑增加了解决问题的难度[1]。借鉴阻力预测系列会议(DPW)的成功经验,美国航空航天学会组织了高升力系列会议(HiliftPW)。通过发布基准研究构型、提供基础计算网格、公布详实的风洞试验结果等组织方式,吸引了世界范围内相关研究机构和CFD工作者的广泛参与[2]。前两届的AIAA高升力会议表明,想要精确的模拟高升力构型的最大升力系数CLmax、俯仰力矩系数以及升阻比仍然存在很多挑战。
AIAA第三届高升力会议(HiliftPW-III)于2017年6月3日在美国科罗拉多州丹佛市举行。会议旨在评估现有的CFD技术(包括网格、数值格式、湍流模型以及高性能计算等)对模拟大展弦比飞机起飞着陆特性方面的能力[2]。会议选取了两个高升力模型:一个是日本太空探索机构JAXA设计的标准模型JSM,另外一个是AIAA阻力预测会议上的通用模型CRM的高升力版本——HLCRM模型。目前HLCRM模型的实验还没有完成,因此选用JSM模型作为研究对象,对高升力构型进行失速特性的模拟研究。
1、JSM模型描述
JSM模型是一个通用航空器模型,代表着可以容纳100位乘客的经济型客机。模型由日本JAXA中心设计,是一个翼身融合体的高升力构型,包括襟翼、缝翼、支撑、吊挂以及滑轨整流罩等众多机构[3]。实验使用的模型是半模,在对称面上有一段60mm长的延伸用作支架,消除风洞侧壁对流动的影响。计算时不考虑这一段支架,使用自由来流计算。JSM模型机翼前缘后掠角33度,前缘缝翼有25度下偏角,后缘襟翼有30度的下偏角,模型展弦比9.42,梢根比0.3。模型平均气动弦长529.2mm,模型半展长2300mm,机身长度4900mm,半模参考面积1123300mm2,力矩中心参考点坐标(2375.7mm,0.0mm,0.0mm)。JSM模型风洞实验于2009年在日本JAXA-LWT1低速风洞完成,风洞实验数据主要包括气动力和力矩、表面压力分布、表面油流、边界层转捩测量、模型变形测量、非定常压力和速度场、气动噪声测量等[4]。
因为JSM模型非常复杂,流动特征难以捕捉,对CFD模拟提出了严峻挑战。从实验油流图来看,JSM模型的基本流动特征包括:(1)后缘襟翼受到主机翼下表面滑轨整流罩影响,导致出现小的分离区;(2)大攻角时翼尖处流动发生分离,从流动分离发生的位置判断可能是因为受到缝翼支撑的影响,并且随着攻角继续增大翼尖的分离也逐渐扩大;(3)失速是由翼根处流动大面积分离引起的。具体如图1所示。
图1JSM模型支撑及整流罩对流动的影响
2、计算模型与网格
网格生成是CFD模拟中非常重要的一个环节。对于模拟复杂流动,同一套模型生成的不同网格很可能会得到不同的模拟结果[5]。而几何模型的处理和简化又是网格生成其中最为关键的一步。CAD几何模型被导入网格生成软件里后,一般要经过修剪、整理与简化才能用于网格生成。对原始几何模型中一些不利于网格生成的部件需要进行合理的修型,修型时应尽可能的保留原有的几何特征,避免因几何模型的简化导致数值模拟结果产生重大误差。
JSM模型的前缘缝翼与主机翼交接处存在剪刀角,这给模型附面层的网格生成带来了困难。针对这种情况,一般需要对几何外形进行简单的处理以便于网格生成。其中,最简单的办法就是直接将缝翼上与主机翼相接处的一小段用平面切掉,简称切缝。这种切缝的方法简单快捷,并且经验表明多数情况下对工程计算结果影响不大。另外一种方法,不做切缝处理,只在缝翼与主机翼交接处做轻微的修型,这也是HiliftPW-III会议上组委会所采用的修型方式。如果交接处夹角是钝角就用一小段平面代替,如果交接处夹角是锐角则用半径很小的球形代替。这种方法对原始模型的改动很小,可以最大限度保留原始几何外形,尽可能的降低几何修型对计算结果的影响。
因为高升力构型复杂,生成结构网格需要的人工工作量太大,所以HiliftPW系列会议在网格生成方面有一个显著特点就是非结构网格占据主导地位。针对两种几何修型方式,使用Pointwise软件生成两套JSM模型非结构网格。其中,Cfix代表不切缝修型方式生成的网格,Pfix代表用平面切缝生成的网格。生成两套网格时严格按照HiliftPW-III组委会制定的网格生成规范进行。
生成网格时附面层第一层网格间距取3.63×10E-6m,壁面y+=2/3,体网格增长比设置为1.16,附面层网格粘性网格层数60层。另外在前缘缝翼的尾迹区进行局部空间网格加密,最终体网格总点数控制在3500万,体网格单元总数约为7600万。
3、计算结果
对切缝(Pfix)与不切缝(Cfix)得到的两套网格,使用WiseCFD-UG软件进行计算。计算状态:来流马赫数Ma=0.172(60m/s),模型总雷诺数Re=193万,来流温度T=306.55K。计算用湍流模型选择SA一方程模型,空间离散格式使用Roe格式。计算攻角选择-10.06°、-5.95°、-4.89°、-2.81°、0.0°、4.36°、8.45°、10.47°、12.52°、14.54°、16.57°、18.58°、19.50°、20.09°、20.58°、21.57°,共计16个计算状态。所有状态均从自由流开始计算,计算15000步,16.57度攻角以下残差收敛4~5个量级,18.58度攻角以上残差收敛3~4个量级,大攻角下计算的力系数取最后1000步的平均值。
两种修型方式得到的力系数计算结果如图2所示。从图中可以看出,在12.52度攻角之前两种修型处理方式的升力系数计算结果非常接近,并且都与实验值吻合的很好。但是从14.54度攻角开始,切缝处理方式的升力系数计算结果明显低于实验值,并且在18.58度攻角之后开始失速,最大升力系数2.671,比实验值2.775要低许多。而不切缝处理方式的升力系数计算结果与实验值吻合的很好,失速攻角在20.09度攻角附近,最大升力系数为2.789,与实验值非常接近。升阻比的计算结果也显示出,大攻角下不切缝处理方式的计算结果与实验值更为接近。
图2升力系数及升阻比结果对比
比较14.54度攻角之后的压力分布计算结果,发现两种修型方式计算得到的压力分布有较大差别。图3展示了18.58度攻角下机翼上靠近翼尖出H-H剖面面压力分布同实验值的对比。从图中可以看出,不切缝处理方式得到的机翼各截面压力分布和实验值吻合的更好,缝翼和主翼上的前缘吸力峰较切缝方式有明显增加。另外,在主机翼外翼H-H截面上,不切缝与切缝的计算结果有很大差别,这是因为不切缝处理方式计算的翼尖分离区比切缝处理方式的分离区要小,更符合真实流动的结果。其他几个截面也存在类似情况,不切缝方式的压力分布计算结果均要好于切缝方式,受篇幅所限不再一一给出。
图3机翼外翼段H-H截面压力分布对比
4、结论
本文在对AIAA第三届高升力会议提供的JSM模型开展数值模拟研究的过程中发现,前缘缝翼与主机翼交接处的几何修型会对计算结果产生较大影响。对使用切缝方式与不切缝方式处理的JSM模型分别生成了两套同等规模的非结构网格,并使用WiseCFD-UG软件对JSM模型Case2c状态开展计算研究。将两套网格得到的计算结果与实验值做了对比分析,可以得到如下结论:(1)几何修型对升力系数计算结果有很大影响,处理几何时应在尽可能保留原始模型的各种关键特征的前提下,对几何模型做一些合理的简化以便于网格生成。(2)JSM模型不切缝处理方式大攻角下升力系数的计算结果和实验值吻合的更好,失速攻角及最大升力系数同实验值更接近,机翼各截面上的压力分布也与实验值更为一致。两种几何修型方式对阻力系数的计算影响不大,阻力系数与俯仰力矩系数的模拟精度有待进一步提高,需要后续开展更深入的研究工作。
姚冰,淮洋,梁益华.几何修型对JSM模型失速特性的影响研究[J].电子测试,2020(08):42-43+33.
分享:
专业导论课程是各个工科专业都要开设的一门必修课程,旨在介绍本专业的研究对象、发展方向、关键技术及课程体系,通过该门课程的教学让学生对专业有一个总体上的了解,并为后续的专业课学习提供系统性指导,同时也需要培养学生树立投身于祖国工业和国防建设的伟大理想,在专业导论课的教学方面,多位教师学者都提出了很多教学改革思路[1,2,3,4,5]。
2024-07-17在网络信息发达、繁杂的时代背景下,青年大学生极易接触到不良的思潮和价值观,存在急功近利的价值取向,缺乏家国情怀,集体责任感以及坚定的社会主义人生理想信念。大学生在校期间的绝大多数时间是在非思政类课程学习中度过。习近平总书记在全国高校思想政治工作会议上强调指出“要用好课堂教学这个主渠道。
2020-12-15腐蚀是普遍存在于国民经济和生产生活中的重要破坏因素,学者们针对材料腐蚀现象进行了众多研究,其中,静态环境下的腐蚀研究已较为丰富。在流动体系中,流体对材料表面及腐蚀产物膜的力学冲刷作用、对腐蚀性离子的传质促进作用以及对界面反应过程的影响使得流动环境下的腐蚀问题更加严重和复杂。
2020-12-03在核反应堆中,控制棒驱动机构是反应堆控制和保护系统中重要的伺服机构,控制棒驱动机构是驱动控制棒组件运动、实现准确调节反应堆内反应性变化的重要部件,控制棒驱动机构具有实现反应堆启动、提升功率、保持功率、负荷跟踪、正常停堆和紧急事故停堆等重要的安全功能。
2020-11-27腐蚀是普遍存在于国民经济和生产生活中的重要破坏因素,学者们针对材料腐蚀现象进行了众多研究,其中,静态环境下的腐蚀研究已较为丰富。在流动体系中,流体对材料表面及腐蚀产物膜的力学冲刷作用、对腐蚀性离子的传质促进作用以及对界面反应过程的影响使得流动环境下的腐蚀问题更加严重和复杂。
2020-11-17房间内污染物来源包括:煤气泄漏、一氧化碳泄漏、室内装修后残留的甲醛等污染物。其中煤气和天然气、CO等污染物通过管道或者开关未关闭泄露至空间中,而室内装修残留物则主要为家具或者地板等使用的材料超出规定指标而不断残留在房间中。本文主要采用CFD的方法进行结构合理性的验证。
2020-11-17为研究一种高减压比组合式空气减压阀的流场特征,基于计算流体动力学(CFD)方法,建立减压阀数值计算模型,针对两级阀体特征参数组合的典型工况进行仿真计算,研究减压阀流场参数分布,分析柱阀间隙和引气位置对流动状态和减压比影响。计算结果表明,组合式减压阀流场状态复杂,第一级柱阀间隙对于流动损失和局部参数分布有重要影响,两级阀分别通过总压损失和流速变化实现减压过程,通过合理的特征参数组合可获得较大减压比及其调节范围。
2020-11-16气体在液体中的运动广泛存在于自然界和工程领域中,如发动机水下排气、水下爆炸引起的气泡运动、石油开采.其中的动力学过程在一个多世纪以来一直受到人们的关注[1,2].数值模拟是研究气体在液体中运动的一种有效方法[3,4].由于液体与气体之间的密度比一般远大于重液体与轻液体之间的密度比,所以气液相界面的捕捉一直是数值模拟的重点与难点.
2020-08-10理论和实践证明,环量控制[1]方法可以通过产生射流在后缘形成Coanda效应[2]来改变环量,从而产生高升力[3]。国内外一些学者针对环量控制参数的影响都进行了研究:Englar[3]、Liu等[4]经研究得出在小攻角、小动量系数下环量控制可产生较大升力的结论;张艳华等[5]通过计算流体力学数值仿真方法(CFD)研究了动量系数和攻角对NASA超临界翼型的升阻特性和附面层分离特性的影响规律。
2020-08-10流体力学是研究流体在静止和运动状态下的运动规律及流体与固壁相互作用力的一门学科,在机械、能源、汽车、化工、航空、航天等诸多领域有着极广泛的应用[1]。该课程是力学、机械等专业本科生的专业基础课程。由于其涉及数学、物理、力学等内容,概念多,公式推导繁琐,理论性强,实际应用广,对学生综合利用知识的要求较高,是普遍公认的教师难教、学生难学的课程之一[2]。
2020-08-10人气:8097
人气:7758
人气:6590
人气:6027
人气:4809
我要评论
期刊名称:力学学报
期刊人气:2858
主管单位:中国科学院
主办单位:中国科学院力学研究所
出版地方:北京
专业分类:科学
国际刊号:0459-1879
国内刊号:11-2062/O3
邮发代号: 2-814
创刊时间:1957年
发行周期:双月刊
期刊开本:大16开
见刊时间:一年半以上
影响因子:0.963
影响因子:1.730
影响因子:0.542
影响因子:0.744
影响因子:1.028
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!