HDMR方法在数值求解详细化学反应机理中的应用研究.pdfVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：094320 HDMR 方法在数值求解详细化学反应机 理中的应用 * 赵真龙，陈 正，陈十一 北京大学 工学院 湍流与复杂系统国家重点实验室 北京 100871 * Tel: Email: cz@pku.edu.cn 摘要：基于详细化学反应机理的数值模拟是当今燃烧领域的重要研究方向。由于燃烧问题涉及众 多的化学反应和广泛的特征时间范围，如何实现燃烧问题的高效数值模拟一直是研究的热点和难 点。本文将HDMR (High Dimensional Model Representation)方法应用于涉及详细化学反应机理的 燃烧数值模拟中，模拟了定压条件下零维均质氢气/空气的着火过程，验证了 HDMR 方法处理详 细化学反应机理的精确性和高效性。研究结果表明通过HDMR 方法，计算速度能够加快一个量级， 而计算精度几乎不受影响。 关键词：HDMR ；详细化学反应；数值模拟；氢气/空气着火 0 前言 在全球能源与环境问题日益严峻的形势下，新型替代燃料和高性能燃烧器的开发和 利用受到了世界各国的高度关注。作为新型替代燃料和高性能燃烧器的重要研究手段， 基于详细化学反应机理的数值模拟是当今燃烧领域的重要研究方向。由于燃烧问题涉及 众多的化学反应[1] [2] 和广泛的特征时间范围 ，如何实现燃烧问题的高效数值模拟一直是 研究的难点。 在涉及详细化学反应机理的燃烧数值模拟中，绝大部分 CPU 时间被消耗在求解化 学反应上。例如，在对一维甲烷/空气球形传播火焰的数值模拟中[3-5] ，求解详细化学反 应机理占 85 ％左右的CPU 时间，而流动只占 15％左右的CPU 时间。因此，高效精确 地求解详细化学反应机理是实现燃烧问题快速数值模拟的关键所在。针对详细化学反应 机理，目前最常用的方法是进行机理简化，将基元反应的数目和组分的种类在给定容忍 误差范围内减到最小，从而降低计算量[1] 。 与传统的反应机理简化方法不同，本文通过 HDMR(High Dimensional Model Representation)方法[6-10]来处理详细化学反应机理。数值求解详细化学反应机理的本质其 实是构建一个自输入变量(当前时间步的温度、压力、各组分的质量分数等)到输出变量 (下一时间步的温度、压力、各组分的质量分数等)之间的映射系统。该系统由一系列刚 性很强的常微分方程控制，而且包含众多变量(详细化学反应机理往往涉及几百甚至几 千个基元反应和中间组分) 。因此，求解刚性常微分方程组的运算量巨大、求解效率很 低。事实上，我们可以将输入变量至输出变量之间的关系视为一个“黑箱子” 的操作过程： 无需知道其中的本质关系如何，而只关心如何建立这种关系使得结果尽可能准确。这种 [6-10] 输入输出关系的可以通过 HDMR 方法来建立 。HDMR 方法作为一个捕捉该输入输 出关系的定量评估以及分析工具，得到了广泛的应用和发展[6-10] 。本文通过HDMR 方法 将求解关于详细化学反应机理的刚性常微分方程组转化为构造最优 HDMR 基函数以及 成员函数，从而极大地加快计算速度，并且能够保证计算精度。 本文首先简要地介绍了HDMR 方法；然后将HDMR 方法应用于模拟零维均质氢气 /空气的自着火过程，验证了HDMR 方法处理详细化学反应机理的精确性和高效性。 [6-10] 1 HDMR 方法 介绍 HDMR 方法的核心