柴油机燃烧计算量研究.pdfVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：084210 柴油机燃烧计算量研究 郭晓平 王占杰 大连理工大学内燃机研究所 大连116023 Tel: 041184708246-8210 E-mail: gxp190@163.com 摘要 研究柴油机燃烧动力学计算量庞大的主要影响因素。组分质量分数太小或反应速率太高都使燃烧 动力学计算迭代步长减小，导致计算量剧增。引起组分质量分数小的关键是流动和扩散携带微量组分 迁移，反应温度是导致反应速率大的主要因素。缸内计算单元数太多是计算量大的又一主要因素。但 单元数与计算量之间是多项式关系，用高级计算机可以实现计算。而组分质量分数与计算量之间的关 系，已超出了依靠编写计算机代码能完成计算的能力。 关键词 柴油机燃烧，化学动力学，数值计算，计算量 0引言 从广泛的意义上讲，“计算”，已经和“理论”与“实验”并列，被普遍认为是人类 认识自然的三大支柱支一。随着计算机软、硬件技术的飞速发展，随着世界对能源与环 保日益尖锐地挑战，为了研究发动机的新型燃烧方式，以及掌握有害物生成机理，人们 需要开展发动机燃烧详细化学动力学计算。虽然计算机的运算速度目前已经达到每秒执 行百万条指令的水平，但是，这一计算速度仍然无法满足发动机燃烧详细动力学计算的 需要。 有文献估计，化学反应计算时间与组分数量的平方成正比。如果真是如此，那么发 动机燃烧计算属于 P 类算法[1]，即多项式时间运行算法，它可以在多项式时间内完成。 P 类算法被认为是可实际执行的算法。但实际上，发动机燃烧详细化学动力学计算目前 处于无法进行的地步[2] 。只将计算量看作组分数量的平方关系似乎过于简单。为此，本 文采用正庚烷氧化详细机理进行柴油机工作过程计算，剖析化学计算步骤，分析导致计 算量大的原因。 1柴油机工作过程数值计算 1.1 柴油机主要参数 4 气门直喷式四冲程压缩燃烧柴油机，缸径 135mm ，行程 150mm，压缩比 16.7， 转速 1500rpm，喷油量 73.5mg/cyc，喷油提前角 20 ℃A ，喷油持续期 18℃A 。喷油器轴 线与气缸轴线重合，喷孔 4 ×0.34mm 均布。喷孔夹角 150°。油束夹角为 15°，燃油喷 雾索特平均直径为 50μm 。喷油速度曲线如图 1 所示。 资助项目：AMD大学计划（2007-10） 图1 喷油速度曲线 图2 计算区域和燃烧室结构 1.2 燃烧模型与工作过程数值计算 简化起见，本文采用美国国家重点实验室研发的正庚烷氧化详细机理(LLNL 机 理)[3]，该机理由 556 个组分 2539 个基元反应式组成。 化学模型改进的 kiva-3[4]软件包用于缸内流动计算。其中湍流采用 RNG k-ε模型,燃 油喷射采用离散液滴模型(DDM)。将 chemkin 软件[5]与 kiva 软件结合起来实现化学与流动 的耦合计算。为了有效模拟柴油机缸内的进气涡流，计算从排气冲程的下止点开始。计算 区域如图 2 所示，包括气缸，燃烧室，进、排气支管，进、排气阀。规定压缩冲程上止点 为360℃A。 计算开始t=0时刻，缸内是温度为400K的空气。进、排气进、出口处压力恒为环境压力 0.10133MPa。喷入气缸的油温为400K。 从喷油开始，柴油机每运行5×10-7秒(约0.005℃A)就输出一次计算中间结果，将这个 时间段称为输出步长δt。 2柴油机燃烧计算 在多维计算中，采用“反应⎯流动”算法，即将流动和化学反应分开交替求解。用 kiva 软件计算缸内流场和其它热力学参数，然后调用 chemkin 气相子程序计算δt时间内 的化学反应。在化学反应计算中，各控制单元的体积、压力和温度保持不变。 调用 ckemkin 程序前，控制单元压力p ，温度 T，组