CN110032784B - 带封严篦齿的高速轴流压气机的低速模化设计方法 - Google Patents
带封严篦齿的高速轴流压气机的低速模化设计方法 Download PDFInfo
- Publication number
- CN110032784B CN110032784B CN201910255611.3A CN201910255611A CN110032784B CN 110032784 B CN110032784 B CN 110032784B CN 201910255611 A CN201910255611 A CN 201910255611A CN 110032784 B CN110032784 B CN 110032784B
- Authority
- CN
- China
- Prior art keywords
- speed
- low
- grid
- compressor
- labyrinth
- 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/15—Vehicle, aircraft or watercraft 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
-
- 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)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
一种带封严篦齿的高速轴流压气机的低速模化设计方法,根据高速轴流压气机的转子与静子叶片叶型模化设计,确定压气机篦齿封严腔低速模化设计方案并分别进行建模和网格划分,并且通过雷诺平均NS方程求解设计点工况的流场数值解,对该流场进行不同方面的对比分析,最后将绝热效率最接近高速原型的低速模化设计方案中的篦齿个数作为低速方案的最优篦齿个数以实现其优化制造。本发明采用低速模化设计,得到低速轴流压气机的篦齿数的最佳个数,为航空发动机轴流压气机篦齿数的设计提供方向。
Description
技术领域
本发明涉及的是一种叶轮机械领域的技术,具体是一种带封严篦齿的高速轴流压气机的低速模化气动设计方法。
背景技术
轴流压气机的气动设计为航空发动机研制中的关键技术,轴流压气机的转子与静子之间多采用篦齿封严腔结构来封气,目的是减少转子与静子之间由于气流逆压梯度作用导致的漏气现象,提高压气机的绝热效率,并减少对压气机动稳定性的不利影响。轴流压气机篦齿封严处因存在径向间隙,从而形成了泄漏流,其与叶片根部表面和环形端壁上的二次流相互掺混,使得该处的流场结构复杂。目前,封严腔泄漏已成为影响航空发动机性能的主要因素之一,进一步研究发现,篦齿封严腔流场特性的准确认识有助于提高航空发动机的性能。
发明内容
本发明针对现有技术存在的上述不足,提出一种带封严篦齿的高速轴流压气机的低速模化气动设计方法,通过该方法设计低速大尺寸模型压气机实验件来开展工作,能够反映出高速原型压气机内部流场的结构特点和流场参数变化趋势,从而对高速原型压气机进行针对性改进设计,实现高速轴流压气机的优化制造。
本发明是通过以下技术方案实现的:
本发明根据高速、低速轴流压气机的转子与静子叶片叶型模化设计,确定压气机篦齿封严腔高速与低速模化设计方案并分别进行建模和网格划分,并且通过雷诺平均NS方程求解设计点工况的流场数值解,对该流场进行不同方面的对比分析,最后将绝热效率最接近高速原型的低速模化设计方案中的篦齿个数作为低速方案的最优篦齿个数以实现其优化制造。
所述的高速、低速轴流压气机的转子与静子叶片叶型模化设计是指:根据相似理论将高速轴流压气机的静子叶片进行放大,再根据流场对静子叶片几何参数进行修改,保证高、低速下静子叶片的气动参数一致。
所述的压气机篦齿封严腔高速与低速模化设计方案包括:压气机篦齿封严腔高速设计方案和压气机篦齿封严腔低速模化设计方案,其中:压气机篦齿封严腔高速设计方案是指:将单个篦齿的几何形状、篦齿径向间隙与叶高的比值、篦齿前后与主流交界面的轴向间隙与静子叶片中径弦长的比值、篦齿腔表面粗糙度;压气机篦齿封严腔低速模化设计方案是指:参照高速原型的篦齿数n,设计篦齿数为n-2、n-1和n三种低速模化方案,第一方案中:篦齿距离左侧篦齿外环为5.8cm、篦齿高度为1.5cm、上方宽度为0.1cm、下方宽度为0.6cm、篦齿间隙为0.08cm、篦齿斜角为72°;第二方案中:篦齿距离左侧篦齿外环为5.8cm、篦齿高度为1.5cm、上方宽度为0.1cm、下方宽度为0.6cm、篦齿间隙为0.08cm、篦齿斜角为72°、两个篦齿之间的距离为0.8cm;第三方案中:篦齿距离左侧篦齿外环距离为3.2cm、篦齿高度1.5cm、上方宽度为0.1cm、下方宽度为0.6cm、篦齿间隙为0.08cm、篦齿斜角为72°、两个篦齿之间的距离为0.8cm。
所述的建模,针对压气机篦齿封严腔低速模化设计方案和光壁方案进行建模。
所述的网格划分是指:主流区采用O4H网格分区,封严腔内采用H型网格,保证主流区与封严腔较好过渡,第一层网格高度设置为3*10-6m,网格的最小正交性保证30°。
所述的雷诺平均NS方程的数学表征为:采用数值模拟方法求解雷诺平均Navier-Stokes方程:其中:为守恒型参数向量, 和分别是无粘矢通量和粘性矢通量,qi为热源项,τij为应力,δij为克罗内克符号,Q为源项, 代表外作用力,Wf代表这些外作用力所做的功,
所述的对该流场进行不同方面的对比分析包括:对比不同低速模化设计方案的篦齿腔内静压云图与流线、静子叶片表面压力系数分布和静子叶片前缘进口俯仰气流角分布、对比高速轴流压气机和低速模化设计方案的低速压气机的绝热效率。
技术效果
与现有技术相比,本发明采用低速模化设计,得到高速轴流压气机的低速模化篦齿数的最佳个数,为航空发动机轴流压气机低速模化篦齿数的设计提供方向。
附图说明
图1为高速轴流压气机篦齿封严腔的结构示意图;
图2为三种篦齿封严腔低速模化设计几何结构,以及篦齿腔内静压云图和流线图;
图中:a为一个篦齿;b为两个篦齿;c为三个篦齿;
图3为三种篦齿数和光壁结构方案的静子叶片压力系数分布图;
图4为三种篦齿数和光壁结构方案的静子叶片前缘进口俯仰气流角分布图。
具体实施方式
如图1所示,本实施例以某航空发动机高速轴流压气机的一级转子与静子的叶片为例,该高速压气机静子叶片的篦齿封严腔采用三个篦齿的几何结构,所有篦齿间隙保持一致,为0.7%相对叶高。
本发明通过以下步骤实现:
1)制定高、低速模化设计方案:在完成高速、低速轴流压气机的转子与静子叶片叶型模化设计后,进行压气机篦齿封严腔的高速、低速模化设计,其中:
高速、低速轴流压气机的转子与静子叶片叶型模化设计是指:根据相似理论将高速压气机叶片进行放大,再根据流场对叶片几何参数进行修改,保证高低速的流场一致。
压气机篦齿封严腔低速模化设计是指:将单个篦齿的几何形状、篦齿径向间隙与叶高的比值、篦齿前后与主流交界面的轴向间隙与静子叶片中径弦长的比值、篦齿腔表面粗糙度全部与高速原型保持一致。
压气机篦齿封严腔低速模化设计是指:参照高速原型的篦齿数,设计了三种低速模化方案,如图2所示,分别为单篦齿方案SS_1tooth、双篦齿方案SS_2tooth和三篦齿方案SS_3tooth,其中:单篦齿方案中:篦齿距离左侧篦齿外环为5.8cm、篦齿高度为1.5cm、上方宽度为0.1cm、下方宽度为0.6cm、篦齿间隙为0.08cm、篦齿斜角为72°;双篦齿方案中:篦齿距离左侧篦齿外环为5.8cm、篦齿高度为1.5cm、上方宽度为0.1cm、下方宽度为0.6cm、篦齿间隙为0.08cm、篦齿斜角为72°、两个篦齿之间的距离为0.8cm;三篦齿方案中:篦齿距离左侧篦齿外环距离为3.2cm、篦齿高度1.5cm、上方宽度为0.1cm、下方宽度为0.6cm、篦齿间隙为0.08cm、篦齿斜角为72°、两个篦齿之间的距离为0.8cm。
2)对设计方案建模、网格划分并计算雷诺平均NS方程,具体步骤包括:选定如图1所示的高速轴流压气机转子与静子叶片,其中静子叶片根部带有篦齿封严腔,作为模化对象的原型,以步骤1)中的三种低速模化方案作为低速模化设计模型,并为了比较篦齿对压气机根部流场的影响,增设了第四个模型,即光壁结构IS,对单篦齿方案SS_1tooth、双篦齿方案SS_2tooth、三篦齿方案SS_3tooth和光壁方案IS建模并进行网格划分,采用数值模拟方法,在设计点工况下,求解雷诺平均NS方程。
所述的网格划分是指:主流区采用O4H网格分区,封严腔内采用H型网格,保证主流区与封严腔较好过渡,第一层网格高度设置为3*10-6m,网格的最小正交性保证30°。
所述的雷诺平均NS方程的数学表征为:采用数值模拟方法求解雷诺平均Navier-Stokes方程:其中:为守恒型参数向量, 和分别是无粘矢通量和粘性矢通量,qi为热源项,τij为应力,δij为克罗内克符号,Q为源项, 代表外作用力,Wf代表这些外作用力所做的功,
3)从不同方面进行计算结果的对比分析:由图2所示的单篦齿方案SS_1tooth、双篦齿方案SS_2tooth和三篦齿方案SS_3tooth模化方案的篦齿腔内静压云图和流线的对比,明确其腔内流场存在差异,该差异具体是指从后气封口到前气封口的静压下降过程差异及封严腔内流场结构的差异。
如图3和图4所示,依次展示了静子叶片表面压力系数和静子叶片前缘进口俯仰气流角分布,其中:叶片表面压力系数 为叶片通道进口总压,pin为叶片通道进口静压,pblade为叶片表面静压。不同篦齿封严结构,导致静子叶片根部2%叶高截面的气动载荷不同,影响到静子叶片前缘的进口俯仰气流角不同,最终导致低速模化压气机的性能发生变化。低速实验压气机模型的几何尺寸、转速、压力、温度等参数与高速压气机不同,但两者的气动特性绝热效率在模化时要保持相同或极其接近,绝热效率其中:为压气机的增压比,为压气机的出口总温,为压气机的进口总温,k为理想空气的绝热指数,得到:高速原型、IS、SS_1tooth、SS_2tooth、SS_3tooth低速模化设计方案的绝热效率计算值分别为:82.97%、89.22%、86.49%、82.70%、85.76%,由此得出,SS_2tooth方案的绝热效率与高速原型最接近,相对误差只有0.33%。因此,带封严篦齿的高速轴流压气机低速模化设计的篦齿数最终选择双篦齿方案SS_2tooth,即篦齿方案中:篦齿距离左侧篦齿外环为5.8cm,篦齿高度为1.5cm,上方宽度为0.1cm,下方宽度为0.6cm,篦齿间隙为0.08cm,篦齿斜角为72°,两个篦齿之间的距离为0.8cm。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。
Claims (7)
1.一种带封严篦齿的高速轴流压气机的低速模化设计方法,其特征在于,根据高速轴流压气机的转子与静子叶片叶型低速模化设计,确定压气机篦齿封严腔低速模化设计方案并分别进行建模和网格划分,并且通过雷诺平均NS方程求解设计点工况的流场数值解,对该流场进行不同方面的对比分析,最后将绝热效率最接近高速原型的低速模化设计方案中的篦齿个数作为低速方案的最优篦齿个数以实现其优化制造;
所述的压气机篦齿封严腔高速与低速模化设计方案包括:压气机篦齿封严腔高速设计方案和压气机篦齿封严腔低速模化设计方案,其中:压气机篦齿封严腔高速设计方案是指:单个篦齿的几何形状、篦齿径向间隙与叶高的比值、篦齿前后与主流交界面的轴向间隙与静子叶片中径弦长的比值、篦齿腔表面粗糙度;压气机篦齿封严腔低速模化设计方案是指:参照高速原型的篦齿数n,设计篦齿数为n-2、n-1和n三种低速模化方案,第一方案中:篦齿距离左侧篦齿外环为5.8cm、篦齿高度为1.5cm、上方宽度为0.1cm、下方宽度为0.6cm、篦齿间隙为0.08cm、篦齿斜角为72°;第二方案中:篦齿距离左侧篦齿外环为5.8cm、篦齿高度为1.5cm、上方宽度为0.1cm、下方宽度为0.6cm、篦齿间隙为0.08cm、篦齿斜角为72°、两个篦齿之间的距离为0.8cm;第三方案中:篦齿距离左侧篦齿外环距离为3.2cm、篦齿高度1.5cm、上方宽度为0.1cm、下方宽度为0.6cm、篦齿间隙为0.08cm、篦齿斜角为72°、两个篦齿之间的距离为0.8cm;
2.根据权利要求1所述的方法,其特征是,所述的高速轴流压气机的转子与静子叶片叶型低速模化设计是指:根据相似理论将高速轴流压气机的静子叶片进行放大,再根据流场对静子叶片几何参数进行修改,保证高、低速下静子叶片的流场一致。
3.根据权利要求1所述的方法,其特征是,所述的建模,针对压气机篦齿封严腔低速模化方案和光壁方案进行建模。
4.根据权利要求3所述的方法,其特征是,所述的网格划分是指:主流区采用O4H网格分区,封严腔内采用H型网格,保证主流区与封严腔较好过渡,第一层网格高度设置为3*10-6m,网格的最小正交性保证30°。
5.根据权利要求1所述的方法,其特征是,所述的对该流场进行不同方面的对比分析包括:对比不同低速模化设计方案的篦齿腔内静压云图与流线、静子叶片表面压力系数分布和静子叶片前缘进口俯仰气流角分布、对比高速轴流压气机和低速模化设计方案的低速压气机的绝热效率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910255611.3A CN110032784B (zh) | 2019-04-01 | 2019-04-01 | 带封严篦齿的高速轴流压气机的低速模化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910255611.3A CN110032784B (zh) | 2019-04-01 | 2019-04-01 | 带封严篦齿的高速轴流压气机的低速模化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110032784A CN110032784A (zh) | 2019-07-19 |
CN110032784B true CN110032784B (zh) | 2023-01-20 |
Family
ID=67237087
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910255611.3A Active CN110032784B (zh) | 2019-04-01 | 2019-04-01 | 带封严篦齿的高速轴流压气机的低速模化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110032784B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111259559B (zh) * | 2020-02-02 | 2023-03-24 | 上海交通大学 | 悬臂静子叶片前加载设计减小损失的流动控制方法 |
CN111255724B (zh) * | 2020-02-02 | 2021-07-16 | 上海交通大学 | 模拟轴流压气机高速旋转轮毂的悬臂静子平面叶栅实验方法 |
CN112945503B (zh) * | 2021-03-09 | 2023-02-10 | 上海交通大学 | 模拟轴流压气机静子叶片根部流动的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105840551A (zh) * | 2016-04-15 | 2016-08-10 | 上海交通大学 | 多工况点高负荷压气机叶片的气动实现方法 |
CN107063612A (zh) * | 2017-04-18 | 2017-08-18 | 东北大学 | 一种整体叶盘动力学相似测试试验台 |
-
2019
- 2019-04-01 CN CN201910255611.3A patent/CN110032784B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105840551A (zh) * | 2016-04-15 | 2016-08-10 | 上海交通大学 | 多工况点高负荷压气机叶片的气动实现方法 |
CN107063612A (zh) * | 2017-04-18 | 2017-08-18 | 东北大学 | 一种整体叶盘动力学相似测试试验台 |
Non-Patent Citations (5)
Title |
---|
Structural analysis of continuous fiber wound composite impellers of a multistage high-speed counter rotating axial compressor for compressing water vapor (R-718) as refrigerant using Finite Element Analysis;Mohit Patil 等;《Materials & design》;20131231;全文 * |
低速相似模化原理在高负荷轴流压气机设计中的应用;于成海;《中国优秀博硕士论文全文数据库 工程科技Ⅱ辑》;20170815;全文 * |
压气机低速模化设计及机匣处理研究;黄天豪;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20150715;全文 * |
压气机转子叶片叶尖流场的低速模化设计;李艳开 等;《航空动力学报》;20151231;全文 * |
高速压气机低速模化相似准则及轴向缝机匣处理流动机理研究;马宁;《中国优秀博硕士论文全文数据库 工程科技Ⅱ辑》;20161115;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110032784A (zh) | 2019-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210209264A1 (en) | Modeling and calculation aerodynamic performances of multi-stage transonic axial compressors | |
Hamakhan et al. | Aerodynamic performance effects of leading-edge geometry in gas-turbine blades | |
CN110032784B (zh) | 带封严篦齿的高速轴流压气机的低速模化设计方法 | |
Korakianitis et al. | Design of high-efficiency turbomachinery blades for energy conversion devices with the three-dimensional prescribed surface curvature distribution blade design (CIRCLE) method | |
Mischo et al. | Flow physics and profiling of recessed blade tips: impact on performance and heat load | |
CN112287580B (zh) | 一种基于全三维数值模拟的轴流压气机喘振边界计算方法 | |
Korakianitis et al. | Two-and three-dimensional prescribed surface curvature distribution blade design (CIRCLE) method for the design of high efficiency turbines, compressors, and isolated airfoils | |
Hopfinger et al. | Preliminary design of a three-stage low-speed research compressor using tandem vanes | |
Yang et al. | Numerical investigation of a cantilevered compressor stator at varying clearance sizes | |
Zhu et al. | Impact of hub gap leakage on stator endwall flow in an axial compressor stage with casing treatment | |
JP2009523211A (ja) | 凹状先端を有するタービンブレード | |
CN111255724B (zh) | 模拟轴流压气机高速旋转轮毂的悬臂静子平面叶栅实验方法 | |
Mischo et al. | Flow physics and profiling of recessed blade tips: Impact on performance and heat load | |
Falla | Numerical investigation of the flow in tandem compressor cascades | |
Uemura et al. | Simulation of unsteady flows through three-stage middle pressure steam turbine in operation | |
McGlumphy et al. | Numerical investigation of tandem airfoils for subsonic axial-flow compressor blades | |
Korakianitis et al. | Two-and three-dimensional prescribed surface curvature distribution blade design (circle) method for the design of high efficiency turbines, compressors, and isolated airfoils | |
Liang et al. | Structure and topology analysis of separated vortex in forward-swept blade | |
Kong et al. | Influence of cavity leakage flow on corner separation in a shrouded stator cascade | |
CN112989720B (zh) | 轴流压气机悬臂静子叶片优化轮毂径向间隙实现方法 | |
Xie et al. | Research On The Application Of Similarity Principle For Multistage Axial Compressor | |
Kulkarni | Development of a Methodology to Estimate Aero-Performance and Aero-Operability Limits of a Multistage Axial Flow Compressor for Use in Preliminary Design | |
Vincekovic et al. | Exploring Topology Optimisation of High Pressure Turbine Blade Tips | |
Korakianitis et al. | Aerodynamic improvements of wind-turbine airfoil geometries with the prescribed surface curvature distribution blade design (CIRCLE) method | |
Abda et al. | Axisymmetric Tip Gap Profiling in a Shroudless 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 |