喷管气体流动的数值模拟和实验研究.pdf

中国工程热物理学会 传热传质学 第叫一届年会 编号：053284 喷管气体流动的数值模拟和实验研究 余敏胡卓焕 。 上海理工大学热工程研究所上海200093 (Tel：021 Email：yuminm@citiz．net) 工况下的流场，并进行了测量实验。实验所得数值根好地吻合了数值模拟的结果。表明喷管通道内气 体流场特性的数值模拟在实际工程也用中是可行的，为喷管的优化设计提供丁有效的技术途径。 关键词：Laval喷管，Fluent，数值模拟。 0引言 喷管这一类部件在动力、航天等领域有着，。泛应用，是许多动力装置(例如燃气轮 机、蒸汽轮机)的重要部件，尤其是燃气轮机的进气道和尾喷管，其性能的优劣对燃气 轮机楚体性能有重大影响。长期以来，人们对其研究的手段常局限于实验研究，由于实 验条件等诸多方面的原因，如燃气轮机，由于工作环境特殊，实验研究很难甚至根本不 可行，而且存在实验代价过大等问题，不能对喷管整体性能及内部流场做很好的研究。 随着数值模拟最近十年的发展，其理论和研究手段逐渐完善，并很好应用于动力、航天 等领域，本文运用数值模拟的理论和工具，研究了喷管在不同工况下的流场，并进行了 实验验证，对喷管这一类部件的设计与开发和优化起到了很好的指导作用。 1数值模拟 1．1控制方程 本文采用任意坐标系卜轴对称非稳态可压缩流体基本控制方程的统一形式【l】，【2】： 专刍p争cy嚣一∥芸卅s+鄙 ㈩ 当。分别为i、u、v、T、k、￡时，则上述方程分别表示为连续方程、两个方向的 动量方程、能量方程、湍流动能及湍流耗散率方程。其中：S(E，q)是(t，q)坐标系下中 的源项；Sp为颓粒群对气相造成的源项。 1．2具体数学计算模型和边界条件的建立 考虑到喷管是一个轴对称系统，可应用对称性得到整个喷管的流动情况。把一个复 杂的三维流场计算问题简化为二维流场计算问题，方便了研究丁作，并减少了计算量。 对上述数值模拟，本文采用CFD软件Fluent进行计算”1。它是用于计算流体流动和传 热问题的程序。对相对复杂的儿何结构网络生存非常有效。并且还可根据计算结果调整 网络，对丁精确求解有较人梯度的流场有很实际的作用。 本文选取较为复杂的一种喷管型式一 Laval喷管作为具体的研究对象，建立的轴截 面数学计算模型如图l所示。 首先，对一给定计算工况下的Laval喷管 网1轴裁而计箩樽犁 进行计算网格的划分。根据边界层原理，考 虑到喷管靠近壁面部分气体流动变化相对比较剧烈，而靠近轴心部分气体流动变化相对 比较平缓的这一情况。为了尽量保证数值模拟的精度，在划分计算模型网格的时候，靠 近壁面部分的网格划分相对比较密，靠近轴心部分的网格划分相对比较稀。 对Laval喷管等比例的轴截面计算模型定义边界条件：左侧为压力一入口侧，右侧为 压力一出口侧，上侧为壁可，下侧为轴心线。这样，就建立了完整的与实验的Laval喷管 等比例的轴截面数学计算模型。 1．3数学计算横型的仿真计算及计算结果 运用Fluent对数学计算模型进行仿真计算。考虑到研究的流体是接近常温的大气， 因此，可以将其作为无粘性的理想气体处理。进口压力设定为大气压，出口压力分别设 行初始化，为了得到更高精度的计算结果，将残差收敛条件设定为1E一6，迭代次数设定 为1000次。在经过1000次的迭代运算后，计算模型的残差都小于1E一5，以保证数值模拟 的精度。经仿真计算一共得到上述不同背压下的流场分布情况，如“r所示。 压力场、速度场、温度场和密度场。图6为背压逐渐降低时，Laval喷管速度场变化 情况，从图6可以清楚的看到：Laval喷管当背压等于计算压力时，它就是计算J．况下 的等熵流动(不计摩阻影响)，当背压火于临界压力时，其喉部(临界截面)首先达到 声速如图6背压为0．09MPa时所示，