凹腔支板火焰稳定器PDF-LES模拟.docxVIP

基金项目：国家重点研发计划(2017YFB0202400，2017YFB0202402)，国家自然科学基金 中国工程热物理学会 燃烧学 学术会议论文 编号：184554 凹腔支板火焰稳定器PDF-LES模拟 王方1.2，窦力1.2，魏观溢1.2，刘邓欢1.2，金捷1.2 (1. 北航能动学院航空发动机数值仿真研究中心，北京 100191 2. 先进航空发动机协同创新中心，北京 100191) (Tel: Email: fwang@buaa.edu.cn) 摘 要 大涡模拟(LES)兼具高精度与合理计算量两个特点，概率密度函数模型(PDF)适用于湍流与复杂化学反应相互耦合的问题。本文在基于PDF-LES的Aero Engine Combustor Simulation Code (AECSC)程序基础上，对凹腔支板火焰稳定器进行数值模拟。使用气相版本对有无凹腔支板结构分别进行三个速度入口条件下的甲烷湍流燃烧模拟，并用两相版本对带凹腔支板结构进行设计工况下煤油喷雾的模拟。结果表明：凹腔结构的火焰稳定性要优于无凹腔结构；凹腔支板结构对于液相燃料的控制能力较气相更强。 关键词 概率密度函数；大涡模拟；凹腔；火焰稳定器；AECSC 0 前言 加力燃烧室的入口条件较主燃烧室更为恶劣，在高温高速低压且氧气含量较低的条件下，维持火焰稳定性有更大的难度。为了得到更好的火焰稳定性能以及减轻发动机总体质量，凹腔结构被采用于火焰稳定器上。将凹腔稳焰技术应用于一体化加力燃烧室，是在原有的一种内突扩加力燃烧室方案的基础上逐渐发展而来的。在此基础上，众多研究者对其进行了改进并提出了一种新的一体化火焰稳定器方案[1]。该方案将涡轮后框架支板改为凹腔火焰稳定器，不仅简化了结构、减少了零件数目，还有效降低总压损失。 随着计算机技术的发展，大涡模拟（LES）方法逐渐成为工程计算的优先选择。而概率密度函数（PDF）输运方程模型在熄火、复燃等有限速率化学反应起决定作用的现象时具有明显的优势。但其在湍流流场的计算上存在困难，考虑到LES在模拟湍流流场方面的诸多优点，将PDF和LES相结合正在逐渐成为新的研究趋势。Popov使用PDF联合LES方法研究Sandia系列钝体火焰，指出PDF-LES方法的结果比原先的计算结果有了较大的改善，尤其是在流场下游区域[2]。这些研究结果表明了将PDF和LES相结合对更好的研究湍流燃烧问题有很大的帮助。 目前PDF-LES方法尚未广泛用于工程计算。本文使用基于PDF-LES方法的AECSC程序对凹腔支板火焰稳定器进行研究，旨在研究该程序对于凹腔支板火焰稳定器的湍流燃烧精确性，同时展示凹腔稳焰特性方面的细节，以对火焰稳定器的改进工作提供参考。 1 程序与模型 1.1 程序简介 Aero Engine Combustor Simulation Code(AECSC)是基于PDF-LES方法的三维亚声速湍流燃烧数值模拟程序，最初版本由帝国理工William Philip Jones教授编写用于计算射流的大涡模拟，取得了很多优异的成果[3,4]。目前北京航空天学数值仿真中心使用并发展了气相和两相版本。 AECSC程序采用模块化编程，拓展性好，计算精度高。程序采用大涡模拟方法求解连续方程和动量方程，采用概率密度函数方法求解能量方程和组分方程。大涡模拟亚格子应力模型使用动力Smagorinsky-Lilly模型，概率密度函数方程使用欧拉随机场解法。蒸发模型使用Abramzon-Sirignano模型，雾化破碎模型使用Claudios破碎模型，气液两相耦合采用欧拉-拉格朗日方法，点火模型使用能量沉积模型。 1.2 数值模型 1.2.1 大涡模拟方法 大涡模拟需要对大尺度湍流运动直接进行计算，同时也需要对小于计算网格的小尺度应力建立模型求解。大涡模拟采用空间滤波函数对N-S方程进行变换。假设?(??, ??)是一个通用的瞬态流动变量，在位置?? = (??1, ??2, ??3)相应的过滤变量定义如下： ? (1) 程序中滤波方式采用箱式滤波，由于程序使用的过滤函数为均匀过滤过程，过滤运算和求导可以交换。在密度的变化的情况下，为了获得需要的变量，需要引入密度加权滤波函数 Favre平均，定义如下： ? (2) 经过Favre平均后的模型方程如下： ? (3) ? (4) 在上述的动量方程中，ρuiuj为未知量，引入τij=ρuiuj-