氢气非预混燃烧的流场输运特性分析.docxVIP

PAGE 6 PAGE 8 中国工程热物理学会 燃烧学 学术会议论文 编号：184324 氢气非预混燃烧的流场输运特性分析 姜婕妤，张俊卿，许晨滨，谭睿 （青岛科技大学，山东省青岛市 266100） (Tel Email:jjieyu2006@163.com) 摘 要 本文采用三维直接数值模拟的方法对氢气非预混燃烧流场进行了精细的预测，并探讨了浮力对氢气燃烧流场输运特性的影响。研究发现，在浮力驱动的燃烧流场中对流是质量、动量及热量输运行为的主要影响因素，而无浮力火焰中扩散主导着流场的输运特性。平均统计结果表明有浮力和无浮力的燃烧流场中都可以捕捉到逆梯度输运现象，且浮力会促进逆梯度输运行为的发生。 关键词 浮力；输运特性；直接数值模拟；非预混燃烧；逆梯度输运 0 前言 氢气以其质量轻、热值高、无污染等优点可广泛应用于现代高科技产业中，如航天器、火箭、汽车发动机等方面，它作为一种替代性的清洁能源有着十分诱人的应用前景。由于氢气具有较强的扩散性，不仅会诱发燃烧不稳定性，还会影响湍流输运行为。氢气燃烧过程中产生的化学组分如H，OH，HO2和H2O2等具有不同的扩散速度，使燃烧过程中遵循优势扩散的行为[1-3]。优势扩散对燃料的混合和化学反应进度产生影响[4-6]。另外，统计结果表明，有热释放的燃烧流场中热扩散会存在逆梯度输运的现象。同样，在动量输运和其他湍流标量的输运中也可以捕捉到逆梯度输运行为。这和传统梯度输运（GD）模型[7]的描述是不一致的。在梯度输运模型中标量通量一般模化为： (1) 其中，ρ，ui，νT，分别表示密度，i方向的速度分量和湍流扩散系数，q表示描述燃烧流场的任意标量，如温度，混合分数，反应进度变量等。为favre平均量，为favre过滤脉动量。但很多研究发现，湍流标量通量与标量梯度可能同为正值或负值，并不符合上述方程所描述的，即在燃烧流场中存在着逆梯度输运现象。 Bray等认为湍流通量与未燃气体和已燃气体速度（uiu 和 uib）直接相关，火焰锋面处已燃气体速度 uib 极有可能大于未燃气速度 uiu，此时对应着逆梯度输运行为[8]。研究发现，流场中梯度输运区域随着湍流强度的增加而增加，随热释放因子 τ 的增加而减少，而逆梯度输运的发生规律则相反[9，10]。另外，通过对输运方程里各组成项的分析发现，逆梯度输运主要受压力梯度驱动的影响[11,12]。Luo和Bray[13]认为，当热释放率足够大时，非预混湍流燃烧中存在局部逆梯度输运现象。逆梯度输运的现象不仅会发生于热量的输运中，对于动量输运和其他标量输运都会发生[14,15]。前人的工作大都是基于预混火焰的流场输运行为，非预混火焰的输运特性和相关修正模型的研究则很少，且浮力多数被归因于压力梯度效应，但实际上浮力在燃烧流场的湍流输运和混合速率中起着重要作用[16]。因此有必要深入探讨氢气非预混火焰湍流输运特性 基金项目：国家自然科学基金项目 (No.; 山东省自然科学基金项目 (ZR2017BEE063)  1 数值方法 针对氢气非预混燃烧问题开展大规模的直接数值模拟研究，如图1所示，燃烧流场的计算区域为Lx = Ly = Lz = 12 D（D为燃料喷嘴直径）的三维空间。分别设置流场有浮力和无浮力两种情况来计算该物理问题，并进行对比分析。模拟中求解的控制方程为三维可压缩的无量纲方程，包括连续性方程、动量守恒方程、能量守恒方程、输运方程以及状态方程,其具体形式如下： (2) (3) (4) (5) (6) (7) 方程组求解时采用的空间导数离散格式为六阶紧致差分格式，它能在一定程度上减少节点间的通讯次数，在保证高精度的同时能够节省大量的计算资源。时间步进格式为三阶显式低存储Runge-Kutta格式，程序中直接嵌入FGM化学反应机制[17]，这种方法是建立于层流火焰面的思想上的，将实际中三维的火焰结构本质上能够看作是由一维的火焰结构组成的，并用来构造相应的子空间。入口处采用特征无反射边界条件 (NSCBC)，出口处和侧边界采用无反射边界条件 (CBC)。为实现程序的高效运行，采用768个处理器核开展并行计算，相应