CN113591416A - 一种适用于汽轮机叶片流热固耦合分析方法 - Google Patents
一种适用于汽轮机叶片流热固耦合分析方法 Download PDFInfo
- Publication number
- CN113591416A CN113591416A CN202110915804.4A CN202110915804A CN113591416A CN 113591416 A CN113591416 A CN 113591416A CN 202110915804 A CN202110915804 A CN 202110915804A CN 113591416 A CN113591416 A CN 113591416A
- Authority
- CN
- China
- Prior art keywords
- solid
- domain
- calculation
- thermal
- fluid
- 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.)
- Pending
Links
- 239000007787 solid Substances 0.000 title claims abstract description 79
- 230000008878 coupling Effects 0.000 title claims abstract description 38
- 238000010168 coupling process Methods 0.000 title claims abstract description 38
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 38
- 238000004458 analytical method Methods 0.000 title claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 85
- 239000012530 fluid Substances 0.000 claims abstract description 51
- 238000000034 method Methods 0.000 claims abstract description 20
- 230000004907 flux Effects 0.000 claims abstract description 14
- 230000005540 biological transmission Effects 0.000 claims abstract description 11
- 239000000463 material Substances 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 229920001971 elastomer Polymers 0.000 claims description 5
- 239000000806 elastomer Substances 0.000 claims description 5
- 239000000853 adhesive Substances 0.000 claims description 4
- 230000001070 adhesive effect Effects 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 claims description 4
- SJRJDHOLKUXKIB-UHFFFAOYSA-N 1-[1-[4-(4-fluorophenyl)-4-oxobutyl]piperidin-4-yl]-3-methylbenzimidazol-2-one Chemical compound O=C1N(C)C2=CC=CC=C2N1C(CC1)CCN1CCCC(=O)C1=CC=C(F)C=C1 SJRJDHOLKUXKIB-UHFFFAOYSA-N 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 239000011343 solid material Substances 0.000 claims description 3
- 230000036760 body temperature Effects 0.000 claims description 2
- 229920001187 thermosetting polymer Polymers 0.000 claims 1
- 238000012821 model calculation Methods 0.000 abstract description 3
- 238000005059 solid analysis Methods 0.000 abstract description 2
- 238000002076 thermal analysis method Methods 0.000 abstract description 2
- 230000035882 stress Effects 0.000 description 25
- 238000004088 simulation Methods 0.000 description 3
- 230000008646 thermal stress Effects 0.000 description 2
- 238000005452 bending Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
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
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Control Of Turbines (AREA)
Abstract
本发明公开了一种适用于汽轮机叶片流热固耦合分析方法,本发明首先在所建立的模型上进行流体计算域与固体计算域间的热‑流耦合计算,为了实现流体域和固体域热通量和温度的守恒传递,二者之间采用热连接面进行数据传递,流体域和固体域内同步求解并实时数据传递。接着在所建立的模型上进行流体计算域与固体计算域间的热‑固耦合计算,其中热‑固分析以传递并施加在固体表面的CFD压力结果作为边界条件,并读入热分析得到的节点温度作为热载荷,同时施加相应的约束条件以及转速,从而实现了模型的综合应力分析,可极大程度的提高模型计算方法的真实性及计算结果可靠性。
Description
技术领域
本发明属于发电领域,具体涉及一种适用于汽轮机叶片流热固耦合分析方法。
背景技术
汽轮机进口参数的不断提高将导致前端叶片承受的热载荷不断增加,叶片长期运行在高温、高压和高转速环境下时,会诱使材料发生热疲劳和高温蠕变。而在诸如启停机、变负荷等瞬态运行过程中,高温部件内部热应力和热变形会进一步加剧并引起低周疲劳,进而影响高温区域叶片的运行可靠性和寿命。由于该区域的流动特性具有极其的三维特征,压力、温度以及速度场分布呈现极不均匀状态,因此准确分析高压高温环境下汽轮机叶片的流动特性、受力特征及应变情况,对高性能叶片开发设计及安全性分析具有非常重要的意义。
目前对高压高温环境下汽轮机叶片的模拟技术通常采用流动、传热及应力场无耦合独立仿真的计算方法,即先通过CFD软件进行流动计算获取气动边界条件后,再将结果单向传递至FEM计算软件进行后续的固体域温度场及应力场计算。这样的简化计算方式使得计算模型的边界条件在一定程度上与真实模型不符,模拟结果将产生不确定的插值误差。
发明内容
本发明的目的在于克服上述不足,提供一种适用于汽轮机叶片流热固耦合分析方法,可极大程度的提高模型计算方法的真实性及计算结果可靠性。
为了达到上述目的,本发明包括以下步骤:
S1,根据汽轮机的流通结构,建立叶片几何模型,并依据实际工况建立流体及固体计算域;
S2,分别对流体计算域及固体计算域进行网格刨分,生成多块结构化网格;
S3,确定流体及固体计算模型的边界条件;
S4,在叶片几何模型上进行流体计算域与固体计算域间的热-流耦合计算,对流体域和固体域之间的传热进行分析,流体域和固体域内同步求解并实时数据传递;
S5,在流体域中通过流动求解器数值求解雷诺时均Navier-Stokes方程组,并引入Boussinesq湍流模型假设使湍流计算雷诺时均Navier-Stokes方程组封闭,通过计算获取流体计算域压力和温度及流速;
S6,在固体域中只对热传导方程进行求解,在热连接面建立流体域和固体域之间的温度分布曲线,并通过热平衡方程进行反复的迭代求解以保证两侧热通量守恒传递,通过计算获取节点温度;
S7,在叶片几何模型上进行流体计算域与固体计算域间的热-固耦合计算,热-固耦合分析以传递并施加在固体表面的流体域压力结果作为边界条件,并读入S4中热-流耦合计算分析得到的节点温度作为热载荷,施加相应的约束条件以及转速,最终对模型进行综合应力分析。
S4中,固体域中对热传导方程进行求解,控制体离散热传导方程为:
S4中,流体域和固体域之间的热连接面采用完全非匹配热连接技术,由FNMB分解算法得到的单元进行热耦合方程求解得到每个单元独立的连接面温度和热通量,而将两侧单元上的面温度通过加权平均计算获得单元体温度。
S4中,热耦合方程如下:
其中,单元面完全未连接时Wi=1,当单元面部分连接时0<Wi<1,当完全连接时Wi=0。
S4中,S5中,通过流动求解器数值求解雷诺时均Navier-Stokes方程组,具体方法如下:
S4中,引入Boussinesq湍流模型假设使湍流计算雷诺时均Navier-Stokes方程组封闭,具体方法如下:
其中,μeff为有效粘性。
S7中,应力具有六个分量,其中三个正应力和三个剪切应力,其矩阵形式为:
S7中,应力的变化矩阵形式为:
S7中,弹性体应变由应力和温度变化两部分引起:
{ε}={εth}+{εel}
其中{εel}为应力引起的应变向量,{εt}为温度导致的应变向量,弹性体应变对于各项同性材料的表达式为:
{εth}=ΔT[a a a 0 0 0]T
式中:a为材料线膨胀系数,ΔT=T-Tref为当地温度与参考温度的温度差。
与现有技术相比,本发明首先在所建立的模型上进行流体计算域与固体计算域间的热-流耦合计算,为了实现流体域和固体域热通量和温度的守恒传递,二者之间采用热连接面进行数据传递,流体域和固体域内同步求解并实时数据传递。接着在所建立的模型上进行流体计算域与固体计算域间的热-固耦合计算,其中热-固分析以传递并施加在固体表面的CFD压力结果作为边界条件,并读入热分析得到的节点温度作为热载荷,同时施加相应的约束条件以及转速,从而实现了模型的综合应力分析,可极大程度的提高模型计算方法的真实性及计算结果可靠性。
附图说明
图1为本发明实施例的计算模型;
图2本发明实施例的叶片流体域网格;
图3本发明实施例的叶片固体域网格;
图4为本发明实施例的叶片温度场计算结果;
图5为本发明实施例的叶片应力场计算结果;其中,(a)为叶片离心应力,(b)为叶片弯扭应力,(c)为叶片热应力,(d)为叶片综合等效应力;
图6为本发明的流程图。
具体实施方式
下面结合附图对本发明做进一步说明。
参见图1至图6,本发明包括以下步骤:
步骤一,参照汽轮机真实通流结构,根据几何图纸通过三维建模软件建立叶片几何模型,并依据实际工况建立流体及固体计算域;
步骤二,采用商用软件分别对流体计算域及固体计算域进行网格刨分,生成多块结构化网格;
步骤三,根据物理实际确定流体及固体计算模型的边界条件;
步骤四,在所建立的模型上进行流体计算域与固体计算域间的热-流耦合计算,对流体域和体域之间的传热进行分析,为了实现流体域和固体域热通量和温度的守恒传递,二者之间采用热连接面进行数据传递,流体域和固体域内同步求解并实时数据传递。
固体传热分析计算中采用的格子中心有限体积法。由于固体内部没有流动现象,其传热方程可以简化为下式,用于求解基于流固耦合面热通量守恒而引起的给定时间步长内固体温度变化:
热连接面采用完全非匹配热连接技术,由FNMB分解算法得到的单元进行热耦合方程求解可以得到每个单元独立的连接面温度和热通量,而将两侧单元上的面温度通过加权平均计算可以进一步获得单元体温度。
式中当单元面完全未连接时Wi=1,当单元面部分连接时0<Wi<1,当完全连接时Wi=0。
步骤五,在流体域中通过流动求解器数值求解雷诺时均Navier-Stokes方程组,并引入Boussinesq湍流模型假设使湍流计算雷诺时均Navier-Stokes方程组封闭,通过计算获取流体计算域压力、温度及流速等重要气动参数;
通过流动求解器数值求解雷诺时均Navier-Stokes方程组,其具体形式如下:
并引入Boussinesq湍流模型假设,使湍流计算雷诺时均Navier-Stokes方程组封闭,其具体形式如下:
式中:μeff为有效粘性
步骤六,在固体域中只对热传导方程进行求解,在热连接面建立流体域和固体域之间的温度分布曲线,并通过热平衡方程进行反复的迭代求解以保证两侧热通量守恒传递,通过计算获取节点温度等边界参数。热传导方程形式如下:
hf(Tw-Tf)=-hs(Tw-Ts)
式中:hf为流体侧导热系数/W·(m·K)-1;hs为固体侧导热系数/W·(m·K)-1;Tw为连接面温度/K;Tf为流体侧第一层网格中心温度/K;Ts为固体侧第一层网格中心温度/K。
固体侧传热系数表达式为:
式中:κ为当地热导率/W·(m·K)-1;Δy为第一层网格中心到壁面距离/m。
而当第一层网格没有在粘性底层内部时需要采用对流热通量,此时流体域和固体域之间的热量传递适合采用结合壁面函数的高雷诺数湍流模型,其传热系数采用下式求解:
在粘性底层中(y+<13.2):
T+=Pr y+
在湍流层中:
步骤七,在所建立的模型上进行流体计算域与固体计算域间的热-固耦合计算,热-固耦合分析以传递并施加在固体表面的流体域压力结果作为边界条件,并读入步骤四热-流耦合计算分析得到的节点温度作为热载荷,进一步施加相应的约束条件以及转速,最终对模型进行综合应力分析。
三维问题中弹性体内任意点的应力具有6个分量:3个正应力和3个剪切应力,其矩阵形式为:
应变的矩阵形式为:
根据热弹性理论,弹性体应变由应力和温度变化两部分引起:
{ε}={εth}+{εel}
其中{εel}为应力引起的应变向量,{εth}为温度导致的应变向量,其对于各项同性材料的表达式为:
{εth}=ΔT[a a a 0 0 0]T
式中:a为材料线膨胀系数/m·K-1;ΔT=T-Tref为当地温度与参考温度的温度差。
Claims (9)
1.一种适用于汽轮机叶片流热固耦合分析方法,其特征在于,包括以下步骤:
S1,根据汽轮机的流通结构,建立叶片几何模型,并依据实际工况建立流体及固体计算域;
S2,分别对流体计算域及固体计算域进行网格刨分,生成多块结构化网格;
S3,确定流体及固体计算模型的边界条件;
S4,在叶片几何模型上进行流体计算域与固体计算域间的热-流耦合计算,对流体域和固体域之间的传热进行分析,流体域和固体域内同步求解并实时数据传递;
S5,在流体域中通过流动求解器数值求解雷诺时均Navier-Stokes方程组,并引入Boussinesq湍流模型假设使湍流计算雷诺时均Navier-Stokes方程组封闭,通过计算获取流体计算域压力和温度及流速;
S6,在固体域中只对热传导方程进行求解,在热连接面建立流体域和固体域之间的温度分布曲线,并通过热平衡方程进行反复的迭代求解以保证两侧热通量守恒传递,通过计算获取节点温度;
S7,在叶片几何模型上进行流体计算域与固体计算域间的热-固耦合计算,热-固耦合分析以传递并施加在固体表面的流体域压力结果作为边界条件,并读入S4中热-流耦合计算分析得到的节点温度作为热载荷,施加相应的约束条件以及转速,最终对模型进行综合应力分析。
3.根据权利要求1所述的一种适用于汽轮机叶片流热固耦合分析方法,其特征在于,S4中,流体域和固体域之间的热连接面采用完全非匹配热连接技术,由FNMB分解算法得到的单元进行热耦合方程求解得到每个单元独立的连接面温度和热通量,而将两侧单元上的面温度通过加权平均计算获得单元体温度。
9.根据权利要求1所述的一种适用于汽轮机叶片流热固耦合分析方法,其特征在于,S7中,弹性体应变由应力和温度变化两部分引起:
{ε}={εth}+{εel}
其中{εel}为应力引起的应变向量,{εt}为温度导致的应变向量,弹性体应变对于各项同性材料的表达式为:
{εth}=ΔT[a a a 0 0 0]T
式中:a为材料线膨胀系数,ΔT=T-Tref为当地温度与参考温度的温度差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110915804.4A CN113591416A (zh) | 2021-08-10 | 2021-08-10 | 一种适用于汽轮机叶片流热固耦合分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110915804.4A CN113591416A (zh) | 2021-08-10 | 2021-08-10 | 一种适用于汽轮机叶片流热固耦合分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113591416A true CN113591416A (zh) | 2021-11-02 |
Family
ID=78256942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110915804.4A Pending CN113591416A (zh) | 2021-08-10 | 2021-08-10 | 一种适用于汽轮机叶片流热固耦合分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113591416A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015215772A (ja) * | 2014-05-12 | 2015-12-03 | 株式会社東芝 | 伝熱シミュレーション装置および伝熱シミュレーション方法 |
US20170235863A1 (en) * | 2014-10-31 | 2017-08-17 | Siemens Aktiengesellschaft | Method for the prediction of turbomachine performances |
CN107506562A (zh) * | 2017-09-29 | 2017-12-22 | 西安科技大学 | 一种水润滑橡胶轴承双向热流固耦合计算方法 |
CN108984920A (zh) * | 2018-07-23 | 2018-12-11 | 安徽江淮纳威司达柴油发动机有限公司 | 一种发动机冷却水套的直接流固耦合传热分析方法 |
CN109063337A (zh) * | 2018-08-03 | 2018-12-21 | 内蒙古工业大学 | 一种永磁风力发电机多场-路耦合模拟计算方法 |
CN109858135A (zh) * | 2019-01-25 | 2019-06-07 | 西安热工研究院有限公司 | 一种汽轮机低压通流区长叶片安全性校核的计算方法 |
CN110032769A (zh) * | 2019-03-18 | 2019-07-19 | 北京工业大学 | 一种通风机叶片流固耦合模型 |
-
2021
- 2021-08-10 CN CN202110915804.4A patent/CN113591416A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015215772A (ja) * | 2014-05-12 | 2015-12-03 | 株式会社東芝 | 伝熱シミュレーション装置および伝熱シミュレーション方法 |
US20170235863A1 (en) * | 2014-10-31 | 2017-08-17 | Siemens Aktiengesellschaft | Method for the prediction of turbomachine performances |
CN107506562A (zh) * | 2017-09-29 | 2017-12-22 | 西安科技大学 | 一种水润滑橡胶轴承双向热流固耦合计算方法 |
CN108984920A (zh) * | 2018-07-23 | 2018-12-11 | 安徽江淮纳威司达柴油发动机有限公司 | 一种发动机冷却水套的直接流固耦合传热分析方法 |
CN109063337A (zh) * | 2018-08-03 | 2018-12-21 | 内蒙古工业大学 | 一种永磁风力发电机多场-路耦合模拟计算方法 |
CN109858135A (zh) * | 2019-01-25 | 2019-06-07 | 西安热工研究院有限公司 | 一种汽轮机低压通流区长叶片安全性校核的计算方法 |
CN110032769A (zh) * | 2019-03-18 | 2019-07-19 | 北京工业大学 | 一种通风机叶片流固耦合模型 |
Non-Patent Citations (1)
Title |
---|
李宇;邹正平;叶建;李维;周志翔;秦代斌;: "二维流热耦合数值模拟程序及方法研究", 燃气涡轮试验与研究, no. 02, 15 May 2007 (2007-05-15), pages 18 - 26 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109858135B (zh) | 一种汽轮机低压通流区长叶片安全性校核的计算方法 | |
Liu et al. | Aeroelastic analysis of a floating offshore wind turbine in platform‐induced surge motion using a fully coupled CFD‐MBD method | |
Lu et al. | Improved decomposed-coordinated kriging modeling strategy for dynamic probabilistic analysis of multicomponent structures | |
CN102799730B (zh) | 一种燃气轮机风扇叶片反扭过程的预估方法 | |
Li et al. | Multivariate ensembles-based hierarchical linkage strategy for system reliability evaluation of aeroengine cooling blades | |
Chi et al. | Coupled aerothermodynamics optimization for the cooling system of a turbine vane | |
Huang et al. | An optimal design for axial-flow fan blade: theoretical and experimental studies | |
Xiaodong et al. | Multidisciplinary design optimization of turbine disks based on ANSYS Workbench platforms | |
Li et al. | The influence of the bound vortex on the aerodynamics of curved wind turbine blades | |
Sørensen et al. | Aerodynamic effects of compressibility for wind turbines at high tip speeds | |
CN116090191A (zh) | 一种综合气流因素下的海上风机的仿真方法及系统 | |
CN113051846B (zh) | 考虑可压缩和热传导效应的壁面第一层网格厚度估算方法 | |
Jafari et al. | Advances in nonlinear stress analysis of a steam cooled gas turbine blade | |
CN113591416A (zh) | 一种适用于汽轮机叶片流热固耦合分析方法 | |
Shi et al. | Influences of wind and rotating speed on the fluid-structure interaction vibration for the offshore wind turbine blade | |
Zhang et al. | Experimental and numerical investigations of aerodynamic loads and 3D flow over non‐rotating MEXICO blades | |
Cao et al. | Study of air compressibility effects on the aerodynamic performance of the IEA-15 MW offshore wind turbine | |
Adibi et al. | Predicting airfoil stalling dynamics using upwind numerical solutions to non-viscous equations | |
CN116167220A (zh) | 一种综合气流因素下的海上风机的仿真方法及系统 | |
Campobasso et al. | Compressible Reynolds-Averaged Navier-Stokes analysis of wind turbine turbulent flows using a Fully-Coupled Low-Speed preconditioned multigrid solver | |
Yu et al. | A coupled CFD-CSD method for predicting HAWT rotor blade performance | |
Liou et al. | Flutter analysis for turbomachinery using Volterra series | |
Birch | Navier-Stokes predictions of transition, loss and heat transfer in a turbine cascade | |
Nadarajah et al. | Multi-objective aerodynamic shape optimization for unsteady viscous flows | |
de la Llave Plata et al. | elsA-Hybrid: an all-in-one structured/unstructured solver for the simulation of internal and external flows. Application to turbomachinery |
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 |