基于FGM的湍流射流扩散火焰超大涡模拟.docxVIP

收稿日期： 修订日期： 作者简介：夏朝阳 博士研究生 计算流体力学与湍流燃烧 基金项目：国家自然科学基金（No，江苏省自然科学基金（No.B和航空动力基金（No.6141B090316）项目资助 中国工程热物理学会 燃烧学 学术会议论文 编号：184219 基于FGM的湍流射流扩散火焰超大涡模拟 夏朝阳，韩省思*，毛军逵 （南京航空航天大学能源与动力学院，南京 210016） ( Email: xshan@) 摘 要 为了研究超大涡模拟(VLES)预测湍流燃烧问题的能力，本文采用VLES结合基于假定概率密度函数的火焰面生成流型（FGM）建表湍流燃烧模型对值班甲烷/空气湍流射流扩散火焰（Sandia Flame D和Flame F）开展数值研究，并与实验结果进行比较。结果表明，VLES-FGM可以较准确地预测出湍流射流扩散火焰中的燃烧过程。 关键词 超大涡模拟；火焰面生成流型；假定概率密度函数；湍流扩散火焰。 0 前言 近些年来，随着计算机能力地不断提高，湍流数值模拟方法迅速发展，成为一类热门的研究方向。在航空和航天工程领域存在着大量的湍流问题，尤其在发动机燃烧室中更是具有复杂的湍流燃烧问题，这对高精度数值模拟方法提出了很大的挑战。 湍流数值模拟方法主要可以分为三大类：雷诺平均方法（RANS），大涡模拟（LES）和直接数值模拟（DNS）。近些年来，大涡模拟由于可以直接求解出湍流大尺度结构，在流动、声学及燃烧领域得到了很多研究应用。大涡模拟需要的网格尺度直接和可解湍流尺度相关，而对于高雷诺数湍流存在很多小尺度运动，特别在近壁面区域需要很高的网格分辨率。因此，高计算资源消耗限制了大涡模拟在高雷诺数工程湍流问题中的应用。 混合RANS-LES方法是近些年来发展迅速的一类方法，它能够结合不同方法的优势，使得大涡模拟应用于高雷诺数工程湍流问题成为可能[1]。研究者们已经提出了很多不同类型的混合方法，其中分离涡方法（DES）[2]被广泛应用与解决工程流动问题，还包括近些年发展的统一湍流模化方法如VLES[3,4]，PANS[5]以及PITM[6]等方法。统一湍流模化方法主要根据局部的网格尺度和湍流尺度实现RANS，LES和DNS之间的转变，通过一种模化框架涵盖三种湍流模拟方式。这种模化方式的主要优势是可以应用于不同湍流模拟，特别是高雷诺数湍流的工程应用，但对高雷诺数化学反应流的研究仍较少。 本文针对美国Sandia国家实验室的射流扩散火焰 [7]，采用VLES方法结合基于beta-PDF的FGM建表湍流燃烧模型[8]开展数值模拟研究，并与实验结果进行比较，以验证VLES-FGM对高雷诺数化学反应流的预测精度。 1 数学物理模型 1.1 湍流模型 本文的湍流流动计算基于标准k-ε湍流模型的超大涡模拟[3]。超大涡模拟是在RANS模拟的基础上，通过分辨率控制函数Fr对湍流粘性μt进行修正，参数Fr是超大涡模拟的核心。Fr的形式如式（1.1）所示，其中Lc，Li和Lk分别为截断长度尺度，积分长度尺度和Kolmogorov长度尺度。 Fr=min (1. SEQ eq 1) Fr的取值根据网格分辨率的不同在0到1.0之间变化，来决定湍流的模化程度，并可以实现在RANS，LES和DNS之间的转变[3,4]。 VLES k-ε模拟的控制方程如式(1.2)-(1.4)所示，其中湍动能k和耗散率ε采用标准k-ε湍流模型的控制方程。湍流粘性系数通过分辨率控制函数进行调整。 D (1. SEQ eq 2) D (1.3) μ (1.4) 其中，D()/Dt为物质导数。 1.2 燃烧模型 本文采用火焰面生成流型(FGM)结合假定概率密度函数分布（beta-PDF）的建表湍流燃烧模型，采用一维层流火焰中的流型代替多维火焰，基本思路和稳态层流火焰面方法（SLFM）相同，都是将复杂、多维的化学反应信息映射到一个或者多个特征变量中以实现降维，这样同时考虑了详细化学反应机理的同时，又减少了求解输运方程的个数，降低计算量。本文针对Sandia Flame D和F采用部分预混FGM，基于特征变量混合物分数和反应进度变量构建部分预混FGM层流表。采用假定PDF方法使得流动和化学反应解耦计算，混合物分数Z和反应进度变量C的beta-PDF函数分别由其平均值和方