考虑冷气组分扩散的气冷涡轮流场的数值模拟第一部分数值模拟方法与计算网格研究.pdfVIP

中国T程热物理学会 热机气动热力学 学术会议论文 编呼：062109 考虑冷气组分扩散的气冷涡轮流场的数 值模拟1 第一部分数值模拟方法与计算网格 王强2王松涛颜培刚冯国泰王仲奇 ·．．· (哈尔滨工业大学458信箱，推进理论与技术研究所) 电话：0451 Email：ehuangwang．w@hit．edu．cn 摘要：推导了自然曲线坐标系下时间平均的多组分扩散方程。并在原程序基础上，增加了对冷气组分 扩散方程的求解，以及多组分对比热和气体常数的影响。然后用改进的程序对某型两级带冷气喷射的 涡轮叶栅进行了全三维N．S方程数值求解。 。 关键词： 涡轮，冷气掺混，组分扩散，变比热 1．前言 随着现代航空工业的迅猛发展，涡轮前燃气温度的不断提高，为了保护涡轮部件， 就必须对其进行有效的冷却川。冷气喷射使得研究人员一方面要考虑冷气在叶片内部冷 气通道的流动与换热12,31，另一方面也必须考虑冷气对流道内流场的影响，本文的研究 研究重点在于第二方面。 目前，考虑冷气掺混的计算方法主要有以下四种： (1)利用冷气掺混分布模型， 在掺混区考虑冷气的影响。这一方法的优点是不必采用过分密集的网格，它适合二维计 算，三维计算冷气分布模型很难确定： (2)冷气孔区与主流区采用不同疏密的网格分 区计算，在交界面处采用非对应的网格计算。这种方式由于主流区可以采用较少的网格， 这种方法是本文计算采用的方法，这种方法可适用丁=单列及多级计算，计算网格比较多， 多级计算要用并行给与并行算法：(3)“面源”法，这种方法的特点是根据守恒定律， 在冷气喷射出口施加一定的以射流方式表示的边界条件。这种做法比方法(2)更适合 多级计算。随着计算工具和计算手段的发展，求解二维NS方程早已成为现实，在全三 维的NS方程的求解程序中，将射流作为流场的一个入口边界的做法被广泛使用14嗡j。这 种方法的缺点是不能细致描述冷气对主流流场结构的影响： (4)上述三种方法均为气 冷涡轮不同的简化方法，它们都没有考虑冷气的组分计算、冷气的组分扩散计算，第四 种方法便是求解多组分方程的计算19,10]。 在前三种冷气掺混的计算方法中，比热通常都是作为常数计算，冷气和主流都是作 为同一种组分进行计算。准确地预测发动机叶片温度，依赖丁．精确的预测发动机真实条 件下的对流换热率…J。如果计算热传导系数的不准确性是±35％，则在典型发动机条件 下将会导致98℃的温度误差ll21。涡轮的工作范围通常很大，尤其在有冷气掺混的情况 下，流场内温度变化更大，所以在计算中考虑变比热是十分必要的。同时，在冷气量比 较大的情况下，比热又受流体的组分影响较大。因此在冷气掺混的计算中，为了得到较 为准确的温度分布，应使用多组分进行计算。但由于考虑了多组分的控制方程其复杂程 度和计算量都大为增加，所以目前的文献中在涡轮叶栅冷气掺混计算中使用多组分方法 进行计算的文章比较少；多组分计算一般运用于化：．亡和环境科学领域[13,141，在空调领域 也有少量文章发表【15】；近年来，在研究不同组分流体混合流动过程时，也有用到多组分 计算的Il61：在CFD软件中，多组分计算也犬多是关于多相流动和化学反应的。文献【9】 。 ’自然科学基金资助项目号：5047028 339 和【10】运川Fluent计算程序，采JHj多组分的守恒方法求解油气比，从而进行变比热与变 气体常数的冷气掺混数值计算。 我们已开发出一套可以很好地计算冷气掺混的程序HIT．3D。该程序可以对冷气掺 混射流流场进行精细模拟，同时稃序还考虑了由丁．濉度引起的变比热，但足对冷气组分 的影响还没有计及。作者在原有程序基础上，增加了对冷气浓度扩散方程的求解，同时 还在原变比热计算基础上，增加了组分浓度对气动热力参数的