CN112632719B - 一种基于一维平均流线法的多级轴流压气机特性校正方法 - Google Patents
一种基于一维平均流线法的多级轴流压气机特性校正方法 Download PDFInfo
- Publication number
- CN112632719B CN112632719B CN202011463281.6A CN202011463281A CN112632719B CN 112632719 B CN112632719 B CN 112632719B CN 202011463281 A CN202011463281 A CN 202011463281A CN 112632719 B CN112632719 B CN 112632719B
- Authority
- CN
- China
- Prior art keywords
- compressor
- scalar
- flow
- point
- streamline
- 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
Images
Classifications
-
- 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
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明的目的在于提供一种基于一维平均流线法的多级轴流压气机特性校正方法,根据多级轴流压气机运行所遵循的物理规律,建立基于一维平均流线法的多级轴流压气机性能分析方法,针对叶栅性能模型发展一种自动校准方法,搜寻压气机压比和效率与特定工作点的实验数据相匹配的叶片尾迹动量厚度系数,建立压气机已知条件范围的完整叶片尾迹动量厚度系数标量数据库;由已知叶片尾迹动量厚度系数标量数据库,求取所求工作点的尾迹动量厚度系数,实现对压气机特性的自动校准。本发明计算速度快,具有通用性,可以对压气机特性进行预测加密和外推,用于进口导叶、静叶可调时的特性预测,或计算没有实验数据或CFD数据转速下的整体性能。
Description
技术领域
本发明涉及的是一种仿真方法,具体地说是燃气轮机仿真方法。
背景技术
一维平均流线法性能求解中通过引入流体粘性作用的影响进行性能分析,因此落后角和损失等经验模型对一维计算精度起到了显著影响。目前公开文献中的叶栅经验模型主要是针对标准叶型。但是随着叶轮机械的发展,通过参数化优化或反设计等方法得到的现代叶型及弯、掠等先进技术被广泛应用于现代高性能压气机中,现有的经验模型对现代叶型的适用性有待进一步检验。因此,如何利用模型自动校准方法对经验模型进行校准成为性能分析的关键。
经验模型通常是测试结果的平均值或统计曲线的拟合,因此不能期望各种经验相关性可以准确地代表每一台压气机。在多级压气机中,不同叶片排所处的工况条件并不相同,经验模型对各叶片排的准确性也存在一定差异,模型偏差的不断累积加剧了性能预测结果对经验模型的敏感性。此外,叶轮机械内流动具有强三维、有旋及非定常等特点,在多级轴流压气机中流动情况尤为复杂,因此经验模型难以普遍适用于各种流动情况。通过修正经验模型实现对压气机特性的自动校准,可以进一步提高性能分析的准确性。
发明内容
本发明的目的在于提供用于进口导叶、静叶可调时的特性预测,或计算没有实验或CFD数据的转速下整体性能的一种基于一维平均流线法的多级轴流压气机特性校正方法。
本发明的目的是这样实现的:
本发明一种基于一维平均流线法的多级轴流压气机特性校正方法,其特征是:
(1)建立一维平均流线法模型
平均流线法基于一维流动假设,在压气机子午流面上指定一条流线,其半径的定义为:
Rt为外径,Rh为内径,Rm为平均半径;该流线称为平均流线;平均流线与各排叶片前后缘的交点即为流场参数的计算点,也称计算站,对于各计算站而言,确定流场收敛的原则是连续性条件得到满足,连续性方程即为平均流线方法的控制方程:
m=ρVmA(1-B)
其中Vm为平均流线上的气流速度,A为通流区域面积,B为堵塞量,ρ为密度,m为质量流量;
由各计算站所处位置的流动参数代表该压气机整体气动布局情况,通过迭代求解,使得连续性条件得到满足,则本站计算完成;根据已知的几何参数和边界条件,通过基本气动关系式对多级压气机逐排求解,并采用落后角及损失等经验模型来引入流体粘性的影响;
根据输入平均流线程序不同的进口温度、压力、转速、质量流量,输出不同的性能分析结果:
[PI,Eff]=f(T0,P0,n,G)
PI为压气机压比,Eff为压气机效率,n为转速,G为质量流量,T0、P0为进口温度与压力;
(2)选定优化变量
由二维低速叶栅损失试验得到的叶片尾迹动量厚度表达式为:
其中为叶型损失中的尾迹动量厚度,K1是常数,K2为综合影响因子,为等效扩压因子;取叶片尾迹动量厚度计算公式中的K1为设计变量,对特定工作点反复运行平均流线模型进行性能分析,通过优化算法搜寻得到与特定工作点压气机压比和效率与实验数据相匹配的最优叶片尾迹动量厚度系数K1;
(3)建立完整叶片尾迹动量厚度系数标量数据库
对不同压气机等转速线取3个特定工作点,包括一个近失速点、一个近堵塞点以及一个中间流量点,对每个工作点,以尾迹动量厚度拟合公式中的系数K1为设计变量,目标函数及约束条件为:
vlb≤K1≤ulb
RMSE为目标函数值,K1最优时,RMSE最小,PI0、Eff0为此工作点下试验数据,vlb为K1下限,ulb为上限;对每个工作点搜寻与特定工作点性能与实验数据相匹配的最优叶片尾迹动量厚度系数K1,建立压气机已知条件范围的完整叶片尾迹动量厚度系数标量数据库;
(4)建立中间标量的插值方法
计算压气机特性线或所需工作点特性时,对落在计算等转速线上的点,依据标量数据库中已知等转速线上的3个特定工作点的标量通过Hermite插值建立中间标量K1,代入平均流线法进行性能分析;对于不在标量数据库中的未知等转速线上的未知工作点,首先依据标量数据库中等转速线上3个特定工作点,通过线性插值建立未知工作点所在等转速线nInter上的3个特定工作点的标量K1,再依据未知工作点所在等转速线nInter上的3个特定工作点的标量K1通过Hermite插值建立未知工作点的中间标量K1,代入平均流线法型性能分析,由此得到性能分析结果,从而实现对压气机特性的自动校准。
本发明的优势在于:将压气机特性的自动校准方法应用于一维平均流线法中的经验损失模型,从而得到较符合试验数据的压气机整体特性参数,可以进一步提高性能分析的准确性。该方法计算速度快、精度较高,具有通用性,可以对压气机特性进行预测加密和外推,用于进口导叶、静叶可调时的特性预测,或计算没有实验或CFD数据的转速下整体性能,具有一定的应用价值。
附图说明
图1为本发明的流程图。
具体实施方式
下面结合附图举例对本发明做更详细地描述:
结合图1,以平均流线与各排叶片前后缘的交点为流场参数的计算点,由各计算站所处位置的流动参数代表该压气机整体气动布局情况,以落后角及损失等经验模型体现粘性、湍流、跨音、非定常等因素的影响,沿轴向分别计算各级(排)叶片出口参数,以此推进计算,最终通过叠加得到整机总体性能与气动布局参数。建立基于一维平均流线法的多级轴流压气机气动分析方法。在已知轴流压气机几何参数的基础上进行一维性能分析。
针对叶栅性能模型发展一种自动校准方法,由于特性和实验特性之间的差值主要是因为所选择的损失模型无法准确计算叶片尾迹动量厚度(也有其他因素的影响,但影响较小),不同的系数使尾迹动量厚度不同,压比与效率特性也不同。对不同压气机等转速线取3个或以上不同质量流量的特定工作点,以尾迹动量厚度拟合公式中的系数为优化变量,对每个工作点反复运行平均流线模型,搜索优化群体中的最优个体,得到与特定工作点压气机压比和效率与实验数据相匹配的最优叶片尾迹动量厚度系数,建立压气机已知条件范围的完整叶片尾迹动量厚度系数标量数据库。
由已知叶片尾迹动量厚度系数标量数据库,建立数据库之外工作点的尾迹动量厚度系数标量。对落在计算等转速线上的点通过Hermite插值建立中间标量;对各等转速线之间的工作点先通过线性插值建立特定工作点标量,再通过Hermite插值建立中间标量。平均流线法进行性能分析,由此得到性能分析结果,实现对压气机特性的自动校准。
本发明具体流程为:
(1)建立一维平均流线法模型
平均流线法基于一维流动假设。在压气机子午流面上指定一条特殊流线,其半径的定义为:
Rt为外径,Rh为内径,Rm为平均半径。
则该流线称为平均流线。平均流线与各排叶片前后缘的交点即为流场参数的计算点,也称计算站。对于各计算站而言,确定流场收敛的原则是连续性条件得到满足,故连续性方程即为平均流线方法的控制方程:
m=ρVmA(1-B) (2)
其中Vm为平均流线上的气流速度,A为通流区域面积,B为堵塞量,ρ为密度,m为质量流量。
由各计算站所处位置的流动参数代表该压气机整体气动布局情况,通过迭代求解,使得连续性条件得到满足,则本站计算完成。根据已知的几何参数和边界条件,通过基本气动关系式对多级压气机逐排求解,并采用落后角及损失等经验模型来引入流体粘性的影响。
根据输入平均流线程序不同的进口温度、压力、转速、质量流量,将输出不同的性能分析结果:
[PI,Eff]=f(T0,P0,n,G) (2)
PI为压气机压比,Eff为压气机效率,n为转速,G为质量流量,T0、P0为进口温度与压力。
(2)优化变量的选定
亚音速压气机中,叶型损失占比较大。研究已经表明,绝大多数叶型损失模型只考虑了尾缘动量损失厚度的作用。由二维低速叶栅损失试验得到的叶片尾迹动量厚度表达式为:
其中为叶型损失中的尾迹动量厚度,K1是常数,一般取值0.0073,K2为综合影响因子,约为1,为等效扩压因子。取叶片尾迹动量厚度计算公式中的K1为设计变量,对某一特定工作点反复运行平均流线模型进行性能分析,通过优化算法搜寻得到与特定工作点压气机压比和效率与实验数据相匹配的最优叶片尾迹动量厚度系数K1。
(3)建立完整叶片尾迹动量厚度系数标量数据库
对不同压气机等转速线取3个特定工作点,包括一个近失速点、一个近堵塞点以及一个中间流量点。对于表现出较强曲率的压气机特性,需要根据所需的精度取更多得工作点。对每个工作点,以尾迹动量厚度拟合公式中的系数K1为设计变量,目标函数及约束条件为:
vlb≤K1≤ulb (6)
RMSE为目标函数值,K1最优时,RMSE最小,PI0、Eff0为此工作点下试验数据,vlb为K1下限,ulb为上限。对每个工作点搜寻与特定工作点性能与实验数据相匹配的最优叶片尾迹动量厚度系数K1,建立压气机已知条件范围的完整叶片尾迹动量厚度系数标量数据库。
(4)建立中间标量的插值方法
计算压气机特性线或所需工作点特性时,以不同压气机等转速线为例:对落在计算等转速线上的点,依据标量数据库中已知等转速线上的3个特定工作点的标量通过Hermite插值建立中间标量K1,代入平均流线法进行性能分析;对于不在标量数据库中的未知等转速线上的未知工作点,首先依据标量数据库中等转速线上3个特定工作点,通过线性插值建立未知工作点所在等转速线nInter上的3个特定工作点的标量K1,再依据未知工作点所在等转速线nInter上的3个特定工作点的标量K1通过Hermite插值建立未知工作点的中间标量K1,代入平均流线法型性能分析。由此得到性能分析结果,从而实现对压气机特性的自动校准。
Claims (1)
1.一种基于一维平均流线法的多级轴流压气机特性校正方法,其特征是:
(1)建立一维平均流线法模型
平均流线法基于一维流动假设,在压气机子午流面上指定一条流线,其半径的定义为:
Rt为外径,Rh为内径,Rm为平均半径;该流线称为平均流线;平均流线与各排叶片前后缘的交点即为流场参数的计算点,也称计算站,对于各计算站而言,确定流场收敛的原则是连续性条件得到满足,连续性方程即为平均流线方法的控制方程:
m=ρVmA(1-B)
其中Vm为平均流线上的气流速度,A为通流区域面积,B为堵塞量,ρ为密度,m为质量流量;
由各计算站所处位置的流动参数代表该压气机整体气动布局情况,通过迭代求解,使得连续性条件得到满足,则本站计算完成;根据已知的几何参数和边界条件,通过基本气动关系式对多级压气机逐排求解,并采用落后角及损失等经验模型来引入流体粘性的影响;
根据输入平均流线程序不同的进口温度、压力、转速、质量流量,输出不同的性能分析结果:
[PI,Eff]=f(T0,P0,n,G)
PI为压气机压比,Eff为压气机效率,n为转速,G为质量流量,T0、P0为进口温度与压力;
(2)选定优化变量
由二维低速叶栅损失试验得到的叶片尾迹动量厚度表达式为:
其中为叶型损失中的尾迹动量厚度,K1是常数,K2为综合影响因子,为等效扩压因子;取叶片尾迹动量厚度计算公式中的K1为设计变量,对特定工作点反复运行平均流线模型进行性能分析,通过优化算法搜寻得到与特定工作点压气机压比和效率与实验数据相匹配的最优叶片尾迹动量厚度系数K1;
(3)建立完整叶片尾迹动量厚度系数标量数据库
对不同压气机等转速线取3个特定工作点,包括一个近失速点、一个近堵塞点以及一个中间流量点,对每个工作点,以尾迹动量厚度拟合公式中的系数K1为设计变量,目标函数及约束条件为:
vlb≤K1≤ulb
RMSE为目标函数值,K1最优时,RMSE最小,PI0、Eff0为此工作点下试验数据,vlb为K1下限,ulb为上限;对每个工作点搜寻与特定工作点性能与实验数据相匹配的最优叶片尾迹动量厚度系数K1,建立压气机已知条件范围的完整叶片尾迹动量厚度系数标量数据库;
(4)建立中间标量的插值方法
计算压气机特性线或所需工作点特性时,对落在计算等转速线上的点,依据标量数据库中已知等转速线上的3个特定工作点的标量通过Hermite插值建立中间标量K1,代入平均流线法进行性能分析;对于不在标量数据库中的未知等转速线上的未知工作点,首先依据标量数据库中等转速线上3个特定工作点,通过线性插值建立未知工作点所在等转速线nInter上的3个特定工作点的标量K1,再依据未知工作点所在等转速线nInter上的3个特定工作点的标量K1通过Hermite插值建立未知工作点的中间标量K1,代入平均流线法型性能分析,由此得到性能分析结果,从而实现对压气机特性的自动校准。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011463281.6A CN112632719B (zh) | 2020-12-11 | 2020-12-11 | 一种基于一维平均流线法的多级轴流压气机特性校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011463281.6A CN112632719B (zh) | 2020-12-11 | 2020-12-11 | 一种基于一维平均流线法的多级轴流压气机特性校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112632719A CN112632719A (zh) | 2021-04-09 |
CN112632719B true CN112632719B (zh) | 2022-05-17 |
Family
ID=75312564
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011463281.6A Active CN112632719B (zh) | 2020-12-11 | 2020-12-11 | 一种基于一维平均流线法的多级轴流压气机特性校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112632719B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114117665B (zh) * | 2021-11-15 | 2024-04-09 | 北京动力机械研究所 | 一种s2流面框架下轴流压气机经验模型标定方法 |
CN115544667B (zh) * | 2022-10-31 | 2024-05-10 | 南京航空航天大学 | 一种基于叶素动量源耦合cfd的等效盘方法 |
CN116663202B (zh) * | 2023-07-27 | 2023-10-17 | 中国航发四川燃气涡轮研究院 | 一种多级轴流压气机性能仿真模型的校核方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105370404A (zh) * | 2014-08-15 | 2016-03-02 | 通用电气公司 | 具有单型低损耗轴承和低密度材料的机械驱动体系 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100350158C (zh) * | 2002-08-23 | 2007-11-21 | 约克国际公司 | 用于探测离心压缩机内的旋转失速的系统和方法 |
-
2020
- 2020-12-11 CN CN202011463281.6A patent/CN112632719B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105370404A (zh) * | 2014-08-15 | 2016-03-02 | 通用电气公司 | 具有单型低损耗轴承和低密度材料的机械驱动体系 |
Non-Patent Citations (7)
Title |
---|
Axial compressor performance modelling with quasi-one-dimensional approach;WHITE N M 等;《Proceedings of Institution of Mechanical Engineers》;20020101;181-193 * |
Axial-flow compressor performance prediction in design and off-design conditions through 1-D and 3-D modeling and experimental study;PEYVAN A 等;《Journal of Applied Fluid Mechanics》;20160731;2149-2160 * |
一种多级轴流压气机特性预估方法探讨;张跃学 等;《航空科学技术》;20110815(第04期);65-67 * |
几何参数变化对离心压气机性能的影响;靳军 等;《科学技术与工程》;20141218;第14卷(第35期);296-302、307 * |
基于时间推进的轴流压气机准二维性能分析;张炯 等;《西北工业大学学报》;20200215;第38卷(第01期);114-120 * |
多级轴流压气机多排可转导/静叶联合调节规律研究;廖吉香 等;《推进技术》;20170111;第38卷(第02期);334-340 * |
流线曲率法在多级跨声速轴流压气机特性预测中的应用;巫骁雄 等;《推进技术》;20170828;第38卷(第10期);2235-2245 * |
Also Published As
Publication number | Publication date |
---|---|
CN112632719A (zh) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112632719B (zh) | 一种基于一维平均流线法的多级轴流压气机特性校正方法 | |
Schnoes et al. | Design optimization of a multi-stage axial compressor using throughflow and a database of optimal airfoils | |
CN111027148B (zh) | 损失落后角模型自动标定及工业轴流压缩机性能计算方法 | |
CN106382253B (zh) | 流量系数0.02管线压缩机模型级及叶轮设计方法 | |
JP2009144716A (ja) | ターボ機械用の多段タービンを設計する方法 | |
CN110175340A (zh) | 一种适用于不确定工况的离心压气机叶轮气动稳健优化设计方法 | |
CN107203364B (zh) | 一种用于压气机全工况特性的预测和辨识方法 | |
CN114254460A (zh) | 基于图卷积神经网络的透平机械气动鲁棒性优化方法 | |
CN107908914B (zh) | 离心压缩机闭式叶轮可加工性判断及中间截面的计算方法 | |
Olivero et al. | Three-dimensional turbulent optimization of vaned diffusers for centrifugal compressors based on metamodel-assisted genetic algorithms | |
CN115017843A (zh) | 一种离心压缩机气动性能优化设计方法 | |
CN116399541A (zh) | 基于深度神经网络的叶栅风洞实验工况参数修正方法 | |
Schmitz et al. | Novel performance prediction of a transonic 4.5 stage compressor | |
Ellbrant et al. | CFD optimization of a transonic compressor using multiobjective GA and metamodels | |
CN117195760A (zh) | 基于径向掺混的轴流风扇或压气机子午面通流计算方法 | |
Lengyel-Kampmann et al. | Generalized optimization of counter-rotating and single-rotating fans | |
Luers et al. | Adjoint-based volumetric shape optimization of turbine blades | |
CN116894298A (zh) | 基于cfd/s2混合维度的多级轴流压气机特性预测方法 | |
CN116595874A (zh) | 叶轮机械性能预测模型参数优化方法及装置、存储介质 | |
CN115906430A (zh) | 轴流压气机篦齿泄漏损失预测方法 | |
CN114720145A (zh) | 一种带整流叶片的低压涡轮性能试验方法 | |
CN111288016B (zh) | 轴流压气机基元叶型造型方法 | |
CN109241585B (zh) | 一种高低压涡轮过渡流道型面反问题设计方法 | |
CN112800554A (zh) | 一种叶片表面粗糙度变化影响压气机稳定性的仿真方法 | |
CN113361028A (zh) | 一种蜗壳二维设计方法 |
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 |