三维典型超声速燃烧室流场数值模拟

超声速燃烧是实现以吸气式冲压发动机为动力进行高超声速飞行的关键技术之一.本文对三维超声速燃烧室流场进行了基于有限体积方法的并行数值模拟,还对JPNal三维典型超声速燃烧室内的激波诱导燃烧流场进行了数值模拟.控制方程为三维多组元带燃烧反应模型的Euler方程,空间离散采用二阶精度VanLeer迎风通量分裂格式,时间推进为考虑了化学反应特征时间的5步Runge-Kutta方法,燃烧反应为氢气/空气十反应模型.为提高计算效率,采用了基于区域分裂技术的并行计算.最后对比分析了单、双气态氢气喷流对燃烧室内流场结构及燃烧特性的影响.     

高超声速侧压式模型进气道不起动特性分析

对某典型高超声速侧压式进气道三维流场进行了数值分析，就壁面压力分布、波系结构和近壁流谱图与实验结果进行了比较，计算反映了流动的基本特征。分析了在不同来流马赫数下的流动特征，随着来流马赫数的减小，激波角增大，压缩波在通道内的反射次数增加，而强度逐渐减弱，总压恢复系数逐渐增大。当马赫数减小到一定数值时，在等直隔离段入口出现喉道截面，进一步减小来流马赫数，流量阻塞，引起压力升高，波系向进口方向移动，波后出现亚声速流场，进气道不起动。同时还发现当不起动现象发生时，由于波后分离包的存在，在进气道的进口前形成向后倾斜的激波而不是正激波。本文还提出了一种确定不起动马赫数的方法。     

氢/空气超声速燃烧研究

Ｈ２／Ａｉｒ在两种不同的燃烧室尺寸、七种燃烧喷注方式下进行了系统的超声速燃烧实验。实验空气的滞止温度在２０００Ｋ左右，滞止压力１～１．４ＭＰａ，总流量２ｋｇ／ｓ，燃烧室进口马赫数２．５，可以模拟飞行Ｍ数为７的超燃冲压发动机中的燃烧工况。新开发的一维超声速燃烧程序ＳＳＣ－１可以估算出燃烧室内的流场参数、燃烧效率和总压损失。计算结果与实验进行了比较，发现较好的一致。实验结果表明，利用垂直喷射，燃烧效率可以超过８０％，同时不引起严重的总压损失。由燃烧室壁面静压分布与燃烧效率的分析发现，燃烧室燃料注射位置应避免过于集中，宜分散按规律分布，使燃烧室静压分布尽量平直以获得高燃烧效率。     

马赫5一级定几何二元混压式进气道反压特性试验研究

对一种设计马赫数为5一级的定几何二元混压式亚燃冲压发动机进气道进行了风洞试验研究,得到了该进气道的反压特性,结果表明：设计状态时,随节流锥堵塞度的增加,进气道出口反压比不断增加,马赫数逐渐下降,总压恢复系数先下降后上升,通道内气流脉动的功率谱密度无明显峰值;节流锥堵塞比为72%时,发生喘振,喘振基频约为48 Hz;随节流锥堵塞比的降低,进气道喘振基频逐渐降低,进气道结束喘振后结尾激波先到达进气道进口处,然后稳定在进气道内收缩段内,随着节流锥堵塞比的进一步降低,结尾激波逐渐进入进气道扩张段。     

某组合发动机进气道抽吸作用分析

对某组合发动机进气道设计马赫数下的三维流场进行了数值模拟。研究了不同位置抽吸对进气道流场及性能的影响,研究结果表明,外压部分抽吸能在较小的抽吸流量下有效改善进气道的起动性能和气动性能,而又不改变进气道的内通道特性;内通道抽吸均化了喉道出口流场,提高了总压恢复系数,但减弱了压缩程度,造成压比、温升比较大幅度地下降,喉道出口马赫数增加,且抽吸流量越大,压比、温升比的降幅和喉道出口马赫数的升幅越大。     

三维高超声速底部喷流干扰流场数值模拟与试验研究

通过数值模拟和风洞试验两种手段对来流马赫数M∞=4、喷流压比Pj/P∞=156.8、不同迎角下的三维高超声速底部喷流干扰流场进行了研究.研究结果表明超声速底部喷流干扰流场结构复杂,有、无喷流时底部流场有很大不同,对气动力系数影响显著;在大喷流压比情况下,喷流干扰使导弹纵向气动力系数下降、压心前移.最后,对数值模拟与风洞试验在结果上的差异进行了分析.     

超声速化学反应流场迎风格式数值模拟及并行计算

基于有限体积迎风格式对超声速燃烧流场进行了的数值模拟.由于超声速燃烧流场绕流的复杂性,要求对多组分Euler/N-S方程求解的数值模拟方法应具有较高的计算精度及效率.本文引用辅助点方法建立了具有空间二阶精度的van Leer迎风矢通量分裂格式,并应用于超声速燃烧流场绕流的数值模拟.化学反应为氢气/空气十反应模型,采用考虑了化学反应特征时间的当地时间步长显式Runge-Kutta时间推进格式.对钝头体模型爆轰现象、后向台阶氢气喷射及二维内外流超声速燃烧流场模型进行了区域分裂技术的并行计算.计算结果与参考文献作了对比,得到了满意的结果.     

高超声速咽式进气道起动特性研究

为了研究高超声速咽式进气道在非设计迎角以及低马赫数下的起动性能,利用流线追踪生成了设计马赫数Ma=7,具有8-7无粘基本流场(即俯仰平面内的斜激波由和自由来流呈8°夹角的斜压缩面产生;偏航平面内的斜激波由和自由来流呈7°夹角的斜压缩面产生)的咽式进气道,并对边界层修正前后的两种咽式进气道进行了数值模拟和高超声速风洞实验.实验观测和记录了各个来流条件下进气道模型唇口的激波系结构,测量了沿进气道模型上下壁面中心线从气流进口到出口的沿程静压分布.结果表明:迎角的增大和来流马赫数的减小都会对迸气道的起动性能造成不利的影响,通过对咽式进气道进行边界层修正,可以提高进气道的总压恢复系数,减小内收缩比,从而扩宽进气道起动的马赫数以及迎角范围,对进气道设计有着积极的作用.     

型面设计马赫数对马赫数分布可控高超声速内收缩进气道的影响

在研究型面设计马赫数Mai对马赫数分布可控基准流场在Ma∞=4.0~7.0的性能影响中发现,降低型面设计马赫数可获得更高的流量系数、高马赫数时的增压比以及更短的长度。基于型面设计马赫数Mai=5.5,6.0和6.5的基准流场分别设计了圆形进口的内收缩进气道,并在Ma∞=4.5~7.0时进行数值模拟。结果表明:基于低型面设计马赫数基准流场设计的进气道具有更好的流量捕获特性和较高的增压比,这与基准流场变化规律基本一致。型面设计马赫数对出口总压恢复系数影响较小。     

Roe格式在多组元燃烧流场数值模拟中的应用

对Roe通量分裂格式进行了详细讨论，并将其推广应用到基于非结构网格的高超声速多组元燃烧流场的数值模拟中。控制方程采用三维多组元Euler方程，重点针对Roe格式开展研究；考虑到燃烧流场求解过程中出现的刚性问题，对化学反应源项采用点隐式方法处理。对超声速条件下钝头体激波诱导燃烧流场进行了数值模拟，比较分析了3种经典化学动力学模型对流场结构的影响，获得与数值模拟、试验数据相符的结果。结果表明发展的方法可以用于多组元燃烧流场的数值分析，且化学动力学模型的选取直接影响激波诱导燃烧流场中诱导区厚度。     

纯净空气来流下的超声速燃烧实验装置及其初步实验结果

采用电阻加热的连续式实验设备,在燃烧室进口气流为高温纯净空气、马赫数Ma=2、总温Tt=1000K,总压Pt=0.8MPa条件下,进行了不同当量油气比的氢和乙烯燃料的超声速燃烧室直连式实验.采用从壁面垂直于主流喷射燃料和以氢作为先锋火焰,实现了乙烯燃料的可靠点火和稳定燃烧.实验测量了燃烧室的壁面压力、空气流量、燃料喷射压力、喷管进口总温等参数,并拍摄了燃烧室出口火焰.本文实验采用的电阻加热设备具有实验介质无污染、稳定运行时间长、工作性能稳定、成本低、操作简单等优点,其主要部件电阻加热器出口的最高温度可达600~1000K,对应的流量为1.5~0.73kg/s、加热器功率为750KW.     

单头部燃烧室流场PIV试验测量

航空发动机燃烧室内流场结构直接影响燃油雾化、油气掺混以及燃烧性能,本文采用粒子图像速度仪(Particle image velocimetry,PIV)对某单头部基准燃烧室内的冷态流场和燃烧流场分别进行了试验测量。在冷态流场试验中,研究了进口空气流量变化对燃烧室内的流场结构、回流区尺寸大小变化的影响规律;而燃烧流场试验测量分别研究了进口空气流量和油气比变化对燃烧流场结构的影响。试验结果表明:由于下壁面中间主燃孔进气射流的强烈影响与挤压,导致旋流器出口处横向截面上的旋转气流不是一个完整的旋流气流;燃烧流场与冷态流场相比,其流场结构基本相似,但中心回流区宽度稍变瘦,随着油气比的增大,中心回流区逐渐变瘦,宽度变窄;随着油气比的增加,轴向速度逐渐变大、回流负速度变大;燃烧流场测量中,在燃烧室头部较好地捕捉到喷嘴喷出的油雾锥上油珠的速度大小。     

不同飞行工况下双模态发动机流动及燃烧特性

为研究煤油燃料矩形截面双模态超燃冲压发动机在不同飞行工况下的流动及燃烧特征,在通过直连式试验验证计算方法的准确性后,对6个不同马赫数及当量比工况进行了三维定常数值模拟,得出了发动机壁面压力、一维质量平均马赫数沿流向的分布规律,分析了各工况下流场中波系结构、释热变化率等特征。研究结果表明:不同工况下发动机明显工作于两类不同的燃烧模态。当发动机处于预燃激波串前传至注油位以前的亚燃模态时,凹槽段波系相对较弱;随着激波串的前移,隔离段中形成明显的分离旋涡结构将燃料卷至上游,部分燃烧在注油位之前已完成;在燃烧室内,分离主要发生于凹槽内部,燃烧释热集中于第一凹槽头部。当发动机处于激波串未前传的超燃模态时,凹槽段波系相对更强,流动参数波动更大,燃烧在注油位以后进行,燃烧室内分离旋涡在流向跨度大,形成从第一凹槽前缘至第二凹槽处的连续流动分离;分离旋涡有助于燃烧向下游传播,因此释热沿流向分布更均匀、更分散。在过渡段诱导流动分离,促使燃烧室内形成大流向跨度的分离旋涡可能有助于燃烧向下游传播,实现分布式释热,避免释热过于集中导致激波串前传。     

蛇形进气道地面工作状态附面层抽吸试验研究

对一种蛇形进气道开展了地面工作状态下的抽吸试验研究,结果表明,在该状态下进气道出口截面的总压恢复系数较低、流场畸变较大。为此,本文采用附面层抽吸技术对其进行了地面抽吸状态下的流场控制试验研究。研究结果表明:(1)地面工作状态下,随着出口马赫数的增加,蛇形进气道出口截面的总压恢复系数不断下降,而稳态周向畸变指数、紊流度和综合畸变指数均上升,稳态径向畸变指数变化不大。本研究的蛇形进气道在出口马赫数为0.45时,总压恢复系数为0.90,综合畸变指数为13.85%,总压恢复较低,畸变较大,超出了一般航空发动机的承受范围。(2)与原型方案的地面抽吸试验结果相比,采用附面层抽吸技术后,进气道出口截面的总压恢复系数得到了提高。在出口马赫数为0.45,相对抽吸量为0.043时总压恢复系数提高了2.6%。     

三维平壁进气道的特性研究

介绍三维进气道和超声速燃烧冲压式喷气发动机 (SCRAMJET)模型一项实验研究的结果。Scramjet模型由进气道和燃烧室组成 ,此模型组件的设计是为了研究流动结构以获得进气道的特性 ,以及研究燃烧室和进气道、燃料 (氢及碳氢燃料 )点火和燃烧的相互影响。这些试验是在下吹式风洞(M =2 ～ 6,Reunit=( 8～ 54) × 1 0 6)和热射式风洞 (M =6及 7 2 ,Reunit =( 1 0 ～ 32 )× 1 0 6,Tt =1 50 0 ～ 2 50 0K)中进行的。     

矩形截面超燃发动机不同燃烧模态下的流场特征

为研究乙烯燃料矩形截面超燃冲压发动机不同燃烧模态下的流动特性，在直连式试验的基础上对冷流和不同当量比的4个状态进行了三维定常数值模拟，比较了试验和计算结果，选择了适用于本构型的模态判别准则，给出了流道内壁面压力、一维平均马赫数的沿程分布规律，分析了各状态下流场中波系结构、流动分离及燃烧的特征。研究结果表明:采用AHL3D对该发动机进行三维计算所得壁面压力与试验壁压吻合良好，试验与计算具有较好的一致性；未注油的冷态情况下流道内形成由多道斜激波与膨胀波组成的反射波系，壁面压力波动较大，波系分布主要受流道结构影响；纯超燃模态时，燃料喷射与主流相互作用使注油位处形成明显激波，压升起点固定在注油位之后，注油位波系对流场结构的影响较大，同时分离结构分布在整个凹槽内；双模态超燃时，流道内主导波系是激波诱导边界层分离形成的斜激波串结构，燃烧室内波系较弱，此时隔离段内激波串前缘后的角区出现分离，凹槽内分离区域减小；双模态亚燃时，随着逆压梯度激波串的前移，隔离段内角区的分离面积不断扩大，凹槽内分离区进一步缩小。发动机处于双模态超燃或双模态亚燃模态时，随着激波串结构的形成与前移，部分燃烧可能在隔离段内完成；而对于纯超燃模态，燃烧仅发生在凹槽与扩张段内，化学反应与高温区的分布相对更集中。     

定几何混压式轴对称超声速进气道的稳定性实验

对一种定几何混压式轴对称超声速进气道的稳定性进行了风洞实验研究,得到了不同状态下进气道的纹影图片、沿程静压分布以及进气道性能随出口马赫数的变化曲线,并进行了分析。结果表明:在超临界工作范围内,随着反压提高结尾激波系前移,进气道出口马赫数下降,总压恢复系数升高,而流量系数保持不变。当工作马赫数大于封口马赫数时,本研究中的混压式进气道亚临界工作状态存在稳定状态,而当马赫数低于封口马赫数时,进气道亚临界状态不稳定。     

Effects of Chemical Reactions in the Hypersonic Reacting Flow around Blunt Bodies

The effects of chemical reactions in the hypersonic reacting flow are investigated using an integrated algorithm considering simultaneously two different reaction mechanisms,i.e.,including the high temperature air nonequilibrium chemical reactions and the H_2-air combustion reactions. The program is validated by the air non-equilibrium flow at Mach number of 25.9 with the RAM C-II configuration and the shock-induced combustion flow at Mach number of 4.512 6 around a sphere,respectively. Furthermore,the mixed reacting flow with the Mach number of 10.0 with an opposing jet of hydrogen is numerically analyzed. The results show that the program is reliable,and the effects of chemical reactions engender in the decrease of peak temperature along characteristic lines,as well as on the surface. The production of water is augmented in the region with high ratio of oxygen to hydrogen and weakened in the area with low ratio of oxygen to hydrogen by the air chemical non-equilibrium effects.