论文部分内容阅读
概率密度函数(probability density function,PDF)能有效地计算湍流与化学反应相互作用,已经成功地模拟了许多实验室尺度的火焰。通常PDF方法在网格中使用大量的颗粒以近似当地物理量的概率密度函数,并在颗粒上计算化学反应。由于化学反应计算的刚性,导致PDF方法模拟时会消耗较多的计算资源。 稀疏拉格朗日过滤密度函数(Sparse-Lagrangian filtered density function,稀疏FDF)方法,期望在保证统计结果可靠的前提下,使用的总颗粒数小于总网格数,从而大幅减小模拟所需的计算量。由于颗粒数量较少,使得稀疏FDF方法与传统FDF方法有许多不同:稀疏FDF方法无法通过网格内的颗粒系综得到瞬态的概率密度函数;稀疏FDF方法求解标量的FDF方程时,需要采用广义多维映射条件(multiple mapping conditioning,MMC)模型保证混合的当地性,同时需要采用全新的混合时间尺度模型;稀疏FDF方法中的密度耦合方法也需要重新设计。 本文针对稀疏FDF方法中的混合时间尺度模型以及反应流中的密度耦合方法进行了深入研究。将现有稀疏FDF方法中的密度耦合方法作了改进,以颗粒上等效焓源项的条件平均值,作为源项直接耦合到大涡模拟(large eddy simulation,LES)求解的等效焓输运方程中。在LES程序基础上开发了能够实现稀疏FDF方法的并行计算程序。并对Sandia火焰D、E、F进行了模拟。 Sandia火焰D、E计算结果与实验结果符合很好,表明改进的密度耦合方法比较合理,可以减少LES结果与实验结果间的偏差,同时也验证了稀疏FDF方法与程序的有效性。但Sandia火焰F的模拟结果与实验结果在局部熄火程度上存在差异,这表明现有的稀疏FDF方法与模型在预测熄火程度高的火焰时还存在有待改进之处。