CN112943668B - 航空轴流压气机复杂进气畸变下动态失速过程预测方法 - Google Patents
航空轴流压气机复杂进气畸变下动态失速过程预测方法 Download PDFInfo
- Publication number
- CN112943668B CN112943668B CN202110208288.1A CN202110208288A CN112943668B CN 112943668 B CN112943668 B CN 112943668B CN 202110208288 A CN202110208288 A CN 202110208288A CN 112943668 B CN112943668 B CN 112943668B
- Authority
- CN
- China
- Prior art keywords
- blade
- force
- circumferential
- compressor
- inlet
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 230000008569 process Effects 0.000 title claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 230000004044 response Effects 0.000 claims description 9
- 230000003068 static effect Effects 0.000 claims description 9
- 238000004445 quantitative analysis Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 5
- 238000013459 approach Methods 0.000 claims description 4
- 238000005316 response function Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 4
- 238000013461 design Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D27/00—Control, e.g. regulation, of pumps, pumping installations or pumping systems specially adapted for elastic fluids
- F04D27/001—Testing thereof; Determination or simulation of flow characteristics; Stall or surge detection, e.g. condition monitoring
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明公开了航空轴流压气机复杂进气畸变下动态失速过程预测方法,该方法通过分布式力源描述叶片固壁效应,基于损失系数及落后角概念有效构建局部叶片基元流道力源与气流参数的动态关联,具备定量描述多级轴流压气机复杂进气畸变下动态失速演化过程的能力。较传统二维计算方法和高阶三维CFD计算方法,该方法可以在兼顾计算效率的同时,有效反映包含旋流及径向/周向空间组合畸变特征在内的复杂进气畸变扰动下航空轴流压气机气动失稳特征,为先进航空发动机高保真度气动稳定性设计提供技术支撑。
Description
技术领域
本发明属于航空发动机数值仿真技术领域,具体为航空轴流压气机复杂进气畸变下动态失速过程预测方法。
背景技术
压气机是构成燃气涡轮发动机压缩系统的基本部件,其内部粘性气体的逆压力梯度决定了其具有气动失稳的特征,直接影响着发动机的工作可靠性。随着附面层摄入型发动机布局、大曲率S弯进气道等先进技术在现代航空飞行器推进系统上的逐步应用,使得发动机进口流场畸变愈加呈现以旋流及径向/周向空间组合畸变为特征的三维非定常性,其对压气机气动性能及稳定裕度的负面影响需要引起足够重视。发展压气机复杂进气畸变下高效可靠的流动稳定性分析模型是研制高性能、高稳定性压气机的重要技术保障。
稳定性分析模型的一种典型建模思路是在CFD的时间推进技术框架内,聚焦压气机内主要尺度流动特征的模化,对叶片效应进行合理降维刻画,降低计算量,使其在满足目前工程应用需求的同时,还能有效反映压气机流动失稳的整个非线性演化过程的典型特征。基于该思想,研究者先后建立了激盘模型及彻体力模型。激盘模型将叶排处理成间断面,对叶片区内流动过程不做求解。与之相比,彻体力模型的理论构成更为完善,其通过分布式力源刻画叶片固壁边界效应,能够有效模拟叶片通道内的主要流动特征。国内外学者开发了各类彻体力模型用于分析均匀进气下航空轴流压气机气动失稳起始过程,但是彻体力模型在进气畸变尤其是空间分布特征更为多样性的复杂畸变问题的应用潜力仍待深入挖掘。Chima曾基于彻体力模型的概念初步构建了一个用于分析压气机进气畸变下动态失速过程的三维非定常数值模型,但受限于其源项建模方法的局限性仅能模拟单级轴流压气机单纯周向总压畸变进气问题(Chima R V.A Three-Dimensional Unsteady CFD Modelof Compressor Stability[R].ASME Paper,GT2006-90040.)。总体而言,彻体力模型是当前航空叶轮机械内流稳定性问题研究的一个有力技术工具。显然,一个完备的叶片力源构建方法是此类模型可靠性的前提保障,尤其是需要对叶片力源与当地气动/几何参数的全面有效关联进行细致考虑,这是有效量化复杂进气畸变影响的理论基础。
发明内容
本发明的主要目的在于:基于彻体力模型建立用于高效分析多级航空轴流压气机复杂进气畸变下的动态失速过程预测方法,该方法在兼顾计算效率的同时,有效反映包含旋流、径向以及周向空间组合畸变特征在内的复杂进气畸变扰动下航空轴流压气机气动失稳特征,为先进航空发动机高保真度气动稳定性设计提供技术支撑。
本明的航空轴流压气机复杂进气畸变下动态失速过程预测方法,包括如下步骤:
步骤1,提取轴流压气机叶片中弧面几何坐标、叶型几何参数及子午流道数据,构建轴流压气机三维计算域,所述的三维计算域包括轴流压气机的叶片区域和无叶区域,且以叶片前尾缘型线区分叶片区域与无叶区域;所述叶型几何参数包括栅距g、弦长c、弯角ζ、最大厚度a和叶片几何角κ,子午流道数据包括轮毂及轮缘型线几何坐标;
步骤2,以三维彻体力模型为基础,针对叶片区域构建含分布式力源的三维非定常可压缩流动控制方程,所述分布式力源是指瞬时叶片力F,具体包括叶片区域中单位体积气体受到的轴向瞬时叶片力Fz、单位体积气体受到的周向瞬时叶片力Fθ及单位体积气体受到的径向瞬时叶片力Fr;
步骤3,构建准稳态叶片力源;
步骤4,通过一阶滞后响应函数对准稳态叶片力进行非定常响应修正获得瞬时叶片力,如下式:
其中,F为瞬时叶片力,t为时间,Ω为转子叶片旋转角速度,θ为周向坐标系,Fss为准稳态叶片力,τ=cm/vm为滞后响应时间常数,cm为叶型子午弦长,vm为子午速度。
步骤5,采用零维节流阀模型确定压气机动态失速过程中的出口压力边界。
步骤6,在计算域进口给定总压、总温及进气方向的径向-周向二维空间分布高自由度描述压气机复杂进气畸变条件;通过调节阀门系数Kt进行数值节流逼近压气机失稳边界,基于时间推进求解获得航空轴流压气机动态失速的时空演化过程。
进一步的,步骤2中,所述含分布式力源的三维非定常可压缩流动控制方程为:
其中,b表示叶型堵塞系数,ρ表示气流密度,vz、vθ和vr分别为轴向、周向及径向气流速度,e*为单位质量气流总能量,p为气流静压,h*为单位质量气流总焓,Ω为转子叶片旋转角速度,z、θ及r分别表示轴向、周向及径向坐标,Fz、Fθ及Fr即为单位体积气体受到的轴向、周向及径向瞬时叶片力。
叶型堵塞系数b定义为:
b=N|θp-θs|/2π
其中N表示叶片数目,θp及θs分别表示叶片压力面和吸力面周向坐标。
其中,ρ为气流密度,T为气流静温,vm为子午速度,△s为叶片基元流道进出口熵增,△m为叶片基元流道长度,γ为气体绝热指数,Ma1、T1 *及分别为转子进口相对马赫数、进口相对总温及出口相对总温或静子进口绝对马赫数、进口绝对总温及出口绝对总温,vrel是相对速度,指压气机相对转动坐标系下的合速度,损失力f平行且反向于相对速度vrel;
其中,Fθ表示单位体积气体受到的周向瞬时叶片力,fθ表示周向损失力,
其中,Ω为转子叶片旋转角速度,vm为子午速度,κ为叶片几何角,δ为落后角,r表示叶片区域径向坐标。
进一步的,步骤5中采用零维节流阀模型确定压气机动态失速过程中的出口压力边界,具体为:
采用基于二次抛物线解析式的零维节流阀模型确定压气机动态失速过程中的出口压力边界:
有益效果:较传统二维计算方法和高阶三维CFD计算方法,该方法可以在兼顾计算效率的同时,有效反映包含旋流、径向/周向空间组合畸变特征在内的复杂进气畸变扰动下航空轴流压气机气动失稳特征,为先进航空发动机高保真度气动稳定性设计提供技术支撑。
附图说明
图1为本发明提供的计算模型理论示意图;
图2为两级轴流压气机径向-周向空间组合总压畸变下模型预测压升特性;
图3为两级轴流压气机径向-周向空间组合总压畸变下预失速小扰动时空演化特征;
图4为两级轴流压气机径向-周向空间组合总压畸变下旋转失速演化过程。
具体实施方式
下面结合附图对模型技术方法做进一步解释。
本明提出的航空轴流压气机复杂进气畸变下动态失速过程预测方法,包括如下步骤:
步骤1,提取轴流压气机叶片中弧面几何坐标、叶型几何参数及子午流道数据,构建轴流压气机三维计算域;所述轴流压气机为多级轴流压气机。所述的三维计算域包括轴流压气机的叶片区域和无叶区域,且以叶片前尾缘型线区分叶片区域与无叶区域,所述叶型几何参数包括栅距g、弦长c、弯角ζ、最大厚度a和叶片几何角κ,子午流道数据包括轮毂及轮缘型线几何坐标。如图1所示。
步骤2,以三维彻体力模型为基础,针对叶片区域构建含分布式力源的三维非定常可压缩流动控制方程,分布式力源用于加载压气机各叶排对气流产生的压升、转向、做功及损失作用,可避免因直接考虑叶片三维几何造型导致的复杂网格划分工作,允许以相对规整的粗网格快速捕获叶轮机械内三维流动的大尺度特征。
所述分布式力源是指瞬时叶片力F,具体包括叶片区域中单位体积气体受到的轴向瞬时叶片力Fz、单位体积气体受到的周向瞬时叶片力Fθ及单位体积气体受到的径向瞬时叶片力Fr。
具体的,叶片区域中含分布式力源的三维非定常可压缩流动控制方程为:
其中,b表示叶型堵塞系数,ρ表示气流密度,vz、vθ和vr分别为轴向、周向及径向气流速度,e*为单位质量气流总能量,p为气流静压,h*为单位质量气流总焓,Ω为转子叶片旋转角速度,z、θ及r分别表示轴向、周向及径向坐标,Fz、Fθ及Fr即为单位体积气体受到的轴向、周向及径向瞬时叶片力。
叶型堵塞系数b定义为:
b=N|θp-θs|/2π
其中N表示叶片数目,θp及θs分别表示叶片压力面和吸力面周向坐标。
步骤3,构建准稳态叶片力源;
步骤3.1,构建损失系数及落后角δ与基元流道进口攻角i、马赫数Ma以及叶型几何参数的定量解析关系,间接完成准稳态叶片力源与流动参数的全面关联,所述流动参数包括基元流道进口攻角i和马赫数Ma,所述叶型几何参包括栅距g、弦长c、弯角ζ、最大厚度a、叶片几何角κ。
其中,ρ为气流密度,T为气流静温,vm为子午速度,△s为叶片基元流道进出口熵增,△m为叶片基元流道长度,γ为气体绝热指数,Ma1、T1 *及分别为转子进口相对马赫数、进口相对总温及出口相对总温或静子进口绝对马赫数、进口绝对总温及出口绝对总温,vrel是相对速度,指压气机相对转动坐标系下的合速度,损失力f平行且反向于相对速度vrel。
其中,Fθ表示单位体积气体受到的周向瞬时叶片力,fθ表示周向损失力,
其中,Ω为转子叶片旋转角速度,vm为子午速度,κ为叶片几何角,δ为落后角,r表示叶片区域径向坐标,
步骤4,为有效量化压气机动态失速过程中叶片力源对于上游扰动的动态滞后响应,通过一阶滞后响应函数对准稳态叶片力进行非定常响应修正获得瞬时叶片力F,所述上游扰动是指由压气机工作的初始条件或边界条件改变产生的流场中的扰动,对准稳态叶片力进行修正获得瞬时叶片力F如下式:
其中,F为瞬时叶片力,t为时间,Ω为转子叶片旋转角速度,θ为周向坐标系,Fss为准稳态叶片力,τ=cm/vm为滞后响应时间常数,cm为叶型子午弦长,vm为子午速度。
步骤5,采用零维节流阀模型确定压气机动态失速过程中的出口压力边界。
采用基于二次抛物线解析式的零维节流阀模型确定压气机动态失速过程中的出口压力边界:
步骤6,在计算域进口给定总压、总温及进气方向的径向-周向二维空间分布高自由度描述压气机复杂进气畸变条件。
通过调节阀门系数Kt进行数值节流逼近压气机失稳边界,基于时间推进求解获得航空轴流压气机动态失速的时空演化过程。
图2给出了基于本发明在一台普通台式机上仅耗时1.5天获得的一台两级轴流压气机径向-周向空间组合复杂总压畸变下的压升特性,数值结果有效体现了来流畸变导致压气机峰值压升系数恶化、失稳流量提前的现象;图3给出了基于本发明获得的该两级轴流压气机复杂总压畸变下预失速小扰动时空演化特征,数值结果有效体现了复杂畸变下预失速小扰动在其周向传播过程中存在的周期性形成、发展、衰减甚至消失的动态演变特征;图4给出了基于本发明获得的两级轴流压气机复杂总压畸变下旋转失速团的演化过程,数值结果有效体现了复杂畸变下形成的完全失速团的周向传播过程。
Claims (5)
1.航空轴流压气机复杂进气畸变下动态失速过程预测方法,其特征在于,
步骤1,提取轴流压气机叶片中弧面几何坐标、叶型几何参数及子午流道数据,构建轴流压气机三维计算域,所述的三维计算域包括轴流压气机的叶片区域和无叶区域,所述叶型几何参数包括栅距g、弦长c、弯角ζ、最大厚度a和叶片几何角κ,子午流道数据包括轮毂及轮缘型线几何坐标;
步骤2,以三维彻体力模型为基础,针对叶片区域构建含分布式力源的三维非定常可压缩流动控制方程,所述分布式力源是指瞬时叶片力F,具体包括叶片区域中单位体积气体受到的轴向瞬时叶片力Fz、单位体积气体受到的周向瞬时叶片力Fθ及单位体积气体受到的径向瞬时叶片力Fr;
步骤3,构建准稳态叶片力源;
步骤3.2,将准稳态叶片力Fss分解为转折力和损失力f,写为:所述损失力f用于模拟叶片附面层流动相关气流损失,由基元流道损失系数确定;所述转折力用于控制叶片区内气流运动方向,由叶片几何角叠加落后角δ确定;
步骤4,通过一阶滞后响应函数对准稳态叶片力进行非定常响应修正获得瞬时叶片力,如下式:
其中,F为瞬时叶片力,t为时间,Ω为转子叶片旋转角速度,θ为周向坐标系,Fss为准稳态叶片力,τ=cm/vm为滞后响应时间常数,cm为叶型子午弦长,vm为子午速度;
步骤5,采用零维节流阀模型确定压气机动态失速过程中的出口压力边界;
步骤6,在计算域进口给定总压、总温及进气方向的径向-周向二维空间分布高自由度描述压气机复杂进气畸变条件;通过调节阀门系数Kt进行数值节流逼近压气机失稳边界,基于时间推进求解获得航空轴流压气机动态失速的时空演化过程。
其中,ρ为气流密度,T为气流静温,vm为子午速度,△s为叶片基元流道进出口熵增,△m为叶片基元流道长度,γ为气体绝热指数,Ma1、T1 *及分别为转子进口相对马赫数、进口相对总温及出口相对总温或静子进口绝对马赫数、进口绝对总温及出口绝对总温,vrel是相对速度,指压气机相对转动坐标系下的合速度,损失力f平行且反向于相对速度vrel;
其中,Fθ表示单位体积气体受到的周向瞬时叶片力,fθ表示周向损失力,
其中,Ω为转子叶片旋转角速度,vm为子午速度,κ为叶片几何角,δ为落后角,r表示叶片区域径向坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110208288.1A CN112943668B (zh) | 2021-02-24 | 2021-02-24 | 航空轴流压气机复杂进气畸变下动态失速过程预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110208288.1A CN112943668B (zh) | 2021-02-24 | 2021-02-24 | 航空轴流压气机复杂进气畸变下动态失速过程预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112943668A CN112943668A (zh) | 2021-06-11 |
CN112943668B true CN112943668B (zh) | 2022-04-22 |
Family
ID=76246009
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110208288.1A Active CN112943668B (zh) | 2021-02-24 | 2021-02-24 | 航空轴流压气机复杂进气畸变下动态失速过程预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112943668B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113720609B (zh) * | 2021-07-07 | 2022-08-05 | 南京航空航天大学 | 基于人工神经网络的可调总压畸变发生器调节指导方法 |
CN114329826B (zh) * | 2021-12-22 | 2024-08-20 | 大连理工大学 | 一种基于生成式对抗网络的平面叶栅定常流动预测方法 |
CN114781203B (zh) * | 2022-03-17 | 2024-08-06 | 西安交通大学 | 一种用于提取叶轮机械当地损失力系数的方法及系统 |
CN114624464B (zh) * | 2022-03-23 | 2022-12-20 | 北京航空航天大学 | 非接触式热膜测试设备及压气机转子部件稳定性判别方法 |
CN114580320B (zh) * | 2022-03-23 | 2024-10-22 | 北京航空航天大学 | 非接触式热膜测试设备、压气机局部失稳判别方法及应用 |
CN117634100B (zh) * | 2024-01-25 | 2024-04-30 | 陕西空天信息技术有限公司 | 多级轴流压气机的子午流道获取方法、装置、设备及介质 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109815624B (zh) * | 2019-02-21 | 2022-01-04 | 南京航空航天大学 | 一种考虑进气总压畸变影响的压气机稳定边界判断方法 |
CN110005628B (zh) * | 2019-03-27 | 2020-03-17 | 南京航空航天大学 | 基于异位方差分析的压气机气动失稳在线辨识方法和系统 |
CN111079232A (zh) * | 2019-11-15 | 2020-04-28 | 南京航空航天大学 | 一种预测旋流畸变进气对航空发动机性能影响的计算方法 |
-
2021
- 2021-02-24 CN CN202110208288.1A patent/CN112943668B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112943668A (zh) | 2021-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112943668B (zh) | 航空轴流压气机复杂进气畸变下动态失速过程预测方法 | |
Plas et al. | Performance of a boundary layer ingesting (BLI) propulsion system | |
Liu et al. | Computational and experimental investigations of performance curve of an axial flow fan using downstream flow resistance method | |
Jang et al. | Optimization of a stator blade using response surface method in a single-stage transonic axial compressor | |
CN111079232A (zh) | 一种预测旋流畸变进气对航空发动机性能影响的计算方法 | |
Song et al. | Aerodynamic and structural multidisciplinary optimization design method of fan rotors based on blade curvature constraints | |
Halder et al. | Coupled CAD-CFD automated optimization for leading and trailing edge of an axial impulse turbine blade | |
Wilkinson | The Design of an axial flow fan for air-cooled heat exchanger applications. | |
Yang et al. | Optimization Design and Experimental Study of Low‐Pressure Axial Fan with Forward‐Skewed Blades | |
Wang et al. | Rapid performance prediction model of axial turbine with coupling one-dimensional inverse design and direct analysis | |
Li et al. | Optimization strategy for a single-stage axisymmetric hub endwall in axial compressor by a modified transonic area rule | |
Templalexis et al. | Development of a 2-D compressor streamline curvature code | |
Su et al. | Adjoint-based geometrically constrained aerodynamic optimization of a transonic compressor stage | |
Ellbrant et al. | CFD optimization of a transonic compressor using multiobjective GA and metamodels | |
Zhiping et al. | A modified streamwise body force model of fan with distorted inflow for rapid propulsion-airframe integrated simulation | |
Joly et al. | Full design of a highly loaded fan by multi-objective optimization of through-flow and high-fidelity aero-mechanical performances | |
Xiang et al. | Aerodynamic modification and optimization of intermediate pressure compressor in marine intercooled recuperated gas turbine | |
Guo et al. | A three-dimensional unsteady through-flow model for rotating stall in axial compressors | |
Zhang et al. | Aerodynamic design of a multi-stage industrial axial compressor | |
Zhang et al. | Numerical investigation and modeling of sweep effects on inlet flow field of axial compressor cascades | |
Li et al. | Development and application of a throughflow method for high-loaded axial flow compressors | |
Guo et al. | Static pressure redistribution mechanism of non-axisymmetric endwall based on radial equilibrium | |
Benini | Advances in aerodynamic design of gas turbines compressors | |
Xie et al. | Development and validation of non-axisymmetric body-force propulsor model | |
Gao et al. | Effects of stator blade camber and surface viscosity on unsteady flow in axial turbine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |