分层旋流火焰特征提取与燃烧不稳定预测.pdf
《分层旋流火焰特征提取与燃烧不稳定预测.pdf》由会员分享,可在线阅读,更多相关《分层旋流火焰特征提取与燃烧不稳定预测.pdf(8页珍藏版)》请在咨信网上搜索。
1、燃烧科学与技术 Journal of Combustion Science and Technology 2023,29(5):):483-490DOI 10.11715/rskxjs.R202308024 收稿日期:2023-03-27 基金项目:国家自然科学基金资助项目(52106128);国家科技重大专项资助项目(J2019-0014-0057).作者简介:张世红(1997),男,博士研究生,.通信作者:韩 啸,男,博士,副研究员,han_.分层旋流火焰特征提取与燃烧不稳定预测 张世红,王欣尧,张 弛,韩 啸,周宇晨(北京航空航天大学航空发动机研究院,北京 100191)摘 要:燃烧振荡
2、是火焰释热脉动与压力脉动相耦合形成的不稳定现象,自发光图像是其研究中较易采集的一类数据本文采用手动提取、无监督和声压级监督神经网络三种特征提取方法,在瞬时图像和时均图像两个数据集上提取了火焰图像的低维特征,并分析了其显著性和物理表征,测试了其声压级预测精度和迁移能力研究结果表明,本文构建的声压级监督网络能够在复现火焰物理含义的基础上实现强鲁棒性、高精度、强迁移性的特征提取和振荡预测 关键词:燃烧振荡;数据集;特征提取;振荡预测 中图分类号:TK11 文献标志码:A 文章编号:1006-8740(2023)05-0483-08 Feature Extraction and Combustion
3、Instabilities Prediction of Stratified Swirling Flame Zhang Shihong,Wang Xinyao,Zhang Chi,Han Xiao,Zhou Yuchen(Research Institute of Aero-Engine,Beihang University,Beijing 100191,China)Abstract:Combustion instabilities are unstable phenomena arising from the coupling of heat release rate oscilla-tio
4、ns and pressure fluctuations.The chemiluminescence image is one of the most accessible data types in experi-mental combustion research.This paper adopted three feature extraction methods:hand-crafted,unsupervised neural network(USN),and sound pressure level supervision neural network(PSN).The featur
5、es were extracted from two data sets:transient image and time-averaged image.Based on the extracted features,their significance and physical characterization were analyzed,and their sound pressure level prediction accuracy and transfer ability were tested.The results show that the proposed PSN in th
6、is paper can achieve feature extraction and instabilities prediction with excellent robustness,accuracy and transfer ability based on the reproduction of the physical mean-ing of flame.Keywords:combustion instabilities;data set;feature extraction;instabilities prediction 贫油预混预蒸发燃烧是现代燃气轮机燃烧室实现低排放运行的一
7、个里程碑技术 贫油预混燃烧通过规避高温区来降低氮氧化物排放,但其释热脉动和压力脉动容易相互耦合诱发燃烧振荡1 为了分析燃烧振荡的诱因和发展,前人应用了一系列数据分析方法 常用的燃烧振荡数据分析方法有频域分析(傅里叶变换、小波变换)、模态分解(POD、DMD)、非线性动力学方法(相空间重构、庞加莱图)、局部瑞利指数 燃烧科学与技术 第 29 卷 第 5 期 484 等方法2 在这些分析方法之外,近年以深度学习为代表的机器学习算法在燃烧振荡问题的分析上也得到了初步发展 机器学习和频域分析、模态分解等都是数据驱动的方法 由于近年来深度学习算法和计算机算力的迅速发展,它在燃烧领域也取得了丰硕的成果,如
8、预测燃烧振荡、火焰形态、异常3-5,修正湍流模型、反应模型、加速数值计算6-10,以及实验数据处理和优化设计等等11-15 仅针对燃烧数据分析而言,卷积神经网络(CNN)和循环神经网络(RNN)可以完美匹配燃烧的图像信号和压力信号 卷积神经网络采用学习到的多个卷积核依次提取火焰的低级到高级特 征16;而循环神经网络采用“记忆单元”选择性遗忘和更新当前状态 近年来,在自然语言处理和图像处理领域 Transformer 及预训练模型发展迅速,但在燃烧领域应用还非常有限17 此外,目前流行的深度学习方法存在着数据获取难度大、解释性不足、泛化性难以保证等问题,这极大限制了机器学习算法在燃气轮机等昂贵的
9、工业产品中的应用 为了缓解这些问题,本文对神经网络提取的特征进行了相关性解释,并对其迁移能力进行了实验验证 本文采用了手动(hand-crafted)、无监督神经网络(unsupervised neural network,USN)、声压级监督神经网络(sound pressure level supervised neural net-work,PSN)3 种方法对瞬态、时均的火焰图像提取了特征,并以手动提取的特征为基准对神经网络特征进行了相关性解释,最后验证了利用瞬态图像数据集预训练的模型在多种实验工况、燃烧室结构的时均图像数据集上的迁移性能 1 实验系统与数据集 本文的数据来源于北航轴向
10、旋流独立分层燃烧器 BASIS(Beihang Axial Swirler Independently Strati-fied burner),其结构如图 1 所示 该燃烧器在常温常压下工作,甲烷燃料在与空气预混后经过主燃级和预燃级两级轴向旋流器进入火焰筒燃烧 火焰筒是一个等直径的石英管为了适应广泛的燃烧研究,BASIS 的旋流器、火焰筒等模块均可灵活更换 该燃烧器的燃油流量和空气流量通过 4 个质量流量控制器调节(两个 Alicat-HM-KM601 和两个Severstar,CS200),其控制精度约为 1%通过调节燃油流量和空气流量即可间接控制预燃级和主燃级的当量比 本文使用的数据有瞬态
11、图像、时均图像及压力信号火焰瞬态图由一个高速 CCD 相机(Photron,Fastcam SA4)采集,采集频率为 5000Hz,采集的图像大小为 768768 像素平方 时均图像使用配套CH 滤镜(4305)nm 的单反相机(NIKON D610)采集,光圈 1.8,曝光时间 1s,ISO100,焦距 50mm 压力信号采用动态压力传感器(PCB,112A22)测量得到,测得的时域离散压力信号经过快速傅里叶变换求得主频幅值,再将主频幅值转化为声压级作为火焰图像的标签 本文高速摄像拍摄的瞬时图像涵盖 112个工况,每个工况拍摄连续 600 张图像;单反拍摄的时均图像涵盖 1171 个工况,每
12、个工况拍摄一张图像 各个工况下的压力信号也被采集 图 1 BASIS燃烧器的结构示意 Fig.1 BASIS burner schematic 为了考察不同结构的燃烧器(套筒张角和火焰筒直径)、不同曝光时间(高速摄像和时均)、不同工况(变空气流量、燃油流量)下特征提取与振荡预测的有效性和迁移性 本文采用了两个差别较大的数据集,分别是基于瞬时图像的预训练数据集和基于时均图像的迁移学习数据集 这两个数据集及其工况条件如表 1 所示 表 1 瞬态和时均图像数据集及其实验工况 Tab.1 Two datasets and their operating conditions 参数 瞬态图像数据 集(预
13、训练)时均图像数据集(迁移学习)预燃级空气流量/(gs-1)2 1.5 主燃级空气流量/(gs-1)10 10.5 预燃级当量比 01.6 01.5 主燃级当量比 01 00.75 套筒张角/()52 28,40,52,70,90火焰筒直径/mm 100 100,140,180 火焰筒长度/mm 800 800 瞬态图像的预训练数据集通过改变预燃级和主燃级的当量比来获取,共计 112 个工况 时均图像的迁移数据集通过改变预燃级和主燃级的当量比、套筒张角、火焰筒直径来获取,共计 1171 个工况 以预训张世红等:分层旋流火焰特征提取与燃烧不稳定预测 燃烧科学与技术 485 练数据集为例,在不同预
14、燃级和主燃级当量比下的声压级和瞬态图像如图 2 所示 图 2 瞬态图像数据集 Fig.2 Transient image dataset 2 数据处理和算法设计 首先,本文对时域压力信号的处理方法是经过傅里叶变换求其主频幅值,再将主频幅值转为声压级作为利用火焰图像预测压力的标签:()5SPFFT20lg AMP(2 10)=Lp(1)本文对图像数据的处理方法是特征提取和声压级预测 特征提取即定义映射()f,对输入的r行c 列像素的图像+r cR,得到特征()=nfRxI,一般?nrc 这样就起到了对高维图像数据的“压缩”作用 本文的特征提取算法分手动(hand-crafted)、无监督(USN
15、)和压力监督(PSN)3 种 手动特征指的是对图像做人工的描述,如火焰图像的平均发光亮度、火焰亮区面积、周长等 这类特征有显著的物理含义,但真正直接有用的特征难以挖掘 神经网络在特征提取任务上比较高效,但往往以可解释性为代价 本文构建的手动特征提取方法(hand-crafted)、无监督特征提取器(USN)、压力监督特征提取器(PSN)、声压级回归器(regressor)的算法框架如图 3所示 2.1 手动特征提取算法 本文手动提取的特征有:表征总体亮度的亮度平均值、亮度方差,表征亮度分布的原点矩、中心矩、Hu 矩,表征火焰形状的火焰亮区周长、火焰亮区面积、等效半径,表征火焰明暗走向的方向梯度
16、频率,表征轮廓复杂性的分形维数等 30 个特征:T30hand,=?pqpqE S mHu G A R NFDx (2)表 2 是这些手动提取特征的数学定义 图 3 手动方法、无监督神经网络 USN、声压级监督神经网络 PSN和回归器 Fig.3 Hand-crafted,unsupervised neural network(USN),sound pressure level supervised neural network(PSN)and regressor 其中中心距和原点矩各有 4 个分量,分别是m00、m01、m10、m11和 00、01、10、11 Hu 矩有 7 个分量,均是对
17、二、三阶中心距的适当组合,使其对平移、旋转和缩放具有不变性 方向梯度频率是对各方向梯度上像素频率的统计,有 9 个分量 这样手动提取的特征共有 30 个分量,包括了图像的亮度、分布、形状、梯度、轮廓等多个方面的表征 需要说明的是,虽然这些特征对燃烧从多个角度去描述,但能直接反映燃烧振荡的物理量并不一定能被找到 这里定义多个类型的手动特征的目的是从多个角度对神经网络提取的特征进行解释,观察神经网络提取特征的侧 重点 2.2 无监督特征提取算法 无监督特征提取的思路是先采用编码器将燃烧图像压缩为低维特征向量,再通过解码器将低维特征向量复原为火焰图像 当解码器能真实复原原图像时则认为中间的特征向量是
18、本文需要的低维无监督特征4USN?x 无监督网络 USN 采用了变分自编码燃烧科学与技术 第 29 卷 第 5 期 486 器的总体框架,其编码器是一个卷积神经网络,用于提取本文关注的低维特征 而解码器部分采用 3 个并行网络:浅层网络,深层网络和 Sin 网络 虽然深层网络只有 6 个卷积层,但仍设计浅层网络和 Sin 网络作为残差连接来防止模型退化 表 2 手动特征及其定义 Tab.2 The hand-crafted features and their definitions 特征 类别 特征 名称 定义 亮度 平均值()()111,=rcijEi jrcII 总体 亮度 亮度 标准差
19、()()()2111,=rcijSi jErcIII 原点矩()()11,=rcpqpqijmi ji jII 中心矩()()()()1000110100,=rcppqijqimmjmmi jII 亮度 分布 Hu 矩 中心距()pqI的组合 火焰亮 区周长()()30=GGirthII 火焰亮 区面积()()()11,301(0)0(0)=rcijAHi jxH xxII 火焰 形状 等效 半径()()()2=RAGIII 明暗 走向 方向梯 度频率()()()111arctan120020=?rcijijNHrcxH xxIII 轮廓 复杂性 分形 维数()()()0lglimlg1=NF
20、DNI 2.3 声压级监督特征提取算法 压力信号监督神经网络采用了一个编码器和一个回归器,先采用编码器将图像压缩为低维特征,回归器仅用来对图像做声压级的预测 压力监督网络PSN 由卷积特征提取器和全连接回归器组成 卷积特征提取器的架构和无监督网络 USN 的编码器相同,用于提取声压级监督下的特征4PSN?x 回归器则由输入层、隐藏层和输出层组成,输入层单元数为4,隐藏层层数为 3,每层隐藏层单元数为 512,输出层单元数为 1,这样即可完成从特征到压力信号的 预测 2.4 回归器 无论手动提取的特征还是神经网络提取的特征,都需要确保其对声压级具有预测能力,因此需要一个回归器来拟合低维特征和声压
21、级 为了公平起见,手工和无监督网络 USN 提取的特征与压力监督网络PSN 提取的特征均采用了相同的回归器架构 即输入层单元数为 4,隐藏层单元数为 512,层数为 3,输出层单元数为 1,激活函数均为 ReLU 手动提取的特征handx和无监督网络提取的特征USNx通过以下步骤合并成 4 个分量同步预测声压级:用主成分分析将手动提取的特征handx降到 4 个分量18;将降维后的手动特征()4handPCAx和无监督网络 USN 提取的特征USNx合并成 8 个分量()4handUSNPCA,xx;将合并的 8 维特征再利用主成分分析降到 4 维()44handUSNPCAPCA,xx;将最
22、终的 4 维数据放入回归器预测声压级 本文采用两次主成分分析降维的原因是让手动提取的特征和无监督网络提取的特征具有相同的重要性,此处将特征整合的结果记为()PCA44handUSN=PCAPCA,xxx 3 结果与讨论 3.1 特征显著性及相关性解释 由于各个特征相对声压级不一定是线性相关的因此本文采用衡量两个变量同步单调性的Spearman 秩相关系数来判断其表征燃烧振荡的显著性,如图 4 所示 相关性分析显示,与声压级正相关性的手动特征主要有等效半径、亮度均值、亮度标准差、亮区面积、亮区周长、部分中心距和原点矩以及 40、280和320方向的像素频率等 负相关性较强的主要是 Hu矩.但在瞬
23、时图像和时均图像上它们和声压级相关性的强弱有所差别 这在一定程度上说明了:发生振荡时火焰趋于大范围、高亮度的燃烧;振荡燃烧时局部点有更大概率在火焰传播的前侧方向亮度增加.无监督网络 USN 和监督网络 PSN 提取的特征与燃烧振荡的声压级也有一定关联 由于采用了声压级监督,PSN 网络提取的 4 个特征(PSN1,PSN2,PSN3,PSN4)和声压级(SPL)存在极端的正相关和极端的负相关 无监督网络 USN 提取的一些特征和声压级的相关性有所降低,但仍然较为显著 为了观察神经网络提取到的特征所包含的信息,在瞬态图像和时均图像上分别以 PSN1 和 USN2 为基准,分析其提取的特征和手动特
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 分层 流火 特征 提取 燃烧 不稳定 预测
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。