燃烧的数值模拟专题

燃烧的数值模拟专题[转自江苏省内燃机学会论坛]

内燃机等应用专业的人从事燃烧数值模拟研究的痛苦、挣扎与解脱

内燃机等应用专业的人从事燃烧数值模拟研究有巨大的痛苦，不象上世纪80年代和90年代处，只要拥有程序或软件，加上会使用，能够反复出论文，让人羡慕。之后，情况发生了很大变化，各种相关程序或软件层出不穷，程序或软件来源渠道广泛，燃烧数值模拟不再具有垄断性，从事实验研究的人也可以用用相关软件，纯粹进行燃烧数值模拟的人发论文越来越难，尤其是高水平论文，“程序等是别人的，你的贡献在那里？有无实验验证，。。。”等魔咒始终笼罩着你。由于内燃机等应用专业的人大多不是正宗流体动力学或相关应用数学等专业出身，改程序或推出程序（作出创新性贡献）谈何容易，我们陷入了痛苦和挣扎。

数年前，当数值模拟的人还有些吃香的时候，我对我的学生说过，你只能用现有软件去算工程上的一些算例，如不尽快开发自已有特色的东西，吃饭饭碗难保，现在这种观点逐步成为现实。

我做过多年KIV A，可以说非常熟悉它，对KIV A进行过大量改造，由于KIV A的复杂性，一动百动，每次改造都化费了大量心血，但产出甚少。

是否是内燃机等应用专业的人不适合从事纯粹的燃烧数值模拟研究？其实，内燃机等应用专业的人有其知识面宽广等优点，我现在的观点是，如果你用某一软件或对一程序进行艰苦的努力，仍不能取得成果上的很好回报，放弃它，人生苦短，我们需要重新定位，如果仍坚定数值模拟方向，是否需要在有关燃烧数值模拟这样一个越来越大的领域中找到自已的一席之地，吃碗“安心饭”，趁年轻时找到它，我等已老矣！

个人意见，仅供参考！

中国科技大学火灾科学国家重点实验室

蒋勇

大家好！

有同学问：

"蒋老师，看了您的这么多帖子，我感觉似乎我们学内燃机的，要想采用某种软件模拟缸内的燃烧过程，特别是想利用chemkin来模拟缸内的燃烧，好像是行不通的，因为两种燃烧的最基本类型就不是一样的，chemkin是零维的，而缸内无疑是多维的，还有燃烧速率的问题。那么，我们要想利用chemkin在柴油机上做点文章，该如何下手呢？毕竟进行小火焰的模拟不是我们的主业。期待蒋老师的解答，谢谢。"

回答：

【1】内燃机燃烧本质上是多维湍流燃烧，但内燃机燃烧数值模拟有零维模拟、准维模拟和多维模拟，不同模拟方式获得的信息不一样，比如零维模拟不可能研究湍流脉动、结构等等对燃烧的影响，但它可以重点研究化学反应动力学等方面的问题，这样当然可以采用CHEMKIN软件，但模拟结果的适用范围要注意，比如你获得了NOX排放，你要注意，这个NOX数据是在某个压力和温度下获得的，是均质层流燃烧背景，但往往够用了，这样的话，你当然不再需要多维计算。再举个例子，如果你对缸内的“VORTEX SHEDDING”等现象感兴趣，二维模拟都不行，你只能进行三维模拟。非常重要的是：研究目标取决于研究手段。

【2】湍流燃烧模型问题。小火焰是湍流燃烧模型的一种。如果你要进行缸内多维燃烧模拟，湍流燃烧模型你就回避不了，如果你采用FLUENT、STAR-CD、PHOENICS等商用软件，你当然可以根据需要，选取里面合适的湍流燃烧模型；如果你要开发相应程序，必须考虑湍

流燃烧模型问题。

【3】燃烧数值模拟研究问题。燃烧学界包括内燃机专业，不能认为小火焰等模型问题与我们无关，只要你搞燃烧数值模拟，它们就是你的主业，这方面我深有教训。我以为，软件应用不等同于进行燃烧数值模拟研究，同时，不能对湍流问题和湍流燃烧问题的进行深入学习和理解，对内燃机燃烧数值模拟也不太可能作到高水平。确实，燃烧数值模拟研究者非常辛苦。

中国科技大学

蒋勇

一名老校友的有关内燃机燃烧过程一名老校友的有关内燃

机燃烧数值模拟研究的一点想法

看了一下本论坛，觉得办得很好，作为江大内燃机老校友，提几个问题，供大家讨论和批

判：

【1】内燃机燃烧过程是湍流燃烧，而湍流燃烧理论和数值模拟的核心困难和问题在哪里。【2】目前商用和通用软件很多，如FLUENT,STAR-CD,CHEMKIN,KIVA等，这些商用和通用软件中的算法和模型并不是目前国际上先进的东西，如果要用这些商用和通用软件，其重点是否应放在对未知或仍不清楚的现象的计算上，而不是对一般现象的重复研究上。【3】目前国际上相应的先进算法和模型有哪些，研究中是否需要实现这些前沿技术，用它或商用（通用软件）来计算同一问题，结论会不会不同。

中国科技大学火灾科学国家重点实验室蒋勇

谈谈多维燃烧计算程序与CHEMKIN的耦合

大家好！

我这里所说的多维燃烧计算程序通称“KIV A”、"FIRE"，等软件。

有许多同学希望计算生物柴油等，所以关心多维燃烧计算程序与CHEMKIN的耦合，其实，我记得在我发过的一个帖子中解释了该问题，为清楚起见，现重新解释如下：

有同学问：“将其动力学文件CHEM.INP、热物性文件THERM.DAT、和输运文件TRAN.DAT“整合”成一个文件然后在“多维燃烧计算程序”中进行计算，不知道这样理解是否

正确？

回答：很不全面。要理解它，要首先知道：CHEMKIN中的燃烧速率计算并没有考虑湍流脉动对燃烧速率的影响，而内燃机中的燃烧无疑是湍流燃烧。所以，如果单纯（或直接）调

用CHEMKIN中的某模块或函数来计算燃烧速率，肯定有问题，除非你的计算网格和时间步长极小，类似直接数值模拟，但需要超级计算机。

所以，你必须面对湍流燃烧模型，现在的湍流燃烧模型很多，你可以选取，在确定湍流燃烧模型后，要编程，编程时可以调用CHEMKIN中的库函数或模块，来实现湍流燃烧功能。举个例子，内燃机中的湍流燃烧计算类似盖楼，CHEMKIN中提供了大量的建筑材料，你必须根据设计（某个湍流燃烧模型），将建筑材料有机地排列好。

顺便说一句，有同学问：湍流燃烧模型应该怎样选择？

这个问题太大，由于不同湍流燃烧模型的特点、应用场合不同，没有统一的答案。

中国科技大学

蒋勇

耦合燃料复杂机理的燃烧模拟

有同学问，我在前些天所说的“内燃机缸内燃烧的Reynolds数较高，模拟计算的量很大，对硬件要求极高”，而有论文的计算采用了几千步的反应机理，为什么能计算呢？非常高兴能

提出问题，现解答如下：

要注意，这些论文一般采用了纯粹的CHEMKIN来计算，维数为零维或一维，并且是层流燃烧计算。我前些天所指的是针对缸内的多维湍流燃烧计算，并且是采用直接模拟或大涡模拟。对于均质各向同性三维非定常的湍流燃烧直接模拟，通过严格的理论分析，所需的计算网格是Reynolds数的3次方，大涡模拟的要求稍低，但其计算量也是十分巨大。但采用RANS方法进行计算，如KIV A，则计算量可大大降低。如果要进行耦合燃料复杂机理的多维湍流燃烧模拟，即使是RANS方法，其计算量也很大，主要原因是：各基元反应的反应特征时间相差巨大（除了氢气等非常简单的燃料），再加上湍流-化学之间的刚性问题。目前已有的一些论文称，进行了耦合燃料复杂机理的多维湍流燃烧模拟，并且仅使用了一般高性能微机，我个人认为，计算结果值得怀疑，因为背离了严格的理论分析结论。说到这里，顺便提及，没有较为坚实的理论和较为深入的对相关机理的认识，搞好工程应用，也

是有难度的。

个人观点，仅供参考！

中国科学技术大学火灾科学国家重点实验室

蒋勇