CN113283198A - 优化压气机机匣处理提高稳定裕度的方法、系统、终端 - Google Patents
优化压气机机匣处理提高稳定裕度的方法、系统、终端 Download PDFInfo
- Publication number
- CN113283198A CN113283198A CN202110649455.6A CN202110649455A CN113283198A CN 113283198 A CN113283198 A CN 113283198A CN 202110649455 A CN202110649455 A CN 202110649455A CN 113283198 A CN113283198 A CN 113283198A
- Authority
- CN
- China
- Prior art keywords
- casing
- prototype
- processing
- treatment
- optimizing
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000004364 calculation method Methods 0.000 claims abstract description 47
- 238000004088 simulation Methods 0.000 claims abstract description 10
- 230000008602 contraction Effects 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims description 129
- 238000005457 optimization Methods 0.000 claims description 27
- 230000003068 static effect Effects 0.000 claims description 25
- 238000013461 design Methods 0.000 claims description 16
- 230000006872 improvement Effects 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 9
- 238000009413 insulation Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 8
- 238000010586 diagram Methods 0.000 description 13
- 230000009471 action Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000000926 separation method Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003534 oscillatory effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/40—Casings; Connections of working fluid
- F04D29/403—Casings; Connections of working fluid especially adapted for elastic fluid pumps
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/66—Combating cavitation, whirls, noise, vibration or the like; Balancing
- F04D29/661—Combating cavitation, whirls, noise, vibration or the like; Balancing especially adapted for elastic fluid pumps
- F04D29/668—Combating cavitation, whirls, noise, vibration or the like; Balancing especially adapted for elastic fluid pumps damping or preventing mechanical vibrations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- Mechanical Engineering (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明属于船舶推进用燃气轮机的压气机技术领域,公开了一种优化压气机机匣处理提高稳定裕度的方法、系统、终端,优化压气机机匣处理提高稳定裕度的方法包括:对原机带轴向斜槽和回流腔机匣及参数改进型机匣进行数值模拟,确定机匣处理的结构形式,通过增大轴向速度,采用整个环形进气,通过调整进气段入口,扩大进气角度,增加进气量;在出气段,采用倾斜收缩喷嘴作为出气口,增大出口速度,优化排气段,计算网格采用IGG/AutoGrid5进行单通道网格的划分,优化压气机机匣处理效果。本发明能够有效优化排气段处理机匣,在低工况下对压气机工作稳定裕度改善效果最好,稳定裕度提高了6.04%。
Description
技术领域
本发明属于压气机技术领域,尤其涉及一种优化压气机机匣处理提高稳定裕度的方法、系统、终端。
背景技术
目前,压气机机匣处理的主要作用是减小失速点的流量,从而增加压气机的喘振裕度,在实际工作中,由于APU引气压力较高或出口负载变化导致压气机工作点逼近喘振边界,可能导致压气机的叶片断裂,甚至APU发动机熄火等严重故障。为了提高压气机的喘振裕度,目前航空发动机主要采用压气机级间放气、进口可调导叶等方法。采用级间放气会使压气机出口压力剧烈波动,并不适用于APU,而采用进口可调导叶结构复杂,会明显增加发动机的重量,且喘振裕度增加有限。
通过上述分析,现有技术存在的问题及缺陷为:现有压气机的机匣处理结构不能适应压气机进气量增大的需求,排气段处理机匣在低工况下喘振裕度不高。
解决以上问题及缺陷的难度为:
本发明的技术难点在于对原有机匣处理结构的改进,之后通过理论计算结合仿真实验的方法进行多次优化修正,得到效果较优的机匣处理结构。
解决以上问题及缺陷的意义为:
通过本发明,通过对现有机匣处理结构的优化改进,可以提高压气机在低工况下的喘振裕度。
发明内容
针对现有技术存在的问题,本发明提供了一种优化压气机机匣处理提高稳定裕度的方法、系统、终端。
本发明是这样实现的,一种优化压气机机匣处理提高稳定裕度的方法,所述优化压气机机匣处理提高稳定裕度的方法包括:
对原机带轴向斜槽和回流腔机匣及参数改进型机匣进行数值模拟,确定参数改进型机匣的结构形式,并通过优化排气段、网格划分优化参数改进型机匣。
所述参数改进型机匣结构:
由喘振机理出发,当冲角增加时:①在压差作用下,气流进入机匣环形腔室并流向级的进口,使叶顶轴向速度增加,抑制叶顶气流分离;②循环流动把进入动叶的主流空气挤向轮毂,使轮毂处轴向速度增加;③腔室将脉动气流能量耗散,抑制了旋转分离的发展。
由前期对原型机匣处理计算,发现排气段对轴向速度的增加量较小,由于喘振发生主要是由于叶顶处轴向速度,流量减少,叶背发生气流分离导致,因此,针对发现问题对其进行改进,在原型机匣的基础上,改进排气段,增大轴向速度,改进结构如图13所示,将环形凹槽进气结构改为整个环形进气,调整进气段入口,扩大进气角度,增加进气量;在出气段,将出口改为倾斜收缩喷嘴,增大出口速度,计算低工况时的扩稳效果。
进一步,优化排气段包括:增大轴向速度,采用整个环形进气,通过调整进气段入口,扩大进气角度,增加进气量;在出气段,采用倾斜收缩喷嘴作为出气口,增大出口速度。
进一步,所述优化网格划分包括:
计算网格采用IGG/AutoGrid5进行单通道网格的划分,将计算区域分为主流区与机匣处理区;
对主流区网格拓扑结构采用O4H型结构,壁面附近的网格沿壁面法向方向按几何级数的规律加密,近壁面第一层网格距离固壁距离取0.001mm,将y+值控制在一定的低雷诺数湍流模型要求的范围之内,边界层外的网格均匀分布;
对机匣处理部分采用ZR技术处理,将机匣进气段设置为环形进气,围绕整个叶轮进行划分。
进一步,所述优化压气机机匣处理提高稳定裕度的方法包括以下步骤:
步骤一,获取原型机匣处理几何模型,建立原型机匣处理平面数值模型;
步骤二,计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分;
步骤三,基于计算的相关数据,计算原型机匣处理的喘振裕度;
步骤四,基于计算的相关数据进行机匣处理的结构优化,并计算结构优化后的机匣处理的喘振裕度;对比改进后机匣处理的喘振裕度与原型机匣处理的喘振裕度确定优化有效性。
进一步,步骤二中,所述计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分包括:
(1)计算域:采用CFD方法对涡轮流场进行数值模拟;
(2)网格划分:计算网格采用可视化IGG/AutoGrid5进行单通道网格的划分,将整个计算区域分为主流区与机匣处理区;所述主流流道利用NUMECA中专门针对叶轮机械部分网格生成模块AutoGrid整体划分网格,主流区网格拓扑结构采用O4H型结构,近壁面第一层网格距离固壁距离取0.001mm,将y+值控制在一定的低雷诺数湍流模型要求的范围之内,壁面附近的网格沿壁面法向方向按几何级数的规律加密,网格数约为该方向网格总数的三分之一,边界层外的网格均匀分布;机匣处理区和转子区之间设置转静子交界面;
(3)设计工况计算:通过设定边界条件、确定收敛标准,对压气机三维流场进行模拟,并基于模拟数值对机匣处理的总压分布、相对马赫数分布、压比流量进行分析。
进一步,所述通过设定边界条件,对压气机三维流场进行模拟包括:
对压气机三维流场进行模拟,进口给定总温,总压和气流角,出口给定平均静压,壁面采用绝热无滑移边界条件;
所述边界条件设置包括:
1)进口边界:流动方向为轴向,设定均匀总压为99300Pa,均匀总温为300K;
2)出口边界:平均静压为430000Pa;
3)固壁:无滑移固壁,绝热;
4)转子转速:7436RPM;
5)工质:空气;
6)转/静交界面处理:主流区各转/静交界面采用Conservation Coupling bypitchwise row,机匣处理结构的转静交界面采用Non reflecting 1D。
进一步,所述收敛标准包括:
总残差降低至一定水平后不再降低;总参数随迭代步数增加不再改变,总参数包括效率、压比、输出功率和扭矩;进出口流量差异不超过0.5%;总参数随迭代步数增加呈周期性波动。
进一步,所述机匣处理的喘振裕度计算公式如下:
本发明另一目的在于提供一种优化压气机机匣处理提高稳定裕度系统,包括:
原型机匣处理平面数值模型构建模块,用于获取原型机匣处理几何模型,建立原型机匣处理平面数值模型;
原型机匣处理的网格划分模块,用于计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分;
原型机匣处理的喘振裕度计算模块,用于基于计算的相关数据,计算原型机匣处理的喘振裕度;
机匣处理的喘振裕度计算模块,用于基于计算的相关数据进行机匣处理的结构优化,并计算结构优化后的机匣处理的喘振裕度;对比改进后机匣处理的喘振裕度与原型机匣处理的喘振裕度确定优化有效性。
本发明另一目的在于提供一种接收用户输入程序存储介质,所存储的计算机程序使电子设备执行所述优化压气机机匣处理提高稳定裕度的方法,包括下列步骤:
步骤一,获取原型机匣处理几何模型,建立原型机匣处理平面数值模型;
步骤二,计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分;
步骤三,基于计算的相关数据,计算原型机匣处理的喘振裕度;
步骤四,基于计算的相关数据进行机匣处理的结构优化,并计算结构优化后的机匣处理的喘振裕度;对比改进后机匣处理的喘振裕度与原型机匣处理的喘振裕度确定优化有效性。
本发明另一目的在于提供一种信息数据处理终端,所述信息数据处理终端包括存储器和处理器,所述存储器存储有计算机程序,所述计算机程序被所述处理器执行时,使得所述处理器执行所述优化压气机机匣处理提高稳定裕度的方法。
结合上述的所有技术方案,本发明所具备的优点及积极效果为:本发明能够有效优化排气段处理机匣,在低工况下对稳定裕度改善效果最好,本发明的方案相对于原始机匣处理将喘振裕度增大了6.04%。
附图说明
图1是本发明实施例提供的优化压气机机匣处理提高稳定裕度的方法流程图。
图2是本发明实施例提供的进出口流量收敛曲线示意图。
图3是本发明实施例提供的全局残差收敛曲线示意图。横坐标为计算迭代步数(无量纲),纵坐标为全局残差(无量纲).
图4是本发明实施例提供的效率收敛曲线示意图。横坐标为计算迭代步数(无量纲),纵坐标为效率(无量纲)。
图5是本发明实施例提供的压比收敛曲线示意图。横坐标为计算迭代步数(无量纲),纵坐标为压比(无量纲)。
图6是本发明实施例提供的通流部分总压分布示意图。
图7是本发明实施例提供的机匣处理截面速度流线分布示意图。
图8是本发明实施例提供的无机匣处理叶顶截面相对马赫数分布示意图。
图9是本发明实施例提供的带机匣处理叶顶截面相对马赫数分布示意图。
图10是本发明实施例提供的工况下原型机匣处理压比-流量特性示意图。横坐标为流量(kg/s),纵坐标为压比(无量纲)。
图11是本发明实施例提供的工况原型机匣处理的效率-流量特性示意图。横坐标为流量(kg/s),纵坐标为效率(无量纲)。
图12是本发明实施例提供的压比-流量特性对比示意图。横坐标为流量(kg/s),纵坐标为压比(无量纲);Original:原始没有机匣处理的结构;Case3:原始机匣处理结构;Case4:本发明的机匣处理结构。
图13是本发明实施例提供的改进型机匣示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
针对现有技术存在的问题,本发明提供了一种优化压气机机匣处理提高稳定裕度的方法,下面结合附图对本发明作详细的描述。
本发明实施例提供的优化压气机机匣处理提高稳定裕度的方法包括:
对原机带轴向斜槽和回流腔机匣及参数改进型机匣进行数值模拟,确定机匣处理的结构形式,通过优化排气段、网格划分优化压气机机匣处理。
本发明实施例提供的优化排气段包括:增大轴向速度,采用整个环形进气,通过调整进气段入口,扩大进气角度,增加进气量;在出气段,采用倾斜收缩喷嘴作为出气口,增大出口速度。
本发明实施例提供的优化网格划分包括:
计算网格采用IGG/AutoGrid5进行单通道网格的划分,将计算区域分为主流区与机匣处理区;
对主流区网格拓扑结构采用O4H型结构,壁面附近的网格沿壁面法向方向按几何级数的规律加密,近壁面第一层网格距离固壁距离取0.001mm,将y+值控制在一定的低雷诺数湍流模型要求的范围之内,边界层外的网格均匀分布;
对机匣处理部分采用ZR技术处理,将机匣进气段设置为环形进气,围绕整个叶轮进行划分。
如图1所示,本发明实施例提供的优化压气机机匣处理提高稳定裕度的方法包括以下步骤:
S101,获取原型机匣处理几何模型,建立原型机匣处理平面数值模型;
S102,计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分;
S103,基于计算的相关数据,计算原型机匣处理的喘振裕度;
S104,基于计算的相关数据进行机匣处理的结构优化,并计算结构优化后的机匣处理的喘振裕度;对比改进后机匣处理的喘振裕度与原型机匣处理的喘振裕度确定优化有效性。
步骤S102中,本发明实施例提供的计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分包括:
(1)计算域:采用CFD方法对涡轮流场进行数值模拟;
(2)网格划分:计算网格采用可视化IGG/AutoGrid5进行单通道网格的划分,将整个计算区域分为主流区与机匣处理区;所述主流流道利用NUMECA中专门针对叶轮机械部分网格生成模块AutoGrid整体划分网格,主流区网格拓扑结构采用O4H型结构,近壁面第一层网格距离固壁距离取0.001mm,将y+值控制在一定的低雷诺数湍流模型要求的范围之内,壁面附近的网格沿壁面法向方向按几何级数的规律加密,网格数约为该方向网格总数的三分之一,边界层外的网格均匀分布;机匣处理区和转子区之间设置转静子交界面;
(3)设计工况计算:通过设定边界条件、确定收敛标准,对压气机三维流场进行模拟,并基于模拟数值对机匣处理的总压分布、相对马赫数分布、压比流量进行分析。
本发明实施例提供的通过设定边界条件,对压气机三维流场进行模拟包括:
对压气机三维流场进行模拟,进口给定总温,总压和气流角,出口给定平均静压,壁面采用绝热无滑移边界条件;
边界条件设置包括:
1)进口边界:流动方向为轴向,设定均匀总压为99300Pa,均匀总温为300K;
2)出口边界:平均静压为430000Pa;
3)固壁:无滑移固壁,绝热;
4)转子转速:7436RPM;
5)工质:空气;
6)转/静交界面处理:主流区各转/静交界面采用Conservation Coupling bypitchwise row,机匣处理结构的转静交界面采用Non reflecting 1D。
本发明实施例提供的收敛标准包括:
总残差降低至一定水平后不再降低;总参数随迭代步数增加不再改变,总参数包括效率、压比、输出功率和扭矩;进出口流量差异不超过0.5%;总参数随迭代步数增加呈周期性波动。
本发明实施例提供的机匣处理的喘振裕度计算公式如下:
本发明还提供一种优化压气机机匣处理提高稳定裕度系统,包括:
原型机匣处理平面数值模型构建模块,用于获取原型机匣处理几何模型,建立原型机匣处理平面数值模型;
原型机匣处理的网格划分模块,用于计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分;
原型机匣处理的喘振裕度计算模块,用于基于计算的相关数据,计算原型机匣处理的喘振裕度;
机匣处理的喘振裕度计算模块,用于基于计算的相关数据进行机匣处理的结构优化,并计算结构优化后的机匣处理的喘振裕度;对比改进后机匣处理的喘振裕度与原型机匣处理的喘振裕度确定优化有效性。
下面结合具体实施例对本发明的技术方案作进一步说明。
实施例1:
1、原型机匣处理数值计算
(1)原型机匣处理几何模型
原型机匣处理模型为带轴向斜槽和回流腔的结构.
(2)建立平面数值模型
原型带轴向斜槽和回流腔机匣处理结构。
2、数值计算
2.1计算域
采用CFD方法对涡轮流场进行数值模拟,为减少计算量,采用周期边界处理方法,计算域只包括一个流道,此外,为保证计算收敛,将出口位置沿Z轴(工质流向)向下游延伸1倍动叶弦长。
2.2网格划分
计算网格采用可视化IGG/AutoGrid5进行单通道网格的划分,整个计算区域分为主流区与机匣处理区。主流流道利用NUMECA中专门针对叶轮机械部分网格生成模块AutoGrid整体划分网格,主流区网格拓扑结构采用O4H型结构。在网格生成过程中,考虑粘性流场计算中低雷诺数湍流模型的应用,捕捉近壁面湍流附面层流动的详细流动特征,加密靠近壁面的网格,近壁面第一层网格距离固壁距离取0.001mm,将y+值控制在一定的低雷诺数湍流模型要求的范围之内,壁面附近(近似于边界层内)的网格沿壁面法向方向按几何级数的规律加密,网格数约为该方向网格总数的三分之一,边界层外的网格均匀分布。计算结果表明,壁面绝大部分区域y+<20,流道区域y+<10。
机匣处理结构采用Autogrid5划分结构化网格,机匣处理区和转子区之间设置转静子交界面,叶轮组三维网格总数为1800万,网格总块数为221,经过网格质量检查,网格参数均在合理取值范围内,检查结果如下:
3、设计工况计算
3.1边界条件
对压气机三维流场进行模拟,进口给定总温,总压和气流角,出口给定平均静压,壁面采用绝热无滑移边界条件,边界条件设置如下:
1)进口边界:流动方向为轴向,设定均匀总压为99300Pa,均匀总温为300K;
2)出口边界:平均静压为430000Pa;
3)固壁:无滑移固壁,绝热;
4)转子转速:7436RPM;
5)工质:空气。
转/静交界面处理:主流区各转/静交界面采用Conservation Coupling bypitchwise row,机匣处理结构的转静交界面采用Non reflecting 1D,为了保证计算的收敛性,将专家参数loccor设置为0。
3.2收敛标准
在三维流场计算中,判断算例收敛的准则为:
总残差降低至一定水平后不再降低;
总参数随迭代步数增加不再改变,总参数包括效率、压比、输出功率和扭矩等;
进出口流量差异不超过0.5%;
总参数随迭代步数增加呈周期性波动。
4、工况计算
4.1边界条件
对低工况运行时,机匣处理进行数值计算,设定边界条件如下:
进口边界:流动方向为轴向,设定均匀总压为99300Pa,均匀总温为300K;出口边界:平均静压为180000Pa;固壁:无滑移固壁,绝热;转子转速:5567RPM;工质:空气。转/静交界面处理:主流区各转/静交界面采用Conservation Coupling by pitchwise row,机匣处理结构的转静交界面采用Non reflecting 1D,为了保证计算的收敛性,将专家参数loccor设置为0。
4.2收敛曲线
本发明给出了经过8000次迭代计算后,全局残差下降到10e-4,效率和压降呈振荡性收敛,且不再发生变化,进出口流量完全匹配,因此可以判断为计算收敛。
4.3计算结果分析
(1)总压分布
低工况下,通流部分总压分布中可以看出,在有机匣处理段,在主流与机匣处理区域部分总压较大,此处产生了总压损失。
进一步分析机匣处理段流场分布,在总压损失区,主流通道与机匣回流腔内部均产生了回转涡流,相比于设计工况,主流区的回转涡流范围较大,此时产生的损失较大,通流部分轴向流速降低。
(2)相对马赫数分布
叶顶截面的相对马赫数分布中从不带机匣处理和有机匣处理,可以看出,低工况下叶顶截面流速明显降低,在机匣回流腔的作用下,带机匣处理叶顶截面相对马赫数相对增大。
其中,图2是本发明实施例提供的进出口流量收敛曲线示意图。横坐标为计算迭代步数(无量纲),纵坐标为流量kg/s。
图3是本发明实施例提供的全局残差收敛曲线示意图。横坐标为计算迭代步数(无量纲),纵坐标为全局残差(无量纲)。
图4是本发明实施例提供的效率收敛曲线示意图。横坐标为计算迭代步数(无量纲),纵坐标为效率(无量纲)。
图5是本发明实施例提供的压比收敛曲线示意图。横坐标为计算迭代步数(无量纲),纵坐标为压比(无量纲)。
图6是本发明实施例提供的通流部分总压分布示意图。
图7是本发明实施例提供的机匣处理截面速度流线分布示意图。
图8是本发明实施例提供的无机匣处理叶顶截面相对马赫数分布示意图。
图9是本发明实施例提供的带机匣处理叶顶截面相对马赫数分布示意图。
(3)低工况下压比流量,效率流量特性
图10是本发明实施例提供的工况下原型机匣处理压比-流量特性示意图。
0.2工况下,计算原型机匣处理的压比流量,效率流量特性图11所示,图12是本发明实施例提供的压比-流量特性对比示意图。
(4)喘振裕度
将计算数据带入稳定裕度改进量计算公式,对原型机匣处理的喘振裕度改进量进行计算,其中,在转速为5567RPM情况下,工作点压比为2.6,流量为50kg/s,机匣处理近失速点压比为2.74,流量为49.35kg/s,得到原型机匣处理的稳定裕度改进量为6.77%。
5、改进型机匣处理的性能
根据对原型机匣处理分析,发现:排气段对轴向速度的增加量较小。
因此,对其进行改进,改进方案:优化排气段,增大轴向速度。
5.1优化排气段
由喘振机理出发,当冲角增加时:①在压差作用下,气流进入机匣环形腔室并流向级的进口,使叶顶轴向速度增加,抑制叶顶气流分离;②循环流动把进入动叶的主流空气挤向轮毂,使轮毂处轴向速度增加;③腔室将脉动气流能量耗散,抑制了旋转分离的发展。
由前期对原型机匣处理计算,发现排气段对轴向速度的增加量较小,由于喘振发生主要是由于叶顶处轴向速度,流量减少,叶背发生气流分离导致,因此,针对发现问题对其进行改进,在原型机匣的基础上,改进排气段,增大轴向速度,改进结构如图所示,将环形凹槽进气结构改为整个环形进气,调整进气段入口,扩大进气角度,增加进气量;在出气段,将出口改为倾斜收缩喷嘴,增大出口速度,计算低工况时的扩稳效果。
具体包括:
1、网格划分
计算网格采用IGG/AutoGrid5进行单通道网格的划分,整个计算区域分为主流区与机匣处理区。主流区网格拓扑结构采用O4H型结构。壁面附近(近似于边界层内)的网格沿壁面法向方向按几何级数的规律加密,近壁面第一层网格距离固壁距离取0.001mm,将y+值控制在一定的低雷诺数湍流模型要求的范围之内,边界层外的网格均匀分布。对机匣处理部分采用ZR技术处理,将机匣进气段设置为环形进气,围绕整个叶轮进行划分。
2、计算及结果分析
计算设定与前期低工况计算类似,设定进口边界:流动方向为轴向,设定均匀总压为99300Pa,均匀总温为300K;出口边界:平均静压为188000Pa;固壁:无滑移固壁,绝热;转子转速:5567RPM;工质:空气。经过8000次迭代计算后达到收敛。
对这种改进型结构与原型机匣及优化周向槽机匣的气动特性进行对比分析。
(1)流场涡系分布
改进型结构机匣处理区流场流动情况所示,在机匣处理回流腔内产生了涡流,尤其在进气段,回转涡流相对于原型机匣作用区域减小,对下游叶栅的影响也相对降低。回转涡系产生了气动损失,为分析气动损失情况,计算熵增分布云图,可以看出,损失主要分布在环形进气与内壁机匣相接区域,但机匣内部整体损失与原型机匣相比明显降低,对下游叶栅造成的损失同样降低。
(2)相对马赫数分布
原型机匣、改进3型机匣与环形进气改进型机匣的叶顶截面马赫数分布所示,环形进气改进型机匣的截面马赫数最大,其流速最快,这主要是由于通过优化进气段与排气出口,增大了进气量同时增大了排气流速,在这种共同作用下,叶顶截面主流场受环形流场的作用,加速了叶顶区的流速。
(3)压比流量,效率流量特性
对原型、改进3型机、环形进气改进型机,0.2工况下的压比流量特性进行计算,压比-流量特性以及效率-流量特性计算所示,可以看出,环形进气改进型机其喘振裕度最大,但在增大了喘振裕度的同时气动效率最低。
将计算数据带入稳定裕度改进量计算公式,对原型机匣处理的喘振裕度改进量进行计算,其中,在转速为5567RPM情况下,工作点压比为2.6,流量为50kg/s,机匣处理近失速点压比为2.77,流量为47.13kg/s,得到稳定裕度改进量为12.9%,环形进气改进型机优化方案相对于原始机匣处理将喘振裕度增大了6.2%。
(4)喘振裕度
计算数据带入稳定裕度改进量计算公式,对原型机匣处理的喘振裕度改进量进行计算,其中,在转速为5567RPM情况下,工作点压比为2.6,流量为50kg/s,机匣处理近失速点压比为2.89,流量为49.13kg/s,得到稳定裕度改进量为13.11%。因此,优化方案相对于原始机匣处理将喘振裕度增大了6.04%。
5.2机匣处理性能对比
通过性能计算结果对比,可以得到如下结论:
优化排气段处理机匣在低工况下对稳定裕度改善效果最好,稳定裕度提高了6.04%。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种优化压气机机匣处理提高稳定裕度的方法,其特征在于,所述优化压气机机匣处理提高稳定裕度的方法包括:
对原机带轴向斜槽和回流腔机匣及参数改进型机匣进行数值模拟,确定参数改进型机匣的结构形式,并通过优化排气段、网格划分优化参数改进型机匣。
2.如权利要求1所述优化压气机机匣处理提高稳定裕度的方法,其特征在于,所述参数改进型机匣包括:
所述优化排气段包括:增大轴向速度,采用整个环形进气,通过调整进气段入口,扩大进气角度,增加进气量;在出气段,采用倾斜收缩喷嘴作为出气口,增大出口速度。
3.如权利要求1所述优化压气机机匣处理提高稳定裕度的方法,其特征在于,所述优化网格划分包括:
计算网格采用IGG/AutoGrid5进行单通道网格的划分,将计算区域分为主流区与机匣处理区;
对主流区网格拓扑结构采用O4H型结构,壁面附近的网格沿壁面法向方向按几何级数的规律加密,近壁面第一层网格距离固壁距离取0.001mm,将y+值控制在一定的低雷诺数湍流模型要求的范围之内,边界层外的网格均匀分布;
对机匣处理部分采用ZR技术处理,将机匣进气段设置为环形进气,围绕整个叶轮进行划分。
4.如权利要求1所述优化压气机机匣处理提高稳定裕度的方法,其特征在于,所述优化压气机机匣处理提高稳定裕度的方法包括以下步骤:
步骤一,获取原型机匣处理几何模型,建立原型机匣处理平面数值模型;
步骤二,计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分;
步骤三,基于计算的相关数据,计算原型机匣处理的喘振裕度;
步骤四,基于计算的相关数据进行机匣处理的结构优化,并计算结构优化后的机匣处理的喘振裕度;对比改进后机匣处理的喘振裕度与原型机匣处理的喘振裕度确定优化有效性。
5.如权利要求4所述优化压气机机匣处理提高稳定裕度的方法,其特征在于,步骤二中,所述计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分包括:
(1)计算域:采用CFD方法对涡轮流场进行数值模拟;
(2)网格划分:计算网格采用可视化IGG/AutoGrid5进行单通道网格的划分,将整个计算区域分为主流区与机匣处理区;所述主流流道利用NUMECA中专门针对叶轮机械部分网格生成模块AutoGrid整体划分网格,主流区网格拓扑结构采用O4H型结构,近壁面第一层网格距离固壁距离取0.001mm,将y+值控制在一定的低雷诺数湍流模型要求的范围之内,壁面附近的网格沿壁面法向方向按几何级数的规律加密,网格数约为该方向网格总数的三分之一,边界层外的网格均匀分布;机匣处理区和转子区之间设置转静子交界面;
(3)设计工况计算:通过设定边界条件、确定收敛标准,对压气机三维流场进行模拟,并基于模拟数值对机匣处理的总压分布、相对马赫数分布、压比流量进行分析。
6.如权利要求5所述优化压气机机匣处理提高稳定裕度的方法,其特征在于,所述通过设定边界条件,对压气机三维流场进行模拟包括:
对压气机三维流场进行模拟,进口给定总温,总压和气流角,出口给定平均静压,壁面采用绝热无滑移边界条件;
所述边界条件设置包括:
1)进口边界:流动方向为轴向,设定均匀总压为99300Pa,均匀总温为300K;
2)出口边界:平均静压为430000Pa;
3)固壁:无滑移固壁,绝热;
4)转子转速:7436RPM;
5)工质:空气;
6)转/静交界面处理:主流区各转/静交界面采用Conservation Coupling bypitchwise row,机匣处理结构的转静交界面采用Non reflecting 1D;
所述收敛标准包括:
总残差降低至一定水平后不再降低;总参数随迭代步数增加不再改变,总参数包括效率、压比、输出功率和扭矩;进出口流量差异不超过0.5%;总参数随迭代步数增加呈周期性波动。
8.一种优化压气机机匣处理提高稳定裕度系统,其特征在于,所述优化压气机机匣处理提高稳定裕度系统包括:
原型机匣处理平面数值模型构建模块,用于获取原型机匣处理几何模型,建立原型机匣处理平面数值模型;
原型机匣处理的网格划分模块,用于计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分;
原型机匣处理的喘振裕度计算模块,用于基于计算的相关数据,计算原型机匣处理的喘振裕度;
机匣处理的喘振裕度计算模块,用于基于计算的相关数据进行机匣处理的结构优化,并计算结构优化后的机匣处理的喘振裕度;对比改进后机匣处理的喘振裕度与原型机匣处理的喘振裕度确定优化有效性。
9.一种接收用户输入程序存储介质,所存储的计算机程序使电子设备执行权利要求1~7任意一项所述优化压气机机匣处理提高稳定裕度的方法,包括下列步骤:
步骤一,获取原型机匣处理几何模型,建立原型机匣处理平面数值模型;
步骤二,计算原型机匣处理的域、设计工况,确定原型机匣处理的网格划分;
步骤三,基于计算的相关数据,计算原型机匣处理的喘振裕度;
步骤四,基于计算的相关数据进行机匣处理的结构优化,并计算结构优化后的机匣处理的喘振裕度;对比改进后机匣处理的喘振裕度与原型机匣处理的喘振裕度确定优化有效性。
10.一种信息数据处理终端,其特征在于,所述信息数据处理终端包括存储器和处理器,所述存储器存储有计算机程序,所述计算机程序被所述处理器执行时,使得所述处理器执行权利要求1~7任意一项所述优化压气机机匣处理提高稳定裕度的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110649455.6A CN113283198A (zh) | 2021-06-10 | 2021-06-10 | 优化压气机机匣处理提高稳定裕度的方法、系统、终端 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110649455.6A CN113283198A (zh) | 2021-06-10 | 2021-06-10 | 优化压气机机匣处理提高稳定裕度的方法、系统、终端 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113283198A true CN113283198A (zh) | 2021-08-20 |
Family
ID=77284206
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110649455.6A Pending CN113283198A (zh) | 2021-06-10 | 2021-06-10 | 优化压气机机匣处理提高稳定裕度的方法、系统、终端 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113283198A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113569363A (zh) * | 2021-09-22 | 2021-10-29 | 中国航发上海商用航空发动机制造有限责任公司 | 多级压气机特性线获取方法和装置 |
CN114279714A (zh) * | 2021-12-27 | 2022-04-05 | 北京航空航天大学 | 高空低雷诺数下航空发动机涡轮试验台、模拟方法及应用 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112160943A (zh) * | 2020-10-13 | 2021-01-01 | 中国人民解放军海军工程大学 | 一种燃气轮机压气机机匣防喘结构及其喘振裕度评价方法 |
CN112685966A (zh) * | 2020-12-22 | 2021-04-20 | 中国船舶重工集团公司第七0三研究所 | 一种船舶燃气轮机压气机自循环式处理机匣设计方法 |
CN112685929A (zh) * | 2020-12-22 | 2021-04-20 | 中国船舶重工集团公司第七0三研究所 | 一种船舶燃气轮机压气机回流腔-扰流片式处理机匣设计方法 |
-
2021
- 2021-06-10 CN CN202110649455.6A patent/CN113283198A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112160943A (zh) * | 2020-10-13 | 2021-01-01 | 中国人民解放军海军工程大学 | 一种燃气轮机压气机机匣防喘结构及其喘振裕度评价方法 |
CN112685966A (zh) * | 2020-12-22 | 2021-04-20 | 中国船舶重工集团公司第七0三研究所 | 一种船舶燃气轮机压气机自循环式处理机匣设计方法 |
CN112685929A (zh) * | 2020-12-22 | 2021-04-20 | 中国船舶重工集团公司第七0三研究所 | 一种船舶燃气轮机压气机回流腔-扰流片式处理机匣设计方法 |
Non-Patent Citations (7)
Title |
---|
李思敏等: "圆弧槽处理机匣对影响跨声速压气机叶根失速先兆发展过程的实验研究", 《推进技术》, no. 12 * |
李霖: "直升机辅助动力装置吞砂试验装置设计及试验", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, 15 September 2019 (2019-09-15), pages 41 * |
王广等: "自循环机匣处理周向覆盖比例和喷嘴轴向位置对轴流压气机扩稳效果的影响", 《推进技术》, no. 06 * |
胡良军;杨策;孙惠民;张继忠: "机匣处理对离心压气机激波与泄漏涡干涉的影响", pages 521 - 526 * |
胡良军;杨策;祁明旭;马朝臣: "离心压气机机匣处理与导叶轮缘放气的扩稳机理分析及比较", pages 138 - 144 * |
邱兴海: "改进处理机匣扩大压气机的稳定工作裕度", 《燃气涡轮试验与研究》, no. 02 * |
马岩;席光;吴广宽;: "拓展小尺寸离心压气机级失速裕度机匣处理方案的数值研究", pages 1123 - 1126 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113569363A (zh) * | 2021-09-22 | 2021-10-29 | 中国航发上海商用航空发动机制造有限责任公司 | 多级压气机特性线获取方法和装置 |
CN113569363B (zh) * | 2021-09-22 | 2022-02-11 | 中国航发上海商用航空发动机制造有限责任公司 | 多级压气机特性线获取方法和装置 |
CN114279714A (zh) * | 2021-12-27 | 2022-04-05 | 北京航空航天大学 | 高空低雷诺数下航空发动机涡轮试验台、模拟方法及应用 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210209264A1 (en) | Modeling and calculation aerodynamic performances of multi-stage transonic axial compressors | |
CN112160943A (zh) | 一种燃气轮机压气机机匣防喘结构及其喘振裕度评价方法 | |
CN113283198A (zh) | 优化压气机机匣处理提高稳定裕度的方法、系统、终端 | |
CN109815624A (zh) | 一种考虑进气总压畸变影响的压气机稳定边界判断方法 | |
CN111079232A (zh) | 一种预测旋流畸变进气对航空发动机性能影响的计算方法 | |
Van den Braembussche et al. | Multidisciplinary multipoint optimization of a transonic turbocharger compressor | |
Li et al. | Aerodynamic optimization of a high-expansion ratio organic radial-inflow turbine | |
Li et al. | Development and application of a throughflow method for high-loaded axial flow compressors | |
Zhang et al. | Aerodynamic performance prediction of transonic axial multistage compressors based on one-dimensional meanline method | |
Kim et al. | Equivalent model for an axial compressor used for aero engines based on 1D and 3D analytical models and performance data | |
Kim et al. | New profile loss model for improved prediction of transonic axial flow compressor performance in choking region | |
Zhang et al. | Internal flow characteristics of a variable mixed flow turbine with partially-rotating vane nozzle | |
Yamagami et al. | Unsteady flow effects in a high-speed multistage axial compressor | |
Biesinger et al. | Speed line computation of a transonic compressor stage with unsteady cfd methods | |
Sieverding et al. | Design of industrial axial compressor blade sections for optimal range and performance | |
Abbasi et al. | Investigation of inlet distortion effects on axial compressor performance based on streamline curvature method | |
Zeng et al. | Evaluation of surge prediction capabilities of body-force model on a high-speed multi-stage axial-radial compressor | |
Rehman et al. | Non-Axisymmetric Endwall Profiling of a Stator Row in the Presence of the Rotor in a High Pressure Turbine | |
Hamzezade et al. | Numerical study of the effect of the tip gap size and using a single circumferential groove on the performance of a multistage compressor | |
Li et al. | Characteristics of inlet guide vane adjustment of multi-stage axial compressor in compressed air energy storage system | |
Ostad et al. | An Investigation on the Effect of Blade Tip Clearance on the Performance of a Single-Stage Axial Compressor | |
Ju et al. | Efficiency variation on a 4-stage low speed research compressor with a redesigned cantilevered stator | |
Tomita et al. | A model for numerical simulation of variable stator axial flow compressors | |
Zhou et al. | Multi-Objective optimization of an IGV for a large axial fan based on NSGA-II | |
Immery et al. | Design of the Compression System of a Geared Turbofan |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210820 |