涡喷加力燃烧室湍流燃烧数值模拟研究.pdfVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：064114 +涡喷加力燃烧室湍流燃烧数值模拟 叶桃红，洪家荣，陈义良 (中国科学技术大学，热科学与能源工程系，安徽，合1t1,．230026) edu．cn；联系电话：0551-3603048 thyef岔ustc 摘要：本文采用无结构网格技术，对具有环形火焰稳定器的三维涡喷加力燃烧室的湍流燃烧过程 进行了数值模拟，运用离散相方法和扩展的EBU化学反应模型得到了合理的温度j手IJ浓度场分布。在此 基础}：，本文比较了一步反应和两步反应机理下的温度场和浓度场，研究．厂火焰稳定器形状对总压恢 复系数和燃烧效率的影响，计算结果与实验值符合较好，对加力燃烧室的设i,l-+I+OI宄提供有利帮助． 关键词：数值模拟；湍流燃烧；加力燃烧室 1 引言 航空发动机加力燃烧室的数值模拟是发动机设计过程中的重要组成部分。由于-一个 完整的加力燃烧室包括火焰稳定器、混合器等多种组件，流场的几何结构及其复杂，划 分网格很困难，所以对整个i维加力燃烧室流场计算很难开展。目前为止，很多研究工 作都足针对仅有一豳或两圈环形稳定器的■维加力燃烧室进行的…叫引。11l‘12|#红儿引等采 用分区曲线坐标模拟稳定器及扩压器几何形状，在二维情况下取得了较好的结果，但曲 线坐标很难应用到具有复杂几何结构的三维加力燃烧室．此外，也有··些针对三维几何 加力燃烧室的研究工作。Ravichandran喁1对具有两圈环形稳定器和…排径向稳定器的加 力燃烧室内的三维冷态流场进行了模拟，同时还对有两个环形稳定器的模型加力燃烧室 内的燃烧场做了研究。章诚¨0J等对内外涵道气流之间布置了波瓣肜混合器，并有两圈环 形稳定器和一周径向稳定器的涡扇加力燃烧室冷态流场进行了研究。义I状9-10的数值 模拟时由于对稳定器、混合器和加力燃烧室扩压器壁面采用台阶近似，叫格划分较粗糙， 对i：t·算结果小町避免带来误差，对流场分析也不够完善。无结构网格技术能够对非常复 杂的几何体进行网格划分，生成计算所需要的网格，而不需要对几何体作近似处理。该 方法非常适合于具有复杂几何结构加力燃烧室的数值模拟，对加力燃烧室的设计和研究 提供技术支持。Unaune…’运用非结构网格和假定的PDF湍流扩散燃烧模型对喝扇加力燃 烧室内的燃烧过程的数值模拟模拟进行了初步探索，但是计算值缺乏与实验数据的比 较，且温度明显偏高． 本文将采用无结构网格技术，对具有三圈环形稳定器的涡喷加力燃烧窒湍流两榭燃 烧进行数值研究。针对涡喷加力燃烧室内部气流速度大、喷油系统复杂、液雾体积率很 小的特点。本文采用欧拉一拉格朗日方法处理两相流动和扩展的EBU湍流燃烧模型，首 先得到了一步反应机理下的加力燃烧室流场。然后进一步研究了／卜同化学反应机理和火 焰稳定器结构对计算结果的影响。 2采用的基本物理模型 2．1两相流动模型 对连续相采用欧拉方法进行处理，湍流模型采用标准的双方程k一占模；．gJ，在欧拉 航空推进技术骏征汁划(AFTD)J二课题(APTD-0502-004) 463 牮标系中气相场诸守恒方程的通用形式为n引 c·， 昙(p“矽)+昙(户渺)+昙(户w≯)=昙(r，罢)+导(r，考]+昙(r一，警)+薯+—_ 其中：矽=l为连续性方程，矽=u。v，w，k，占。h，m。，分别对应子乏个速度分量、 湍流动能、湍流动能耗散率，焓及组分的质量分数：％为颗粒相和连续相的质量、动 量和能量的相瓦作用源项，通过随机的轨道模型计算。 下进行，赋了；每个离散相颗粒一定的初始尺寸、温度、位置和速度，通过对Newton第 二定律建立的微分方程进行积分得到颗粒的瞬时速度。在对速度的微分表达式进行积分 就得到离散相颗粒的轨迹。颗粒的湍流扩散应用随机颗粒轨道模型进行模拟。 2．2湍流燃烧模型 航空煤油用cx以表示，一步和两步反应分别为： ． 对于一步反应： cxq+(x+％)q_x