论文部分内容阅读
超声速燃烧冲压发动机因为能够实现高超声速飞行(5<Ma<15)成为近年来的研究热点。由于超声速燃烧实验测量手段复杂、实验设备昂贵,计算流体力学方法成为研究超声速燃烧冲压发动机的重要手段。超燃冲压发动机燃烧室内存在复杂的激波/边界层、激波/火焰的相互作用以及自点火、局部熄火、再点燃、火焰驻定等非定常的燃烧过程,发展超声速燃烧模型和数值方法对研究超燃冲压发动机内的流动和燃烧特征非常重要。本文主要开展以下几个方面的研究: (1)发展基于化学热力学建表的可压缩湍流燃烧模型;(2)基于发展的可压缩火焰面进度变量方法对超声速支板射流DLR燃烧室的燃烧特征分别开展RANS和LES数值研究;(3)采用LES方法结合可压缩修正的自点火燃烧模型分别研究高焓值超声速横侧射流Gamba燃烧室的混合特性和燃烧特性。超声速流场中的可压缩效应如压缩/膨胀过程产生的压力、密度和温度的变化,对化学反应过程产生重要影响。因此,将低马赫数下的化学热力学建表方法应用到超声速可压缩流时需要考虑高速流体的可压缩性。本文对化学热力学建表模型的可压缩修正方法进行了深入研究。本文引入了温度和压力修正方法。在可压缩化学热力学建表方法中,通过直接求解能量方程得到温度值,可以耦合部分可压缩性影响;通过分析不同压力下层流化学热力学数据表,给出反应进度变量源项的压力修正系数,对反应进度变量源项进行压力修正。此方法的优点是能够在不增加化学热力学数据表的大小的基础上考虑流场的压力变化,并可以推广应用到采用反应进度变量建表的模型中;超声速横侧射流燃烧算例中,在以上温度与压力修正的基础上,提出了修正自点火模型中初始温度的方法,考虑燃料射流近场由于压缩/膨胀引起温度的不均匀分布。采用可压缩火焰面进度变量方法对超声速支板射流DLR燃烧室三维流场和燃烧场开展RANS和LES数值研究。相比与RANS方法,LES方法得到的冷态场中压力分布、波系分布、速度分布以及燃烧场中波系分布、速度分布、温度分布结果均与实验值符合更好。LES方法可以预测超声速中的大尺度湍流拟序结构以及非定常的流动及燃烧过程,能够捕捉激波与边界层的相互作用、湍流混合、火焰驻定以及熄火再燃等问题;LES结果表明,钝体两侧剪切层内形成了熄火再燃的不稳定火焰,在燃烧室中心回流泡内形成了稳定的富燃部分预混火焰;在LES框架下,分别采用的ε与β两种概率密度分布对反应进度封闭。对比结果发现,假定概率密度的β分布计算结果与实验符合较好,假定概率密度的δ分布没有考虑反应进度变量的亚格子脉动,得到平均温度偏高。因此,反应进度变量的亚格子脉动对预测超声速燃烧非常重要。采用LES方法研究了超声速横侧射流Gamba燃烧室内的大尺度序结构和混合特性。采用混合分数概率密度函数分析燃料射流喷口上游回流区、射流近场以及射流远场的混合特征,并探讨射流动量通量比对流场结构、射流穿透深度、标量分布以及混合效率的影响。结果表明,低射流动量通量比时,反向旋转涡对CVP结构与壁面边界层相互作用更强烈,燃料的质量分数沿着横向方向扩散的更快,混合效率更高。考虑到超声速来流温度高于燃料的自点火温度,采用修正的自点火燃烧模型研究了Gamba燃烧室燃烧场内的流场结构、射流穿透深度、涡结构以及燃烧特征。由于反应放热,燃烧场中,反向旋转涡对CVP与尾迹反向旋转涡对TCVP结构变大、射流穿透深度提高。计算得到燃烧场三个燃烧反应区域:射流入口上游的点火点、射流剪切层以及近壁面射流尾迹区。近壁面射流尾迹区火焰稳定主要由尾迹反向旋转涡对TCVP结构控制。