模型燃烧室污染物生成的大涡模拟.pdfVIP

江苏省工程热物理学会第一届学术会议 编号； T3—35 模型燃烧室污染物生成的大涡模拟 颜应文赵坚行 张靖周 刘勇 (南京航空航天大学能源与动力学院，南京，210016) 摘要：在i维任意曲线坐标系下采用欧拉一拉格朗日两相大涡模拟方法对模型燃烧室气液两相瞬态 喷雾燃烧过程进行数值模拟，并采用业网格污染物生成模璎模拟燃烧室污染物生成过程。计算结果与 PⅣ测量的瞬态速度场以及出口温度和污染物浓度分布实验数据进行比较，表明在j维贴体坐标系下 采用欧拉一拉格朗日两相大涡模拟方法，数值模拟模氆燃烧室两相喷雾燃烧流场及污染物生成过程， 所采用的哑网格污染物生成模型可以用于燃烧宅气液两相喷雾燃烧流场中污染物排放的预估。 关键词：欧拉一拉格朗臼两相大涡模拟，k方程亚网格尺度模型，亚网格EBU燃烧模型，亚网格污染 物生成模型，随机离散模型，PIV 近几十年来军用和民用航空燃气涡轮发动机的迅速发展，发动机性能参数日益提 高，燃气温度越来越高，发动机排放到大气中的污染物日益增多。为减少航空发动机污 染物的排放量，要了解航空发动机污染物的生成机理。本文采用数值计算与试验测量相 结合的方法对某航空发动机燃烧室研究污染物主要成分CO和NO的生成机理。 目前数值模拟方法…根据湍流数值模拟精细程度不同，通常有直接模拟(DNS)、大 数的湍流流动，RANS方法理论成熟，但对湍流及燃烧过程的细微结构了解不够，因此 本文采用大涡模拟方法【2刈来模拟某航空发动机燃烧室气液两相喷雾燃烧过程及相应污 染物生成。 本文的研究对象为带二级扩压器、火焰筒及内外冷却通道的矩形燃烧室整体流场， 结构形状非常复杂，为了避免外形与坐标轴线不一致而带来的计算误差，因此在任意曲 线坐标系下生成三维贴体网格”’引，并进行相应的流场计算。本文数值分析了气液两相 燃烧过程，通过数值分析结果与PIV测量的速度矢量图以及燃烧室出口截面温度和污染 物浓度分布实验数据比较，表明本文采用的各种亚网格尺度模型可以较好的模拟航空发 动机中污染物的排放过程，为设计高性能、低排放污染物的航空发动机提供设计依据。 1．研究对象 为了加工方便，本文把环形燃烧室简化为矩形燃烧室，如图l所示，该燃烧室包括 二级扩压器、内外冷却通道以及火焰筒。火焰筒上下壁丽上各开有l排主燃孔、l排掺 混孔和8排气膜冷却孔，主燃孔为3个，掺混孔为2个，气膜孔是由许多小孔组成。数 值模拟实际燃烧实验工况，其进口气流参数为：温度500K，压力0．13Mpa，流量为 0．28吲s，航空煤油从喷嘴中喷出，与双级旋流进来的空气混合、燃烧，进口油气比为 国家自然科学基金资助项目(编 614 001518。 圉l模型燃烧宦结构小意陶 2．气相大涡模拟基本控制方程 对于两相喷雾燃烧流场，通过滤波后，气相人涡模拟基本控制方程在贴体坐标系下 通用形式为…o】： z％怕^)]+ ，等+善(∥q+南(，r嘞+毒(，脚)=南[争(“吒+m (1) 对于奉文所研究的三维两相喷雾燃烧流场，采用k方程亚网格尺度模型来模拟湍流 亚网格粘性，亚网格EBu燃烧模型来模拟化学反应速率。上式中通用变量中分别为速 度“、v和w，亚网格动能七昝，焓^，组分儿，六通量辐射模型中的i个辐射通量R，、 E、见。‰代表气液两相之问相互作用的源项，。吼为各方程源项，各方程的源项 见文献【10]， 对于能量方程中的化学反应速度瓦，本文采用亚阚格EBu燃烧模型”川估算： 吃=面n(以。，吒。) (2) 2I‘elp(一E，月丁)。 其中吒=√2∥“％。五)一(鸭舱【燃料D，‰=4。p 有效粘性以=∥+一，，其中∥为层流粘性，H为亚阿格尺度粘性， 投靠的定义见文献[8～lO