圆盘钝体燃烧室速度场的大涡模拟和实验验证.pdfVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：084308 圆盘钝体燃烧室速度场的 大涡模拟和实验验证 郭培卿 臧述升 葛冰 （上海交通大学机械与动力工程学院叶轮机械研究所，上海市 200240 ） （Tel: 021E-mail: kimireggie@sjtu.edu.cn ） 摘要：在悉尼/Sandia 钝体的基础上，针对几何形状更为复杂的圆盘凸台钝体，对燃烧室进行了大涡模 拟，并利用 PIV 实验得到相应工况下的流场，以验证数值模拟的准确性。两种钝体几何形状相似但凸 台的存在导致流场分布的巨大差异。通过与PIV 实验结果以及不同精度网格下模拟结果的对比，成功 地将大涡模拟推广到流场更为复杂的凸台钝体流场。诸如空气在钝体回流边界处产生的切应力层等在 大涡模拟中出现的湍流流动的脉动特征，在实验中得到了验证，并且通过粗细两种网格的比较，表明 计算精度上具有进一步提高的空间。扩大了的计算区域基本消除了凸台边界对流场的影响，另一方面， 在流动特别复杂的区域，当前的网格精度已无法满足大涡模拟的要求，说明从简单形式到凸台钝体过 程中需要在网格精度、亚网格模型的选择等方面作进一步深入的研究。 关键词：大涡模拟；PIV ；圆盘钝体；瞬态流场 0 前言 目前，CFD 计算中常用的湍流模拟方法包括雷诺时均模型（RANS ）、直接数值模拟 （DNS ）以及大涡模拟（LES ）。与DNS 模拟方法相比，LES 的求解范围更大，对大量 存在的高 Re 数湍流及湍流燃烧现象的工程实际中进行有效地求解。另一方面， RANS 模拟方法对特性上具有本质区别的大涡与小涡使用通用的湍流模型，必然忽略和丢失反 映湍流特征的重要信息，而 LES 则对受边界条件和流动类型影响较大的大涡和基本各向 同性的小涡分别进行直接求解和模型求解，因此求解结果更为精确。此外，LES 有助于 揭示湍流结构及其与化学反应之间的相互作用机制，并且由于 RANS 仍然为工程应用上 的主要模拟手段，而 LES 可以用于检验 RANS 中的封闭模型的合理性问题，加上随着 计算机容量的进一步增加，LES 在工程计算方面具有双重重要的意义。 湍流燃烧大涡模拟近年来在国内外取得了显著的研究进展，已经有在一定条件下模 拟工程中的高雷诺数湍流流动和燃烧的实例，比如 D. Bissieres 等[1]对工业燃气轮机燃 烧室进行大涡模拟来验证和预测流场信息、Samer Abdel- Wahab 等[2]将大涡模拟应用于 旋转肋片管道换热研究，J. Janicka 等[3]对大涡模拟在热态流场中的应用作了阶段性总 结。国内方面，赵坚行等[4]对航空用加力燃烧室的冷态及热态流场进行了二维大涡模拟， 周力行等[5]提出的二阶矩亚网格燃烧模型并得到了实验验证。 由于在冷态条件下可以排除化学反应对流场的影响，同时又可以研究燃烧室内部流 场结构以及涡团运动等湍流特性，因此对冷态流场的计算和测量具有非常重要的意义。 本文利用大涡模拟方法对常温（300K 左右）下钝体燃烧室的速度场进行了数值模拟， 同时与PIV 实验进行对比，观察流场内部诸如涡的破碎等细致结构，同时了解所需要的 网格精度，为今后开展钝体燃烧室燃烧数值模拟以及亚网格燃烧模型的研究奠定基础。 1 大涡模拟 大涡模拟从湍流流动研究中发展起来，其基本思想为：在流场的大、小尺度结构之 间选择一滤波宽度对控制方程进行滤波，从而将方程中所有变量分成大、小两种尺度量。 对大尺度量进行直接数值模拟，对小尺度量则采用亚格子模型进行计算。因此，大涡模 拟针对大尺度结构得到的是真实的变量分布，而对小尺度结构而言，由于其具有各向同 性这一重要特点，对其采用统一的亚网格模型仍然是合理的。总体上来看，只要对控制 方程进行的滤波符合物理情况，那么大涡模拟的结果将是十分精确和真实的。 1.1 大涡模拟控制方程 直角坐标系下通过 Favre 滤波函数滤波，得到的不可压连续性方程