一种混合多尺度燃烧动力学模拟方法.pdfVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：104259 一种混合多尺度燃烧动力学模拟方法 1* 2 3 2 苟小龙 ，孙文廷 ，陈 正 ，琚诒光 （1. 重庆大学动力工程学院，重庆400044；3. 北京大学工学院，北京100871； 2. 普林斯顿大学机械及航天工程系，新泽西08544，美国） (* Tel:023 Email:simgxl@) 摘要：本文介绍了一种高效模拟涉及详细化学反应机理燃烧过程的混合多时间尺度模拟方法 (Hybrid Multi Time Scale, HMTS) 。在模拟过程中通过动态地计算各个组分的特性时间，并依据其特征时间将不同基元组 分和不同基元反应划分为不同的计算组，各个计算组采用不同的时间步长进行显式计算，且将最小特征时 间的一个或几个组采用隐式算法进行计算。基于该方法采用详细化学机理，对癸烷/空气、甲烷/空气的自然 着火过程进行模拟，并将模拟结果与欧拉法、隐式积分算法进行了比较。结果证明该方法能够精确地模拟 组分的质量份额和温度等的动态过程，能够得到精度的着火迟延时间。同时，与隐式积分算法相比，HMTS 方法大约能提高一个数量级的计算速度。因此该方法的应用将有助于实现燃烧问题的高效精确数值模拟。 关键词：混合多时间尺度模拟；详细化学反应机理；甲烷；癸烷 0. 引 言 随着计算机技术的发展以及相关高效的计算方法的出现，数值模拟得到了越来越广泛的 应用。目前数值模拟方法已经成为反应流分析的最重要的手段之一，大大促进了先进发动机 和推进技术的发展[1, 2] 。然而由于燃烧过程中反应流存在的多物理尺度和时间尺度问题给数 值模拟带来巨大的挑战，成为目前反应流模拟的巨大障碍。首先，由于碳氢燃料的详细动力 学机理常常包含成百上千种组分和数千种基元反应，所以反应流的模拟常常需要巨大的计算 量。Chen[3] 的研究证明表明，对于甲烷（36种组分、217种基元反应）的反应流模拟中，大 约80%的CPU时间花费在燃烧过程的模拟中。其次，反应流中各种组分时间尺度之间具有巨 大的差异，该强刚性问题给计算同样带来极大的挑战。采用如欧拉等显式算法，需要采用各 种组分中的最小时间步长进行计算，将会导致难以容忍的计算量。 为了解决基于详细反应机理的反应流计算中的问题，研究人员进行了不懈的努力。为了 降低详细机理带来的基元组分数目过大的问题，一些机理简化方法应运而生[4, 5, 6] ，并取得了 很好的简化计算的效果。然而不容乐观的是，对于常见碳氢燃料的精确模拟常需要大量的组 分和基元反应，尤其对低温燃烧过程的模拟，大大限制了机理简化的程度。为了降低燃烧过 程模型的刚性，目前常采用基于特征值分析的计算奇异摄动法（Computational Singular Perturbation, CSP ）以判别准静态（Quasi-steady State Species, QSS ）组分的方法得到了广泛的 应用。然而在温度、压力等初始工况偏离预定参数时，准静态假设的有效性难以保证，尤其 在低温过程中。 针对反应流计算过程中的强刚性问题，一些高效的隐式计算方法[7,8]得到了广泛的应用。 Curtiss和Hirschfelder[9]证明这些隐式的计算方法尤其是先后差分法（Backward Difference Formulas, BDF) ，其性能均好于显式法。近年来，各种基于Gear算法的隐式求解算法得到了 快速发展[10] ，其中，VODE[11]算法的应用最为广泛，该方法已经应用于商业软件CHEMKIN[12] 。 虽然如VODE等算法既适合非刚性问题，也适合求解刚性问题，但是对于非稳定的燃烧问题 的求解，其计算量常常依然难以承受。同时，由于隐式算法的采用，将不可避免地丢失特征 时间小的组分的中间过程，如果该组分对于整个过程的分析具有重要的作用，那么这些隐式 计算方法将不再适用。为此，有必要开发一种能够适用于多组分、多反应