CN114880764B - 部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法 - Google Patents

部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法 Download PDF

Info

Publication number
CN114880764B
CN114880764B CN202210413920.0A CN202210413920A CN114880764B CN 114880764 B CN114880764 B CN 114880764B CN 202210413920 A CN202210413920 A CN 202210413920A CN 114880764 B CN114880764 B CN 114880764B
Authority
CN
China
Prior art keywords
turbine
domain
temperature
coupling
solid
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
Application number
CN202210413920.0A
Other languages
English (en)
Other versions
CN114880764A (zh
Inventor
罗凯
智若阳
秦侃
王瀚伟
党建军
黄闯
李代金
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202210413920.0A priority Critical patent/CN114880764B/zh
Publication of CN114880764A publication Critical patent/CN114880764A/zh
Application granted granted Critical
Publication of CN114880764B publication Critical patent/CN114880764B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • F01D15/00—Adaptations of machines or engines for special use; Combinations of engines with devices driven thereby
    • F01D15/02—Adaptations for driving vehicles, e.g. locomotives
    • 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
    • F01D25/08—Cooling; Heating; Heat-insulation
    • F01D25/12—Cooling
    • 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • 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
    • 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 Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Mechanical Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Turbines (AREA)
  • Turbine Rotor Nozzle Sealing (AREA)

Abstract

本发明公开了部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，通过建模、对模型的网格划分，数据模拟分析，得到受轮盘变形影响的涡轮机气动性能和应力应变，本发明的方法适用于水下部分进气轴向冲动式涡轮发动机，为喷水冷却涡轮机的性能提供了直观的数据支撑，能够较好的解决水下涡轮机部分进气结构对涡轮机气热弹耦合仿真造成的影响，并能判断在喷水冷却下涡轮盘的可靠性。

Description

部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法
技术领域
本发明属于工程热力学技术领域，具体涉及部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法。
背景技术
目前，服役于世界各国海军的潜艇及水面舰艇性能大为改善，特别是航速有了较大的提高，为了适应未来海战的需要，用于反潜反舰的鱼雷急需增大其航速、航程及航深，因此迫切需要提高鱼雷动力推进***的性能。采用涡轮机的热动力***具有比功率大、体积小、使用维护成本低和寿命周期长等优点，因此在动力推进***功率达到上千千瓦的情况下，鱼雷动力推进***发动机的首选方案是涡轮机。
因为涡轮机尺寸受限于水下航行器的空间，所以为了提高涡轮机的性能，其转速通常会达到上万转，现有应用于水下航行器的涡轮机转速可到50000rpm，轮盘根部接近材料屈服强度极限，难以继续提高转速以增加涡轮机效率和功率。同时，提高涡轮前燃气入口温度是提供较大可用焓降，提高涡轮机有效功率的手段之一，但因燃气温度的提高，涡轮盘热负荷过大，在大转速运转下易发生结构的龟裂变形，影响涡轮机安全与稳定工作，如今水下航行器涡轮机入口温度达到1100℃，显然涡轮机入口温度已经超过了金属材料所能承受的极限，所以需对涡轮盘采取一定的冷却措施来确保涡轮机的安全稳定运行。现有涡轮机常用的冷却方式为气体冷却且常用在航空涡轮发动机中，其需要在叶片内部布置复杂的冷却通道，与航空发动机不同的是，水下涡轮机的短叶片以及尺寸小的特点使其无法采用复杂的冷却结构，同时在水下环境中缺少足够持续的冷却气体，而且液态水相比于空气具有更加优越的冷却性能，因此可以使用充足的海水作为冷却介质。与纯气动涡轮机相比，喷水冷却的燃气涡轮机内部流动包含复杂的多相流动、多物理场耦合现象，包括超音速高温燃气与低温液态水的复杂能量和动量交换现象，涡轮盘及叶片与掺混的燃气、液态水及水蒸气的气热弹耦合现象，轮盘的温度和热应力分布直接受流体域的影响。开展具有高速旋转轮盘、超音速流动以及高温高压特点的涡轮机试验时非常困难的，此时仿真方法就成为了一种不可或缺的研究手段。
近年来，基于气热耦合方法进行涡轮叶片流动和传热特性的研究已经在国内外叶轮机械领域被广泛采用，而现有的气热弹耦合仿真方法主要集中在全周进气涡轮机。与全周进气涡轮机不同的是，具有部分进气结构的水下涡轮机的流动过程具有强瞬态性，若要准确使用仿真方法预测轮盘温度必须在计算中考虑旋转效应，现有的气热弹耦合仿真计算方法在求解过程中不适用于水下部分进气涡轮机，其计算结果往往未考虑旋转效应，与实际情况偏差很大，而对于部分进气涡轮机的温度场和应力场求解正确与否直接决定涡轮机的安全稳定运行，因此需要一种高精度、低计算成本的数值仿真方法模拟喷水冷却下部分进气涡轮机的气热弹耦合过程，以获得喷水冷却涡轮机的性能。
发明内容
本发明的目的是提供部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，适用于水下部分进气轴向冲动式涡轮发动机，为喷水冷却涡轮机的性能提供了直观的数据支撑，能够较好的解决水下涡轮机部分进气结构对涡轮机气热弹耦合仿真造成的影响，并能判断在喷水冷却下涡轮盘的可靠性。
本发明所采用的技术方案是，部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，具体按照以下步骤实施：
步骤1、建立水下涡轮机流体域和固体域的三维模型；
步骤2、利用多重分块结构化网格方法对三维模型进行网格划分；
步骤3、采用流场仿真计算软件，建立具有超音速流场、多相流、传质相变的流场数值模拟仿真计算方法，设置喷水冷却涡轮机流体域耦合面的温度初值，对流体域进行数值模拟，得到流体域耦合面上的热流密度和近壁面温度；
步骤4、考虑轮盘的旋转效应，将各个耦合面上的热流密度和近壁面温度沿周向进行平均后计算得到耦合面上的对流换热系数；
步骤5、建立固体域导热模型，将各个耦合面上的对流换热系数和近壁面温度作为固体域的边界条件对固体域进行数值模拟，得到固体域耦合面的壁面温度；
步骤6、将固体域耦合面的壁面温度进行松弛后，作为流体域的耦合面壁面温度，再次对流体域进行数值模拟，执行步骤3～5，得到新的固体域耦合面的壁面温度，将新的固体域耦合面的壁面温度与上一循环步骤5中固体域耦合面的壁面温度相比，判断壁面温度是否收敛，若收敛，则结束循环，若不收敛，则返回步骤3，直到壁面温度收敛，得到喷水冷却涡轮机涡轮盘的温度分布；
步骤7、由最后一次循环中流体域数值模拟的结果得到气动载荷，将温度载荷、气动载荷及离心载荷加载到涡轮盘上对轮盘进行有限元计算分析，得到喷水冷却涡轮机涡轮盘的总应力和总应变；
步骤8、根据喷水冷却涡轮机涡轮盘的总应力和总应变建立变形后的涡轮机流体域和固体域模型，将变形后的涡轮机流体域和固体域模型作为三维模型，执行步骤2～7，得到受轮盘变形影响的喷水冷却涡轮机涡轮盘的总应力和总应变。
本发明的特点还在于：
步骤1具体过程为：根据实际工程中的水下涡轮机三维模型参数，利用三维建模软件进行实体几何造型的构造，得到固体域实体模型，并反解出流体域实体模型，形成水下涡轮机流体域和固体域的三维模型。
三维模型的流体域分为喷管流道域、转子叶栅域、叶顶间隙域和轴向间隙域，在喷管流道域、轴向间隙流域采用O型网格进行加密，转子叶栅域处则采用H型block进行加密，喷管流域在靠近喉部位置沿流动方向进行加密。
步骤3具体过程为：对喷水冷却涡轮机的流体域进行仿真，其中流体域的入口为压力入口，总温总压为实际进气总温总压，出口为压力出口，静压为实际出口静压，轴向和径向喷水口设置为质量流量入口，喷入的冷却水的流量和温度由实际工程决定，采用3D定常压力基求解器，双精度，使用COUPLED算法，燃气和水蒸气为理想可压气体并打开组分输运模型考虑质量和热量的扩散，液态水为不可压液体，液体和气体之间采用VOF模型进行仿真，湍流模型采用k-ωSST模型，与固体域相耦合的壁面设置一个初始的壁面温度，对其他壁面设置无滑移绝热边界条件，采用MRF动参考系模型，通过流体仿真软件计算求解出流体域耦合面上的热流密度和近壁面温度。
步骤4具体过程为：考虑水下涡轮机部分进气结构和轮盘的旋转效应对轮盘温度分布的影响，通过流体仿真软件读取各个耦合面上的热流密度和近壁面温度，通过MATLAB实现各个耦合面上的热流密度和近壁面温度沿周向进行平均；
对涡轮盘前、后轮盘面和围带与流体域的耦合面采用周向平均的方法对热流密度和近壁面温度进行平均，对涡轮盘叶片与流体域的耦合面采用相对位置平均的方法对热流密度和近壁面温度进行平均；
根据平均后的耦合面上各个节点处的热流密度和近壁面温度，由以下公式计算耦合面上各个节点处的对流换热系数，
式中，h为对流换热系数，q为热流密度，Tnearwall为近壁面温度，第k次迭代中即本次迭代中耦合面的壁面温度，在这里为设置的初始壁面温度。
步骤5具体过程为：建立基于常物性、无内热源、稳态的固体域导热模型，即Laplace方程：通过材料手册获得固体域材料的密度、比热和热导率，将各个耦合面上的对流换热系数和近壁面温度作为固体域的边界条件，结合固体域导热模型对固体域进行数值模拟，得到固体域耦合面的壁面温度。
步骤6具体过程为：通过以下公式将固体域耦合面的壁面温度进行松弛：
式中为第k+1次迭代中流体域耦合面的壁面温度，/>为第k次迭代得到的固体域耦合面的壁面温度，β为松弛因子。
步骤6中判断壁面温度是否收敛的具体过程为：设置收敛阈值，将新的固体域耦合面的壁面温度与上一循环步骤5中固体域耦合面的壁面温度相减，若得到的差小于收敛阈值，则收敛；否则不收敛。
步骤7具体过程为：通过材料手册获得固体域材料的密度、杨氏模量、泊松比和平均线膨胀系数，将涡轮盘的轴面设置为固定面，将涡轮盘离心载荷加载到涡轮盘上对轮盘进行有限元计算得到涡轮盘的机械应力和机械应变；将喷水冷却涡轮机涡轮盘的温度分布即涡轮盘的温度载荷加载到涡轮盘上，对轮盘进行有限元计算得到轮盘的热应力和热应变；将离心载荷、温度载荷和气动载荷同时加载到涡轮盘上对轮盘进行有限元计算得到喷水冷却涡轮机的涡轮盘的总应力和总应变。
本发明有益效果是：
本发明部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，根据仿真结果数据能够得到喷水冷却涡轮机的性能，适用于水下部分进气轴向冲动式涡轮发动机，为喷水冷却涡轮机的性能提供了直观的数据支撑，能够较好的解决水下涡轮机部分进气结构对涡轮机气热弹耦合仿真造成的影响，并能判断在喷水冷却下涡轮盘的可靠性。
附图说明
图1是本发明部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法流程图；
图2是本发明中的水下涡轮机流体域模型图；
图3为本发明中的水下涡轮机固体域模型图。
图中，1-燃气入口，2-转子，3-喷管，4-叶顶间隙，5-轴向冷却水入口，6-径向冷却水入口，7-轴向间隙，8-叶片，9-轮盘面，10-围带，11-轴面，12-端面。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，如图1所示，具体按照以下步骤实施：
步骤1、建立水下涡轮机流体域和固体域的三维模型；具体过程为：根据实际工程中的水下涡轮机三维模型参数，利用三维建模软件进行实体几何造型的构造，得到固体域实体模型，并反解出流体域实体模型，形成水下涡轮机流体域和固体域的三维模型。
为了方便，三维模型的流体域分为喷管流道域、转子叶栅域、叶顶间隙域和轴向间隙域，流体域模型如图2所示，流体域模型由转子2、喷管3、叶顶间隙4和轴向间隙7这四部分组成，轴向间隙7为固体域中轮盘面9和涡轮机机壳之间的间隙，轴向间隙7是由于喷管3出口的燃气泄漏造成的，轴向间隙7中部开设轴向冷却水入口5，叶顶间隙4为固体域中围带10和涡轮机机壳之间的间隙，叶顶间隙7是由于进入转子2的燃气泄漏造成的，叶顶间隙4上部开设径向冷却水入口6，水下涡轮机部分进气的结构使得在一部分圆弧上装设喷管3，喷管3前端为燃气入口1，喷管出口正对着转子2。
水下涡轮机固体域模型如图3所示，包括轮盘面9，多个叶片8均匀安装在轮盘的边缘，多个叶片8外套接围带10，轮盘后的中部有一轮毂，轮毂周向表面为轴面11，轮毂需要固定的面为端面12，除轴面11和端面12以外，其余面均为固体域的耦合面，流体域中与固体域耦合面相对应的面为流体域的耦合面。
步骤2、利用多重分块结构化网格方法对三维模型进行网格划分，全流体域采用六面体结构化网格，近壁面处需要边界层，且对回转型O型网格更容易实现对回转型流域的加密，能更好的捕捉壁面附近处流动的细节，所以在喷管流域、轴向间隙流域需要采用O型网格进行加密，而转子叶栅处则采用H型block。为了准确获取喉部处的流动信息，在靠近喉部位置沿流动方向进行了加密，由于围带处泄露由喷管喷出的高速气体从而产生叶顶间隙损失，为准确获取这部分损失，围带处则采用C型block以准确对应齿型密封，以实现叶顶间隙处的准确加密。固体域与流体域之间也应做到网格相互对应以提高耦合计算的收敛性，应保证固体域与流体域网格节点误差在0.01mm以内。
步骤3、采用流场仿真计算软件，建立具有超音速流场、多相流、传质相变的流场数值模拟仿真计算方法，设置喷水冷却涡轮机流体域耦合面的温度初值，对流体域进行数值模拟，得到流体域耦合面上的热流密度和近壁面温度；具体过程为：对喷水冷却涡轮机的流体域进行仿真，其中流体域的入口为压力入口，总温总压为实际进气总温总压，出口为压力出口，静压为实际出口静压，轴向和径向喷水口设置为质量流量入口，喷入的冷却水的流量和温度由实际工程决定，采用3D定常压力基求解器，双精度，使用COUPLED算法，燃气和水蒸气为理想可压气体并打开组分输运模型考虑质量和热量的扩散，液态水为不可压液体，液体和气体之间采用VOF模型进行仿真，湍流模型采用k-ωSST模型，与固体域相耦合的壁面设置一个初始的壁面温度，对其他壁面设置无滑移绝热边界条件，采用MRF动参考系模型，通过流体仿真软件计算求解出流体域耦合面上的热流密度和近壁面温度。
步骤4、考虑水下涡轮机部分进气结构和轮盘的旋转效应对轮盘温度分布的影响，通过流体仿真软件读取各个耦合面上的热流密度和近壁面温度，通过MATLAB实现各个耦合面上的热流密度和近壁面温度沿周向进行平均；
由于步骤3仿真过程中转子旋转时数据传输应用混合平面或者周向平均边界条件，故前后轮盘面和围带与流体域的耦合面采用周向平均的方法对热流密度和近壁面温度进行平均，因每个转子叶片的网格节点的相对位置一样，故叶片与流体域的耦合面采用相对位置平均的方法对热流密度和近壁面温度进行平均。
根据平均后的耦合面上各个节点处的热流密度和近壁面温度，由以下公式计算耦合面上各个节点处的对流换热系数，
式中，h为对流换热系数，q为热流密度，Tnearwall为近壁面温度，第k次迭代中即本次迭代中耦合面的壁面温度，在这里为设置的初始壁面温度。
步骤5、建立固体域导热模型，将各个耦合面上的对流换热系数和近壁面温度作为固体域的边界条件对固体域进行数值模拟，得到固体域耦合面的壁面温度；具体过程为：建立基于常物性、无内热源、稳态的固体域导热模型，即Laplace方程：通过材料手册获得固体域材料的密度、比热和热导率，将各个耦合面上的对流换热系数和近壁面温度作为固体域的边界条件，结合固体域导热模型对固体域进行数值模拟，得到固体域耦合面的壁面温度。
步骤6、通过以下公式将固体域耦合面的壁面温度进行松弛：
式中为第k+1次迭代中流体域耦合面的壁面温度，/>为第k次迭代得到的固体域耦合面的壁面温度，β为松弛因子。
松弛后得到的温度作为流体域的耦合面壁面温度，再次对流体域进行数值模拟，执行步骤3～5，得到新的固体域耦合面的壁面温度；设置收敛阈值，将新的固体域耦合面的壁面温度与上一循环步骤5中固体域耦合面的壁面温度相减，若得到的差小于收敛阈值，则收敛；否则不收敛。
若收敛，则结束循环，若不收敛，则返回步骤3，直到壁面温度收敛，得到喷水冷却涡轮机涡轮盘的温度分布。
步骤7、由最后一次循环中流体域数值模拟的结果得到气动载荷，喷水冷却涡轮机涡轮盘的温度分布即为温度载荷，将温度载荷、气动载荷及离心载荷加载到涡轮盘上对轮盘进行有限元计算分析，得到喷水冷却涡轮机涡轮盘的总应力和总应变；具体过程为：通过材料手册获得固体域材料的密度、杨氏模量、泊松比和平均线膨胀系数，将涡轮盘的轴面设置为固定面，将涡轮盘离心载荷加载到涡轮盘上对轮盘进行有限元计算得到涡轮盘的机械应力和机械应变；将喷水冷却涡轮机涡轮盘的温度分布即涡轮盘的温度载荷加载到涡轮盘上，对轮盘进行有限元计算得到轮盘的热应力和热应变；将离心载荷、温度载荷和气动载荷同时加载到涡轮盘上对轮盘进行有限元计算得到喷水冷却涡轮机的涡轮盘的总应力和总应变。
步骤8、根据喷水冷却涡轮机涡轮盘的总应力和总应变建立变形后的涡轮机流体域和固体域模型，将变形后的涡轮机流体域和固体域模型作为三维模型，执行步骤2～7，得到受轮盘变形影响的喷水冷却涡轮机涡轮盘的总应力和总应变。
本发明中建立水下部分进气涡轮机的三维模型，目的是建立出与实际中的水下部分进气涡轮机相同的模型，为解决水下涡轮机部分进气结构对涡轮机气热弹耦合仿真造成的影响提供真实有效的模型，以准确计算出部分进气涡轮机气热弹耦合的结果；利用多重分块结构化网格方法对三维模型进行网格划分的目的是利用有限体积法求解控制方程，能够准确的求解出气热弹耦合中涡轮机流体域耦合面的热流密度和近壁面温度以及固体域的温度分布；建立流场数值仿真计算方法是为了更好的模拟出流体域的流场情况，能够准确的求解出气热弹耦合中涡轮机流体域的流场结果，为气热弹耦合提供精确的流体域耦合面的热流密度和近壁面温度；将各个耦合面上的热流密度和近壁面温度沿周向进行平均，是为了解决部分进气结构对涡轮机气热弹耦合仿真造成的影响，沿周向进行平均是考虑到了轮盘的旋转效应，使涡轮盘的温度分布更加均匀，使其不受到部分进气结构的影响。计算得到耦合面上的对流换热系数是为气热弹耦合中固体域的温度计算提供边界条件，使得计算出的固体域温度分布更加准确；建立固体域导热模型，设置固体域的边界条件是为了能够准确的求解出气热弹耦合中固体域的温度分布，得到固体域耦合面的壁面温度；将固体域耦合面的壁面温度进行松弛是为下一循环中流体域耦合面的壁面温度提供数值，执行步骤3～5直到固体域的温度收敛，至此完成了涡轮机气热弹耦合中的气热耦合的过程，得出了流体域的流场结果和固体域的温度分布。将离心载荷、由流体域的流场结果得到的气动载荷以及固体域的温度分布即温度载荷加载到涡轮盘进行有限元计算分析，是为了得到喷水冷却涡轮机涡轮盘的总应力和总应变，为涡轮机气热弹耦合提供涡轮盘的变形；根据涡轮盘的总应力和总应变建立受变形影响的涡轮机流体域和固体域的模型后，执行步骤2～7后，得到受轮盘变形影响的喷水冷却涡轮机涡轮盘的总应力和总应变，至此完成了整个涡轮机气热弹耦合的过程，根据此时的涡轮盘的总应力和总应变可判断喷水冷却下涡轮盘的可靠性。
通过上述方式，本发明部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，根据仿真结果数据能够得到喷水冷却涡轮机的性能，适用于水下部分进气轴向冲动式涡轮发动机，为喷水冷却涡轮机的性能提供了直观的数据支撑，能够较好的解决水下涡轮机部分进气结构对涡轮机气热弹耦合仿真造成的影响，并能判断在喷水冷却下涡轮盘的可靠性。

Claims (9)

1.部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，其特征在于，具体按照以下步骤实施：
步骤1、建立水下涡轮机流体域和固体域的三维模型；
步骤2、利用多重分块结构化网格方法对三维模型进行网格划分；
步骤3、采用流场仿真计算软件，建立具有超音速流场、多相流、传质相变的流场数值模拟仿真计算方法，设置喷水冷却涡轮机流体域耦合面的温度初值，对流体域进行数值模拟，得到流体域耦合面上的热流密度和近壁面温度；
步骤4、考虑轮盘的旋转效应，将各个耦合面上的热流密度和近壁面温度沿周向进行平均后计算得到耦合面上的对流换热系数；
步骤5、建立固体域导热模型，将各个耦合面上的对流换热系数和近壁面温度作为固体域的边界条件对固体域进行数值模拟，得到固体域耦合面的壁面温度；
步骤6、将固体域耦合面的壁面温度进行松弛后，作为流体域的耦合面壁面温度，再次对流体域进行数值模拟，执行步骤3～5，得到新的固体域耦合面的壁面温度，将新的固体域耦合面的壁面温度与上一循环步骤5中固体域耦合面的壁面温度相比，判断壁面温度是否收敛，若收敛，则结束循环，若不收敛，则返回步骤3，直到壁面温度收敛，得到喷水冷却涡轮机涡轮盘的温度分布；
步骤7、由最后一次循环中流体域数值模拟的结果得到气动载荷，将温度载荷、气动载荷及离心载荷加载到涡轮盘上对轮盘进行有限元计算分析，得到喷水冷却涡轮机涡轮盘的总应力和总应变；
步骤8、根据喷水冷却涡轮机涡轮盘的总应力和总应变建立变形后的涡轮机流体域和固体域模型(9)，将变形后的涡轮机流体域和固体域模型(9)作为三维模型，执行步骤2～7，得到受轮盘变形影响的喷水冷却涡轮机涡轮盘的总应力和总应变。
2.根据权利要求1所述的用于部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，其特征在于，步骤1具体过程为：根据实际工程中的水下涡轮机三维模型参数，利用三维建模软件进行实体几何造型的构造，得到固体域实体模型，并反解出流体域实体模型，形成水下涡轮机流体域和固体域的三维模型。
3.根据权利要求1所述的用于部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，其特征在于，所述三维模型的流体域分为喷管流道域、转子叶栅域、叶顶间隙域和轴向间隙域，在喷管流道域、轴向间隙流域采用O型网格进行加密，转子叶栅域处则采用H型block进行加密，喷管流域在靠近喉部位置沿流动方向进行加密。
4.根据权利要求1所述的用于部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，其特征在于，步骤3具体过程为：对喷水冷却涡轮机的流体域进行仿真，其中流体域的入口为压力入口，总温总压为实际进气总温总压，出口为压力出口，静压为实际出口静压，轴向和径向喷水口设置为质量流量入口，喷入的冷却水的流量和温度由实际工程决定，采用3D定常压力基求解器，双精度，使用COUPLED算法，燃气和水蒸气为理想可压气体并打开组分输运模型考虑质量和热量的扩散，液态水为不可压液体，液体和气体之间采用VOF模型进行仿真，湍流模型采用k-ωSST模型，与固体域相耦合的壁面设置一个初始的壁面温度，对其他壁面设置无滑移绝热边界条件，采用MRF动参考系模型，通过流体仿真软件计算求解出流体域耦合面上的热流密度和近壁面温度。
5.根据权利要求4所述的用于部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，其特征在于，步骤4具体过程为：考虑水下涡轮机部分进气结构和轮盘的旋转效应对轮盘温度分布的影响，通过流体仿真软件读取各个耦合面上的热流密度和近壁面温度，通过MATLAB实现各个耦合面上的热流密度和近壁面温度沿周向进行平均；
对涡轮盘前轮、后轮盘面和围带与流体域的耦合面采用周向平均的方法对热流密度和近壁面温度进行平均，对涡轮盘叶片与流体域的耦合面采用相对位置平均的方法对热流密度和近壁面温度进行平均；
根据平均后的耦合面上各个节点处的热流密度和近壁面温度，由以下公式计算耦合面上各个节点处的对流换热系数，
式中，h为对流换热系数，q为热流密度，Tnearwall为近壁面温度，第k次迭代中即本次迭代中耦合面的壁面温度，在这里为设置的初始壁面温度。
6.根据权利要求5所述的用于部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，其特征在于，步骤5具体过程为：建立基于常物性、无内热源、稳态的固体域导热模型，即Laplace方程：通过材料手册获得固体域材料的密度、比热和热导率，将各个耦合面上的对流换热系数和近壁面温度作为固体域的边界条件，结合固体域导热模型对固体域进行数值模拟，得到固体域耦合面的壁面温度。
7.根据权利要求6所述的用于部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，其特征在于，步骤6具体过程为：通过以下公式将固体域耦合面的壁面温度进行松弛：
式中为第k+1次迭代中流体域耦合面的壁面温度，/>为第k次迭代得到的固体域耦合面的壁面温度，β为松弛因子。
8.根据权利要求6所述的用于部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，其特征在于，步骤6中所述判断壁面温度是否收敛的具体过程为：设置收敛阈值，将新的固体域耦合面的壁面温度与上一循环步骤5中固体域耦合面的壁面温度相减，若得到的差小于收敛阈值，则收敛；否则不收敛。
9.根据权利要求1所述的用于部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法，其特征在于，步骤7具体过程为：通过材料手册获得固体域材料的密度、杨氏模量、泊松比和平均线膨胀系数，将涡轮盘的轴面设置为固定面，将涡轮盘离心载荷加载到涡轮盘上对轮盘进行有限元计算得到涡轮盘的机械应力和机械应变；将喷水冷却涡轮机涡轮盘的温度分布即涡轮盘的温度载荷加载到涡轮盘上，对轮盘进行有限元计算得到轮盘的热应力和热应变；将离心载荷、温度载荷和气动载荷同时加载到涡轮盘上对轮盘进行有限元计算得到喷水冷却涡轮机的涡轮盘的总应力和总应变。
CN202210413920.0A 2022-04-15 2022-04-15 部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法 Active CN114880764B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210413920.0A CN114880764B (zh) 2022-04-15 2022-04-15 部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210413920.0A CN114880764B (zh) 2022-04-15 2022-04-15 部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法

Publications (2)

Publication Number Publication Date
CN114880764A CN114880764A (zh) 2022-08-09
CN114880764B true CN114880764B (zh) 2024-03-19

Family

ID=82672263

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210413920.0A Active CN114880764B (zh) 2022-04-15 2022-04-15 部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法

Country Status (1)

Country Link
CN (1) CN114880764B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115494892B (zh) * 2022-08-31 2023-06-30 中国航发四川燃气涡轮研究院 高空模拟试车台进气环境模拟***解耦控制方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4338780A (en) * 1977-12-02 1982-07-13 Hitachi, Ltd. Method of cooling a gas turbine blade and apparatus therefor
CN108138584A (zh) * 2015-09-30 2018-06-08 西门子能源有限公司 燃烧器涡轮机壳体后集气室的螺旋冷却
CN108984920A (zh) * 2018-07-23 2018-12-11 安徽江淮纳威司达柴油发动机有限公司 一种发动机冷却水套的直接流固耦合传热分析方法
CN111460581A (zh) * 2020-04-01 2020-07-28 三一重机有限公司 一种发动机温度分布的计算方法、装置及存储介质
CN113217193A (zh) * 2021-05-14 2021-08-06 西北工业大学 一种涡轮机轮盘喷水冷却结构
CN113408180A (zh) * 2021-07-17 2021-09-17 北京工业大学 一种火箭发射场用涡街流量计流固耦合数值分析方法
CN113669160A (zh) * 2021-08-06 2021-11-19 西北工业大学 一种涡轮机燃气喷水冷却装置
WO2022011723A1 (zh) * 2020-07-17 2022-01-20 天华化工机械及自动化研究设计院有限公司 一种基于ansys的高温碳化炉多场耦合应力分布模拟方法
WO2022011724A1 (zh) * 2020-07-17 2022-01-20 天华化工机械及自动化研究设计院有限公司 一种基于workbench的低温碳化炉多耦合场应力分布的模拟方法
CN114117961A (zh) * 2021-11-30 2022-03-01 山东大学 一种基于计算流体动力学的涡轮机械优化方法及***

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4338780A (en) * 1977-12-02 1982-07-13 Hitachi, Ltd. Method of cooling a gas turbine blade and apparatus therefor
CN108138584A (zh) * 2015-09-30 2018-06-08 西门子能源有限公司 燃烧器涡轮机壳体后集气室的螺旋冷却
CN108984920A (zh) * 2018-07-23 2018-12-11 安徽江淮纳威司达柴油发动机有限公司 一种发动机冷却水套的直接流固耦合传热分析方法
CN111460581A (zh) * 2020-04-01 2020-07-28 三一重机有限公司 一种发动机温度分布的计算方法、装置及存储介质
WO2022011723A1 (zh) * 2020-07-17 2022-01-20 天华化工机械及自动化研究设计院有限公司 一种基于ansys的高温碳化炉多场耦合应力分布模拟方法
WO2022011724A1 (zh) * 2020-07-17 2022-01-20 天华化工机械及自动化研究设计院有限公司 一种基于workbench的低温碳化炉多耦合场应力分布的模拟方法
CN113217193A (zh) * 2021-05-14 2021-08-06 西北工业大学 一种涡轮机轮盘喷水冷却结构
CN113408180A (zh) * 2021-07-17 2021-09-17 北京工业大学 一种火箭发射场用涡街流量计流固耦合数值分析方法
CN113669160A (zh) * 2021-08-06 2021-11-19 西北工业大学 一种涡轮机燃气喷水冷却装置
CN114117961A (zh) * 2021-11-30 2022-03-01 山东大学 一种基于计算流体动力学的涡轮机械优化方法及***

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
多级涡轮气热耦合仿真性能及强度对比研究;张鹏飞;贾少锋;许开富;朱东华;任众;;火箭推进;20200215(01);全文 *
水下小型特种燃气轮机轮盘热-流耦合仿真分析;肖炎彬;史小锋;伊寅;封启玺;伊进宝;赵海涛;;水下无人***学报;20190415(02);全文 *

Also Published As

Publication number Publication date
CN114880764A (zh) 2022-08-09

Similar Documents

Publication Publication Date Title
Han Fundamental gas turbine heat transfer
Bunker Axial turbine blade tips: function, design, and durability
CN114880764B (zh) 部分进气结构的喷水冷却轴向涡轮机气热弹耦合计算方法
Zhang et al. Design of an air-cooled radial turbine: Part 1—Computational modelling
Zou et al. A new similarity method for turbomachinery with different working media
Xia et al. A novel design of cooling air supply system with dual row pre-swirl nozzles
Baines et al. Computer aided design of mixed flow turbines for turbochargers
Sadasivan et al. Numerical simulation of diffuser of a gas turbine using the actuator disc model
CN113609797B (zh) 一种基于cfd的动叶端壁复合射流下气膜冷却特性仿真方法
Jinsong et al. Study on heat transfer performance of shaft seal in HTR-PM circulator
Amano et al. Study of flow through a stationary ribbed channel for blade cooling
Lv et al. Study on flow characteristics of leakage in the last stage inner blade tip clearance of steam turbine
Čížek et al. On CFD investigation of radial clearance of labyrinth seals of a turbine engine
Rogalev et al. The design approach for supercritical CO2 gas turbine
CN112149327B (zh) 一种基于仿生的迷宫密封流-固-热多场分析及承载隔热性能表征方法
Gaitanis et al. Towards real time transient mGT performance assessment: Effective prediction using accurate component modelling techniques
Wang et al. Fluid-thermal analysis of partial admission axial turbines
Senoo et al. Development and verification tests of titanium 50-inch and 60-inch last stage blades for steam turbines
Bohn et al. Influence of open and closed shroud cavities on the flowfield in a 2-Stage turbine with shrouded blades
Niu et al. Gas-dynamics design of reversible turbine for marine gas turbine engine
Ma et al. Axial hole air-entraining cooling technique of radial-flow turbines
Choi et al. Local heat/mass transfer measurements in a shrouded rotor-stator cavity
Pei et al. Heat transfer and aerodynamics of complex shroud leakage flows in a low-pressure turbine
CN108875154A (zh) 一种气冷涡轮低热应力气热设计方法
Wang et al. Fluid–Thermal–Structural Analysis of Partial Admission Axial Impulse Turbines With Liquid Jet Impingement Cooling

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