[发明专利]一种模拟固体燃料热化学转化的数值计算方法

说明书

本发明提供了一种模拟固体燃料热化学转化的数值计算方法，属于热化学领域，更具体地，该方法具体为：构建包含固体燃料和气体环境在内的几何模型；针对构建的几何模型建立控制方程；定义固体燃料的化学反应动力学模型，并根据其确定控制方程的源项，从而获得实际控制方程；基于努森扩散机理，对实际控制方程中不同气体的扩散通量进行修正，以获得修正后的实际控制方程；根据构建的几何模型和修正后的实际控制方程对热化学转化过程进行计算。本发明能够避免因可燃气体和氧化剂或气化剂向较高温度区域的迅速集中而导致计算模型局部出现严重超温崩溃，从而弥补现有固体燃料热化学转化数值计算方法的缺陷。

技术领域

本发明属于热化学领域，更具体地，涉及一种模拟固体燃料热化学转化的数值计算方法。

背景技术

在热化学转化领域，数值计算是一种能精确快速描述热化学转化过程的研究方法，但由于固体燃料的热化学转化过程十分复杂，运用数值计算研究固体燃料的热化学转化过程时，要充分考虑固体燃料的特点。

燃料热化学转化包括氧化剂或气化剂到达燃料表面和氧化剂或气化剂与燃料发生反应两个过程，这两个过程进行的快慢是决定固体热化学转化速率两个主要因素。氧化剂移动至燃料表面的过程指的是氧化剂的传质过程，依靠氧化剂的流动和扩散共同实现。对于氧化剂在气体中的传质，流动和扩散同样重要，然而，对于固体的热化学转化，氧化剂不仅需要到达固体表面，还需要穿过固体表面的惰性层与固体内部的活性组分接触才可能完成热化学转化，也就是说，固体的热化学转化还必须考虑氧化剂或气化剂在固体内部的传质过程，且固体的尺寸越大，固体内部传质过程的影响也越大。在固体内部，气体的流动阻力通常很大，氧化剂或气化剂在固体内部难以自由流动，氧化剂或气化剂的传质只能依靠扩散。要想对固体燃料的热化学转化进行模拟计算，不仅要考虑气体在气体中的扩散过程，还必须明确氧化剂或气化剂在固体内部的扩散能力。例如，对于煤颗粒，在颗粒表面的活性物质消耗殆尽之后，颗粒内部燃料被不可燃的紧实灰层包裹，灰层中，氧化剂或气化剂的流动阻力大，流动困难，气体的传质严重受阻，导致内部燃料无法接触氧化剂或气化剂，热化学转化速率明显下降。

因此，想要高效地模拟固体燃料热化学转化过程，就需要对传质这一物理过程进行简洁且准确的数学描述。但是，目前用于模拟固体燃料热化学转化过程的数值模型通常忽略了大量的物理过程，尤其是对传质过程的数学描述过于简略，使得计算结果难以真实可信。CN110727906A公开了一种基于多组分热解气体的固体可燃物燃烧特性求解方法，该方法给出了固体燃烧过程中热解气的燃烧过程，较全面地给出了固体挥发分分解燃烧的反应动力学过程，但该发明并没有考虑固体燃料较缓慢的传质特性对固体燃料燃烧的影响，这使得其对固体燃料燃烧过程的模拟是必然失真的，在失真的气体浓度场中，可燃气体和氧化剂或气化剂易向较高温度区域集中，迅速燃烧放热，导致计算燃烧数值计算模型局部出现严重超温，最终计算结果脱离实际，失去参考价值。综上所述，现有技术缺乏一种简洁的、准确的、鲁棒性高的处理固体燃料热化学转化的数值模拟方法。

发明内容

针对现有技术的缺陷，本发明的目的在于提供一种模拟固体燃料热化学转化的数值计算方法，旨在解决现有的固体燃料燃烧过程模拟失真的问题。

为实现上述目的，本发明提供了一种模拟固体燃料热化学转化的数值计算方法，该方法包括如下步骤：

S1构建包含固体燃料和气体环境在内的几何模型；

S2针对步骤S1构建的几何模型建立控制方程；

S3定义所述固体燃料的化学反应动力学模型，并根据其确定步骤S2中控制方程的源项，从而获得实际控制方程；

S4基于努森扩散机理，利用下式对所述实际控制方程中不同气体的扩散通量进行修正，以获得修正后的实际控制方程，