CN112163691A - 基于标量输运方程的气膜冷却二维有效度预测方法及系统 - Google Patents
基于标量输运方程的气膜冷却二维有效度预测方法及系统 Download PDFInfo
- Publication number
- CN112163691A CN112163691A CN202010861994.1A CN202010861994A CN112163691A CN 112163691 A CN112163691 A CN 112163691A CN 202010861994 A CN202010861994 A CN 202010861994A CN 112163691 A CN112163691 A CN 112163691A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- effectiveness
- equation
- scalar
- prediction
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- 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)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Business, Economics & Management (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Pure & Applied Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Computational Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Game Theory and Decision Science (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Evolutionary Computation (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Algebra (AREA)
- Quality & Reliability (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明公开了一种基于标量输运方程的气膜冷却二维有效度预测方法及系统,其中,该方法包括:构建有效度的二维标量输运方程,获取气膜孔几何参数和主流来流气动参数作为其输入参数;基于边界层理论,通过量级分析将二维标量输运方程抛物化;基于空间推进方法求解单孔排的二维有效度分布,引入二维Sellers叠加模型,将单孔排的二维有效度预测模型推广到多孔排预测,得到全覆盖气膜有效度的二维分布结果。该方法基于标量输运方程将一维有效度预测方法推广到二维,实现了气膜有效度的二维高精度预测。
Description
技术领域
本发明涉及燃气轮机气膜冷却技术领域,特别涉及一种基于标量输运方程的气膜冷却二维有效度预测方法及系统。
背景技术
燃气轮机气膜冷却是一种广泛采用的冷却技术。如图1所示,从压气机侧抽气,不经燃烧室,在透平部件的表面通过离散的气膜孔注入高温的主流边界层内部,形成冷却气膜,隔绝部件表面与高温燃气,防止透平部件的烧蚀。
气膜冷却区别于传统的槽缝冷却,具有全三维的流动结构,尤其是在透平端区,存在着复杂的冷气迁移机制,气膜冷却有效度在壁面上的分布表现为显著的二维特性。
但现有的气膜冷却有效度的预测方法,是基于横向平均意义上的一维有效度预测方法。针对单孔排气膜冷却,基于数学方法拟合流向方向上的横向平均气膜有效度分布,并通过一维Sellers叠加模型,将单孔排的结果推广至多孔排。
但现有一维气膜有效度预测方法难以满足三维气膜冷却的设计需求,如:
(1)气膜冷却在壁面上表现出二维的有效度分布特征,现有的一维预测方法难以精确定位部件表面的局部超温现象;
(2)针对多孔排气膜冷却,孔排间在周向方向上存在顺列、错列的关系,不同的周向位错也影响孔排下游有效度的实际分布规律。然而,现有的一维预测方法针对不同的孔排位错关系,预测的结果均相同。
因此,亟待一种针对全覆盖气膜冷却的二维有效度预测方法,解决燃气轮机气膜冷却有效度的性能预测问题。
发明内容
本发明提供一种基于标量输运方程的气膜冷却二维有效度预测方法及系统,以用于解决现有的一维预测方法难以精确定位部件表面的局部超温现象、针对不同的孔排位错关系,预测的结果均相同,不能实现全覆盖气膜冷却的高效、可靠预测的技术问题。
本发明一方面提供一种基于标量输运方程的气膜冷却二维有效度预测方法,包括:
构建控制方程,获取气膜孔几何参数和主流来流气动参数作为所述控制方程的输入参数,其中,所述控制方程为二维标量输运方程;
基于边界层理论,通过量级分析将所述二维标量输运方程抛物化;以及
基于空间推进方法求解单孔排的二维有效度分布,引入二维Sellers叠加模型,将单孔排的二维有效度预测模型推广到多孔排预测,得到全覆盖气膜有效度的二维分布结果。
本发明另一方面提供一种基于标量输运方程的气膜冷却二维有效度预测系统,包括:
构建模块,用于构建控制方程,获取气膜孔几何参数和主流来流气动参数作为所述控制方程的输入参数,其中,所述控制方程为二维标量输运方程;
抛物化模块,用于基于边界层理论,通过量级分析将所述二维标量输运方程抛物化;
预测模块,用于基于抛物化的二维标量输运方程,预测得到单孔排气膜有效度的二维分布,基于二维Sellers叠加公式,将单孔排结果推广得到多孔排的全覆盖气膜有效度的二维分布结果。
根据本发明实施例提出的基于标量输运方程的气膜冷却二维有效度预测方法及系统,基于标量输运方程将一维气膜冷却有效度预测方法拓展到了二维,实现了全覆盖气膜冷却的高效、可靠预测,能够广泛应用于重型燃气轮机、航空发动机气膜冷却的性能预测和优化设计。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为气膜冷却原理图;
图2为根据本发明一个实施例的基于标量输运方程的气膜冷却二维有效度预测方法流程图;
图3为根据本发明一个实施例的基于标量输运方程的气膜冷却二维有效度预测方法具体执行流程图;
图4为根据本发明一个实施例的基于标量输运方程的气膜冷却二维有效度预测系统的结构示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
下面参照附图描述根据本发明实施例提出的基于标量输运方程的气膜冷却二维有效度预测方法及系统,首先将参照附图描述根据本发明实施例提出的基于标量输运方程的气膜冷却二维有效度预测方法。
图2是本发明一个实施例的基于标量输运方程的气膜冷却二维有效度预测方法。
如图2所示,该基于标量输运方程的气膜冷却二维有效度预测方法包括以下步骤:
在步骤S201中,构建控制方程,获取气膜孔几何参数和主流来流气动参数作为控制方程的输入参数,其中,控制方程为二维标量输运方程。
具体地,如图3所示,控制方程为气膜有效度的二维标量输运方程:
输入参数为气膜孔几何参数和主流来流气动参数,作为控制方程的初值条件和边界条件,其中,η*为归一化的气膜有效度,定义为:其中,为横向平均的气膜有效度分布,t为时间,Uj为输运速度,x为笛卡尔坐标流向方向上的坐标分量,y为笛卡尔坐标横向方向上的坐标分量,xj为笛卡尔坐标系的三个坐标分量,κ为分子扩散系数。
进一步地,基于涡粘假设,将右端第二项由于湍流脉动所造成的标量的扩散写成关于平均标量场的梯度形式:
其中,κt为涡扩散系数,因此标量输运方程可改写成统一形式,综合考虑分子粘性和湍流脉动造成的标量的扩散。
在步骤S202中,基于边界层理论,通过量级分析将二维标量输运方程抛物化。
需要说明的是,本发明实施例同时考虑在射流和主流交互作用下的对流和扩散效应,预先对二维标量输运方程进行量级分析。具体地,如图2所示,扩散效应与主流来流湍流度有关,综合考虑分子扩散和湍流扩散,引入有效扩散系数κeff模化两者的综合作用效果,可以将有效扩散系数表示为来流湍流度Tu和射流速比VR的函数:
在边界层假设下,对流项与扩散项同量级,通过壁面的极限流线模化气膜孔出口肾型涡的卷吸作用:
其中,κeff,0为特征的有效扩散系数,Tu为主流进口湍流度,VR为射流速度比,Uy为横向方向的输运速度,α为壁面极限流线与轴向方向的夹角。另外,公式中变量上带~表示通过主流来流速度U∞,气膜孔孔径D进行无量纲化之后的无量纲量。
进一步地,通过量级分析后,将标量输运方程抛物化,采用空间推进方法求解,大大节省了计算成本,控制方程的最终形式(即抛物化后的二维标量输运方程)为:
在步骤S203中,基于空间推进方法求解单孔排的二维有效度分布,引入二维Sellers叠加模型,将单孔排的二维有效度预测模型推广到多孔排预测,得到全覆盖气膜有效度的二维分布结果。
也就是说,如图3所示,将控制方程的最终形式引入二维Sellers叠加模型,将单孔排的二维有效度预测模型推广到多孔排预测,二维多孔排气膜有效度的预测方法为:
其中,m为多排孔multi-row的缩写,i为第i排孔,n为共计n排孔,s为单排孔single-row的缩写;
再根据多孔排预测模型预测全覆盖气膜有效度的二维分布结果。
本发明实施例提出的基于标量输运方程的气膜冷却二维有效度预测方法,至少实现了如下有益的技术效果:
1)预测精度高,本发明实施例基于实验数据拟合得到的二维分布结果,在吹风比M≤3.0时,有效度横向分布的拟合优度R2满足在0.98以上;2)气膜有效度的二维模型,本发明实施例能够精确预测透平叶片表面气膜有效度的二维分布,在初步设计过程中,精确定位可能发生局部烧蚀的区域;3)计算代价低,本发明实施例通过对控制方程抛物化,能够基于空间推进的方法进行数值求解,多孔排有效度二维预测的计算时间能够有效地控制在0.1s内;4)能够考虑孔排间位错的影响,本发明实施例有效克服现有一维Sellers公式的横向均匀分布假设条件,消除孔排间位错带来的预测误差,进而有效提高多孔排的预测精度。
综上,根据本发明实施例提出的基于标量输运方程的气膜冷却二维有效度预测方法,基于标量输运方程将一维气膜冷却有效度预测方法拓展到了二维,实现了全覆盖气膜冷却的高效、可靠预测,能够广泛应用于重型燃气轮机、航空发动机气膜冷却的性能预测和优化设计。
其次参照附图描述根据本发明实施例提出的基于标量输运方程的气膜冷却二维有效度预测系统。
图4是本发明一个实施例的基于标量输运方程的气膜冷却二维有效度预测系统的结构示意图。
如图4所示,该系统10包括:构建模块100、抛物化模块200和预测模块300。
其中,构建模块100,用于构建控制方程,获取气膜孔几何参数和主流来流气动参数作为控制方程的输入参数,其中,控制方程为二维标量输运方程。抛物化模块200,用于基于边界层理论,通过量级分析将二维标量输运方程抛物化。预测模块300,用于基于抛物化的二维标量输运方程,预测得到单孔排气膜有效度的二维分布,基于二维Sellers叠加公式,将单孔排结果推广得到多孔排的全覆盖气膜有效度的二维分布结果。
根据本发明实施例提出的基于标量输运方程的气膜冷却二维有效度预测系统,基于标量输运方程将一维气膜冷却有效度预测方法拓展到了二维,实现了全覆盖气膜冷却的高效、可靠预测,能够广泛应用于重型燃气轮机、航空发动机气膜冷却的性能预测和优化设计。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (6)
1.一种基于标量输运方程的气膜冷却二维有效度预测方法,其特征在于,包括以下步骤:
构建控制方程,获取气膜孔几何参数和主流来流气动参数作为所述控制方程的输入参数,其中,所述控制方程为二维标量输运方程;
基于边界层理论,通过量级分析将所述二维标量输运方程抛物化;以及
基于空间推进方法求解单孔排的二维有效度分布,引入二维Sellers叠加模型,将单孔排的二维有效度预测模型推广到多孔排预测,得到全覆盖气膜有效度的二维分布结果。
6.一种基于标量输运方程的气膜冷却二维有效度预测系统,其特征在于,包括:
构建模块,用于构建控制方程,获取气膜孔几何参数和主流来流气动参数作为所述控制方程的输入参数,其中,所述控制方程为二维标量输运方程;
抛物化模块,用于基于边界层理论,通过量级分析将所述二维标量输运方程抛物化;
预测模块,用于基于抛物化的二维标量输运方程,预测得到单孔排气膜有效度的二维分布,基于二维Sellers叠加公式,将单孔排结果推广得到多孔排的全覆盖气膜有效度的二维分布结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010861994.1A CN112163691B (zh) | 2020-08-25 | 2020-08-25 | 基于标量输运方程的气膜冷却二维有效度预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010861994.1A CN112163691B (zh) | 2020-08-25 | 2020-08-25 | 基于标量输运方程的气膜冷却二维有效度预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112163691A true CN112163691A (zh) | 2021-01-01 |
CN112163691B CN112163691B (zh) | 2022-05-17 |
Family
ID=73859797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010861994.1A Active CN112163691B (zh) | 2020-08-25 | 2020-08-25 | 基于标量输运方程的气膜冷却二维有效度预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112163691B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112683743A (zh) * | 2021-03-12 | 2021-04-20 | 杭州晟视科技有限公司 | 血小板状态分析方法、装置、电子设备和存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108416075A (zh) * | 2017-08-29 | 2018-08-17 | 沈阳工业大学 | 基于cfd技术的风力机气动计算模型建模方法 |
CN110032756A (zh) * | 2019-02-27 | 2019-07-19 | 上海交通大学 | 基于流函数分数坐标系变换的流动边界层数值分析方法 |
-
2020
- 2020-08-25 CN CN202010861994.1A patent/CN112163691B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108416075A (zh) * | 2017-08-29 | 2018-08-17 | 沈阳工业大学 | 基于cfd技术的风力机气动计算模型建模方法 |
CN110032756A (zh) * | 2019-02-27 | 2019-07-19 | 上海交通大学 | 基于流函数分数坐标系变换的流动边界层数值分析方法 |
Non-Patent Citations (2)
Title |
---|
LANG WANG: "EFFECT OF ROW SPACING ON THE ACCURACY OF FILM COOLING SUPERPOSITION METHOD", 《PROCEEDINGS OF ASME TURBO EXPO 2018》 * |
陈子聿,苏欣荣,袁新: "气冷透平端区射流与主流的交互作用研究", 《风机技术》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112683743A (zh) * | 2021-03-12 | 2021-04-20 | 杭州晟视科技有限公司 | 血小板状态分析方法、装置、电子设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN112163691B (zh) | 2022-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liou et al. | Fluid flow in a 180 deg sharp turning duct with different divider thicknesses | |
CN110083968B (zh) | 基于修正气封泄露量影响数值模型的压气机特性预测方法 | |
JP2019535987A (ja) | 構造的熱交換器 | |
CN112163691B (zh) | 基于标量输运方程的气膜冷却二维有效度预测方法及系统 | |
Jain et al. | CFD modeling of a radiator axial fan for air flow distribution | |
Noman Danish et al. | Performance evaluation of tandem bladed centrifugal compressor | |
CN113095552A (zh) | 一种海洋平台泄漏可燃气云团体积预测方法及系统 | |
CN115982892B (zh) | 叶片设计方法、叶片及相关设备 | |
CN114297807B (zh) | 一种防冰系统笛形管结构参数优化方法 | |
Lin et al. | Effective boundary conditions and numerical method for flow characteristics of aeroengine compressor at high Mach flight | |
CN114282323B (zh) | 一种涡轮叶片层板冷却结构的流量分配预测方法 | |
Madadi et al. | Quasi-3D inverse design of S-shaped diffuser with specified cross-section distribution; super-ellipse, egg-shaped, and ellipse | |
CN113901602A (zh) | 基于分区相似的多孔隔热屏壁面温度计算方法 | |
Liao et al. | Effects of curvature on the film cooling effectiveness of double-jet film cooling | |
Kim et al. | Shape optimization of a bended film-cooling hole to enhance cooling effectiveness | |
Il’inkov et al. | Algorithm of numerical checking calculation of turbine cooled blades of gas turbine engines using experimental data on heat transfer and resistance | |
CN113722900B (zh) | 一种航空换热器非设计点性能设计与分析方法 | |
Abbasi et al. | Investigation of inlet distortion effects on axial compressor performance based on streamline curvature method | |
Moghimi EsfandAbadi et al. | Design and Analysis of Supersonic Inlet for Ramjet Engines: Aerodynamic Considerations and Performance Optimization | |
Sewall et al. | Large eddy simulation of flow and heat transfer in the developing flow region of a rotating gas turbine blade internal cooling duct with coriolis and buoyancy forces | |
CN113065292B (zh) | 用于一体化涡轮过渡段非定常流动计算的边界模型方法 | |
CN117688698B (zh) | 涡轮叶片冷却结构多学科设计方法及装置 | |
Kucuk | Application of response surface methodology to optimize aerodynamic performance of NACA inlet | |
Kurzrock et al. | Transonic flow around rotor blade elements | |
CN115773891A (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 |