跨声速压气机非定常流动及损失的数值研究
高负荷跨音速轴流压气机的叶型优化设计
能。
压气机的全三维数值优化技 术集三维 流场求解 、 优化 技
术 为 一 体 , 过 对 原 型 机 进 行 参 数 化 而 扰 动 几 何 构 型 , 从 通 以
中寻优。 目前较常用 的方法是结合神 经网络 ( 响应面模 型的
一
种) 与遗传算法 的优 化方法。众多国外学者 已在该领域 进
第2 卷 第7 9 期
文 章 编 号 :06 94 (0 2 0 - 0 5 0 1 0 - 3 8 2 1 )7 07 - 5
计
算
机
仿
真
21年7 02 月
荷 跨 音 速 轴 流压 气 机 的 叶型 优 化 设 计
李相 君 , 武利 , 皓光 , 楚 张 米 攀
( 西北工业大学 动力与 能源学 院, 陕西 西安 7 07 102)
定 范 围 内扰 动 , 可 以 得 到 一 系 列 新 的 叶 型 , 就 以此 构 建 叶 型
面处则使用 周向平 均法进 行数 值传 递。计算 区域 包括进 口
导 叶、 转子及静子 排 。另外 为提 高计 算效 率 , 采用 了多重 网
格 法 、 部 时 间 步 长 和 隐 式 残 差 光 顺 三 种 措 施 加 速 计 算 收 局
一
7 — 5
表 1 计 算 网格 数 量分 布
压气机跨音速叶栅叶型的数值模拟
求解器采用耦合求解器隐式方案 , 根据湍流模型
收 稿 日期 :0 8— 6—0 20 0 3 作 者 简介 : 张广 (9 3一 , , 龙 江 哈 尔滨 人 , 读 硕 士 , 18 )男 黑 在 主要 研
对弯度. 最大相对厚度 e不变的基础上 , 厂 _ 、 将最大 弯度 位置 、 最大 厚度 位 置 向后 缘移 动 , 过 这 通
摘
要: 以抑制叶型吸力面激波强度 、 附面层分离 , 降低叶型损失为 目标 , 对某型 压气机跨 音速 叶
栅 叶型进行气 动优化 。运用数值模拟 的方法对 比研究优化后 的叶栅 同原型叶栅在相 同工况下 的 气动性能 , 概述了气动性能产生变化 的机 理。与优化 前对 比发现 , 优化后 的叶 型具 有低损 失 , T
采用 迎风 格式 的有 限体积 法对 二维 定 常 Ⅳ 一 s方程 进行 数 值 求 解 , 维 流 场 计 算 区域 采 用 H 二
型 网格 , 网格 数 为 9 ( 距 方 向 ) ×3 0 流 动 方 0栅 8(
2 气 动 优 化 思 路
借 鉴高 亚音 速叶栅 以及 超音 速 叶栅 叶型特点 对 该跨 音速 叶 型进 行 优 化设 计 , 索 验 证 叶 型最 探
马赫数 。叶栅 前后 延伸 部分 均采用 周期 性边 界条 件 。叶盆和 叶背均 采用 绝热 的无 滑移壁 面条 件 。
研究 方 向。本 文通 过 求 解 可 压缩 雷 诺 平 均 N—s
方程 , 数值 模拟 了压 气机 叶栅跨 音速 流 动情况 。
1 计 算方 法 与边 界 条 件
作 范 围大 的特 点 。
关键 词 : 跨音速压气机 ; 叶型 ; 气动优化 ; 数值计算
跨声速轴流压气机非设计点特性计算
t u ba nn t e p n w s a rd n mi p r mee s n c a a t r t u v s t d s n n o - e i p i t. T e h s o ti i g h s a — ie eo y a c a a tr a d h rc ei i c re a e i a d f d sg sc g n o ns h c mp r o i e tr s l n iae h tt e l s n e it n a g e mo e s efc ie fr t e c mp e s r p roma c o ai n w t t s e u t i d c ts t a h o s a d d va i n l d li f t o h o r so e fr n e s h s o e v
Nu ei a i l to fa c ran t p fta s nc a ilfo c mpr so sc rid o tu i te mlne c r aur t d, m rc lsmu ai n o et i y e o r n o i x a— w o l e s r i a re u sng sra i u v t e meho
摘要 : 基于公开发表 的研究成果 , 发展了一种适用于跨声速轴流压气机 的损失和落后 角模型 , 采用流线 曲率法对某
单级跨声速轴流压气机进行了数值模拟 , 到了展 向参数分 布和 全工况下 的性能 曲线 。通 过与实验值 的 比较 , 得 验 证 了该方法和模型在跨声速轴流压气机非设计点性能预测 中的有效性 。
prdito e c in. Ke r y wo ds:t e m ln ur at e;t a on c;de i to a sr a i e c v ur r ns i v a i n ngl ls e; o s
跨声速轴流压气机转子Rotor37周向槽机匣处理的数值研究
跨声速轴流压气机转子Rotor37周向槽机匣处理的数值研究何文博;史磊【摘要】以跨声速压气机转子Rotor37为研究对象,采用商业软件NUMECA数值研究了单槽处理机匣的轴向位置对于压气机性能及内部流场的影响.周向槽处理机匣的宽度为3 mm,深度为10倍叶尖间隙,即3.56 mm,起始位置分别位于轮缘机匣尖部型面的10%、20%、30%、40%、50%相对弦长处.数值计算结果表明:原始光壁压气机转子的失速原因为叶尖泄漏流动引发的低速区对于尖部叶片通道的堵塞,其稳定工作裕度为14.74%.采取的周向槽机匣处理能够改变转子叶尖流动堵塞状况.当机匣处理起始位置位于30%相对弦长时,压气机转子稳定工作裕度的提升量最大,相比原始压气机转子的稳定裕度提高了1.86%.【期刊名称】《科学技术与工程》【年(卷),期】2018(018)028【总页数】6页(P164-169)【关键词】跨声速压气机;周向槽机匣处理;叶尖泄漏流动;稳定工作裕度【作者】何文博;史磊【作者单位】中国民航大学天津市民用航空器适航与维修重点实验室,天津300300;中国民航大学天津市民用航空器适航与维修重点实验室,天津300300【正文语种】中文【中图分类】V231.3近些年来,航空发动机的推重比在升高,压气机的轴向尺寸在不断缩短,级数在逐渐减少,平均级负荷水平一直在提高。
然而负荷水平的提高使得压气机的稳定工作范围减小,因此有必要拓宽压气机的稳定工作边界;其中机匣处理即作为一种有效的扩稳技术手段。
人们对机匣处理的扩稳认识始于20世纪60年代初,此后美国NASA、法国国家、航空航天研究院、英国剑桥大学、日本东京大学等纷纷开展了对机匣处理技术的研究[1]。
进入70年代,大量的实验研究工作侧重于槽类和缝类机匣处理的设计方法、扩稳效果及结构尺寸的优化,Lewis研究中心针对机匣处理就做了大量的工作,得出了机匣处理可以使压气机裕度改进,但同时又使压气机效率降低的结论[2,3]。
吸附式跨声速压气机参数化设计研究
1 设 计 和 分 析 方 法
1 I数 值 方 法 .
罗建枫 一朱俊 强 卢新 根 , , ,
( .中国科 学院工程 热物 理研 究所 ; . 国科 学院研 究 生院 北京 10 8 ) 1 2 中 0 0 0
摘 要 : 文 借 助 S/ 2两 类 流 面 迭 代 和 s 流 面 正反 问 题 混 合 计 算 方 法 , 叶 片 吸 力 面 边 界 层 吸 气 为手 段 , 本 1¥ l 以
详 细研 究 了边界 层抽 吸 的吸气 量 和吸气 位置 对跨声 速压 气机 叶栅 气 动 性 能 的影 响 , 果 表 明 叶 栅来 流 结
Ma 和方 向一定 时 , 吸气 位 置和 吸气 量是相 互 关联 的 关键 参数 , 想的 吸 气 位 置应 该 是 在 靠 近 激 波后 边 理
迭代 计算 , 中无 粘 计 算 部分 使 用 牛顿 法 求 解 欧 拉 其 方程组 , 粘性 影响 由边界层 积分方 程计算 得 出的边 界
中 图分 类 号 : K 7 . T 44 8 1
提高转 子 的转速 或增 加 叶片 的气动 负荷 是提 高 轴 流压气 机级 压 比的 有 效途 径 , 而转 子 转 速 的 提 然 高 除受 到叶 片所 能 承受 最 大 应 力 的 限 制外 , 使 激 还 波本 身 以及激 波 所 诱 发 的 一 系列 损 失 增 加 。因 此 ,
讨 , 文将 从 以上两 方面 展开பைடு நூலகம்详细 研究 。 本
压气机失速信号分析
压气机失速信号分析作者:姜椿阳来源:《科技传播》2016年第11期摘要以压气机失速信号分析为目标,选用MATLAB作为信号分析平台,采用时域分析法、频域分析法进行信号分析,并考虑在得到压气机叶片参数和转数等信息后,采用滤波的方法做更深入的分析。
关键词压气机;失速信息;信号分析中图分类号文献标识码 A 文章编号 1674-6708(2016)164-0229-02在本课题中,由于未给出压气机转子叶片数量、压气机转速等信息,所以,这些信息我们必须通过信号分析得到,此外,我们还需得到失速征兆以及由扰动产生到进入失速所经历的时间,同时,在可能的情况下,我们还应判断出失速的基本类型,包括辨析出是突变失速还是渐进失速、全展向失速还是部分展向失速等,以及失速团的个数。
为了得到以上信息,我选用MATLAB作为信号分析的平台,主要采用了时域分析法、频域分析法进行信号分析,并考虑在得到压气机叶片参数和转数等信息后,采用滤波的方法做更深入的分析。
1 信号分析时域分析:信号1:105800A.dat如图1所示,即为近设计点工况时域图,横坐标按采集点个数均分,由于采样频率恒定为12K,采样时间为5s,所以采集点个数也就表征了时间。
由上图可以看出,近工况状态下,前缘、50%弦长和尾缘的压强信号基本保持稳定,只是在小范围内有脉动,将测得参数按按时间平均,并设前缘静压平均量为P1,50%弦长处P2,尾缘P3,则有P1=2.544,P3=3.013,于是得近设计点工况增压比n1=P1/P3=1.184。
如图2所示,压气机的近失速点实际即为压气机叶片气流分离点位于50%弦长处时,此时,微小扰动即可导致压气机进入失速状态。
按上一数据的处理方法,我们得到近失速点压比为n2=1.156,较近设计点工况略小,这是由于流动的分离导致压气机效率降低造成的/但是当进入失速状态时(近失速状态),我们可以定性的分析出尾缘处静压与前缘处静压在数值上呈现近似的负相关关系,这是由失速团的性质决定的,并且,节流阀开启后,压气机增压绝对值突然上升又剧烈下降,最后,成周期性变化,由此可以初步判定该失速类型为突变型失速。
基于CFD
第8卷㊀第6期2023年11月气体物理PHYSICSOFGASESVol.8㊀No.6Nov.2023㊀㊀DOI:10.19527/j.cnki.2096 ̄1642.1088基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析李泳德ꎬ㊀郭㊀力ꎬ㊀季㊀辰(中国航天空气动力技术研究院ꎬ北京100074)CharacterizationofTransonicAerodynamicDampingofRocketsBasedonCFD/CSDCouplingLIYong ̄deꎬ㊀GUOLiꎬ㊀JIChen(ChinaAcademyofAerospaceAerodynamicsꎬBeijing100074ꎬChina)摘㊀要:随着新型大推力火箭的发展ꎬ弯曲模态频率的不断降低ꎬ以及流动分离和跨声速飞行时产生的激波震荡等因素ꎬ其在跨声速飞行过程中更容易出现非定常振动发散ꎮ文章以某带助推的运载火箭模型为研究对象ꎬ通过数值计算获取火箭强迫振动时的气动阻尼ꎬ并对影响火箭气动阻尼的因素进行了分析ꎮ包括结构节点位置㊁振动振幅大小㊁脉动压力等ꎮ研究表明:助推主要起到增大气动阻尼的作用ꎻ前节点主要影响收缩段的气动阻尼ꎻ振动振幅大小和脉动压力对气动阻尼的影响可忽略不计ꎮ关键词:气动阻尼ꎻ数值计算ꎻ跨声速ꎻ气动弹性ꎻ运载火箭㊀㊀㊀收稿日期:2023 ̄09 ̄25ꎻ修回日期:2023 ̄10 ̄23第一作者简介:李泳德(1995 ̄)㊀男ꎬ工学硕士ꎬ助理工程师ꎬ主要研究方向为气动弹性分析ꎮE ̄mail:562064169@qq.com通信作者简介:季辰(1982 ̄)㊀男ꎬ工学博士ꎬ研究员ꎬ主要研究方向为气动弹性力学ꎮE ̄mail:jichen167@hotmail.com中图分类号:V475.1㊀㊀文献标志码:AAbstract:Withthedevelopmentofnewhigh ̄thrustrocketsꎬthedecreasingfrequencyofthebendingmodesoftherocketꎬaswellasthefactorssuchasflowseparationandshockoscillationsgeneratedduringtransonicflightmakeitmorepronetonon ̄constantvibration.Inthispaperꎬalaunchvehiclemodelwithboostwastakenastheresearchobjectꎬandtheaerody ̄namicdampingoftherocketduringforcedvibrationwasobtainedthroughnumericalcalculation.Thefactorsaffectingtheaerodynamicdampingoftherocketwereanalyzedꎬincludingthepositionofstructuralnodesꎬthemagnitudeofvibrationam ̄plitudeꎬpulsatingpressureandsoon.Thestudyshowsthattheboostmainlyplaystheroleofincreasingaerodynamicdamp ̄ingandthefrontnodemainlyaffectstheaerodynamicdampingofthecontractionsection.Thevibrationamplitudesizeandthepulsatingpressurehaveanegligibleeffectontheaerodynamicdamping.Keywords:aerodynamicdampingꎻnumericalcalculationꎻtransonicꎻaeroelasticityꎻlaunchvehicle引㊀言通常情况下人们认为气动力对火箭的振动起到阻尼作用ꎬ即气动阻尼为正值ꎮ然而随着大推力火箭发展ꎬ火箭的长细比逐渐加大ꎬ导致弯曲刚度越来越小ꎬ同时为了满足有效载荷的外形要求ꎬ火箭头部整流罩尺寸不断加大ꎬ后续箱体的直径却保持不变ꎬ形成了典型的锤头体外形ꎮ国内外大量的火箭研制经验表明[1 ̄9]ꎬ对于此类锤头体外形火箭的气动设计ꎬ必须要进行动态气动载荷与动态气弹稳定性分析ꎬ否则设计的疏忽可能会导致火箭结构出现毁灭性的破坏进而导致发射失败ꎮ目前常用的衡量气弹稳定性的方法是通过风洞试验来获取气动阻尼系数ꎮ早在1963年ꎬ美国国家航空航天局Ames研究中心(NASAAmesRe ̄searchCenter)采用半刚性模型开展试验研究[10]ꎬ获取火箭头部的气动阻尼来评估其稳定性ꎬ但这只能用来模拟火箭弯曲振型前节点之前部分的结构动力学特性ꎮ直到兰利研究中心(NASALangleyResearchCenter)开发了全弹性模型气动阻尼试验气体物理2023年㊀第8卷技术ꎬ其可以模拟整体的结构动力学特性以及气动外形ꎬ并应用于多款运载火箭研制[11 ̄15]ꎮ国内ꎬ中国航天空气动力技术研究院对气动阻尼问题开展过较多的研究[16 ̄20]ꎬ从模型设计方法㊁模型制作工艺㊁试验机构设计和数据处理等诸多方面ꎬ逐步改进实现了从半刚性模型到全弹性模型的过渡ꎬ并在多个型号上得到验证ꎮ然而通过风洞试验研究气动弹性问题ꎬ技术难度大ꎬ试验成本高ꎬ同时几乎不可能开展全尺寸试验ꎮ因此通过数值计算的方法开展相关研究是另一种重要的手段ꎮ刘子强等[21]实现了通过数值计算确定气动阻尼系数的技术和方法ꎬ并与试验结果进行对比ꎬ证实了该方法的可靠性ꎮ冉景洪等[22]通过模态数据结合准定常理论的方法分析了减阻杆加后体这一弹性结构的气动阻尼ꎬ结果表明减阻杆造成的分离流会对后体的气动阻尼系数产生影响ꎮ朱剑等[23]针对新一代捆绑式运载火箭发展了非结构网格下的气动阻尼计算方法ꎬ并分析了攻角㊁Mach数等参数对气动阻尼的影响ꎮ本文在之前的计算方法[23]的基础上采用IDDES模型ꎬ考虑脉动压力的影响ꎬ通过强迫振动的方式ꎬ针对捆绑式运载火箭的某一特定模态进行数值计算仿真ꎬ研究前节点位置ꎬ振动振幅ꎬ脉动压力等参数对气动阻尼的影响规律ꎮ1㊀计算方法图1为本文所用的捆绑式运载火箭的计算模型ꎬ是典型的锤头体结构ꎮ在跨声速阶段ꎬ其头部会产生激波造成激波边界层干扰ꎬ而在锤头体外形的过渡段会出现气流分离ꎮ为探究各部分气动阻尼的变化ꎬ将整个箭体分为头部㊁过渡段㊁弹身3个部分ꎮ图1㊀表面网格及区域划分Fig.1㊀Surfacegridandregiondivision1.1㊀流场仿真模型本文分别用Reynolds平均法(Reynolds ̄averagedNavier ̄StokesꎬRANS)和改进的延迟分离涡模拟(improveddelayeddetached ̄eddysimulationꎬID ̄DES)[24 ̄25]进行计算ꎬ在RANS方程中ꎬ将变量分为平均值和波动值两部分ꎬ对于速度分量有ui=ui+uᶄi其中ꎬi=1ꎬ2ꎬ3ꎬui和uᶄi分别代表平均量和波动量ꎬ对于压强和其他标量也采用类似的形式ꎬ将这种形式代入连续性方程和动量方程中ꎬ并写成张量形式∂ρ∂t+∂∂xi(ρui)=0(1)∂∂t(ρui)+∂∂xj(ρuiuj)=∂p∂xi+∂∂xjμ∂ui∂xj+∂uj∂xi-23δij∂uk∂xkæèçöø÷éëêêùûúú+∂∂xj(-ρuᶄiuᶄj)(2)其中ꎬiꎬjꎬk可分别取1ꎬ2ꎬ3ꎻρ是密度ꎻt是时间ꎻ当i=j时δij取0ꎬ否则取1ꎮ式(1)㊁(2)是RANS方程ꎬ由方程可知RANS方法将湍流脉动对平均流动的作用模化为Reynolds应力项即-ρuᶄiuᶄjꎬ之后采用湍流模型进行封闭ꎬ本文采用的湍流模型为SSTk ̄ω模型ꎬ其输运方程为∂∂t(ρk)+∂∂xi(ρkui)=∂∂xjΓk∂k∂xjæèçöø÷+Gk-Yk∂∂t(ρω)+∂∂xi(ρωui)=∂∂xjΓω∂ω∂xjæèçöø÷+Gω-Yω其中ꎬk和ω分别代表湍流动能和湍流耗散率ꎬΓk和Γω分别代表k和ω的有效扩散系数ꎬGk和Gω分别代表k和ω的生成率ꎬYk和Yω分别代表k和ω的耗散率ꎮ因此RANS方法只能计算大尺度的平均流动ꎬ本文采用IDDES方法计算脉动压力对气动阻尼的影响ꎮIDDES方法是由分离涡模拟(detached ̄eddysimulationꎬDES)方法改进而来ꎬ其本质思想与DES方法相同ꎬ是想以网格尺度和模型中的特征尺度隐式划分RANS和大涡模拟(large ̄eddysimulationꎬLES)区域ꎬ使其既能处理RANS方法无法得到的脉动场ꎬ也能降低LES方法在模拟高Reynolds数流动时所需的计算资源ꎮ区别在于当边界层较厚或者分离区域较窄时ꎬDES方法会出现如模型应力损耗(modeledstressdepletionꎬMSD)ꎬ网格诱导分离(grid ̄inducedseparationꎬGIS)以及对数层不匹配(logarithmic ̄layermismatchꎬLLM)问题[24]ꎬ而IDDES模型通过改良计算区域划分ꎬ结合延迟分离涡模拟(delayeddetached ̄eddysimulationꎬDDES)和03第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析壁面模型大涡模拟(wall ̄modeledlarge ̄eddysimula ̄tionꎬWMLES)ꎬ定义新的长度尺度解决了这些问题ꎬ具体公式详见文献[25]ꎮ流场网格如图2㊁图3所示ꎬ边界层采用棱柱层结构ꎬ并调整第1层网格高度使得y+小于1ꎬ远场部分采用六面体结构网格ꎬ与边界层的过渡层采用非结构网格ꎮ整体网格单元数量为4.2ˑ106ꎮ图2㊀y方向截面网格示意图Fig.2㊀Schematicdiagramofcross ̄sectionalgridinthey ̄direction图3㊀x方向截面网格示意图Fig.3㊀Schematicdiagramofcross ̄sectionalgridinthex ̄direction物面边界条件为无滑移壁面条件ꎬ远场采用压力远场边界条件ꎬ湍流模型采用SSTk ̄ω模型ꎬ采用密度基求解ꎬ气体黏性采用Sutherland定律ꎬ空间离散采用2阶迎风格式ꎬ对流通量采用Roe格式ꎮ1.2 结构分析模型结构与流场耦合分析过程中ꎬ结构部分可以采用模态方法描述ꎮ结构模态可以通过有限元方法与结构模态试验方法获得ꎮ本文采用有限元分析结果获得的模态ꎬ图4所示为结构的前3阶模态ꎬ本文只分析计算结果中气动阻尼最小的第2阶模态ꎮ(a)f=1.200Hz(b)f=2.460Hz(c)f=2.957Hz图4㊀结构的前3阶模态Fig.4㊀Firstthreemodesofthestructure由于火箭结构外形简单ꎬ一般不考虑其扭转影响ꎬ因此可以将其简化为简单的梁模型ꎬ这样就可以给出其模态振动方程q㊆i+2biωiq˙i+ω2iqi=fi(3)式中ꎬqi为第i阶模态的广义位移ꎬbi为第i阶模态的结构阻尼系数ꎬωi为第i阶模态的固有频率ꎬ13气体物理2023年㊀第8卷fi为第i阶模态下质量归一化的广义气动力ꎮ若将fi按照Taylor展开并略去高阶项ꎬ可以将其转化为气动阻尼项与气动刚度项的形式ꎬ则式(3)可写为q㊆i+2(bi+Bi)ωiq˙i+(Ki+1)ω2iqi=0(4)式中ꎬBi为气动阻尼系数ꎬKi为气动刚度系数ꎬ研究表明[26]ꎬ气动刚度相对于结构刚度为小量可以忽略不计ꎬ而在计算中结构阻尼往往设置为0ꎬ因此气动阻尼可以直接反映其气弹稳定性ꎮ1.3㊀气动阻尼分析原理气动阻尼的分析可以采用强迫振动或者自由振动的方式进行ꎬ这两种方法获得的时域数据不同ꎬ提取气动阻尼的方式也不同ꎮ强迫振动方法初始演化过程较短ꎬ因此计算量较小ꎬ同时能够分析某一种振动形式的气动阻尼ꎬ明确该振动形式是收敛还是发散ꎮ分析过程中能够获得不同部位与部件的气动阻尼ꎮ但是对于多模态相互作用引起的发散(例如颤振)较难预测ꎮ自由振动方法需要一定的自由演化时间才能够对时域数据进行分析ꎬ不过自由振动方法能够获得最能够吸收能量的模态及其振动频率ꎮ对于本研究所关注的问题ꎬ气动载荷对结构振动的过程中气动阻尼的影响较大ꎬ而对气动刚度与气动质量影响较小ꎬ即结构的固有振动频率受到来流的影响较小ꎬ其稳定性问题主要由气动阻尼的正㊁负引起ꎬ所以采用强迫振动方法分析ꎮ强迫振动下结构做简谐模态振动qi(t)=Asin(ωit)式中ꎬA表示振动的振幅ꎬ将其代入计算气动力的公式中[21]并做正交积分可得Bi=ʏl0Bx(x)dx=-1MiAω2iTʏl0ʏt0+Tt0G(xꎬt)cos(ωit)dtdx(5)式中ꎬMi为第i阶模态的模态质量ꎬT为整数倍周期ꎬG为广义气动力ꎮ根据式(5)便可以得到局部或分区域的气动阻尼ꎮ1.4㊀耦合计算流程首先进行模态分析ꎬ以确定结构的模态频率与振型ꎬ用以设计强迫振动的频率和振幅ꎮ非定常流场计算前先进行定常流场计算ꎬ来加快非定常计算的演化速度并增强收敛性ꎬ结构节点位移通过径向基函数(RBF)插值方法[27]映射到气动网格节点上ꎬ来进行网格的变形ꎬ这里径向基函数选用WendlandC2ꎬ如下所示φ(x)=(1-x)4(4x+1)最后将计算出来的广义力提取出来ꎬ截取演化完毕的整数倍周期ꎬ进行气动阻尼计算ꎮ耦合计算流程图如图5所示ꎮ图5㊀耦合计算流程图Fig.5㊀Flowchartofcoupledcalculation2㊀结果分析与讨论2.1㊀流场分析结果计算的来流Mach数范围为0.7~1.2ꎮ其中中截面的压力分布如图6所示ꎮ可以看出在头部出现了膨胀波以及跨声速激波ꎬ在过渡段存在流动分离ꎬ随着Mach数的增大ꎬ头部低压区域逐渐扩张ꎬ并且能明显看到ꎬ在流动再附的位置产生了再附激波ꎮ(a)Ma=0.7023第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析(b)Ma=0.75(c)Ma=0.80(d)Ma=0.85(e)Ma=0.88(f)Ma=0.90(g)Ma=0.92(h)Ma=0.96(i)Ma=0.9833气体物理2023年㊀第8卷(j)Ma=1.00(k)Ma=1.05(l)Ma=1.10图6㊀不同Mach数下的中截面压力分布Fig.6㊀PressuredistributioninthemiddlesectionatdifferentMachnumbers2.2 气动阻尼分布通过上述流场分析ꎬ可以看出火箭不同部位流动结构并不相同ꎬ在头部与箭身上ꎬ流动主要为附着流动ꎬ而在过渡段会出现较为复杂的波系结构以及流动分离ꎮ针对不同的流动结构随流向站位x的变化ꎬ设该位置上广义力与广义位移的相位差为φ(x)ꎬ并且简谐振动没有引入其他模态的广义力ꎬ则广义力的表达式为G(xꎬt)=Fgen sin[ωt+φ(x)]+F0(6)其中ꎬFgen为广义力的振动幅度ꎬF0为广义力的常数偏移量ꎮ将式(6)代入到式(5)中得到B(x)=-FgenMAω2Tʏt0+Tt0sin[ωt+φ(x)]cos(ωt)dt其中ꎬ广义力的常数偏移量F0的积分为0ꎬ因此省略ꎮ通过将等式中的正弦函数部分进行和差化积得到B(x)=-FgenMAω2Tʏt0+Tt0sin(ωt)cos[φ(x)]cos(ωt)dt+[ʏt0+Tt0sin[φ(x)]cos(ωt)cos(ωt)dt](7)式(7)中第1部分在整个周期中的积分为0ꎬ只有第2部分保留ꎬ因此得到B(x)=-Fgensin[φ(x)]MAω2Tʏt0+Tt0cos2(ωt)dt(8)式(8)中积分部分恒为正值ꎬ决定整个气动阻尼的部分只有相位角φ(x)的正弦值sin[φ(x)]ꎬ为了能够更加直观地获得相位角与气动阻尼B之间的关系ꎬ须将符号转化为对应的正弦函数转角ꎬ根据正弦关系ꎬ此转角为πꎬ因此得到B(x)=-Fgen(x)sin[φ(x)+π]MAω2Tʏt0+Tt0cos2(ωt)dt(9)图7为气动阻尼变化曲线ꎬ可以看出随着Mach数的增大ꎬ整体气动阻尼先增大后减少ꎬ在Mach数为0.98时达到最大值ꎬ过渡段与箭体的气动阻尼变化趋势与整体基本相同ꎬ而头部区域则不同ꎬ是随着Mach数的增大一直增大ꎬ只是增长速率变缓ꎮ图7㊀有助推时气动阻尼变化曲线Fig.7㊀Aerodynamicdampingchangecurvewithboost根据式(9)ꎬ得到相位角与气动阻尼B之间的关系为:当φ(x)ɪ(-πꎬ0)时ꎬ相位角滞后ꎬ气动阻尼B为负值ꎻ当φ(x)ɪ(0ꎬπ)ꎬ相位角提前ꎬ43第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析气动阻尼B为正值ꎻ为当φ(x)=0时ꎬ无相位角差别ꎬ气动阻尼B为0ꎮ在过渡段上ꎬ复杂的波系结构以及流动分离ꎬ使得气动力与结构位移之间会出现较为明显的迟滞现象ꎬ从而导致相位角φ(x)ɪ(-πꎬ0)ꎬ由此在过渡段上产生了负的气动阻尼ꎮ计算过程中的广义力与广义位移随时间变化曲线如图8所示ꎬ可以看出所有工况计算结果都表现良好ꎬ需要注意的是在非定常计算初期ꎬ演化的不完全导致广义力存在一些突变异常的结果ꎬ计算气动阻尼时须剔除ꎬ选择后面演化完全的周期ꎮ本文计算了9个周期ꎬ剔除了第1个周期出现的错误结果ꎬ采用后8个周期进行气动阻尼分析ꎮ强迫运动振幅为芯级直径的0.5%ꎮ(a)Ma=0.70㊀㊀㊀(b)Ma=0.75(c)Ma=0.80㊀㊀㊀(d)Ma=0.85(e)Ma=0.88㊀㊀㊀(f)Ma=0.9053气体物理2023年㊀第8卷(g)Ma=0.92㊀㊀㊀(h)Ma=0.96(i)Ma=0.98㊀㊀㊀(j)Ma=1.00(k)Ma=1.05㊀㊀㊀(l)Ma=1.10图8㊀不同工况下的广义力与广义位移随时间变化曲线Fig.8㊀Timedependentcurvesofgeneralizedforceandgeneralizeddisplacementunderdifferentoperatingconditions2.3㊀气动阻尼影响因素2.3.1㊀有无助推对气动阻尼的影响捆绑式运载火箭相比于传统的运载火箭ꎬ最大的区别就是在尾部四周捆绑了助推器ꎬ使得其流场特性变得复杂ꎬ因此须分析其对气动阻尼的影响ꎮ图7㊁图9分别为有无助推时气动阻尼变化曲线ꎬ可以看出随着Mach数的增大整体气动阻尼先增大后减少ꎬ在Mach数为0.98时达到最大值ꎬ过63第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析渡段与箭体的气动阻尼变化趋势与整体基本相同ꎬ而头部区域则不同ꎬ是随着Mach数的增大一直增大ꎬ只是增长速率变缓ꎮ对比两个图可知ꎬ助推主要起增大气动阻尼的作用ꎮ还可以看出有无助推情况下头部的气动阻尼变化很小ꎬ意味着在箭体尾部施加控制很难影响到头部的气动阻尼ꎬ特别是在超声速流场中ꎮ图9㊀无助推时气动阻尼变化曲线Fig.9㊀Aerodynamicdampingchangecurvewithoutboost2.3.2㊀前节点位置影响为了考察前节点位置变化对气动阻尼的影响ꎬ在保持振动频率不变㊁头部最大振型位置与振幅不变的条件下移动前节点ꎬ变化后的振型如图10所示ꎮ(a)Frontnodeafterthetransitionregion(b)Frontnodeinthetransitionregion(c)Frontnodebeforethetransitionregion图10㊀前节点变化后的振型Fig.10㊀Vibrationmodeafterthechangeofformernode根据对计算结果的分析分别获得了不同前节点位置的整体气动阻尼对比与过渡段气动阻尼对比ꎬ如图11㊁图12所示ꎬ可以看出前节点位置的改变并没有影响整体气动阻尼随Mach数增大而增大的趋势ꎬ且前节点在过渡段上与过渡段前的整体气动阻尼相差不大ꎬ而前节点在过渡段后的整体气动阻尼要高于另两种情况ꎬ因此过渡段与头部放在同一侧有助于提高气动阻尼ꎮ过渡段的气动阻尼会随着前节点的变化发生剧烈改变ꎬ前节点在过渡段前后随Mach数增大的变化规律相反ꎬ节点前后的振动相位变化导致不同节点位置过渡段的振动相位不同ꎬ进而导致气动阻尼发生变化ꎮ图11㊀不同节点位置的整体气动阻尼Fig.11㊀Overallaerodynamicdampingatdifferentnodepositions图12㊀不同节点位置的过渡段气动阻尼Fig.12㊀Aerodynamicdampingofthetransitionregionatdifferentnodepositions2.3.3㊀强迫振动振幅大小对气动阻尼的影响为了考察强迫振动振幅大小对气动阻尼的影响ꎬ在保证流场结构不发生改变的前提下ꎬ振动振幅分别为原来的一半和两倍ꎬ根据工程经验ꎬ如果振幅超过芯级直径的5%ꎬ则须考虑流场结构改变所造成的影响ꎮ图13㊁图14分别为不同振幅下的整体与头部气动阻尼ꎮ73气体物理2023年㊀第8卷图13㊀不同振幅下整体气动阻尼Fig.13㊀Overallaerodynamicdampingatdifferentamplitudes图14㊀不同振幅下头部气动阻尼Fig.14㊀Aerodynamicdampingoftheheadregionatdifferentamplitudes可以发现改变振幅无论是对整体气动阻尼还是头部气动阻尼来说变化都很小ꎬ这意味着气动阻尼的大小主要取决于气动力与结构振动的相位差ꎬ不依赖于振动幅度的大小ꎮ2.3.4㊀脉动压力对气动阻尼的影响为了模拟出脉动压力的影响ꎬ采用IDDES方法对火箭气动阻尼进行计算ꎬ计算来流Mach数为0.92ꎬ计算过程中的广义力与广义位移如图15所示ꎬ相较于图8可以看出广义力随时间变化曲线并不光滑ꎬ脉动压力的存在导致广义力由多个频率叠加而成ꎮ由于第2阶模态的频率为2.46Hzꎬ而由分离流㊁激波振荡等引起的脉动压力频率往往远大于此频率ꎬ因此这里选择3.5Hz为分界ꎬ将高于3.5Hz的部分视为由抖振脉动压力引起的广义力ꎬ低于3.5Hz的部分视为强迫振动引起的广义力ꎬ通过低通滤波把高于3.5Hz的广义力滤掉ꎬ可以获得由强迫振动引起的广义力与广义位移变化曲线ꎬ如图16所示ꎬ通过此广义力计算的气动阻尼为2.08ɢꎮ同样地ꎬ进行高通滤波将低于3.5Hz的广义力滤掉ꎬ可以获得由抖振脉动压力引起的气动阻尼为(2.94ˑ10-3)ɢꎬ由此得到脉动压力引起的气动阻尼变化为0.14%ꎬ可以忽略不计ꎮ同时使用RANS方法计算的气动阻尼为2.07ɢꎬ与IDDES的计算结果相比误差约为(2.94ˑ10-3+2.08-2.07)/2.07ʈ0.48%ꎬ这说明针对气动阻尼的模拟ꎬ抖振引起的脉动压力对气动阻尼的计算结果影响很小ꎬ起主要作用的还是广义力的变化ꎬ该变化由强迫振动引起的结构边界变化所导致ꎮ图15㊀基于IDDES的广义力与广义位移变化曲线Fig.15㊀VariationcuresofgeneralizedforceandgeneralizeddisplacementbasedonIDDES图16㊀滤波后的广义力与广义位移变化曲线Fig.16㊀Variationcuresofgeneralizedforceandgeneralizeddisplacementvariationcurveafterfiltering3㊀结论本文通过数值计算方法研究了火箭的气动阻尼特性ꎮ根据流动特征分析与理论推导ꎬ发现火箭过渡段几何外形的收缩导致该区域出现复杂的分离与激波结构ꎬ从而造成了气动力相对于结构振动83第6期李泳德ꎬ等:基于CFD/CSD耦合的火箭跨声速气动阻尼特性分析相位的滞后ꎬ导致了该区域为气动负阻尼ꎬ即气动不稳定性的主要来源ꎮ在此机理的基础上ꎬ分析了前节点位置㊁振动振幅㊁脉动压力等因素对气动阻尼的影响规律ꎮ可以得出以下结论:1)助推增加了正阻尼区域的面积ꎬ从而相对于没有助推的构型起到了增加气动阻尼的作用ꎮ2)前节点位置的改变对过渡段气动阻尼影响很大ꎬ节点前后的振动方向相反ꎬ导致节点在过渡段前后的气动阻尼变化规律也截然相反ꎬ将过渡段与头部区域放在节点的同一侧有助于增加气动阻尼ꎮ3)在不改变流场结构的前提下ꎬ改变振动的振幅ꎬ气动力也会产生相应幅度的变化ꎬ因此结构振幅对气动阻尼的影响可忽略不计ꎮ4)高频部分的广义力对气动阻尼的贡献很小ꎬ即结构振动引起的广义力变化对气动阻尼起主要作用ꎬ而脉动压力对计算气动阻尼影响不大ꎬ可忽略不计ꎮ参考文献(References)[1]㊀CoeCF.Steadyandfluctuatingpressuresattransonicspeedsontwospace ̄vehiclepayloadshape[R].NASATMX ̄503ꎬ1961.[2]ColeSRJrꎬHenningTLꎬRaineyAG.NASAspaceve ̄hicledesigncriteria[R].NASASP ̄8001(REV)ꎬ1964. [3]EricssonLEꎬRedingJP.Analysisofflowseparationeffectsonthedynamicsofalargespacebooster[J].Jour ̄nalofSpacecraftandRocketsꎬ1965ꎬ2(4):481 ̄490. [4]RedingJPꎬEricssonLE.Effectofaeroelasticconsidera ̄tionsonseasat ̄Apayloadshrouddesign[J].JournalofSpacecraftandRocketsꎬ1981ꎬ18(3):241 ̄247. [5]程镇煌.宇航飞行器跨音速气动弹性问题探讨[J].上海航天ꎬ1997(6):16 ̄21.ChengZH.Discussionontransonicpneumaticelasticityofspacecraft[J].AerospaceShanghaiꎬ1997(6):16 ̄21(inChinese).[6]倪嘉敏.我国运载火箭气动设计回顾[C].近代空气动力学研讨会论文集ꎬ2005.NiJM.ReviewofaerodynamicdesignofChinaᶄslaunchvehicle[C].Proceedingsofthemodernaerodynamicssymposiumꎬ2005(inChinese).[7]臧涛成ꎬ胡焕性.大长细比弹箭弹性效应研究综述[J].弹道学报ꎬ1999ꎬ11(3):89 ̄93ꎬ96.ZangTCꎬHuHX.Areviewofgreatslendernessratioprojectileelasticeffectresearch[J].JournalofBallisticsꎬ1999ꎬ11(3):89 ̄93ꎬ96(inChinese). [8]吴志刚ꎬ杨超.细长体弹箭的气动弹性问题与研究方法[C].第九届(2005年)全国空气弹性学术交流会论文ꎬ2005.WuZGꎬYangC.Aeroelasticityproblemsandresearchmethodsofslenderbodymissiles[C].9thNationalSym ̄posiumonAeroelasticityꎬ2005(inChinese). [9]张贺ꎬ黄晓鹏.弹性细长旋转弹箭运动稳定性问题的研究进展[C].第九届(2005年)全国空气弹性学术交流会论文ꎬ2005.ZhangHꎬHuangXP.Researchprogressonthestabilityofelasticslenderrotatingprojectiles[C].9thNationalSymposiumonAeroelasticityꎬ2005(inChinese). [10]ColeHAJr.Dynamicresponseofhammerheadlaunchvehiclestotransonicbuffeting[R].NASATND ̄1982ꎬ1963.[11]BartelsREꎬWiesemanCDꎬMineckRE.ComputationalaeroelasticanalysisoftheAreslaunchvehicleduringas ̄cent[R].AIAA2010 ̄4374ꎬ2010.[12]AzevedoJLF.Aeroelasticanalysisoflaunchvehiclesintransonicflight[J].JournalofSpacecraftandRocketsꎬ1989ꎬ26(1):14 ̄23.[13]SinclairAꎬFlowersG.Low ̄orderaeroelasticmodeloflaunch ̄vehicledynamics[R].AIAA2010 ̄7725ꎬ2010. [14]DotsonKW.Transientcouplingoflaunchvehiclebendingresponseswithaerodynamicflowstatevariations[J].JournalofSpacecraftandRocketsꎬ2001ꎬ38(1):97 ̄104.[15]ColeSRꎬHenningTL.Buffetresponseofahammerheadlaunchvehiclewind ̄tunnelmodel[J].JournalofSpacecraftandRocketsꎬ1992ꎬ29(3):379 ̄385.[16]崔尔杰.流固耦合力学研究与应用进展[C].钱学森科学贡献暨学术思想研讨会论文集ꎬ2001.CuiEJ.Researchandapplicationprogressoffluid ̄struc ̄tureinteractionmechanics[C].ProceedingsofSeminarofQianXuesenScientificContributionsandAcademicThoughtsꎬ2001(inChinese).[17]冯明溪ꎬ王志安.火箭跨音速动导数和抖振实验[J].宇航学报ꎬ1987(1):55 ̄63.FengMXꎬWangZA.Experimentsoftransonicderiva ̄tivesandbuffetingofrocket[J].JournalofAstronauticsꎬ1987(1):55 ̄63(inChinese).[18]白葵ꎬ冯明溪.弹性模型实验技术[J].流体力学实验与测量ꎬ1999ꎬ13(1):38 ̄42.BaiKꎬFengMX.Aeroelasticmodelandthebuffetex ̄perimentaltechnique[J].ExperimentsandMeasurementsinFluidMechanicsꎬ1999ꎬ13(1):38 ̄42(inChinese). [19]JiCꎬRanJHꎬLiFꎬetal.Theaerodynamicdamping93气体物理2023年㊀第8卷testofelasticlaunchvehiclemodelintransonicflow[C].Proceedingsofthe64thInternationalAstronauticalCon ̄gressꎬ2013.[20]季辰ꎬ吴彦森ꎬ何岗ꎬ等.运载火箭气动阻尼风洞试验研究[C].第十二届全国空气弹性学术交流会论文集ꎬ2011.JiCꎬWuYSꎬHeGꎬetal.Experimentalstudyonaero ̄dynamicdampingwindtunneloflaunchvehicle[C].Pro ̄ceedingsofthe12thNationalAeroelasticityConferenceꎬ2011(inChinese).[21]刘子强ꎬ白葵ꎬ毛国良ꎬ等.锤头体弹性振动跨音速气动阻尼系数的确定[J].宇航学报ꎬ2002ꎬ23(6):1 ̄7.LiuZQꎬBaiKꎬMaoGLꎬetal.Thedeterminationofaerodynamicdampingonhammerheadlaunchvehiclesattransonicspeeds[J].JournalofAstronauticsꎬ2002ꎬ23(6):1 ̄7(inChinese).[22]冉景洪ꎬ刘子强ꎬ胡静ꎬ等.减阻杆气动阻尼研究[J].力学学报ꎬ2014ꎬ46(4):636 ̄641.RanJHꎬLiuZQꎬHuJꎬetal.Researchofaero ̄dampingforbluntwithspike[J].ChineseJournalofThe ̄oreticalandAppliedMechanicsꎬ2014ꎬ46(4):636 ̄641(inChinese).[23]朱剑ꎬ冉景洪ꎬ吴彦森ꎬ等.捆绑式运载火箭的气动阻尼数值计算方法[C].第十三届全国空气弹性学术交流会论文集.哈尔滨:中国力学学会ꎬ中国空气动力学会ꎬ2013.ZhuJꎬRanJHꎬWuYSꎬetal.Numericalcalculationmethodforaerodynamicdampingofbundlelaunchvehicles[C].Proceedingsofthe13thNationalAeroelasticityCon ̄ference.Harbin:ChineseSocietyofTheoreticalandAp ̄pliedMechanicsꎬChinaAeromechanicsSocietyꎬ2013(inChinese).[24]SpalartPR.Detached ̄eddysimulation[J].AnnualReviewofFluidMechanicsꎬ2009ꎬ41:181 ̄202. [25]GritskevichMSꎬGarbarukAVꎬSchützeJꎬetal.Devel ̄opmentofDDESandIDDESformulationsforthek ̄ωshearstresstransportmodel[J].FlowꎬTurbulenceandCombustionꎬ2012ꎬ88(3):431 ̄449.[26]季辰ꎬ吴彦森ꎬ侯英昱ꎬ等.捆绑式运载火箭跨声速气动阻尼特性试验研究[J].实验流体力学ꎬ2020ꎬ34(6):24 ̄31.JiCꎬWuYSꎬHouYYꎬetal.Experimentalstudyofaerodynamicdampingcharacteristicsofalaunchvehiclewithboostersintransonicflow[J].JournalofExperimentsinFluidMechanicsꎬ2020ꎬ34(6):24 ̄31(inChinese). [27]AllenCꎬRendallTCS.UnifiedapproachtoCFD ̄CSDinterpolationandmeshmotionusingradialbasisfunctions[R].AIAA2007 ̄3804ꎬ2007.04。
2023 年湍流与噪声和 CFD 方法暑期高级讲习班 会议手册说明书
2023年湍流与噪声和CFD方法暑期高级讲习班2023 Advanced Summer Program on Turbulence,Noise and CFD Methods会议手册时间:2023年7月28至8月5日地点:香港科技大学主办单位:中国空气动力学会承办单位:香港科技大学(HKUST)上海大学南方科技大学复旦大学中国空气动力学会CFD专委会中国空气动力学会低跨超专委会上海市应用数学和力学研究所上海市力学信息学前沿科学基地上海市能源工程力学重点实验室粤港澳数据驱动下的流体力学与工程应用联合实验室中国航空学会航空声学分会协办单位:《空气动力学学报》《实验流体力学》《Advances in Aerodynamics》二零二三年七月二十六日2023年湍流与噪声和CFD方法暑期高级讲习班为了促进流体力学与空气动力学的发展、推动学术交流与合作、培育培养优秀人才,助力解决流体力学与空气动力学等相关领域“卡脖子”技术,经中国空气动力学会批准,2023年湍流与噪声和CFD 方法暑期高级讲习班将于2023年7月28日至8月5日在香港科技大学(HKUST)举行。
会议邀请内地与香港地区在湍流、噪声和CFD方法等方面的专家学者、青年学者为讲习班授课。
现诚邀内地与港澳台地区研究生、工程师、相关领域专家学者以及高年级本科生参会。
本次讲习班由中国空气动力学会主办,香港科技大学(HKUST)、上海大学、南方科技大学、复旦大学、中国空气动力学会CFD专委会、中国空气动力学会低跨超专委会、上海市应用数学和力学研究所、上海市力学信息学前沿科学基地、上海市能源工程力学重点实验室、粤港澳数据驱动下的流体力学与工程应用联合实验室等单位承办。
本次讲习班采用线上线下同时进行的方式,其中线上使用腾讯会议App进行直播,会议号码:964-8147-9182,也可直接扫描下面的二维码参会:2023年湍流与噪声和CFD方法暑期高级讲习班专家报告日程安排报告安排以专家自选日程排列,不分先后次序,后续如有变动以最终表格为准。
【国家自然科学基金】_跨声速压气机_基金支持热词逐年推荐_【万方软件创新助手】_20140730
2011年 科研热词 推荐指数 序号 压气机 2 1 内燃机 2 2 高精度weno(weighted essential 1 non-oscillatory)格式 3 面积律 1 4 非定常作用 1 5 非均匀 1 6 间隙泄漏涡 1 7 轴流压气机 1 8 转子 1 9 跨音速压气机 1 跨声速 1 跨声压气机 1 计算效率 1 端壁 1 离心压气机 1 激波 1 涡轮增压器 1 涡轮增压 1 流动稳定性 1 流动控制 1 泄漏 1 气动负荷 1 机匣处理 1 旋涡 1 新型高分辨率格式 1 数值模拟 1 数值仿真 1 掠、非定常 1 弯 1 子午流道 1 奇异分析 1 失速 1 多分辨分析 1 复杂流场计算 1 堵塞 1 周向局部处理机匣 1 叶尖间隙 1 二维小波 1 三维小波 1
推荐指数 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2010年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39
2014年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
2014年 科研热词 离心压气机 数值模拟 压气机 非定常流动 非定常 静叶 进气畸变 轴流压气机 跨声速轴流压气机 跨声速压气机转子 试验设计 角区分离 角区 自循环机匣处理 统计学分析 端壁融合 稳定裕度 激波结构 泄漏涡 模态分析 旋转不稳定 数值计算 抽吸环槽 总压畸变范围 总压畸变 失速机理 叶顶泄漏涡 叶顶喷气 叶身 叶尖泄漏流 前缘 倒圆 推荐指数 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
阻尼台对轴流跨声速压气机性能的数值研究
密 。处理 机匣 采用H型 网格 。为 了准确模 拟叶尖 间
隙 区的流动 ,在间 隙区沿径 向均匀分 布 l个节 点 , 7
并对转 子顶 部间 隙区采用 了“ 蝶型 网格” ,获得 正交 性 良好 的网格 。图 1 给出 了带 阻尼 台的单 级轴 流 、2
压 气机 三维 图 、 片通 道计算 网格简 图。8 周 向 ) 叶 4(
22 压 气 机 内部 流场 对 比分 析 .
图 4 出了约 6 %处 叶展处 流 面上 的相对 马 给 5
赫数 分布 云图 , 中 6 其 5%叶 高的 流面靠 近阻尼 台
的上壁 面 。从 图 4中可观察 到 ,带 阻尼 台后 ,动 叶 吸力 面附近 的低相对 马赫数 区域较无 阻尼 台时的范 围大得 多 ,特 别是靠 近尾缘处 时 ,这 意味该处 的流 动损失 大 。同时还观察 到带阻尼 台后 ,通道激波 位
第2 5卷 第 2期 21 0 0年 3月
海 军 航 空 工 程 学 院 学 报
J u n l f v l r n u ia n  ̄ n u e l i e st o r a Na a o a t l d As o a t a v r i o Ae c a i Un y
、 12 No2 , .5 o . Ma . 0 O r2 l
内部流场进行模拟 ,着重讨论了压气机转子 中带阻尼台部分 区域及端壁 区域 的流场 。数值计算结果指 出阻尼台 的引入降低了该压气机 的性 能。流场分析表明 ,带阻尼台后增大 了动叶 中部通道 ( 阻尼台附近 )及静叶靠 近机 匣的区域 的流动损失 ,但对动叶顶部流场 、动叶吸力面附近的分离线位置和形状影响很小 。 关键词 :跨声速风扇 ;阻尼台;数 值模 拟 总性 能 比较 .
跨声速轴流压气机不同工况下转子激波与导流叶片尾迹相互干涉研究
关 键 词 : 声 速 轴 流 压 气 机 ; 口导 流 叶 片 ; 波 一 迹 干涉 ; 跨 进 激 尾 数值 模 拟 ;
Iv sia in o n e a t n ewe n R t rSh c n GV W a e f n e t t n It r ci s b t e o o o ka d I g o o k or Tr n o i xaI o Co p e s ri fe e tCo dt n a s nc A i w m r s o Di r n n io s Fl n i
s g rno i a i o o pesrw scn utd b he ies n i osfo t e t sn xa f w c m rso a od c y tre dm ni a vs u w a a c l l e ol c l c cl in s toe Fn ̄ u ec T e c aat ii c  ̄ s o e t no i ai o l t a ua o o wl ie m ca h h rc r t u e t r snc x a f w f ' es c f h a ll c m rso e ot nd T e ntay c c l i o o pesr it n o W S o pesr w r ba e . h used a u o f cm rso ne a f w a e i l a n t rl l
摘要: 应用 3维黏 性流动计算软 件 Fn/ u e a对 某型 1 ieN m c , . 5级跨声速 轴流压 气机
进行 内部流场和全工况特性数值模拟 , 到该压 气机特性 曲线。 得 在最高效率点和近失速
点, 2种质量流 量下 , 在 对压气机 内部 流场进行非定 常对比计算 , 分析 了在质量 流量变
【国家自然科学基金】_叶尖泄漏流_基金支持热词逐年推荐_【万方软件创新助手】_20140801
2012年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
பைடு நூலகம்
科研热词 高负荷 非定常相互作用 间隙流 轴流压气机 涡轮端区 泄漏涡 泄漏流 敏感度 小波分析 叶尖间隙 叶尖二次涡 动态压力测量 前掠 位势场 二次流
推荐指数 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2008年 序号 1 2 3 4 5 6 7 8 9 10 11
科研热词 数值模拟 叶尖间隙高度 非定常流动 被动控制 肋条叶尖 离心压缩机 涡轮叶栅 实验研究 叶尖间隙流 叶尖泄漏涡 clocking效应
推荐指数 2 2 1 1 1 1 1 1 1 1 1
2009年 序号 1 2 3 4
科研热词 涡系结构 泄漏流动 对转涡轮 叶尖凹槽
推荐指数 3 2 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2014年 序号 1 2 3 4 5 6 7 8 9 10
2014年 科研热词 非定常流动 非定常数值模拟 间隙泄漏流 跨声速压气机转子 流场 模态分析 旋转不稳定 失速 叶尖泄漏流 压气机 推荐指数 1 1 1 1 1 1 1 1 1 1
推荐指数 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2011年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
2011年 科研热词 涡轮效率 涡轮 泄漏流 周向载荷 叶尖间隙 逆向涡流器 逆向涡流发生器 跨声压气机 泄漏涡 泄漏 机匣端壁引气 引气量 引气位置 堵塞 压气机 推荐指数 2 2 2 2 2 1 1 1 1 1 1 1 1 1 1
跨声速机翼突然失速现象的CFD研究
第十六届全国计算流体力学会议跨声速机翼突然失速现象的 CFD 研究吴佳莉,马海,卜忱(中航工业空气动力研究院,哈尔滨市第 88 号信箱 150001) 摘 要:20 世纪 90 年代中后期,F/A-18E/F 在跨声速机动飞行中多次出现了非指令横向运动引起的“掉翼尖”现象,使飞机的横向稳定性和操纵性下降,对飞行安全构成极大的威胁。
本文基于中航工业空气动力研究院自主研 发的 UNSMB 非结构网格计算平台对带有前缘锯齿的某先进布局飞机基本构型开展了跨声速机翼突然失速 (AWS) 特性计算研究, 分析了引起 AWS 现象的复杂流动机理, 达到了初步预测机翼突然失速特性的能力。
研究结果显示,1阶段原型机进行试飞时,技术人员发现该飞机在马赫数 0.7~0.95,高度 3000~12000 米,迎角 6°~12°的跨 声速飞行条件下,一侧的机翼会发生突然失速(AWS) ,诱发另一侧机翼下沉,伴随非指令性的横向运动, 称之为“掉翼尖”现象,飞机的横向稳定性和操纵性下降,对飞行安全构成极大的威胁[1]。
这个问题后来由 波音公司相继通过改变前缘襟翼偏转程序和在机翼折叠机构上增开了一个多孔门得到初步解决。
AWS 问题 一度使 F/A-18E/F 的研制试飞工作陷于停滞,因此,在 20 世纪 90 年代末美国空军启动了“机翼突然失速” (AWS)专项研究项目,以求借助 CFD 计算、 风洞试验、 地面模拟和飞行试验等多种手段, 了解 F/A-18E/F 和其他飞机布局 AWS 现象发生的机理,并发展用于评估和预测机翼突然失速的计算品质因素。
AWS 现象的实验研究主要是通过对跨声速风洞静态/动态试验得到的力、力矩、压力和 PSP 图像进行 研究得到。
McMillin[2]等人在 LaRC 16ft 跨音速风洞对 8%缩比的预生产型 F/A-18E 模型的 AWS 现象进行了 系统研究(见图 1) 。
转子叶尖间隙形状对跨声速轴流压气机性能影响机理分析
收稿日期:2021-05-14作者简介:张成烽(1996),男,硕士。
引用格式:张成烽,张国臣,徐志晖,等.转子叶尖间隙形状对跨声速轴流压气机性能影响机理分析[J].航空发动机,2023,49(3):66-74.ZHANG Chengfeng ,ZHANG Guochen ,XU Zhihui ,et al.Analysis of effect mechanism of rotor tip clearance shapes on transonic axial compressor[J].Aeroengine ,2023,49(3):66-74.转子叶尖间隙形状对跨声速轴流压气机性能影响机理分析张成烽,张国臣,徐志晖,孙丹,刘鹏程(沈阳航空航天大学航空发动机学院,沈阳110136)摘要:为研究转子不同叶尖间隙形状对跨声速轴流压气机性能的影响机理,分别对平行式叶尖间隙进行渐变式和阶梯式改型优化,并利用商业软件NUMECA 进行数值模拟。
结果表明:对平行式叶尖间隙进行渐变式和阶梯式改型优化后,压气机性能有较大提高,改型后0.204-0.408的各类间隙压气机性能优于0.408-0.204间隙的。
相比较平行间隙PTC 0.204-0.204,渐变式TTC 0.204-0.408和阶梯式STC 0.204-0.408的失速裕度分别提高1.12%和1.61%,峰值效率基本不变,同时近失速工况下的总压比和效率也略有提高。
叶尖泄漏涡得到抑制,间隙处的流体低速区明显减小,流动损失减小,流场得到较大改善。
转子通道的总压比在85%叶高处明显提高,分别提高了1.01%和3.13%。
阶梯式叶尖间隙压气机的静子通道40%叶高处总压损失系数减小达75.4%。
对平行式叶尖间隙进行改型处理能够有效提高压气机性能,且阶梯式叶尖间隙比渐变式的对压气机性能提高的效果更加显著。
关键词:叶尖间隙;叶尖泄漏流;转子;失速裕度;间隙形状;总压比;总压损失系数;跨声速轴流压气机中图分类号:V231.1文献标识码:Adoi :10.13477/ki.aeroengine.2023.03.009Analysis of Effect Mechanism of Rotor Tip Clearance Shapes on Transonic Axial CompressorZHANG Cheng-feng ,ZHANG Guo-chen ,XU Zhi-hui ,SUN Dan ,LIU Peng-cheng (School of Aero Engine ,Shenyang Aerospace University ,Shenyang 110136,China )Abstract :In order to study the effects of different rotor tip clearance shapes on the performance of a transonic axial compressor,the Parallel Tip Clearances (PTC)were modified by Tapered Tip Clearances (TTC)or Stepped Tip Clearances (STC),and the commercial soft⁃ware NUMECA was used for numerical simulation.The results show that the performance of the modified compressor is significantly im⁃proved.After the modifications,the performance of the modified compressor with 0.204-0.408tip clearances is better than that with 0.408-0.204tip pared with the compressors with PTC 0.204-0.204,stall margins of the compressors with TTC 0.204-0.408and STC 0.204-0.408is increased by 1.12%and 1.61%,respectively,while the peak efficiencies remain essentially unchanged.In addition,the total pressure ratios and efficiencies also increase slightly under near-stall conditions.The tip leakage vortices are suppressed,and the tip low-speed flow areas are significantly reduced,resulting in reduced flow loss and greatly improved flow fields.The total pressure ratios at 85%span-height of the rotor passage increase significantly,and they increase by 1.01%and 3.13%respectively.The total pressure losscoefficient at 40%span-height of the stator passage with stepped blade tip clearances decreases by 75.4%.It can be seen that the modifica⁃tion of the parallel tip clearance can effectively improve the performance of the compressor,and the stepped shape has a more prominent ef⁃fect on the compressor performance than the tapered shape.Key words :tip clearance;tip leakage flow;rotor;stall margin;clearance shapes;total pressure ratio;total pressure loss coefficient;transonic axial compressor航空发动机Aeroengine0引言压气机的内部流场复杂多变,尤其在转子叶尖间隙区域流动更为复杂,存在端壁附面层分离、叶片表面附面层分离、间隙泄漏涡、二次流和多次泄漏流等[1-2]。
压气机的内部损失
压气机的内部损失
压气机内部损失主要包括以下几种:
1. 摩擦损失:压气机内部的零部件之间存在摩擦力,导致能量的损失。
摩擦损失会导致压气机效率降低,因此需要采取措施减少摩擦损失,例如使用润滑剂、优化零部件设计等。
2. 泄漏损失:压气机内部的气体可能会从高压区域泄漏到低压区域,导致能量的损失。
泄漏损失会导致压气机效率降低,因此需要采取措施减少泄漏损失,例如提高密封性能、优化零部件设计等。
3. 气流分离损失:压气机内部的气流可能会发生分离,导致能量的损失。
气流分离损失会导致压气机效率降低,因此需要采取措施减少气流分离损失,例如优化叶片设计、提高进口气流速度等。
4. 二次流损失:压气机内部的气流可能会形成二次流,导致能量的损失。
二次流损失会导致压气机效率降低,因此需要采取措施减少二次流损失,例如优化叶片设计、提高进口气流速度等。
5. 冲击损失:压气机内部的气流可能会发生冲击,导致能量的损失。
冲击损失会导致压气机效率降低,因此需要采取措施减少冲击损失,例如优化叶片设计、提高进口气流速度等。
总之,压气机内部损失是影响压气机效率的重要因素,需要采取措施减少这些损失,提高压气机的性能。
跨声速轴流压气机特性预测的损失模型研究
20 0 7年第 3 3卷第 4期
吴
虎 孙
娜 杨金广
(. 1西北工业大学 , 陕西西安 707 ;. 102 2北京动力机械研究所 , 北京
10 7 ) 004
摘要 : 基于公开发表的研究成果 , 完善 了 1 种新 的跨声速轴流压气机 总压损失及落后角预测模 型 , 并发展 了相
Ab ta t s r c :Ba e n t ep bi h d rs ac c i v me t ,t ett l r s u els n e i t n a ge p e it g mo - s d o h u l e e e rh a h e e n s h oa e s r o sa d d vai n l r d ci d s p o n es o e a s me a i o c mp era e i mv d d t e o —d sg e o ma c n l ss me o ft e l fa n w t n o x a f w o os r mp e I a h f r ll n e i n p r r n e a ay i t d o f h h c re p n i g ta s n c a i o c mp e s r i d v lp d,a d t e c re p n i g c mp tr smu ain p o r m s o rs o d n rn o i x a f w o rs o s e eo e ll n or s o d n o ue i lt rg a i h o b i .T e n me ia i l ai n ft e d sg d o —d sg r r a c r r r e n t e t n o i x a o ul t h u r l s c mu t s o e in a f o h n e i p f m n e a ep f m d i r s n c a i f w n e o e o h a ll e mp ' o .T ec mp rs n b t e ac l t n r s t n e t s t i dc t s a e mo e d a ay i t o s o i  ̄ r h o a io ewe n c ua o e u sa d ts u s n i ae t d l l i l e r l h h t t n n a l ssme d h
NACA0012翼型跨声速强迫运动非定常气动力模型
第41卷第11期2020年11月哈㊀尔㊀滨㊀工㊀程㊀大㊀学㊀学㊀报Journal of Harbin Engineering UniversityVol.41ɴ.11Nov.2020NACA0012翼型跨声速强迫运动非定常气动力模型张庆1,2,叶正寅3(1.西安航空学院飞行器学院,陕西西安710077;2.南洋理工大学机械与航空工程学院,新加坡639798;3.西北工业大学航空学院,陕西西安710072)摘㊀要:传统的一阶线性叠加的气动力模型不再适用于现代高机动性飞行器的非定常气动力建模,为了考察更高阶的气动力模型对非定常迟滞效应模拟的适用程度,本文采用自主发展的求解器,分别计算了NACA0012翼型在跨声速来流条件下做单自由度强迫沉浮㊁俯仰以及沉浮/俯仰两自由度耦合运动的非定常气动力的变化规律㊂然后在Etkin 气动力模型的基础上,探讨了不同类型的高阶的气动导数在非定常气动力建模中的作用㊂研究结果表明:将Etkin 气动力模型中升力和俯仰力矩对迎角的导数项由一阶拓展至二阶就可以较为精确地重构出翼型在强迫运动各阶段的非定常升力和俯仰力矩㊂关键词:跨声速;气动导数;气动力建模;沉浮运动;俯仰运动;耦合运动;Etkin 模型;非定常气动力DOI :10.11990/jheu.201903018网络出版地址:http :// /kcms /detail /23.1390.u.20201028.1424.010.html 中图分类号:V211.41㊀文献标志码:A㊀文章编号:1006-7043(2020)11-1683-06Unsteady aerodynamic model of NACA 0012associated with forcedoscillations and translations in transonic flightZHANG Qing 1,2,YE Zhengyin 3(1.School of Aircraft,Xiᶄan Aeronautical University,Xiᶄan 710077,China;2.School of Mechanical &Aerospace Engineering,Nan-yang Technological University,Singapore 639798,Republic of Singapore;3.School of Aeronautics,Northwestern Polytechnical Uni-versity,Xiᶄan 710072,China)Abstract :It is unsuitable to model unsteady aerodynamics for high-agility modern aircraft by the traditional first-or-der linear superposition theory.For this study,in order to investigate the applicability of a higher order aerodynam-ic model to the simulation of hysteresis effects,the unsteady time histories of aerodynamics for NACA0012associat-ed with single-freedom forced motions,plunging and pitching,and coupled plunging and pitching motion under transonic conditions were investigated computationally based on in-house codes.The effects of various aerodynamic derivatives on aerodynamic models are discussed based on the Etkin aerodynamic model.Final results indicated that unsteady lift and pitching moment in forced single or coupled motions could be accurately regenerated if the Etkin model is expanded to the second order derivative of the angle of attack with respect to time.Keywords :transonic;aerodynamic derivative;aerodynamic model;plunging;pitching;coupled motion;Etkin model;unsteady aerodynamics收稿日期:2019-03-06.网络出版日期:2020-10-28.基金项目:国家自然科学基金重点项目(11732013);校级科研基金项目(2018KY1226).作者简介:张庆,男,讲师,博士;叶正寅,男,教授,博士生导师.通信作者:张庆,E-mail:zhangqing2220@.㊀㊀气动导数作为描述飞行器机动飞行和受扰动时气动特性的关键性气动参数,在飞行器气动性能㊁控制系统和总体设计中扮演着非常重要的作用[1-4]㊂在传统的飞行动力学相关问题的研究中,气动力的数据往往基于小扰动线性叠加原理计算出来,在这种准定常假设情况下,气动力仅仅表示为瞬时飞行状态参数的函数,并且可以以一种简单的解析函数关系式表示出来[2-5]㊂但是,现代飞行器的飞行包线普遍向大迎角区域扩展,在大迎角下飞机机动飞行产生的三维非定常分离流和涡流使得空气动力呈现高度非线性特性,气动力和力矩不仅依赖于瞬时迎角㊁侧滑角㊁姿态角等参数,而且与它们的时间历程有关,因此原来使用的低阶线性叠加模型将不再适用[5-6]㊂同时,由于机动飞行状态涵盖了较大的迎角㊁侧滑角㊁角速率的变化范围,如果采用风洞实验或是数值计算模拟,其时间成本和经济成本都难以接受[7-10]㊂哈㊀尔㊀滨㊀工㊀程㊀大㊀学㊀学㊀报第41卷因此有必要建立起较大飞行包线内普适性较好的的非定常气动力模型[1,4]㊂Etkin模型是目前动导数求解时最常用的一种非定常气动力模型,Etkin模型物理意义明确,考虑了时间历史效应对气动导数的影响[3]㊂但是,在非定常气动力建模时,该模型中的各项气动导数对不同运动形式的非定常气动力的影响规律和适用程度尚不清楚㊂为此,本文结合Etkin气动力模型,研究了气动力关于迎角的一阶和二阶导数在气动力模型的作用,希望能精确地重构出翼型单自由度或是耦合强迫运动过程中的非定常气动力,为未来发展高效的㊁可靠的气动力模型提供参考数据㊂1㊀强迫运动非定常气动力模型本文的计算采用课题组自己开发的柔性体动力学问题求解软件GMFlow[11-13],其中流场求解部分采用基于SA模型的有限体积法[13],强迫运动时的网格变形方法为弹簧网格变形方法[14-16]㊂为了验证求解方法的正确性,首先计算了标准算例NACA0012翼型强迫俯仰运动的非定常气动力变化情况,将计算结果与文献中的计算结果和实验结果对比,对比结果见文献[13]㊂俯仰运动的运动规律可以描述为[15]:α(t)=α0+A sin(ωt)=α0+A sin(2πft)(1)式中:α0是初始位置处的迎角;A是简谐振动的振幅;ω是简谐振动的圆频率;f是简谐振动的频率㊂本文定义减缩频率为:k=ωC2Vɕ(2)式中C是翼型的弦长㊂在本文中,强迫运动时自由来流的马赫数为0.755,翼型弦长为1.0m,强迫运动的减缩频率为0.0814㊂俯仰运动的初始迎角为0.016ʎ,俯仰振幅为2.51ʎ㊂图1(a)是强迫俯仰运动时的升力系数和关于1/4弦点的俯仰力矩系数随时间的变化曲线,图中计算了3个周期的气动力,由图可知,在第1个计算周期的初始阶段,计算的结果收敛性较差,这主要是由于定常计算的步数不足㊂从第2个周期开始,力和力矩系数已经达到了较好的谐振性,可以认为计算结果已经收敛㊂因此,为了减小计算量,本文的所有强迫运动过程都只计算了3个运动周期㊂图1(b)是翼型强迫沉浮运动时的力和力矩系数变化情况,其运动规律为:z(t)=z0+z m sin(ωt)(3)式中:z0=0是初始位置处的纵向位移;z m=0.1m 是沉浮运动的振幅㊂考虑洗流影响,在沉浮运动的任一时刻,瞬时迎角为:α(t)=α0-ωz m cos(ωt)/Vɕ(4)㊀㊀图1(c)是翼型强迫俯仰/沉浮耦合运动时的力和力矩系数变化情况,其运动规律为式(1)和式(3)叠加㊂对比图1可知,虽然耦合运动形式是俯仰和沉浮运动的叠加,但是耦合运动的气动力和力矩并不等于俯仰运动和沉浮运动的简单叠加,这也说明了翼型强迫运动时气动力的非线性迟滞特性比较复杂,并不是简单的线性叠加关系㊂图1㊀不同运动过程升力和力矩系数随时间变化Fig.1㊀History of lift/moment coefficients in different mo-tions㊃4861㊃第11期张庆,等:NACA0012翼型跨声速强迫运动非定常气动力模型1.1㊀一阶气动模型根据Etkin 气动力模型[2-3],强迫运动过程中的非定常气动力可以表示为:ΔC j =C j -C j 0=C jαΔα+C j ̇αC 2V ɕ()Δ̇α+C jq C 2V ɕ()Δq (5)式中C j 0是平衡位置处的力系数或是力矩系数㊂由于式(5)中C j ̇α和C jq 的量纲相同,都是空气动力系数对角度随时间一阶变化率的导数,所以在本文中称式(5)为一阶气动力模型㊂根据强迫俯仰运动时运动规律可知:̇α(t )=q =Aωcos(ωt )(6)㊀㊀俯仰运动的非定常气动力可以表示为:ΔC j =A sin(ωt )㊃C jα+kA (cos(ωt )-1)(C j ̇α+C jq )(7)㊀㊀此处需要注意,由于ΔC j 是相对于初始位置的变化量,因此右侧是(cos(ωt )-1)而不是cos(ωt )㊂所以:C j ̇α+C jq =ʏT n +1T nΔC j cos(ωt )d t()kAπω()(8)㊀㊀根据强迫沉浮运动时运动规律可知:̇α(t )=ω2z m sin(ωt )/V ɕ(9)㊀㊀沉浮运动的非定常气动力可以表示为:ΔC j =C jα-ωz m cos(ωt )V ɕ+C j ̇αω2z m V ɕC2V ɕ()sin(ωt )(10)㊀㊀所以:C j ̇α=ʏT n +1T nΔC jsin(ωt )d t ()ω2z m πC2V 2ɕ()(11)㊀㊀将式(8)和式(11)相减就可以得到单独的C j ̇α和C jq ,根据图1(a)和图1(b)的结果,可以求得一阶气动力模型的各个导数值,升力系数对̇α和q 的导数值分别为-38.9457和6.6747,力矩系数对̇α和q 的导数值分别为-2.0595和-1.3344㊂1.2㊀二阶气动模型根据Etkin 气动力模型[2-3],非定常气动力可以表示为:ΔC j =C jαΔα+C j ̇αD 2V ɕ()Δ̇α+C jα㊃㊃D 2V ɕ()2Δα㊃㊃+C jqD 2V ɕ()Δq +C j ̇qD 2Vɕ()2Δ̇q (12)㊀㊀由于式中C jα㊃㊃和C j ̇q 的量纲相同,都是空气动力系数对角度随时间二阶变化率的导数,所以在本文中称式(12)为二阶气动力模型㊂与1.1节类似,由俯仰运动的气动力变化规律可以得到:C j ̇q=C jα/k 2-ʏT n +1T n ΔC j sin(ωt )d t ()k 2A πω()C j ̇α+C jq =ʏT n +1TnΔC j cos(ωt )d t ()kA πω()ìîíïïïïïï(13)㊀㊀由沉浮运动的气动力变化规律可以得到:C j ̇α=ʏT n +1T nΔC jsin(ωt )d t ()ω2z m πC 2V ɕ2()C jα㊃㊃=ʏT n +1T nΔC j cos(ωt )d t()(z m πV ɕ()+C jα)/k 2ìîíïïïïïï(14)㊀㊀式(13)减去式(14)就可以得到单独的C j ̇α和C jq ,根据图1(a)和图1(b)的结果,可以求得一阶和二阶气动力模型的各个气动导数值,一阶导数与上节完全相同,升力系数对α㊃㊃和̇q的导数值分别为578.5118和-38.2752,力矩系数对α㊃㊃和̇q的导数值分别为18.8996和12.9044㊂式(10)㊁(11)与式(12)~(14)相比,一阶气动导数完全一样,这也间接说明传统上忽略高阶导数的做法对低阶气动导数的求解结果并没有影响,这是传统上普遍采用Etkin 气动力模型进行小迎角㊁小扰动飞行包线范围内动态稳定性分析的重要原因㊂2㊀气动力建模结果比较为了定量考察这些气动力模型对强迫运动过程非定常迟滞效应模拟的适用程度,本节对比了这些气动力模型的计算结果与直接采用CFD 进行计算得到的结果㊂图2分别是采用一阶和二阶Etkin 气动力模型计算得到的强迫俯仰运动㊁强迫沉浮运动以及耦合运动的气动力与采用CFD 方法得到的气动力迟滞曲线的对比图㊂由图2(a)可知,对于强迫俯仰运动,采用二阶气动导数得到的升力系数与CFD 计算值完全重合,而采用一阶气动导数得到的升力系数误差随着迎角的增加而增大,在最大迎角位置比CFD 计算值大50%㊂对于俯仰力矩系数,一阶模型的误差较大,而二阶模型的结果与CFD 计算值虽然不像升力系数那样完全重合,但是吻合程度也较好㊂表1是俯仰运动不同位置处的不同变量对该时刻非定常气动力的贡献情况,需要注意的是,强迫俯仰运动时̇α与q 的数值相等,α㊃㊃与̇q的数值相等㊂由表1可知,在俯仰运动1/4周期时,到达抬头最大位置处,此时α㊃㊃对应的非定常升力和力矩分别占总的非定常升力和力矩的-37.98%和-45.72%,̇q对应的非定常升力和力矩分别占总的非定常升力和力矩㊃5861㊃哈㊀尔㊀滨㊀工㊀程㊀大㊀学㊀学㊀报第41卷的2.51%和-31.21%㊂由于忽略了α㊃㊃与̇q的影响,所以一阶气动力模型对应的升力和力矩系数与CFD 计算值差别较大,见图2(a)㊂在俯仰运动3/4周期时,到达低头最大位置处㊂此时α㊃㊃与̇q这2项对非定常特性的贡献也比较大,所以与1/4周期时的情况类似,一阶模型的计算结果误差较大㊂而在1/2周期时,翼型低头经过初始位置,由于Δα㊁Δα㊃㊃以及̇q的数值刚好为0,此时非定常气动力主要由̇α和q 项产生,所以一阶气动力模型就能较为准确地重现出非定常气动力㊂在一个周期时,翼型抬头经过初始位置,此时的非定常气动力贡献情况与1/2周期时相反,此时由于Δ̇α和Δq 的数值刚好为0,此时非定常气动力主要由α㊁α㊃㊃以及̇q项产生㊂这些分析结果与图2(a)的结论一致,说明在俯仰运动过程中,二阶模型才能更准确地再现出非定常气动力和力矩㊂图2㊀不同运动过程升力和力矩系数迟滞曲线Fig.2㊀Comparison of lift /moment coefficients indifferent motions表1㊀俯仰运动不同位置非定常气动力分布情况Table 1㊀Percentage distributions at different time in thepitching motion %周期α̇αα㊃㊃q ̇q 0.25T 0.50T 0.75T 1.00TΔC L 109.4431.41-37.98-5.38 2.51ΔC m 76.0861.20-45.7239.65-31.21ΔC L0120.680-20.680ΔC m 060.68039.320ΔC L 228.27-65.52-79.2211.23 5.24ΔC m -74.8160.1844.9538.9930.69ΔC L 147.950-51.350 3.40ΔC m -8935.8505369.583666.27㊀㊀图2(b)是采用一阶和二阶Etkin 气动力模型计算得到的强迫沉浮运动的气动力与采用CFD 方法得到的气动力迟滞曲线的对比图㊂由图2(b)可知,对于强迫沉浮运动,采用二阶气动导数得到的升力系数与CFD 计算值几乎重合,而采用一阶气动导数得到的升力系数误差较大,在最大纵向位移位置比CFD 计算值大90%㊂对于俯仰力矩系数,一阶模型的误差较大,而二阶模型的结果与CFD 计算值虽然不像升力系数那样完全重合,但是吻合程度也较好㊂表2是沉浮运动不同位置处迎角的各阶导数对非定常气动力的贡献情况,由于强迫沉浮运动时没有俯仰角速度,所以在表2中没有出现q 和̇q 项对应的非定常气动力㊂由表2可知,在沉浮运动1/4周期时,到达纵向最大位置处,此时α㊃㊃对应的非定常升力和力矩分别占总的非定常升力和力矩的-94.84%和148.24%㊂由于忽略了α㊃㊃的影响,所以一阶气动力模型对应的升力和力矩系数与CFD 计算值差别较大,见图2(b)㊂在沉浮运动3/4周期时,到达纵向最小位置处㊂此时α㊃㊃对非定常特性的贡献也比较大,其数值分别为-36.92%和-49.93%,所以与1/4周期时的情况类似,一阶模型的计算结果误差较大㊂而在1/2周期时,翼型下沉经过初始位置,由于Δ̇α的数值刚好为0,此时非定常气动力主要由α和α㊃㊃产生㊂在一个周期时,翼型回到初始位置,此时的非定常气动力贡献情况与1/2周期时刚好相反,由于Δα和Δα㊃㊃的数值为0,所以此时非定常气动力主要由Δ̇α产生㊂这与图2(b)的结论一致,说明在沉浮运动过程中,二阶模型才能更准确地再现出非定常气动力和力矩㊂图2(c)是采用一阶和二阶Etkin 气动力模型计算得到的强迫俯仰/沉浮耦合运动的气动力与采用CFD 方法得到的气动力迟滞曲线的对比图㊂由图2(c)可知,对于强迫耦合运动,采用二阶气动导数得到的升力系数与CFD 计算值完全重合,而采用一阶气动导数得到的升力系数误差较大,在最大纵向位移位置比CFD 计算值大145%㊂对于俯仰力矩系数,一阶模型的误差较大,而二阶模型的结果与CFD 计算值虽然不像升力系数那样完全重合,但是㊃6861㊃第11期张庆,等:NACA0012翼型跨声速强迫运动非定常气动力模型吻合程度也较好㊂表2㊀沉浮运动不同位置非定常气动力分布情况Table 2㊀Percentage distributions at different time in theplunging motion %周期α̇αα㊃㊃0.25T 0.50T 0.75T 1.00TΔC L273.27-78.43-94.84ΔC m -246.70198.46148.24ΔC L153.150-53.15ΔC m 250.570-150.57ΔC L 106.3930.53-36.92ΔC m 83.0966.84-49.93ΔC L 0100.000ΔC m 0100.000㊀㊀表3是耦合运动不同位置处各导数对非定常气动力的贡献情况㊂由表3可知,在耦合运动1/4周期时,到达纵向最大位置处,此时也处于抬头的最大位置,此时α㊃㊃对应的非定常升力和力矩分别占总的非定常升力和力矩的-45.35%和-70.82%,̇q 对应的非定常升力和力矩分别占总的非定常升力和力矩的2.19%和-35.25%㊂由于忽略了α㊃㊃与̇q 这2项的影响,所以一阶气动力模型对应的升力和力矩系数与CFD 计算值差别较大,见图2(c)㊂在耦合运动3/4周期时,到达纵向最小位置处,此时也处于低头的最大位置㊂此时α㊃㊃对非定常特性的贡献也比较大,其数值分别为-245.66%和21.26%,̇q 对非定常特性的贡献分别为25.87%和23.00%,所以与1/4周期时的情况类似,一阶模型的计算结果误差较大㊂而在1/2周期时,翼型下沉经过初始位置,由于Δ̇q的数值刚好为0,对非定常气动力没有贡献㊂在一个周期时,翼型回到初始位置,此时Δq 的数值刚好为0,对非定常气动力没有贡献㊂虽然耦合运动气动力并不是单独运动的简单叠加,但是通过对经典Etkin 气动力模型的二阶延拓,能准确地再现出沉浮/俯仰耦合运动过程的非定常气动特性㊂表3㊀耦合运动不同位置非定常气动力分布情况Table 3㊀Percentage distributions at different time in thecoupled pitching /plunging motion %周期α̇αα㊃㊃q̇q0.25T 0.50T 0.75T 1.00TΔC L 130.6617.18-45.35-4.68 2.19ΔC m 117.8643.43-70.8244.78-35.25ΔC L77.3459.74-26.84-10.240ΔC m 25.2154.58-15.1535.360ΔC L 707.86-443.48-245.6655.4125.87ΔC m -35.2261.8521.1629.2123.00ΔC L 172.08-16.19-59.930 4.04ΔC m -349.6492.22210.81146.613㊀结论1)不论是强迫俯仰运动㊁沉浮运动,还是俯仰/沉浮耦合运动,将气动导数拓展至迎角和俯仰角的二阶导数,都可以十分精确地重现出强迫运动过程中的非定常升力变化情况㊂2)由于俯仰力矩的迟滞曲线并不是简单的椭圆形,二阶模型计算出的强迫运动过程的俯仰力矩与CFD 计算值的吻合程度不像升力那么好,说明俯仰力矩的模型要比升力更加复杂㊂3)俯仰/沉浮耦合运动的非定常气动力并不是俯仰运动和沉浮运动的简单叠加,说明精确的气动力建模还需要深入考虑其他变量的影响㊂本文的研究结果表明,Etkin 气动力模型对于非线性较强的气动力建模仍然具有较好的适用性,但是,对于三维流动以及接近失速迎角情况下的非定常气动力的建模,需要更加深入地讨论马赫数㊁减缩频率㊁更高阶导数以及交叉导数在非定常气动力模型中的作用㊂参考文献:[1]杨磊,叶正寅.倾转涵道倾转过渡阶段的非定常气动力[J].航空动力学报,2015,30(1):155-163.YANG Lei,YE Zhengyin.Unsteady aerodynamic force oftilt ducted fan during transition period[J].Journal of aero-space power,2015,30(1):155-163.[2]张庆,叶正寅.基于气动导数的类X -37B 飞行器纵向稳定性分析[J].北京航空航天大学学报,2020,46(1):77-85.ZHANG Qing,YE Zhengyin.Longitudinal stability analysis for X -37B like trans-atmospheric orbital test vehicle based on aerodynamic derivatives[J].Journal of Beijing Universi-ty of Aeronautics and Astronautics,2020,46(1):77-85.[3]ETKIN B.Dynamics of atmospheric flight[M].Mineola:Dover Publications,2012.[4]袁先旭,陈琦,谢昱飞,等.动导数数值预测中的相关问题[J].航空学报,2016,37(8):2385-2394.YUAN Xianxu,CHEN Qi,XIE Yufei,et al.Problems innumerical prediction of dynamic stability derivatives [J].Acta aeronautica et astronautica sinica,2016,37(8):2385-2394.[5]和争春,钱炜祺,朱国林,等.飞行器跨声速区俯仰力矩系数建模方法研究[J].空气动力学学报,2005,23(4):470-475.HE Zhengchun,QIAN Weiqi,ZHU Guolin,et al.Re-search on the modeling of pitching moment coefficient in transonic condition for flight vehicle[J].Acta aerodynami-ca sinica,2005,23(4):470-475.[6]GHOREYSHI M,CUMMINGS R M.Challenges in the aer-odynamics modeling of an oscillating and translating airfoil at large incidence angles[J].Aerospace science and tech-nology,2013,28(1):176-190.[7]NELSON R C,PELLETIER A.The unsteady aerodynamicsof slender wings and aircraft undergoing large amplitude㊃7861㊃哈㊀尔㊀滨㊀工㊀程㊀大㊀学㊀学㊀报第41卷maneuvers[J].Progress in aerospace sciences,2003,39 (2/3):185-248.[8]MARZOCCA P,LIBRESCU L,SILVA W A.Aeroelastic response of nonlinear wing sections using a functional series technique[J].AIAA journal,2002,40(5):813-824.[9]SILVA W A.Application of nonlinear systems theory to transonic unsteady aerodynamic responses[J].Journal of aircraft,1993,30(5):660-668.[10]SILVA W A.Simultaneous excitation of multiple-input/multiple-output CFD-based unsteady aerodynamic systems [J].Journal of aircraft,2008,45(4):1267-1274.[11]ZHANG Qing,HUA Ruhao,YE Zhengyin.Experimentaland computational investigation of novel vertical tail buffet suppression method for high sweep delta wing[J].Science China technological sciences,2015,58(1):147-157.[12]ZHANG Qing,YE Zhengyin.Novel method based on in-flatable bump for vertical tail buffeting suppression[J].Journal of aircraft,2015,52(1):367-371. [13]张庆.高速再入飞行器动力学问题研究[D].西安:西北工业大学,2018.ZHANG Qing.Research on flight dynamics of high-veloci-ty reentry vehicles[D].Xiᶄan:Northwestern Polytechni-cal University,2018.[14]RAUSCH R D,BATINA J T,YANG H T Y.Spatial ad-aptation of unstructured meshes for unsteady aerodynamic flow computations[J].AIAA journal,1992,30(5): 1243-1251.[15]BATINA J T.Unsteady Euler airfoil solutions usingunstructured dynamic meshes[J].AIAA journal,1990, 28(8):1381-1388.[16]KANDIL O A,CHUANG H A.Unsteady transonic airfoilcomputation using implicit Euler scheme on body-fixed grid[J].AIAA journal,1989,27(8):1031-1037.本文引用格式:张庆,叶正寅.NACA0012翼型跨声速强迫运动非定常气动力模型[J].哈尔滨工程大学学报,2020,41(11):1683-1688. ZHANG Qing,YE Zhengyin.Unsteady aerodynamic model of NACA0012associated with forced oscillations and translations in transonic flight[J].Jour-nal of Harbin Engineering University,2020,41(11):1683-1688.㊃8861㊃。
燃气轮机高压涡轮动静叶非定常干涉数值研究
计算 的边界条件都取 为设计 工况下 的边 界条件 , 给定静
叶进 口总温沿 径 向的变 化 曲线 , 口总 压21000 a 轴 向 进 0 0 P ,
为这部分正好是静 叶尾 缘激波 的压 力面 分支在 吸力 面上 的 反射位置 , 之激波后 流动 为亚音 速 流动 , 以该位 置的 流 加 所
燃气轮机是高技术的集成 , 是清 洁发电和动力 系统的关 键设备… 。燃 气轮机能够带 动多学科 和领 域研发 、 计和制 设
造能力的提升与发 展, 被誉 为制造业 “ 皇冠上的明珠 ” 。燃气
电站 以其高效 、 环保 的特点 , 迅速扩展 到国 内各大城市 , 成为 传统火电机组在市 区内的唯一替代品 。
跨声 涡轮 高压 级中 , 具有跨声速 流场的静叶 型面静压 波动范
式 中, Ⅳ为计算所取 的叶片通 道数 , 为栅距 , 下标 1 2分别 、
表示上下游 叶栅 。本文对 所采 用 的叶片 数据 动静 叶片数 进 行了调整 , 并保 持节 弦 比不 变 , 选取 静 叶与动 叶片数 目之 比
O
8
2 计算结果分析
在非定常条件下 , 涡轮中尾迹和势流周 期性 的干扰引起 叶片型面压力周期性 的变化 , 如图 2所示 为静叶及其 下游 的
采用 4阶 R n e—K t ug ut 时 间推进 迭代 求解 , a法 全多 重 网格 法、 局部时 间步法 和隐式 残差 光顺加 速 收敛 , 流模型 采用 湍
期 内每个 物理时 间步 的参 数进 行分 析 即可 。由于计算 中动 静 叶采用 了成 比例的叶片数 , 因此 静叶 的一 个周期 内对应动
部 的波动幅度较其它部分要强很 多, 明下 游动 叶顶 部区域 说
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
跨声速压气机非定常流动及损失的数值研究压气机性能的进一步的提高很大程度上依赖于我们对其内部非定常流动机理及相应损失机制的认识。
数值模拟的方法在这方面发挥着日益重要的作用,但是还面临很多问题。
本文将在计算压气机内部复杂非定常流动的数值方法和合理衡量压气机的损失等方面开展工作,并基于这些工作,对跨声速压气机级Rotor35/Stator37内的损失进行分析。
在数值方法方面,构造了以下三个不同程度上反应多维效应的数值格式,(1)将精度和稳健性都已得到广泛认可的二维旋转迎风格式推广为三维。
(2)通过在计算数值通量时近似求解完整的流体力学方程来体现多维效应的影响,本文构造了多维输运格式。
对于粘性流动,该格式无粘通量的计算能够考虑粘性通量的影响。
(3)基于线化多维欧拉方程的次特征线理论构造了考虑无限多方向影响的局部演化迦辽金格式。
这三种方法的共同特点是相对于已有的考虑多维效应的格式更加简单、高效,且可以直接推广到非结构网格的计算。
本文将三维旋转迎风格式格式应用到自行开发的多级叶轮机械模拟程序中。
数值验证表明多维格式在模拟跨声速压气机内复杂流动方面有一定优势。
在衡量压气机内损失方面,本文基于热力学定律,在有用功分析的框架提出了一种相比现有损失评估方法更适于非定常流场局部损失分析的损失评估参数-损失功,其强度量-损失功强度便于用来将流场结构与损失联系起来。
这两个参数物理意义明确,能够被用来方便的、定量的研究局部区域内由于不可逆过程导致的熵产引起的损失。
相比于先前同类方法,这两个参数能够考虑具体流动结构对于局部熵增的影响,从而能够更准确的衡量局部自身熵产的大小。
这两个参数为压气机的设计、优化和改型提供了一个有效工具。