变来流下翼型动态失速的协同射流控制数值模拟.pdf
《变来流下翼型动态失速的协同射流控制数值模拟.pdf》由会员分享,可在线阅读,更多相关《变来流下翼型动态失速的协同射流控制数值模拟.pdf(11页珍藏版)》请在咨信网上搜索。
1、文章编号:0258-1825(2023)09-0059-11变来流下翼型动态失速的协同射流控制数值模拟贾天昊,高超*,许和勇,徐泽阳(西北工业大学翼型、叶栅空气动力学国家级重点实验室,西安710072)摘要:针对直升机前飞时旋翼在变来流下出现动态失速的问题,发展了基于协同射流的翼型动态失速控制方法。选取 NACA0012 翼型为研究对象,基于转捩 SST 湍流模型求解非定常雷诺平均 Navier-Stokes 方程,开展不同参数下协同射流控制翼型动态失速的数值模拟。研究结果表明,协同射流能够有效抑制变来流条件下的翼型动态失速。在变来流下,射流流道对翼型原始气动特性产生不利影响,功率系数的增长速
2、度快于射流动量系数的增加,协同射流存在具有较好控制效果的最佳工作区间。协同射流通过与主流掺混来加速涡系演化,以抑制动态失速,通过增强弦向气流的动能以克服逆压梯度,从而抑制流动分离和促进流动再附着。在马赫数 0.283、减缩频率 0.151、前进比 0.25 的条件下,协同射流使翼型升力提高、阻力下降、负俯仰力矩峰降低、流动再附着提前,翼型气动特性得到明显改善。关键词:动态失速;非定常;协同射流;变来流;前进比中图分类号:TP69;O358;TJ011.+5文献标识码:Adoi:10.7638/kqdlxxb-2022.0156Numerical simulation on dynamic st
3、all control of airfoil based onco-flow jet under variable free streamJIATianhao,GAOChao*,XUHeyong,XUZeyang(National Key Laboratory of Science and Technology on Aerodynamic Design and Research,Northwestern Polytechnical University,Xian710072,China)Abstract:Aimingattheproblemofdynamicstallofhelicopter
4、rotorsundervariablefreestreamduringforwardflight,adynamicstallcontrolmethodofairfoilsbasedontheco-flowjet(CFJ)isdeveloped.TakingNACA0012airfoilastheresearchobject,theunsteadyReynolds-averagedNavier-Stokesequationsaresolvedbased on the transition SST turbulence model,and numerical simulations of the
5、dynamic stall of airfoilscontrolledbyCFJunderdifferentparametersarecarriedout.TheresultsshowthatCFJcaneffectivelysuppressthedynamicstallofairfoils.Undertheconditionofvariablefreestream,thejetchannelhasanegativeimpactontheoriginalaerodynamiccharacteristicsoftheairfoil,andthepowercoefficientincreasesm
6、uchfasterthanthejetmomentumcoefficient.ThereisanoptimalworkingconditionforCFJtoachievethebestcontroleffect.CFJacceleratestheevolutionofmulti-vortexbymixingwiththemainstreamtosuppressthedynamicstall,andenhancesthekineticenergyofthechordwiseflowtoovercometheadversepressuregradienttosuppressflowseparat
7、ionandpromoteflowreattachment.UndertheconditionofMachnumber0.283,reducingfrequency0.151andadvanceratio0.25,CFJcanincreasetheliftoftheairfoil,reducethedragandthenegativepitchingmomentpeak,and advance the reattachment of the flow.The overall aerodynamic characteristics of the airfoil areobviouslyimpro
8、vedbytheCFJcontrol.Keywords:dynamicstall;unsteadiness;co-flowjet(CFJ);variablefreestream;advanceratio收稿日期:2022-09-22;修订日期:2022-10-10;录用日期:2022-11-02;网络出版时间:2022-11-14基金项目:国家自然科学基金(12172299,11972306)作者简介:贾天昊(1998),男,四川绵阳人,硕士研究生,研究方向:动态失速及相关流动控制研究.E-mail:通信作者:高超*(1960),男,陕西榆林人,教授,研究方向:复杂流动控制及应用,跨声速空气动
9、力学.E-mail:引用格式:贾天昊,高超,许和勇,等.变来流下翼型动态失速的协同射流控制数值模拟J.空气动力学学报,2023,41(9):5969.JIATH,GAOC,XUHY,etal.Numericalsimulationondynamicstallcontrolofairfoilbasedonco-flowjetundervariablefreestreamJ.ActaAerodynamicaSinica,2023,41(9):5969(inChinese).doi:10.7638/kqdlxxb-2022.0156第41卷第9期空气动力学学报Vol.41,No.92023年9月AC
10、TA AERODYNAMICA SINICASept.,2023 0 引言动态失速是以形成动态失速涡(dynamicstallvortex,DSV)和发生大尺度流动分离为基本特征的复杂流动现象1。直升机旋翼在复杂的跨声速流场中进行周期变距运动,其后行桨叶在大攻角下易发生动态失速。动态失速发生时,气动力出现明显的非线性迟滞,出现升力损失、阻力突增和气动中心位置波动。研究翼型动态失速问题是提高其气动特性的关键途径之一。近年来,针对翼型俯仰振荡引起的动态失速开展了大量实验和数值模拟。McCroskey2选取多种翼型对定常来流下的动态失速问题开展了实验研究。Benton 等3通过大涡模拟方法研究了NA
11、CA0012翼型在高雷诺数下的动态失速问题。目前大部分的翼型动态失速研究都是在定常来流下开展,不符合前飞状态时旋翼的相对来流马赫数周期变化的特征。因此,基于变来流的翼型动态失速研究更能反映其真实的气动特性4-5。Hird 等6使用变来流跨声速风洞研究了 SSC-A09 翼型在变来流下俯仰振荡的气动特性,发现随时间变化的来流速度对非定常气动载荷具有较大影响。对于动态失速现象产生的不利影响,一般采用流动控制技术来改善翼型的气动性能。流动控制技术分为被动流动控制技术和主动流动控制技术。典型的被动流动控制技术有涡流发生器7-8、后缘襟翼9和格尼襟翼10等。这些控制装置能够有效地提高设计点的气动性能,但
12、是难免会影响非设计点的气动性能。主动流动控制技术是通过向流场中注入能量以改善翼型的气动特性,例如边界层吸气可以抑制前缘分离泡11-12、等离子体激励器能够电离局部空气以加速边界层流动13-15、合成射流与合成双射流激励器可以产生涡流和等离子体射流作用于主流16-18。协同射流(co-flowjet,CFJ)是主动控制技术新发展的基于吹/吸气的控制方法19-21。CFJ 翼型的实验和数值模拟表明其具有较高的能量输出水平,能够显著地提高升力并减小阻力22-23。CFJ 的实验最初是依靠外部气源实现的,这种方法虽然方便但不符合工程应用的要求。Yang 和 Zha 等24改进了内置气源的设计并进行了实
13、验,不仅验证了该设计的有效性,还发现 CFJ翼型具有高升力的特性和增强流动抗逆压梯度的特点。该技术基于内置气泵,具有零质量流的优势,加上控制效果显著,被认为具有较大的发展潜力。协同射流方法已被应用于翼型动态失速的控制中,但相关研究都是基于定常来流条件25。因此,在变来流条件下开展协同射流控制翼型动态失速的数值模拟研究,对于改善旋翼动态失速问题具有现实意义。针对上述情况,本文开展了变来流下翼型动态失速的协同射流控制研究,建立了基于风扇(简称为FAN)边界的协同射流数值模拟方法,研究了变来流下翼型动态失速的气动特性,对比了 CFJ 翼型的射流流道对原始翼型气动特性的影响,评估了协同射流方法的控制效
14、率,分析了不同参数对 CFJ 控制翼型动态失速的影响及控制机理。1 计算方法 1.1 CFJ 翼型设计选择 NACA0012 翼型为基准翼型,并在此基础上进行外形修改,得到的协同射流翼型命名为 CFJ0012。协同射流方法工作原理如图 1 所示。InjectionFAN(PUMP)Suction图 1 CFJ0012 翼型工作原理Fig.1 Working principle of the CFJ0012 airfoil翼型弦长为 c,吹气口距离前缘 8.3%c、高度为0.91%c,吸气口距离前缘 83.3%c、高度为 1.94%c。为了使射流更好地附着于翼型表面,在吹/吸气口之间将翼型上表面
15、设计了 0.6%c 的微量下沉,形成了上表面的射流流道。此外,使用无厚度的 FAN 边界模拟协同射流时,翼型内部具有气流通道的设计。FAN 边界是用于模拟内置气泵,高度为 4%c。该边界条件能使气流获得设定的压升 P 以形成 CFJ,使之更贴合实际的工程应用。协同射流开启时,内置气泵输送气流,从前缘吹气口将加速后的气流吹出。射流与主流发生掺混作用,加快边界层的流动。在后缘吸气口处等质量流的空气被吸入,自动实现了零质量流条件。1.2 数值模拟方法Ret数值模拟采用 Fluent 软件,基于有限体积法,选择速度压力耦合算法。远场为 Riemann 不变量的无反射边界条件,壁面为无滑移壁面边界条件。
16、流场求解的控制方程为非定常雷诺平均 Navier-Stokes 方程26。采用转捩 SST 湍流模型,该湍流模型耦合了k-SST 湍流模型和转捩模型,可以对流场不同区域的流动特征进行正确预测。k-SST 模型既有k-模型预测边界层内层流动的优势,又有 k-模型模拟自由剪切层流动的优势,还引入剪切应力输运方程以改善对逆压梯度流动的预测性能。Menter 等27-2860空气动力学学报第41卷研究发现,转捩 SST 湍流模型预测跨声速流动的准确性较高,不仅能较好地预测逆压梯度流动,还能较好地预测流动分离区域。计算网格采用 O 型拓扑的结构网格,对翼型表面、射流流道和 FAN 边界处的网格进行加密处
17、理。翼型俯仰振荡时网格整体刚性运动,周期性旋转中心为 0.25c 处。翼型俯仰振荡可描述为简谐运动,其瞬时攻角表示为:=0+msin(t)(1)式中:0为平均攻角;m为振荡幅度;为翼型俯仰角频率。减缩频率 k 描述来流的非定常效应,k 值与非定常特性正相关。减缩频率表示为:k=c2V(2)假设翼型对应桨叶径向 70.75%处的剖面,桨尖马赫数为 0.4,基准马赫数为 0.283。来流马赫数与桨尖马赫数之比为前进比,定常来流下=0,变来流下 选为 0.25。翼型相对来流马赫数5表示为:Ma=0.283+0.4sin(t+)(3)其中 为来流马赫数与瞬时攻角之间的相位差。直升机前飞状态下前行桨叶流
18、速大而后行桨叶流速小,因此设定=180以减小升力的不平衡。不同 P 形成不同强弱的射流。射流强弱程度可由射流动量系数衡量,其定义如下:C=mVjet12V2S(4)m式中:为质量流量;S 为参考面积。主动流动控制的能耗直接影响技术理论向工程应用的转化。气泵输送气流以克服总压损失,功率与质量流量和总压比相关。功率系数21的定义为:Pc=mcpT2(P1P2)1112V3S(5)式中:cp为定压比热容;T2为吸口总温;为气泵的效率;P1为喷口总压;P2为吸口总压;为比热比。2 计算方法验证 2.1 计算方法验证选择转捩 SST、k-SST、k-三种湍流模型对NACA0012 翼型俯仰运动进行数值模
19、拟,与实验结果29的对比如图 2 所示。实验模型的弦长为 0.601m,实验状态参数为:0=14.91,m=9.88,k=0.151,Ma=0.283,Re=3.45106。网格单元数量为 10 万,时间步长为 0.13ms。0510/()(a)升力系数(b)阻力系数(c)俯仰力矩系数1520250510/()15202500.200.20.40.60.81.21.00.80.60.20.4CD00.22.5ExpTransition SSTk-SSTk-ExpTransition SSTk-SSTk-ExpTransition SSTk-SSTk-2.01.01.5CL0.50510/()1
20、52025Cm图 2 不同湍流模型预测的气动力迟滞回线对比Fig.2 Comparison of the aerodynamic hysteresis loops predictedby different turbulence models-Ret如图 2 所示,不同湍流模型的预测的气动力与实验值大致吻合,但在翼型下俯阶段差别较大。翼型处于深失速状态时,k-湍流模型的计算值与实验值存在较大的差距,对流动分离后的演化过程预测不准确。计算结果不合理的原因可能是 k-湍流模型对自由来流中湍动能的变化过度敏感。转捩 SST 湍流模型依靠模型对转捩过程的正确捕捉,相比 k-SST 湍流模型更接近实验值
21、。转捩 SST 湍流模型对第9期贾天昊等:变来流下翼型动态失速的协同射流控制数值模拟61该流动趋势预测的准确性表明,本文基于转捩 SST湍流模型计算得到气动力结果是可靠的,开展数值模拟研究是可行的。为验证计算方法在变来流条件下开展数值模拟的有效性,选择王清等5的计算结果与 SSC-A09 翼型的实验结果进行对比,如图 3 所示。数值模拟的状态参数为:来流马赫数 Ma=0.4+0.08cos(t13.3),瞬时攻角=8.513cos(t),减缩频率为 0.05。图 3中升力系数的变化趋势与王清等5的计算结果显示出较好的一致性,但在大攻角下计算值与实验值存在一定的偏差。在实验中变来流风洞的气流速度
22、难以精准控制,会导致来流速度的波形发生畸变5,还会导致气流的湍流度增大。而数值模拟中来流速度的变化规律受到准确的控制,因此可能导致计算值与实验值存在偏差。但气动力变化趋势基本吻合,表明了计算方法具有可靠性与有效性。5CL2.01.51.00.500.50510152025/()ExpCFD by Wang5Present CFD图 3 变来流下翼型升力系数曲线对比Fig.3 Comparison of airfoil lift coefficient curvesunder variable free stream 2.2 CFJ 验证本小节采用具有实验结果的 CFJ6415 翼型进行CFJ
23、数值模拟方法验证。DANO 等30在 NACA6415翼型的基础上设计了 CFJ6415 翼型,并进行了实验验证。实验模型的弦长 30.48mm,展长 59.06mm。实验状态参数为:来流马赫数 0.03,温度 288.15K,标准大气压力 101.325kPa,质量流率 0.06kg/s。实验的协同射流系统由独立控制的高压气源和低压真空气罐组成。为验证 FAN 边界在 CFJ 数值模拟中的可靠性与准确性,本文在相同状态参数下开展了数值模拟。计算结果与实验值的对比如图 4 所示,15攻角下的流场如图 5 所示。可见,计算值与实验值的升力系数具有较好的一致性。攻角大于 15时计算的阻力系数偏低,
24、功率系数略低于实验值,但变化趋势与实验值吻合。结果表明,基于 FAN 边界的数值模拟方法能够较好地模拟 CFJ 翼型的气动特性,为 CFJ 的流场研究提供了较为可靠的方法。0510152025/()43210Exp(CL)Exp(CD)Exp(PC)CFD(CL)CFD(CD)CFD(PC)CL,CD,PC图 4 升力、阻力和功率系数的数值结果与实验值的对比Fig.4 Comparison of the lift,drag and power coefficientsbetween the numerical results and the experimental data0 Velocit
25、y122436图 5 15攻角下的流场Fig.5 Flow field at the angle of attack =15 3 结果与分析本节基于以上 CFJ 模拟方法,对 CFJ0012 翼型在变来流条件下的俯仰运动绕流进行数值模拟,研究射流流道的影响,评估射流功耗特性,分析 CFJ 控制动态失速效果以及流动机理。3.1 CFJ 翼型射流流道的影响气泵关机时的 P 设置为 0,该状态下的翼型命名为 CFJ-off。选取=0 和=0.25,对 CFJ-off 翼型开展数值模拟,并与基准翼型对比。不同方位角 下的气动系数如图 6 所示。选取气流再附着过程中的=30、基准翼型的失速时刻=243和
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 流下 动态 失速 协同 射流 控制 数值 模拟
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。