基于颗粒质量实现考虑差异扩散的混合模型.docxVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：184048 基于颗粒质量实现考虑差异扩散的混合模型 杨天威1,2，周华1*，任祝寅1,2 （1 清华大学燃烧能源中心，北京 100084 2清华大学航天航空学院，北京 100084） (Email: zhouhua12@mails.tsinghua.edu.cn) 摘 要 在输运概率密度函数方法的湍流燃烧模拟中，考虑差异扩散效应对于提升预测精度十分重要，尤其是在火焰面燃烧状态下各组分长度尺度和扩散系数的明显差异会导致显著的差异扩散效应。差异扩散效应对于混合过程的影响体现为不同组分有着不同的混合时间尺度。当前被广泛应用的混合模型对所有组分使用单一混合时间尺度，无法考虑差异扩散效应，且如果直接使用多种混合时间尺度，则会违背可实现性条件β=1NsYβ(n)=1。本文提出一种通用的基于颗粒质量的混合模型实现方法，能够以多种不同的混合时间尺度作为输入参数，并且自然地满足了可实现性条件。用该思想改进了原始的IEM和修正Curl模型以考虑差异扩散效应，得到改进的IEM-DD与MC-DD模型，并从理论上推导它们的均值守恒性，在假定“颗粒质量不变”下推证了方差的指数衰减特性，这正是混合模型预期达到的组分统计值变化特性。采用封闭的惰性颗粒系统对IEM-DD与MC-DD模型进行数值实验，验证了均值守恒性，并且验证了每种 关键词 输运概率密度函数方法；混合模型；差异扩散； 平均交换相互作用(IEM)模型；修正Curl模型 0 前言 湍流燃烧是一种重要的能源利用方式，内燃机、燃气轮机和航空发动机等设备中化学反应都发生在湍流流动中。随着计算能力的极大提升，数值模拟方法在科学研究及工业应用中均发挥着越来越重要的作用。作为湍流燃烧数值模拟的一种方法，输运概率密度函数方法(TPDF)因为精确封闭了非线性化学反应源项，所以对于湍流-化学反应强烈相互作用的湍流燃烧问题能够实现准确预测。研究者曾用TPDF方法成功预测了湍流非预混火焰中的局部熄火/再燃现象和污染物排放[1,2]。然而，TPDF方法中组分的分子扩散项是未封闭的，需要使用混合模型进行建模[3,4]。 每种混合模型都指定了一种混合模式用以描述混合过程，并以给定的标量混合时间尺度进行混合。一些混合模型因为实现简易且性能良好，已被广泛应用于TPDF方法中，例如平均交换相互作用模型(IEM) [3]，修正Curl模型(MC)[4-6]，欧氏距离最小生成树模型(EMST)[7]和二项朗之万模型[8]。组分标量混合时间尺度{τβ, β=1,…,Ns}是混合模型的重要输入参数，τβ表征物质β的方差衰减的特征时间尺度，因此τβ的定义使得dσβ2dt=-σβ2τβ，即在τβ给定时，方差σβ2的衰减需满足指数形式，这是混合模型需要满足的重要条件。Subramaniam等人[7] 湍流燃烧过程中差异扩散现象有时会十分明显。对于湍流预混火焰，即使在高雷诺数条件下，反应-扩散平衡仍然可能存在，尤其是在火焰面燃烧状态下。由于组分标量的梯度受反应-扩散平衡影响，不同组分的长度尺度会存在差异。与此同时，不同组分的扩散系数也彼此不同。因此不同组分的混合时间尺度是不同的[7]。Richardson等人[9,10]通过研究湍流预混本生火焰的直接数值模拟中不同物质的混合速率，发现中间物质的混合速率可能高达主要产物的10倍。因此混合模型被期望能够对不同物质以不同的时间尺度进行混合，而当前的混合模型对所有物质仅使用单一时间尺度[11]，因此无法考虑由于长度尺度及扩散系数不同导致的差异扩散现象[7]。模型在使用不同混合时间尺度时会很容易遇到可实现性问题。Richardson等人[10]曾对IEM及EMST模型引入后续修正以考虑差异扩散并同时满足可实现性条件。但这种方法的修正过程较为复杂。在大涡模拟(LES)框架下，采用均值漂移(mean drift)方法处理组分标量因分子扩散导致的空间输运[12]能够解析过滤尺度以上的差异扩散效应，但对于过滤尺度以下的差异扩散效应仍需要改进混合模型加以考虑。 本文基于颗粒质量实现了考虑差异扩散效应的混合模型，能够以多种不同的混合时间尺度作为输入参数。这种实现方法避免了可实现性困难，不需引入后续修正，具有很好的鲁棒性，满足均值守恒性并且方差的衰减能够十分近似地满足指数形式。利用基于颗粒质量的思想，原始IEM和MC模型被改进以考虑差异扩散效应，分别记为IEM-DD