CN112883575A - 一种考虑表面粗糙度的叶轮机械边界层转捩模型修正方法 - Google Patents
一种考虑表面粗糙度的叶轮机械边界层转捩模型修正方法 Download PDFInfo
- Publication number
- CN112883575A CN112883575A CN202110211561.6A CN202110211561A CN112883575A CN 112883575 A CN112883575 A CN 112883575A CN 202110211561 A CN202110211561 A CN 202110211561A CN 112883575 A CN112883575 A CN 112883575A
- Authority
- CN
- China
- Prior art keywords
- transition
- boundary layer
- roughness
- blade
- reynolds number
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明涉及燃气涡轮发动机叶轮机械高精度数值仿真技术领域,涉及一种适用于真实叶轮机械叶片表面边界层计算的兼顾叶片表面粗糙度的转捩模型修正方法,结合表面粗糙度影响叶轮机械叶片表面边界层转捩的特点,基于现有适用于光滑表面叶片的γ‑Reθ转捩模型进行修正,通过对光滑表面的转捩雷诺数Rθt并利用锥形截面模化粗糙度模型进行粗糙度修正,针对叶轮机械叶片边界层转捩高精度数值仿真,以准确模拟燃气涡轮发动机使用寿命期间叶片表面粗糙度对叶片表面边界层转捩和传热等特征的影响。本发明具有计算简单,易于编程实现,通用性好,计算精度高等优点。
Description
技术领域
本发明涉及燃气涡轮发动机叶轮机械高精度数值仿真技术领域,涉及一种适用于叶轮机械叶片表面边界层转捩模型修正方法,具体涉及一种适用于真实叶轮机械叶片表面边界层计算的兼顾叶片表面粗糙度的转捩模型修正方法。
背景技术
众所周知,叶轮机械叶片表面粗糙度随燃气涡轮发动机服役时间而发生巨大变化。叶片表面粗糙度对边界层层流、转捩和湍流等流动状态具有显著影响。数据显示,随着粗糙度高度的逐渐增加,边界层转捩点位置在叶片吸力面和压力表面上都向前移动,进而影响叶轮机械气动热效率,并对燃气涡轮发动机整机效率和全寿命周期总成本产生重大影响。叶片表面粗糙度不仅影响气动损失,还会增加局部热负载,进而导致叶片寿命缩短。因此,对设计和研发人员来说,精确捕捉叶片表面粗糙度对边界层转捩和传热等特征影响,以便在叶轮机械设计阶段就考虑燃气涡轮发动机使用寿命期间叶片表面粗糙度对气动损失和传热特性影响,对支撑高性能长寿命燃气涡轮发动机研制具有重要意义。
叶轮机械叶片表面边界层经历层流、转捩和湍流三种状态,即便在高来流湍流度情况下,层流、转捩始终存在。因此,在叶轮机械雷诺时均数值计算(RANS)时,需要采用专门的转捩模型来预测叶片表面边界层转捩。目前普遍采用的Langtry和Menter的γ-Reθ转捩模型,该模型由两个输运方程组成,即转捩雷诺数(transition Reynolds number)和数值间歇因子(numerical intermittency)输运方程,以及三个相关函数。γ-Reθ转捩模型是针对光滑叶片试验数据发展起来的,在光滑叶轮机械叶片表面边界层转捩数值仿真上具有相当高准确性,而针对具有表面粗糙度的叶轮机械叶片表面边界层转捩模拟显然不适合。
发明内容
根据以上问题,结合表面粗糙度影响叶轮机械叶片表面边界层转捩的特点,针对叶轮机械叶片边界层转捩高精度数值仿真,本发明提出了一种考虑表面粗糙度的叶轮机械叶片表面边界层转捩模型修正方法,以准确模拟燃气涡轮发动机使用寿命期间叶片表面粗糙度对叶片表面边界层转捩和传热等特征的影响。
本发明为实现其技术目的所采用的技术方案为:
一种考虑表面粗糙度的叶轮机械叶片表面边界层转捩模型修正方法,所述修正方法基于光滑表面叶片的γ-Reθ边界层转捩模型进行修正,所述基于光滑表面叶片的γ-Reθ边界层转捩模型,包括用于评估边界层状态的间歇因子γ输运方程和转捩雷诺数输运方程,其中,
所述间歇因子γ输运方程如式(1)所示:
式中,γ为间歇因子,Pγ为间歇因子产生项,Eγ为消散源项,ρ为密度,t为时间,Uj为速度向量,xj为空间坐标向量,μ为层流黏性系数,μt为湍流黏性系数,σf为间歇方程系数,为转捩雷诺数,Pθt为转捩动量厚度雷诺数产生项,σθt为扩散系数,
并且其中,间歇因子产生项Pγ、转捩动量厚度雷诺数产生项Pθt定义如下:
Pγ=Flenhthca1ρS(γFonset)0.5(1-ce1γ)
式中,Flength为转捩长度控制项,Fonset为转捩位置控制项,ca1为间歇方程系数,S为应变率幅值,ce1为间歇方程系数,cθt为源项幅值控制系数,tP为时间尺度,Rθt为光滑表面转捩雷诺数,Fθt为自由流与边界层切换系数,
其特征在于,
对所述基于光滑表面叶片的γ-Reθ边界层转捩模型进行修正时,利用锥形截面模化粗糙度模型对光滑表面的转捩雷诺数Rθt进行粗糙度修正,具体为:
fTu=max[0.9;1.61+1.15exp(-Tu)]
在所述基于光滑表面叶片的γ-Reθ边界层转捩模型中,在层流边界层中将间歇因子γ设置为0,在转捩流动时将间歇因子γ从0逐渐升高到1,在湍流边界层中将间歇因子γ设置为1。
同现有技术相比,本发明的考虑表面粗糙度的叶轮机械边界层转捩模型修正方法具有以下特点:1)计算简单,易于编程实现:基于现有边界层转捩模型进行粗糙度修正;2)通用性好:该边界层转捩模型修正方法不仅适用于压气机、涡轮叶片转捩预测,也适用于气冷涡轮叶片传热特性预测;3)计算精度高:基于锥形截面粗糙度模型构建了叶轮机械叶片表面常见的粗糙度分布规律,同时考虑了随机粗糙度分布的影响,且修正模型是通过大量实验获得,计算精度高。
附图说明
图1为锥形截面粗糙度模型示意图,图中,d为锥形粗糙度粒直径,Aw为粗糙度粒展开面积,AF为粗糙度粒投影面积,VR为为粗糙度粒体积,AS为单个粗糙度粒作用平面面积,t1为粗糙度粒周向间距,t2为粗糙度粒流向间距,k为粗糙度高度。
图2为基于现有光滑表面转捩模型进行数值计算结果(CFD)与实验结果(EXP)对比示意图,工况条件为雷诺数R=250×103,湍流度Tu=4.6%。
图3为基于本发明的考虑粗糙度后的边界层转捩修正模型进行数值计算结果(CFD)与实验结果(EXP)对比示意图,工况条件为雷诺数R=250×103,湍流度Tu=4.6%。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参照图1~3,本发明的考虑表面粗糙度的叶轮机械叶片表面边界层转捩模型修正方法,基于现有适用于光滑表面叶片的γ-Reθ转捩模型进行修正,结合表面粗糙度影响叶轮机械叶片边界层转捩的特点,针对叶轮机械叶片边界层转捩高精度数值仿真,可准确模拟燃气涡轮发动机使用寿命期间叶片表面粗糙度对叶片表面边界层转捩和传热等特征影响。
现有适用于光滑表面叶片的γ-Reθ转捩模型需两个额外输运方程,其中第一个是间歇因子γ输运方程,主要作用是评估边界层状态,在层流边界层中将γ设置为0,在转捩流动时将其从0逐渐升高到1,然后在其他湍流边界层将其设置为1。方程式如下:
间歇因子产生项Pγ定义,
Pγ=Flenhthca1ρS(γFonset)0.5(1-ce1γ)
该间歇因子产生项决定边界层转捩,其中Flength和Fonset函数分别确定转捩的长度和位置。Fonset涉及了转捩促发的条件,它是基于Rθt相关函数,而该函数取决于另一个转捩雷诺数输运方式,其方程式如下
上述方程式(1)和(2)组成了叶轮机械光滑表面转捩模型,各方程式中,γ为间歇因子,Pγ为间歇因子产生项,Eγ为消散源项,ρ为密度,t为时间,Uj为速度向量,xj为空间坐标向量,μ为层流黏性系数,μt为湍流黏性系数,σf为间歇方程系数,为转捩雷诺数,Pθt为转捩动量厚度雷诺数产生项,σθt为扩散系数,Flength为转捩长度控制项,Fonset为转捩位置控制项,ca1为间歇方程系数,S为应变率幅值,ce1为间歇方程系数,cθt为源项幅值控制系数,tP为时间尺度,Rθt为光滑表面转捩雷诺数,Fθt为自由流与边界层切换系数。
fTu=max[0.9;1.61+1.15exp(-Tu)]
其中,湍流度Tu用百分比表示,δ1是边界层位移厚度。Rθt,smooth为光滑表面转捩雷诺数,k为粗糙度高度,δ1为边界层位移厚度,fTu为有效湍流度,Tu为湍流度,用百分比表示,fΛ为粗糙度形状函数,ΛR为间距参数,且hm为平均粗糙度。
图2和3分别对比了在雷诺数R=250×103,湍流度Tu=4.6%的工况条件下,粗糙度涡轮叶片传热特性努塞尔数试验结果与现有光滑表面转捩模型和本发明的粗糙度转捩修正模型数值计算结果。结果表明,采用本发明的粗糙度转捩修正模型的计算结果与试验一致。
通过上述实施例,完全有效地实现了本发明的目的。该领域的技术人员可以理解本发明包括但不限于附图和以上具体实施方式中描述的内容。虽然本发明已就目前认为最为实用且优选的实施例进行说明,但应知道,本发明并不限于所公开的实施例,任何不偏离本发明的功能和结构原理的修改都将包括在权利要求书的范围中。
Claims (2)
1.一种考虑表面粗糙度的叶轮机械叶片表面边界层转捩模型修正方法,所述修正方法基于光滑表面叶片的γ-Reθ边界层转捩模型进行修正,所述基于光滑表面叶片的γ-Reθ边界层转捩模型,包括用于评估边界层状态的间歇因子γ输运方程和转捩雷诺数输运方程,其中,
所述间歇因子γ输运方程如式(1)所示:
式中,γ为间歇因子,Pγ为间歇因子产生项,Eγ为消散源项,ρ为密度,t为时间,Uj为速度向量,xj为空间坐标向量,μ为层流黏性系数,μt为湍流黏性系数,σf为间歇方程系数,为转捩雷诺数,Pθt为转捩动量厚度雷诺数产生项,σθt为扩散系数,
并且其中,间歇因子产生项Pγ、转捩动量厚度雷诺数产生项Pθt定义如下:
Pγ=Flengthca1ρS(γFonset)0.5(1-ce1γ)
式中,Flength为转捩长度控制项,Fonset为转捩位置控制项,ca1为间歇方程系数,S为应变率幅值,ce1为间歇方程系数,cθt为源项幅值控制系数,tP为时间尺度,Rθt为光滑表面转捩雷诺数,Fθt为自由流与边界层切换系数,
其特征在于,
对所述基于光滑表面叶片的γ-Reθ边界层转捩模型进行修正时,利用锥形截面模化粗糙度模型对光滑表面的转捩雷诺数Rθt进行粗糙度修正,具体为:
fTu=max[0.9;1.61+1.15exp(-Tu)]
2.根据权利要求1所述的修正方法,其特征在于,在所述基于光滑表面叶片的γ-Reθ边界层转捩模型中,在层流边界层中将间歇因子γ设置为0,在转捩流动时将间歇因子γ从0逐渐升高到1,在湍流边界层中将间歇因子γ设置为1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110211561.6A CN112883575B (zh) | 2021-02-25 | 2021-02-25 | 一种考虑表面粗糙度的叶轮机械边界层转捩模型修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110211561.6A CN112883575B (zh) | 2021-02-25 | 2021-02-25 | 一种考虑表面粗糙度的叶轮机械边界层转捩模型修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112883575A true CN112883575A (zh) | 2021-06-01 |
CN112883575B CN112883575B (zh) | 2023-10-03 |
Family
ID=76054520
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110211561.6A Active CN112883575B (zh) | 2021-02-25 | 2021-02-25 | 一种考虑表面粗糙度的叶轮机械边界层转捩模型修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112883575B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113959728A (zh) * | 2021-10-20 | 2022-01-21 | 中国科学院工程热物理研究所 | 基于温降热成像的叶轮机械叶片表面边界层转捩测试方法 |
CN114004404A (zh) * | 2021-11-03 | 2022-02-01 | 中南大学 | 一种边界层分布式粗糙元诱导转捩预测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108287054A (zh) * | 2017-12-25 | 2018-07-17 | 中国航天空气动力技术研究院 | 一种飞行条件下的转捩雷诺数获取方法 |
CN109033525A (zh) * | 2018-06-27 | 2018-12-18 | 浙江大学 | 一种基于简化三方程转捩模型的高超声速转捩预测方法 |
CN110702356A (zh) * | 2019-10-12 | 2020-01-17 | 空气动力学国家重点实验室 | 一种考虑表面粗糙度效应的高超声速横流转捩预测方法 |
-
2021
- 2021-02-25 CN CN202110211561.6A patent/CN112883575B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108287054A (zh) * | 2017-12-25 | 2018-07-17 | 中国航天空气动力技术研究院 | 一种飞行条件下的转捩雷诺数获取方法 |
CN109033525A (zh) * | 2018-06-27 | 2018-12-18 | 浙江大学 | 一种基于简化三方程转捩模型的高超声速转捩预测方法 |
CN110702356A (zh) * | 2019-10-12 | 2020-01-17 | 空气动力学国家重点实验室 | 一种考虑表面粗糙度效应的高超声速横流转捩预测方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113959728A (zh) * | 2021-10-20 | 2022-01-21 | 中国科学院工程热物理研究所 | 基于温降热成像的叶轮机械叶片表面边界层转捩测试方法 |
CN113959728B (zh) * | 2021-10-20 | 2023-12-22 | 中国科学院工程热物理研究所 | 基于温降热成像的叶轮机械叶片表面边界层转捩测试方法 |
CN114004404A (zh) * | 2021-11-03 | 2022-02-01 | 中南大学 | 一种边界层分布式粗糙元诱导转捩预测方法 |
CN114004404B (zh) * | 2021-11-03 | 2024-06-07 | 中南大学 | 一种边界层分布式粗糙元诱导转捩预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112883575B (zh) | 2023-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210209264A1 (en) | Modeling and calculation aerodynamic performances of multi-stage transonic axial compressors | |
Mayle | The 1991 IGTI scholar lecture: the role of laminar-turbulent transition in gas turbine engines | |
CN112883575A (zh) | 一种考虑表面粗糙度的叶轮机械边界层转捩模型修正方法 | |
Wagner et al. | Inlet boundary layer effects in an axial compressor rotor: part I—blade-to-blade effects | |
CN109815624A (zh) | 一种考虑进气总压畸变影响的压气机稳定边界判断方法 | |
Mahmood et al. | Flow characteristics of an optimized axial compressor rotor using smooth design parameters | |
Diener et al. | Multi-Disciplinary optimization of a mixed-flow compressor impeller | |
Chen et al. | Influence of local surface roughness of rotor blade on performance of an axial compressor stage | |
Zheng et al. | Effects of Reynolds number on the performance of a high pressure-ratio turbocharger compressor | |
Zangeneh et al. | Optimization of 6.2: 1 pressure ratio centrifugal compressor impeller by 3D inverse design | |
Piggush et al. | Flow measurements in a first stage nozzle cascade having endwall contouring, leakage and assembly features | |
Feng et al. | Aerothermodynamic design and numerical simulation of radial inflow turbine impeller for a 100kW microturbine | |
Dickens et al. | The design of highly loaded axial compressors | |
Zhao et al. | A novel clocking effect between inlet bend and volute in an automotive turbocharging system | |
An et al. | Internal Flow Phenomena of Two-Way Contra-Rotating Axial Flow Pump-Turbine in Pump Mode under Variable Speed | |
Cheng et al. | Influence of surface roughness on a highly loaded axial compressor stage performance at low Reynolds number | |
Flegel et al. | Aerodynamic effects of high turbulence intensity on a variable-speed power-turbine blade with large incidence and Reynolds number variations | |
Li et al. | Development and application of a throughflow method for high-loaded axial flow compressors | |
Fei et al. | Application of new empirical models based on mathematical statistics in the through-flow analysis | |
Zamiri et al. | Influence of the inclined leading edge diffuser vanes on the aerodynamic performance of a transonic centrifugal compressor | |
Ding et al. | Optimization design of agricultural fans based on skewed-swept blade technology | |
Mårtensson et al. | Design of a sub-scale fan for a boundary layer ingestion test with by-pass flow | |
Segawa et al. | Development of a highly loaded rotor blade for steam turbines (2nd Report, three-dimensional stage performance verifications) | |
Abbasi et al. | Investigation of inlet distortion effects on axial compressor performance based on streamline curvature method | |
Li et al. | Surrogate model on the extension operation range of an isolated centrifugal fan |
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 |