正庚烷氧化机理的简化与加速计算.pdfVIP

中国内燃机学会燃烧节能净化分会2010 年学术年会 CSICE2010-069 正庚烷氧化机理的简化与加速计算 陈 正 北京大学 工学院 湍流与复杂系统国家重点实验室 北京 100871 摘要：详细化学反应机理的引入会给化学反应流的数值模拟带来巨大的困难：一方面，由于不同组分对应 不同的特征时间，详细化学反应机理会导致化学反应流涉及到广泛的时间尺度；另一方面，随着燃料所含 碳原子数目的增加，其详细化学反应机理中所含的组分数目与基元化学反应数目会呈指数增长，这直接导 致计算量的急剧增加。为了解决这两方面的困难，本文以正庚烷氧化机理为例，通过化学反应机理简化(PFA 方法)与加速算法(PM 方法)实现了在确保计算精度的条件下极大程度地提高计算效率。 关键词：正庚烷；简化机理；PFA 方法；PM 方法 基于详细化学反应动力学的燃烧数值 (561 个组分和 2539 个基元化学反应)[5] 。通 模拟有助于深入认识基础燃烧现象、揭示其 过模拟不同条件下的均质着火过程以及均 中的物理化学机理，是开发新型替代燃料和 匀搅拌器中的着火与熄火过程，PFA 首先确 高性能燃烧器的重要研究手段[1]。然而，详 定各个组分之间的关联系数[6]。接下来，根 细化学反应机理的引入会给化学反应流的 据给的关联阈值，删除次要组分及其相关反 数值模拟带来巨大的困难：一方面，由于不 应，获得了包含 66 个组分和 389 个基元反 同组分对应不同的特征时间，详细化学反应 应的简化机理。该简化机理的适用范围为： 机理会导致化学反应流涉及到广泛的时间 当量比位于0.5 到 2.0 之间，温度位于 1000K 尺度[2]；另一方面，随着燃料所含碳原子数 到 2000K 之间，位力间于 1atm 到 30atm 之 目的增加，其详细化学反应机理中所含的组 间。 分数目与基元化学反应数目会呈指数增长， 为了验证简化机理的准确性，本文研究 这直接导致计算量的急剧增加[3]。在关于甲 了不同条件下的正庚烷/ 空气均质着火过 烷/ 空气预混火焰传播的数值模拟中，耗费 程。对于一个均质系统而言，温度 T 以及各 在求解详细化学反应机理(53个组分，325个 组分的质量分数 Yk (k = 1, 2,…n, n 为组分的 基元反应)上的CPU 时间占了总体的95%左 种类数) 的控制方程为： 右[4]；对于大分子碳氢燃料，更多的计算时 dY / dt ω W / ρ (1) k k k 间将会耗费在求解详细化学反应机理上。因 n dT / dt −Σ (h ω W ) /(ρC ) (2) k 1 k k k P 此，如何有效地处理好广泛时间尺度与复杂 其中t, ρ, C , ω, W 分别为时间，密度，定压 P k k 化学反应机理是实现对化学反应流高效精 热容，组分k 的净生成率，组分k 的分子量。 确数值模拟的关键所在。