气固两相燃烧双流体大涡模拟的模型.docVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：184058 气固两相燃烧双流体大涡模拟的模型 周力行 (清华大学航天航空学院工程力学系, 北京100084) (Tel: 010 E-mail: zhoulx@mail.tsinghus.edu.cn) 摘要：目前大多数两相燃烧的大涡模拟（large-eddy simulation，LES) 是欧拉-拉氏模拟，这比欧拉-欧拉或双流体模拟耗时大得多。本文提出气固两相燃烧双流体大涡模拟的数学物理模型，包括双流体框架内大涡模拟过滤的控制方程，两相亚网格（sub-grid scale，SGS)应力模型和亚网格气体燃烧模型。对亚网格应力本文建议用作者提出的两相亚网格动能方程模型。对挥发份和CO燃烧建议用作者提出的二阶矩（second-order moment，SOM) 亚网格燃烧模型。对颗粒的热解挥发和焦炭燃烧模型也在双流体模型的框架内加以应用。 关键词：双流体模型；大涡模拟；气固两相燃烧 0 引言 近年来两相燃烧的大涡模拟 (large-eddy simulation，LES) 引起了越来越多的注意，这是因为它可以给出瞬态流动和火焰结构，比雷诺平均模拟（Reynolds-averaged Navier-Stokes simulation，RANS模拟) 得到更准确的统计结果，同时又比直接数值模拟（direct numerical simulation, DNS）的计算量小得多，能用于解决实际工程问题。目前绝大多数的LES 用欧拉-拉氏（Eulerian-Lagrangian，E-L)模拟 [1-4]。从RANS 模拟的经验中已经了解到颗粒的拉氏模拟：(1)难以给出复杂几何形状的三维空间内各处的颗粒速度、温度和浓度的详细连续分布；(2) 常常低估颗粒的弥散 [5]; (3) 颗粒难以进入或离开边角回流区; (4) 要提高计算精度须计算大量的轨道，这就要有很大的计算量。另一方面，欧拉-欧拉（Eulerian-Eulerian，E-E）或者双流体模拟在两相流动燃烧的RANS 模拟中已经得到成功的发展和应用 [6]。它可以比E-L 模拟用少得多的计算量给出详细的颗粒信息。因此有必要研究双流体大涡模拟。 无论是E-L 还是E-E的 LES ，其模拟的合理性和精度取决于: (1) 足够细的网格; (2)比RANS 模拟用更精确的差分格式; (3) 合适的亚网格（sub-grid scale，SGS) 应力和燃烧模型。在E-L 模拟中，往往对气体用单相流动的SGS 应力模型，没有颗粒的SGS应力模型。最普遍应用的是Smagorinsky 涡粘模型 [7], Germano 动态涡粘模型 [8] 和Kim 的SGS 动能方程模型 [9]， 其中不考虑颗粒对气体SGS 应力的影响。Yuu [10] 和H. Zhou [11] 分别提出了有颗粒作用的气体SGS 应力模型。Poinsot 等 [12] 和 Simonin等[13] 报导了两相燃烧的双流体大涡模拟，其中对颗粒也采用了仿效气体Smagorinsky ________________________________ 基金项目：国家自然科学基金，批准号的模型，没有考虑两相之间的相互作用。本文中建议采用我们自己提出的两相SGS 动能方程模型 [14] ，用于两相燃烧双流体大涡模拟中。 对LES 中的气体燃烧模型，常常是在预混燃烧中用G 方程或线性涡模型，在非预混燃烧中用层流火焰面模型。然而这些模型只适用于特定类型或特定结构的火焰。过滤的概率密度方程（filtered PDF equation，FDF)模型的适用性比较广，但是在LES 中的计算量非常大。本文作者提出了二阶矩（second-order moment，SOM) 燃烧模型 [15], 已经成功地用于单相气体预混和非预混燃烧以及两相燃烧的E-L LES [16]，计算量不大。因此本文中将此模型用于气固两相燃烧的双流体LES。 本文将给出气固两相燃烧双流体大涡模拟的数学物理模型，包括双流体框架下特有的过滤后的控制方程，两相亚网格应力模型，亚网格气体燃烧模型，颗粒热解挥发和焦炭燃烧模型等。 气固两相燃烧双流体大涡模拟过滤的控制方程 取白噪音过滤，即计算网格内的体积平均，可以得到气固两相燃烧的双流体大涡模拟过滤后的气体连续和气体组分及气体和颗粒的、动量和能量方程: