H2N2O2射流扩散火焰的直接数值模拟.pdfVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：094034 H2/N2/O2 射流扩散火焰的直接数值模拟 易富兴[1] 李德波[2] 卢树强[2] 樊建人* （1 浙江大学热能工程研究所 能源清洁利用国家重大实验室 浙江杭州 310027） E-mail：spkaka@，Tel: 0571 摘要： 直接数值模拟(DNS)是与雷诺时均模拟(RANS)、大涡模拟(LES)并列的三大模拟方法之一，而射流 扩散火焰在燃烧理论与实验研究中都扮演着十分重要的角色，本文采用完全可压缩的 Naive-Stokes 方程，对喷口直径为 D=1mm，Re=2000 的射流扩散火焰进行了直接数值模拟，分析了冷态时 H2，O2 的混 合；以及燃烧后的统计特性，主要包括速度场，密度，温度以及主要组分与混合分数的分布，并与实 验结果进行对比分析。 关键词： 射流扩散火焰 DNS 无反射边界条件 0 前言 如摘要所述，到目前为止总共有三类流动模拟方法，分别为雷诺时均法(RANS),大 涡模拟(LES),直接模拟(DNS)。引言的目的不在于比较这三类方法在流动模拟中的优劣， 只粗略介绍一下这三类方法在射流燃烧模拟中的一些应用：RANS 是一种成熟而且为工程 界所应用的一种方法，在湍射流扩散火焰的模拟中，主要采用的是 k-ε模型或雷诺应力 模型模拟流动，采用 PDF 模型或层流小火焰模型模拟燃烧，但模型对混合问题(Mixing) 模拟较差，与实验相比差距明显；LES 在最近十几年取得了飞速发展，如果检索最近几 年在 SCI 上收录的有关大涡模拟的文章，可以发现文章数近似几何级数增长[1] ，在 射 流 燃烧领域，大涡模拟也同样表现非凡，H.Pitsch[2]曾对 Sandia flame D(甲烷/空气射流， 导引火焰，Re10,000)进行过非常成功的大涡模拟；L.Selle[3]等人曾对一个小型化的旋 流燃烧器进行过比较成功的大涡模拟；但大涡模拟在射流燃烧领域的研究还有很多地方 不是很完善，特别是亚网格模型的验证非常困难。直接模拟(DNS)由于能给出所有的统 计信息而被视为最精确的模拟方法，但由于对计算机要求非常之高，已发表的直接模拟 三维射流燃烧的文章非常之少，Y.Mizobuchi[4]等人首次将 DNS 应用于射流扩散火焰领 域，但由于其网格数与格式的精度相去 DNS 甚远，所以也一直为人所诟病；T.Lu[5]等人 在 2008 年的 AIAA J 上发表了对 H2 湍射流抬升火焰(Lifted flame)的 DNS 研究结果，研 究所需网格数与格式的精度基本上达到了 DNS 的要求，为了完成这一研究，总共使用了 大约 10 亿左右的网格。 本文借助于本所现有条件，对 H2/N2 射流扩散火焰进行了直接数值模拟，模拟的管 内直径为 1mm，Re 大约为 2000，总共采用了 48 个 CPU，900 万左右网格。离散格式采用 了中心 8 阶显式中心差分格式，具体的计算方法与结果分析将在下面讨论。 1 计算方法 基金项目： 国家自然科学基金资助项目（Nos 计算中采用的方程为完全可压缩 NS 方程，详细地推导过程可见文献[6]，具体形式下： ?ρ ?ρu + i = 0 (1) ?t ?xi 2 N ?ρu j ?ρui u j ?p ? τ ij (2) + = ? + + ρ∑Yk f k, j ?t ?xi ?x j ?xi k =1 ?ρY ? k (3