滞止点催化燃烧的一维和二维数值分析研究.pdfVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：094154 滞止点催化燃烧的一维和二维数值分析* 徐侃，刘明侯，邢丹，王亚青，刘东，陆游 ( 中国科学技术大学，合肥，230027) (TEL: 0551-3603127, Email: mhliu@ustc.edu.cn) 摘要：采用FLUENT耦合CHEMKIN表面反应方法实现二维轴对称滞止点甲烷催化燃烧的数值模拟，并与基于 Spin程序的一维计算结果进行比较。研究结果表明，二维模拟时由于考虑径向扩散，相同拉伸率时，二维 的火焰前沿更靠近滞止点，且最大组分浓度、最大温度等值均小于一维结果；二者轴向速度趋势一致，变 化幅度相当，但由于火焰前沿位置变化，导致流场整体平移；各个拉伸率时催化滞止面上的径向温度梯度 都很小，但空间反应区低拉伸率时径向温度梯度很大。高拉伸率时，一维简化模拟可以得到较好的轴向温 度和组分浓度分布。 关键词：滞止点，催化燃烧 1 引言 催化燃烧是实现贫燃料燃烧最有效的手段之一。采用催化燃烧技术，甲烷燃烧可在更 低的燃料/空气(1-5 ％)范围内和更低温度下进行，在保证燃料转化率的情况下，可降低燃烧 温度，减少CO、NOx 的产生，起到保护环境和提高资源利用率的作用。一直以来，催化燃 烧的研究主要是采用实验手段优化燃烧器设计、提高催化转化效率和寻找合适的催化剂等 方面。Wood[1] 曾指出进一步实验研究和模型建立将是多孔介质催化燃烧的重点。 美国 Sandia 国家实验室建立了催化反应动力学计算理论，开发出已被广泛应用的气相 和催化反应动力学模拟软件CHEMKIN[2] 。运用DETCHEM 软件，Redenius[3]等建立了基于 化学反应动力学的表面催化反应数值模拟方法，计算结果与实验结果符合很好。钟北京等 人[4]联合使用FLUENT 和DETCHEM ，分析微通道中CH4/Air 预混气体催化燃烧时温度、 当量比对甲烷转化率，转化速度的影响。 滞止点反应器常被用来研究催化燃烧。由于滞止点反应器的中心轴线上的组分和温度 可以简化为一维模型，减小了建模和计算的难度，使数值方法更多集中于解决化学反应的 刚性问题。然而，一维假设有一定的限制条件。为更好指导实际情况下的实验设计，必须 开发更合理的数值模型，或者提出一维假定的适用条件。以该反应器为平台，Li[5]数值分析 了滞止点催化燃烧的熄火特性。Yuan[6]等把一维和二维计算结果进行比较，分析二维长宽比 对氢气催化燃烧的影响。本研究在FLUENT 中耦合CHEMKIN 表面反应实现二维轴对称滞 止点甲烷催化燃烧的数值模拟，并与一维程序Spin[7] 的结果进行比较，分析不同拉伸率(SR) 时滞止点附近的火焰稳定特性和中心轴线上组分浓度和温度的一维特性。 2 数值模型和算法 2.1 物理模型和计算网格 ―――――――――――― *资助项目：国家863 专题资助项目(No. 2007AA05Z236) 图1 滞止点流动示意图 滞止点甲烷催化燃烧几何模型如图 1 所示。甲烷/空气的预混气从直径为D 的喷嘴射出 撞击其下方 L 处的圆盘，并在圆盘中心(0,0)处形成滞止点。圆盘表面带有催化剂涂层，用 于催化燃烧；从混合气入口到圆盘的整个反应区域称为空间反应区域；圆盘表面附近为表 面催化反应区域。预混气体入口条件不同，火焰会稳定在不同的地方。理论上，该流场结 构具有轴对称特性。在射流中心轴线上，虽然流动有二维特性，但组分浓度和温度等都只 与喷嘴到圆盘的距离 x 有关，因此，无限大圆盘假设条件下该问题具有一维特性。但是， 实际圆盘直径并非无穷大，因此，需要考虑圆盘半径W 和离燃料入口的距离L 。二维轴对 称模型能很好的模拟有限 情况下的反应流过程。图 中给出了二维结构化网格，包含圆 W 1 盘壁面附近的催化反应区，沿圆盘径向有60 个结点，轴向110 个结点。只考虑催化反应计 算时，基于 Spin 一维程序自适应网格加密后轴向网格结点均小于 100，且二维计算中加密