CN112983560B - 一种涡轮气冷叶片温度场快速评估方法 - Google Patents
一种涡轮气冷叶片温度场快速评估方法 Download PDFInfo
- Publication number
- CN112983560B CN112983560B CN202110394939.0A CN202110394939A CN112983560B CN 112983560 B CN112983560 B CN 112983560B CN 202110394939 A CN202110394939 A CN 202110394939A CN 112983560 B CN112983560 B CN 112983560B
- Authority
- CN
- China
- Prior art keywords
- function
- blade
- temperature
- turbine air
- cooling
- 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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D5/00—Blades; Blade-carrying members; Heating, heat-insulating, cooling or antivibration means on the blades or the members
- F01D5/12—Blades
- F01D5/14—Form or construction
- F01D5/18—Hollow blades, i.e. blades with cooling or heating channels or cavities; Heating, heat-insulating or cooling means on blades
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D21/00—Shutting-down of machines or engines, e.g. in emergency; Regulating, controlling, or safety means not otherwise provided for
- F01D21/003—Arrangements for testing or measuring
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01D—NON-POSITIVE DISPLACEMENT MACHINES OR ENGINES, e.g. STEAM TURBINES
- F01D25/00—Component parts, details, or accessories, not provided for in, or of interest apart from, other groups
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本申请属于涡轮气冷叶片领域,特别涉及一种涡轮气冷叶片温度场快速评估方法。包括:步骤一、构建带有隔热涂层的涡轮气冷叶片导热模型;步骤二、对综合冷效θ进行解耦,获取表征外部冷却影响的第一函数以及表征内部冷却影响的第二函数,并将所述第一函数与所述第二函数通过共用横坐标x进行关联;步骤三、分别对所述第一函数、所述第二函数通过等值线求交点进行求解,并通过共用的横坐标x对所述第一函数以及所述第二函数进行耦合求解;步骤四、求解叶片外表面壁温。本申请的涡轮气冷叶片温度场快速评估方法,需要的初始参数少,计算速度快,能够实现多参数图形化显示。
Description
技术领域
本申请属于涡轮气冷叶片领域,特别涉及一种涡轮气冷叶片温度场快速评估方法。
背景技术
涡轮叶片的外表面壁温是衡量叶片工作状态的重要指标,是叶片设计过程中的关键参数。但由于叶片外部流动换热非常复杂、内部冷却形式多样,导致影响冷却叶片外表面壁温的因素繁多,包括气膜冷却有效性、叶片壁厚、热障涂层壁厚、叶片导热系数、热障涂层导热系数、燃气冷气温度以及对流换热系数等,且各个因素相互影响,呈强耦合状态,难以在冷却方案设计阶段准确把握每种因素的影响规律。
现有的三维CFD流热热耦合计算虽然可以得到较为精确的叶片外表面壁温,但是该方法不仅需要给出冷却结构的详细结构参数,比如气膜孔排数、气膜孔孔径、气膜孔间距等,而且需要详细的冷气参数,比如冷气温度、冷气压力、冷气质量流量等,这些参数在方案设计阶段只有一个大概的方向,没法达到非常详细的程度。并且,三维CFD耦合的每次计算都需要处理几何模型、网格划分和求解器设置这三个步骤,每一个步骤都需要手动做大量的设置,工作量巨大。三维CFD耦合计算花费的时间非常长,一个算例从开始到结束往往需要半个月甚至一个月的时间,如果修改多个参数进行对比,会花费几个月的时间,反馈周期漫长。在三维CFD耦合计算中,每一个参数的更改都需要重新建模、重新分网、重新计算,修改非常繁琐,不够灵活。另外,三维CFD计算出来的结果是一个综合的结果,呈现的是各个因素耦合作用的结果,没法体现单独一个因素在整个方案中所占的比重以及单一因素的动态影响。
因此,希望有一种技术方案来克服或至少减轻现有技术的至少一个上述缺陷。
发明内容
本申请的目的是提供了一种涡轮气冷叶片温度场快速评估方法,以解决现有技术存在的至少一个问题。
本申请的技术方案是:
一种涡轮气冷叶片温度场快速评估方法,包括:
步骤一、构建带有隔热涂层的涡轮气冷叶片导热模型,包括:
综合冷效θ为:
气膜冷却有效性η为:
叶片毕涡数Bim为:
涂层毕涡数Bit为:
由一维导热原理,可得综合冷效θ为:
其中,δ1为涂层厚度,λ1为涂层导热系数,δ2为金属壁厚度,λ2为金属壁导热系数,Tg为燃气温度,Tce为冷气温度,Tw2为叶片外表面壁温,Taw为叶片绝热壁温,hg为燃气侧对流换热系数,hc为冷气侧对流换热系数;
步骤二、对综合冷效θ进行解耦,获取表征外部冷却影响的第一函数以及表征内部冷却影响的第二函数,并将所述第一函数与所述第二函数通过共用横坐标x进行关联,其中,
所述第一函数为:
f1(θ,η,x)=0
所述第二函数为:
f2(Bim,Bit,hg/hc,x)=0;
步骤三、分别对所述第一函数、所述第二函数通过等值线求交点进行求解,并通过共用的横坐标x对所述第一函数以及所述第二函数进行耦合求解;
步骤四、求解叶片外表面壁温Tw2:
Tw2=Tg-θ*(Tg-Tce)。
可选地,步骤一中,
由:
可得:
由:
可得:
因此,有:
设通过壁面的热流密度为q,由一维导热原理可得:
代入综合冷效θ,可得:
可选地,步骤二中,
构建的共用横坐标x的表达式为:
则所述第一函数为:
θ-η+(η-1)*x=0
所述第二函数为:
可选地,步骤二中,将所述第一函数与所述第二函数通过共用横坐标x关联到同一二维平面上。
可选地,所述第一函数在二维平面的上半区进行显示,所述第二函数在二维平面的下半区进行显示。
发明至少存在以下有益技术效果:
本申请的涡轮气冷叶片温度场快速评估方法,只需要提供9个初始参数就能求出一维温度评估所需要的5个集总参数,进而对叶片表面的温度场进行计算,不需要知道冷却结构的详细参数,非常适用于叶片冷却方案的设计阶段;计算速度快,一个冷却方案的评估在10分钟以内;通过多参数图形化显示,将各个因素的相互影响关系以图形化方式直观的表现出来,物理意义更加明晰。
附图说明
图1是本申请一个实施方式的带有隔热涂层的涡轮气冷叶片导热模型示意图;
图2是本申请一个实施方式的涡轮气冷叶片温度场快速评估方法的多参数图形显示原理图;
图3是本申请一个实施方式的涡轮气冷叶片温度场快速评估方法的第一函数解耦计算过程示意图;
图4是本申请一个实施方式的涡轮气冷叶片温度场快速评估方法的第二函数解耦计算过程示意图;
图5是本申请一个实施方式的涡轮气冷叶片温度场快速评估方法的耦合求解计算过程示意图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本申请一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。下面结合附图对本申请的实施例进行详细说明。
在本申请的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请保护范围的限制。
下面结合附图1至图5对本申请做进一步详细说明。
本申请提供了一种涡轮气冷叶片温度场快速评估方法,包括以下步骤:
步骤一、如图1所示,构建带有隔热涂层的涡轮气冷叶片导热模型,包括:
综合冷效θ为:
气膜冷却有效性η为:
叶片毕涡数Bim为:
涂层毕涡数Bit为:
由一维导热原理,可得综合冷效θ为:
其中,δ1为涂层厚度,λ1为涂层导热系数,δ2为金属壁厚度,λ2为金属壁导热系数,Tg为燃气温度,Tce为冷气温度,Tw2为叶片外表面壁温,Taw为叶片绝热壁温,hg为燃气侧对流换热系数,hc为冷气侧对流换热系数。
本申请中通过公式(5)对叶片外表面壁温进行快速评估。
公式(5)的具体推到过程如下:
由综合冷效的定义公式(1)可得:
从气膜有效冷却有效性定义式(2)可得:
将公式(b)代入公式(a),有:
设通过壁面的热流密度为q,由一维导热原理可得:
将公式(d)、(e)代入公式(c),可得综合冷效θ为:
由公式(5)可知,影响综合冷效的因素有4个,分别是气膜冷却有效性、金属叶片毕涡数、涂层毕涡数和内外流对流换热系数之比,加上综合冷效,公式中共包含了5个变量。
在一个二维平面上最多表达3个变量间的相互关系,分别为横坐标、纵坐标以及等值线,无法清晰的表达5个变量的相互关系。为解决该问题,本申请通过构建共用横坐标轴将两个二维平面图拼接在一起的方法,实现5个变量间相互关系的清晰表达,如图2所示,图2中A和A’即为同一个状态点。
步骤二、对综合冷效θ进行解耦,获取表征外部冷却影响的第一函数以及表征内部冷却影响的第二函数,并将所述第一函数与所述第二函数通过共用横坐标x进行关联。
将公式(5)拆分成两个子公式,一个子公式体现外部冷却的影响,即综合冷效θ与气膜冷效η间的函数关系f1(θ,η,x)=0;另一个子公式体现内部冷却的影响,即叶片毕涡数Bim,涂层毕涡数Bit及冷热对流换热系数之比hg/hc间的函数关系f2(Bim,Bit,hg/hc,x)=0;两个子公式通过构建的共用横坐标x联系在一起,将第一函数与第二函数通过共用横坐标x关联到同一二维平面上,第一函数f1在图形上半区进行显示,第二函数f2在图形下半区进行显示。如图所示,图3为函数f1上半区图示,图4为函数f2下半区图示。
具体的,
构建的共用横坐标x的表达式为:
第一函数为:
f1(θ,η,x)=0
即:
θ-η+(η-1)*x=0 (7)
第二函数为:
f2(Bim,Bit,hg/hc,x)=0;
即:
步骤三、在对应的上下半区分别对第一函数、第二函数通过等值线求交点进行求解,并通过共用的横坐标x对第一函数以及第二函数进行耦合求解,如图5所示。
步骤四、求解叶片外表面壁温Tw2:
Tw2=Tg-θ*(Tg-Tce) (9)
通过公式(5)可以求得综合冷效θ,最终解算出叶片外表面壁温Tw2。
在本申请的一个实施方式中,在模型中输入参数:δ1=0.5mm,λ1=2W/(mK);δ2=10mm,λ2=50W/(mK);Tg=1000K,Tce=100K;hg=5000W(m2K),hc=500W(m2K);η=0.5。
求解出Bim=1,Bit=1.25,hg/hc=10,代入公式(8)可得共用x坐标值x=0.1698;
将η、x值代入公式(7)可得综合冷效θ=0.4151;
最后,将综合冷效θ、燃气温度Tg、冷气温度Tce代入公式(1),可得叶片外表面壁面温为Tw2=626.41K。
本申请的涡轮气冷叶片温度场快速评估方法,利用相似定理,将多个参数组合在一起形成集总参数,将影响叶片温度场的诸多参数简化为5个集总参数,即综合冷效、气膜冷效、叶片毕涡数、涂层毕涡数、冷热对流换热系数比;基于一维导热原理,推导出综合冷效、气膜冷效、叶片毕涡数、涂层毕涡数、冷热对流换热系数比这5个集总参数的函数关系式,用该函数关系式快速评估叶片外表面壁温;实现多参数图形化显示,将5个集总参数显示在一张图上,直观的体现各个因素对整体方案的影响。
本申请需要的初始参数少,只需要提供9个初始参数就能求出一维温度评估所需要的5个集总参数,进而对叶片表面的温度场进行计算,不需要知道冷却结构的详细参数,非常适用于叶片冷却方案的设计阶段;计算速度快,一个冷却方案的评估在10分钟以内;通过多参数图形化显示,将各个因素的相互影响关系以图形化方式直观的表现出来,并且将内部冷却和外部冷却分别放在上下半区显示,物理意义更加明晰。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (1)
1.一种涡轮气冷叶片温度场快速评估方法,其特征在于,包括:
步骤一、构建带有隔热涂层的涡轮气冷叶片导热模型,包括:
综合冷效θ为:
气膜冷却有效性η为:
叶片毕涡数Bim为:
涂层毕涡数Bit为:
由一维导热原理,可得综合冷效θ为:
其中,δ1为涂层厚度,λ1为涂层导热系数,δ2为金属壁厚度,λ2为金属壁导热系数,Tg为燃气温度,Tce为冷气温度,Tw2为叶片外表面壁温,Taw为叶片绝热壁温,hg为燃气侧对流换热系数,hc为冷气侧对流换热系数;
步骤二、对综合冷效θ进行解耦,获取表征外部冷却影响的第一函数以及表征内部冷却影响的第二函数,并将所述第一函数与所述第二函数通过共用横坐标x进行关联,其中,
所述第一函数为:
f1(θ,η,x)=0
所述第二函数为:
f2(Bim,Bit,hg/hc,x)=0;
步骤三、分别对所述第一函数、所述第二函数通过等值线求交点进行求解,并通过共用的横坐标x对所述第一函数以及所述第二函数进行耦合求解;
步骤四、求解叶片外表面壁温Tw2:
Tw2=Tg-θ*(Tg-Tce);
步骤一中,
由:
可得:
由:
可得:
因此,有:
设通过壁面的热流密度为q,由一维导热原理可得:
代入综合冷效θ,可得:
步骤二中,
构建的共用横坐标x的表达式为:
则所述第一函数为:
θ-η+(η-1)*x=0
所述第二函数为:
步骤二中,将所述第一函数与所述第二函数通过共用横坐标x关联到同一二维平面上;
所述第一函数在二维平面的上半区进行显示,所述第二函数在二维平面的下半区进行显示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110394939.0A CN112983560B (zh) | 2021-04-13 | 2021-04-13 | 一种涡轮气冷叶片温度场快速评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110394939.0A CN112983560B (zh) | 2021-04-13 | 2021-04-13 | 一种涡轮气冷叶片温度场快速评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112983560A CN112983560A (zh) | 2021-06-18 |
CN112983560B true CN112983560B (zh) | 2022-11-22 |
Family
ID=76338244
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110394939.0A Active CN112983560B (zh) | 2021-04-13 | 2021-04-13 | 一种涡轮气冷叶片温度场快速评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112983560B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114741817B (zh) * | 2022-05-07 | 2023-05-05 | 中国航发沈阳发动机研究所 | 一种带热障涂层的涡轮冷却叶片设计方法 |
CN115688327A (zh) * | 2022-12-27 | 2023-02-03 | 中国航发沈阳发动机研究所 | 一种带分区变厚度热障涂层的涡轮叶片温度场计算方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005069236A (ja) * | 2004-12-10 | 2005-03-17 | Toshiba Corp | タービン冷却翼 |
CN108009340A (zh) * | 2017-11-28 | 2018-05-08 | 中国航发沈阳发动机研究所 | 一种涡轮冷却叶片三维稳态温度场的计算方法 |
KR102067103B1 (ko) * | 2018-04-05 | 2020-01-16 | 두산중공업 주식회사 | 터빈의 온도 추정 시스템 |
CN109446592B (zh) * | 2018-10-09 | 2020-02-07 | 湘潭大学 | 一种涡轮叶片热障涂层应用效果的评价方法 |
CN112507586A (zh) * | 2020-12-02 | 2021-03-16 | 中国航发沈阳发动机研究所 | 一种涡轮气冷叶片二维温度、强度的快速评估方法 |
-
2021
- 2021-04-13 CN CN202110394939.0A patent/CN112983560B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112983560A (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112983560B (zh) | 一种涡轮气冷叶片温度场快速评估方法 | |
Atkinsona et al. | Two-and three-dimensional numerical models of flow and heat transfer over louvred fin arrays in compact heat exchangers | |
Van Zante et al. | Recommendations for achieving accurate numerical simulation of tip clearance flows in transonic compressor rotors | |
Khoshvaght Aliabadi et al. | 3D-CFD simulation and neural network model for the j and f factors of the wavy fin-and-flat tube heat exchangers | |
Lotfi et al. | Thermo-hydraulic performance enhancement of finned elliptical tube heat exchangers by utilizing innovative dimple turbulators | |
Kong et al. | Impacts of geometric structures on thermo-flow performances of plate fin-tube bundles | |
Tang et al. | Parameters study and analysis of turbulent flow and heat transfer enhancement in narrow channel with discrete grooved structures | |
Gowardhan et al. | Evaluation of a fast response pressure solver for flow around an isolated cube | |
Välikangas et al. | Conjugated heat transfer simulation of a fin-and-tube heat exchanger | |
Ahmed et al. | Three-Dimensional transient heat conduction equation solution for accurate determination of heat transfer coefficient | |
Miao et al. | Experimental and numerical analyses of thermal-hydraulic characteristics of aluminium flying-wing fins | |
Taler | Direct and inverse heat transfer problems in dynamics of plate fin and tube heat exchangers | |
CN105302964A (zh) | 一种用于芯片结构的热分布分析方法 | |
CN113051846B (zh) | 考虑可压缩和热传导效应的壁面第一层网格厚度估算方法 | |
He et al. | Cross flow and heat transfer of hollow-fiber tube banks with complex distribution patterns and various baffle designs | |
Kirsch et al. | Heat transfer measurements of oblong pins | |
Chen et al. | Shape optimization of fluid cooling channel based on Darcy reduced-order isogeometric analysis | |
Gramespacher et al. | The influence of element thermal conductivity, shape, and density on heat transfer in a rough wall turbulent boundary layer with strong pressure gradients | |
Djemel et al. | Three-dimensional numerical study of the performance of the intercooler equipped with vortex generators | |
Li et al. | A one-dimensional analytical method for turbine blade preliminary cooling design | |
Gomes et al. | The Concept of Adiabatic Heat Transfer Coefficient and Its Application to Turbomachinery | |
Satish et al. | Comparison of predictions obtained on an exhaust manifold analysis using conformal and indirect mapped interface | |
Wang | Investigation on heat transfer in circular fins with H-shaped cavity using finite element method and constructal design | |
Marcinkowski et al. | Calculating the Efficiency of Complex-Shaped Fins. Energies 2021, 14, 577 | |
Kuntysh et al. | Influence of the angle of inclination of round-finned tubes in a staggered tube bundle on the free convective heat exchange between it and an unbounded air space |
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 |