基于网格化Damkohler数的燃烧室贫熄边界预测方法探究.docxVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：2018.05.15.0042 基于网格化Damk?hler数的燃烧室贫熄边界预测方法探究 王中豪1,2，张军华1,2，邓爱明1，胡斌1,2基金项目：国家重点研发计划（2016YFB0901402）；国家自然科学基金 基金项目：国家重点研发计划（2016YFB0901402）；国家自然科学基金*通讯作者：胡斌（1982-），男，博士，副研究员，研究领域为航空发动机中的燃烧。Email: hubin@iet.cn 1中国科学院工程热物理研究所，北京 100190 2中国科学院大学，北京 100049 3中国科学院轻型动力重点实验室，北京 100190 摘 要 本文发展了一种基于网格化Damk?hler (Da) 数的贫熄预测方法。将基于OH基浓度定义的反应区均匀划分为尺度一致的子区域，以降低反应区组分、温度不均匀对预测精度的影响。针对整个反应区及反应区内的每个子区域分别计算Da数。分别采用新鲜混气在反应区的实际停留时间和能够触发反应的最短停留时间定义流动时间和化学反应时间。局部Da数计算结果表明，不同流量下的近贫熄工况，燃烧核心集中在回流区上游滞止点附近，设计工况下的燃烧核心却集中在回流区下游滞止点附近。整体Da数计算结果表明，对于空气来流流量从0.6kg/s增加到1.2kg/s的贫熄工况，整体Da数维持在7.9上下波动，平均波动范围在5%左右，初步证明了网格化Damk?hler方法应用于贫熄边界预测的可能性。 关键词 旋流稳燃；贫熄预测；Damk?hler数；反应区 0 前言 因贫油熄火导致的飞行安全问题一直被研究者所重视，在航发燃烧室设计阶段，需要贫熄性能评估和燃烧室结构改进的多次循环，而发展一套能够高效高精度预测航发燃烧室贫熄边界的方法，对于缩短燃烧室设计周期，降低设计成本具有重要意义[1]。 从早期Longwell等人提出均匀搅拌反应器模型[2]，Zukoski等人提出回流区点火模型开始[3]，基于能量平衡和时间平衡的半经验模型就在大量贫熄实验的基础上不断被发展完善[4-9]。随着计算机技术的发展，大涡模拟方法被广泛应用于分析熄火前的火焰行为，部分研究者将大涡模拟方法直接应用于燃烧室的贫熄边界预测。通过大涡模拟方法对燃烧室贫熄边界进行预测，一般是通过再现贫油熄火的瞬态过程实现的，而基于不同的燃烧类型及运行工况恰当地选择燃烧模型，对于能否直接通过数值模拟准确再现贫熄过程有决定性作用。 Kim等人应用涡破碎模型针对钝体稳燃预混火焰进行大涡模拟[10]，在当量比0.52时模拟出熄火，但是在油气比为化学恰当比时模拟的钝体后回流区长度却与实验有一定偏差，Kim等人将此偏差归因于该工况下发生的强烈的非平衡化学反应。Gokulakrishnan等人采用涡耗散概念模型和层流燃烧模型对钝体稳燃预混火焰进行了大涡模拟[11]，在当量比为0.45时，采用涡耗散概念模型能够模拟出全局熄火，与实验观察到的结果基本一致，而采用层流燃烧模型的模拟中燃烧依然维持稳定，与实验结果不符。Gokulakrishnan等人认为这是由于涡耗散概念模型假定燃烧发生在微尺度涡团中，而层流燃烧模型在亚格子尺度下假定标量不发生脉动，前者更好地模拟了湍流和化学反应的相互作用。Black等人采用概率密度输运模型对低排放旋流雾化燃烧进行了大涡模拟[12]。在实验获得的贫熄当量比下，数值模拟结果中燃烧依然稳定，未能获得预期的贫熄结果。Black等人认为是流场中局部富油而产生的高温区影响了预测结果的准确性，强调了精确模拟燃油轨迹及蒸发对于雾化燃烧贫熄预测的重要性。此外还有Garmory，[13]Zhang[14]等人采用条件矩方法对Sandia DF火焰进行的大涡模拟。模拟能够再现近贫熄工况频繁发生的局部熄火和再点火过程，Zhang等人通过数值模拟获得了不同来流流速下的贫熄边界[14]，与实验结果误差在25%到60%之间。Menon等人采用线性涡模型对旋流预混火焰进行了大涡模拟[15]。模拟结果展示出了火焰拉伸对火焰传播的影响，但是由于一些尺度的涡没有准确求解，未能捕捉到准确的贫熄边界。 在当前计算机发展水平下，尽管能够对燃烧室中的流场进行大涡模拟，其计算成本依然很高,并不十分适用于工程计算。在Black等人的工作中[12]，尽管只采用了一步总包反应，每一个工况仍然需要在16900MHz的计算机上计算20天。考虑到大涡模拟仍然需要较高昂的计算成本，将贫油熄火半经验模型