湍流-化学作用的喷雾燃烧模拟 | 基于OpenFOAM的FGM模型实现与分析

引言

喷雾燃烧是内燃机研究领域中一个重要且富有挑战性的课题。本文重点讨论柴油喷雾燃烧，其特点是高温非预混燃烧。为了加深对内燃机的理解以便更好地对其进行设计，必须考虑详细的化学机理和TCI（turbulence-chemistry interaction）效应。准确地模拟非预混喷雾自点火和氧化过程以及污染物排放，特别是多环芳烃物种的演化过程，详细的化学计算至关重要。

许多TCI模型已被应用于喷雾火焰的建模。例如，输运概率密度函数(TPDF)方法、代表性交互火焰面(RIF)、火焰面/进度变量(FPV)模型、火焰面生成流形(FGM) 和建表火焰面模型(TFM)。

在这些湍流燃烧模型中，基于火焰面思想的模型具有计算效率高的特点，因此可以使用详细的化学反应动力学。火焰面方法的基本思想是，多维湍流火焰可以看作是嵌入在湍流流场中的被拉伸的一维层流火焰(称为火焰面)的集合。引入混合分数Z以消除非线性化学反应源项求解的困难。由此，化学可以在混合分数坐标下求解，然后映射到流场。基于火焰面的模型与化学建表方法相结合，通过将3D-CFD和层流火焰面计算解耦，降低了计算成本。这使得火焰面模型能够使用复杂化学反应机理，且计算成本相对较低。此外，基于火焰面的模型能够通过预设概率密度函数（PDF）有效地解释TCI现象。只当特征化学时间尺度比混合时间尺度短时，火焰面假设才是有效的，就像在大多数相关条件下类似柴油的燃烧一样。

本文使用FGM燃烧模型对正十二烷燃料的ECN sprayA进行RANS模拟。此外，由于传统观点认为高温非预混燃烧受限于混合过程，其进度变量的方差很大程度上依赖于混合物的形成速度，因此进度变量的方差经常被忽略。本研究考虑了进度变量的方差，类似于预混系统中进度变量的处理。

本研究的目的是为了增进对自动点火过程的了解，并揭示混合分数的变化和进度变量对自动点火过程和火焰结构的影响。

1、数值方法

1.1 气象模拟

FGM (Flamelet Generated Manifolds) 模型与火焰面方法具有相同的思想，即多维火焰可被视为一维火焰的集合。FGM模型的特征还在于存储和检索过程。除了混合分数之外，本研究引入了另一个控制变量，即反应进度变量来考虑不稳定过程，这可以捕捉喷雾燃烧中的自燃现象。本研究中使用的FGM模型可以总结为如下步骤：

(1) 一维火焰面原型的计算。

(2) 一维火焰面的解到控制变量空间的变换。

(3) 采用预设PDF法进行PDF积分。

(4) FGM表的存储。

(5) 根据控制变量及其脉动的三维控制方程求解得到数据库的索引，从FGM数据库中检索热化学量。

1.2 湍流化学相互作用

采用Bilger的公式计算混合分数:

上面Y和W分别是质量分数和分子质量。下标C、H、O分别表示碳、氢、氧元素，下标1、2分别表示纯燃料和纯氧化剂。

2、配置和计算设计

<