CE/SE方法在非定常爆轰计算中的应用

建立了非定常爆轰波的理论模型。根据CE/SE方法的基本思想，推导出适合于求解二维N-S方程的长方形网格的CE/SE方法。对于包含刚性反应源项的守恒方程组，将其分为两部分，其中刚性反应源项单独求解。本文以爆轰管内装填C3H8/O2预混气体为例，数值模拟管内爆轰波的传播过程以及爆轰波脱离爆轰管后流场的图谱。计算结果显示，爆轰波在管内传播时，其波系极其复杂。当爆轰波离开出口后，迅速熄灭为非反应的激波。管外不仅有激波、膨胀波的传播，而且还存在激波与涡的相互作用。数值实践表明，该CE／SE方法可以有效捕获爆轰波与激波等强间断。     

氢氧气体爆轰光谱的研究

在以前用摄谱仪研究激波管内氢氧气体爆轰光谱的基础上，用光学多道分析仪收集了不同初始压强条件３００￣８００ｎｍ光谱范围氢氧气体爆轰的低分辨和高分辨光谱。     

激波聚焦火焰聚心起爆器的数值研究

针对短距离内的激波聚焦起爆问题,设计了一种激波聚焦火焰聚心起爆结构,通过数值模拟方法对该结构内部流场变化进行研究。结果表明:加速区是火焰失稳加速的重要区域;聚焦区能量汇聚造成局部爆炸是燃烧转爆轰加速的关键因素,且该区域内由于几何结构形成的卷吸涡能够增强聚焦并阻碍激波回传入加速区;稳燃区主要实现爆轰波从过驱状态到稳定爆轰状态的过渡区。确定中心轴线上温度、压力波动是由于聚焦腔聚焦作用造成的,而爆轰波锋面温度压力变化是由于三波点运动造成。      

基于OpenFOAM的三维H2/Air连续旋转爆轰流场数值模拟

为进一步研究旋转爆轰流场特征,基于开源计算流体动力学软件OpenFOAM,采用9组分19步的基元化学反应模型,对H2/Air连续旋转爆轰流场进行了三维数值模拟,得到了旋转爆轰波稳定传播时燃烧室内部流场的详细结构,研究了燃烧室头部激波的传播特性,分析了旋转爆轰燃烧室的压力增益性能。结果表明:旋转爆轰波后的第一道反射激波在由燃烧室外壁面向内壁面传播过程中反射激波的高度增加并在靠近内壁面附近与滑移线交汇形成局部高温高压区域;旋转爆轰波在外壁面位置处相位约落后于内壁面0.003rad～0.15rad,其相位差随燃烧室曲率差的增加而增大;燃烧室头部反射激波数目受到曲率差和进气总压的影响,燃烧室曲率差增大,反射激波数目减少,进气总压增大,反射激波数目增多;燃烧室压力增益保持在0.3以上,在进气总压一定的条件下,压力增益随着燃烧室曲率差的增大有增加的趋势。研究结果揭示了三维旋转爆轰流场的精细结构和燃烧室头部激波的传播规律。     

高焓激波风洞爆轰驱动技术研究

激波风洞爆轰驱动技术利用引爆可燃混合气体快速释放的化学能产生强激波,压缩激波管的试验气体,提供产生超高速流动所需的试验气源,是近十几年来发展成功的激波风洞强驱动方法.本文分布介绍了反向爆轰驱动、正向爆轰驱动和反向爆轰膨胀驱动模式,分析了应用这些驱动技术产生的高焓、高雷诺数、高超声速流动的气源特点,探讨了不同驱动模式影响激波风洞性能的关键因素.并重点介绍了反向爆轰膨胀驱动模式,分析了影响缝合条件的参数以及二次波现象.应用这些爆轰驱动技术,研制了能够产生总焓为1000K～8000K,具有较长试验时间的高品质超高速气流.为开展高超声速气动实验研究奠定了良好的基础.     

高温区碰并诱导气相爆轰二次起爆的数值研究

针对心形管内气相爆轰二次起爆过程进行了数值研究。采用二阶附加半隐的龙格-库塔法和5阶WENO格式求解二维欧拉方程,采用基元反应来描述爆轰化学反应过程。结果表明:气相爆轰在心形管中传播,由于壁面几何结构的作用,爆轰波熄灭后会发生二次起爆。该二次起爆由高温区碰并所诱导,且在管道不同位置诱导了两处二次起爆。爆轰流场中存在未反应气囊、射流、爆轰波马赫反射、激波与火焰作用及其所诱导的不稳定性等现象。计算得到的二次起爆过程胞格结构演变与实验基本一致。     

驻定斜爆轰波起爆条件与结构的数值模拟研究

为了探究驻定斜爆轰波的起爆条件以及来流马赫数和楔面角对斜爆轰波形态结构的共同影响,基于二维欧拉方程和单步化学反应模型,对驻定斜爆轰波进行数值模拟研究。结果发现,斜爆轰波的起爆和结构受到来流马赫数和楔面角的显著影响,来流马赫数越小,形成斜爆轰波的楔面角条件越苛刻。随着来流马赫数和楔面角增加,波后混合气体温度升高,爆轰反应变得更加剧烈,使得斜激波诱导区域变短,未反应气团减少,接触间断面更靠近楔面,下游出现的涡结构也减少。     

旋转爆轰流场的数值模拟

基于带化学反应的二维Euler方程,采用氢气-空气的9组分19步基元反应简化模型,对充有当量比的氢气-空气预混气和空气的环形旋转爆轰流场,从点火燃烧到发展成旋转爆轰的过程进行了数值模拟。根据数值结果分析了波后流场中爆轰产物受激波、高温和离心力等作用而挤向外壁,形成有利于充入燃料,实现持续稳定旋转爆轰的流场特征。还讨论流场中爆轰波、激波与间断面和内外壁面反射或折射,从而形成多个激波相交的波系特征。为认识和理解旋转爆轰流场,开展旋转爆轰的实验研究等具有现实的指导意义。     

高马赫数来流条件下斜激波与流向涡对相互作用

针对斜激波入射流向涡对形成的复杂流动干扰现象,在来流马赫数6条件下,利用后掠斜坡式涡流发生器生成流向涡对,采用数值模拟结合激波风洞实验,研究了斜激波与流向涡对相互作用引起的激波变形和旋涡演变。结果表明,斜激波与流向涡对相互作用的剧烈程度大致分为强、中、弱三种。在强相互作用下,激波面向上游凸起的程度最剧烈,呈现类“T”字型特征,其对称面处接近于一道正激波,波后出现回流区和局部流动滞止,使得流向涡对发生破碎;在中等相互作用下,对称面处存在局部正激波,但没有出现回流区;在弱相互作用下,不再出现局部正激波。进一步地分析发现,斜激波与流向涡对相互作用产生的流向涡对破碎现象,不符合经典的激波作用下单个流向涡破碎理论的预测。采用流向涡对对称面处的气流参数,对经典的流向涡破碎理论进行修正后,能够预测本文流向涡对在斜激波作用下是否发生破碎。本文为掌握激波作用下流向涡对破碎现象的规律,提供了重要参考。     

预混超声速气流斜激波诱导脱体爆轰研究

为了研究高静温预混超声速气流中爆轰的机理问题,利用自主设计的连续式加热器来产生马赫数为2.6,静温大于700 K的氢/空气预混超声速气流,并采用高速纹影技术来观测斜激波诱导脱体爆轰的直接起爆和发展的动态过程.研究表明,实验得到的起爆前斜激波和起爆后脱体爆轰波的角度与理论分析结果非常一致.并发现由激波诱导爆燃向激波诱导脱体爆轰的转变过程非常快;而且对于当前初温较高当量比较低的混合物,会出现爆轰波突然熄灭与重新起爆的现象.当前研究为爆轰燃烧在高超声速推进中的应用提供了重要参考.      

脉冲爆震发动机外流场数值模拟及实验

通过实验和数值方法,对脉冲爆震发动机(Pulse Detonation Engine,PDE)外流场进行了可视化研究。实验中采用YA-16高速阴影系统,拍摄了爆震外流场的时序阴影照片。对高温火团引发的氢气-氧气-氮气爆震过程和管内、外流场分布进行轴对称数值模拟,其中考虑了包含19个基元反应和9种组份的H2-O2-N2详细化学反应动力学机理。根据计算结果,由计算光学获得爆震外流场的时序计算阴影图。数值计算结果和实验结果基本一致,均形象地描述了管外流场的变化。根据实验和数值计算结果,详细地讨论了涡环和悬吊激波产生的动力学机理及其发展变化过程。     

某典型二元高超声速进气道内二次流分析

对某典型二元高超声速进气道三维流场进行了数值分析，将固壁面压力分布的计算结果与试验结果进行了比较，二者吻合较好，表明数值模拟方法正确，结果可信。分析了高超声速进气道流场的波系结构以及进气道(含隔离段)内二次流特征。给出了外压缩波系下的角涡的形成和发展规律，在侧壁与压缩面相交角区形成为逆时针方向角涡，角涡的形成、发展与近壁激波／侧壁附面层干扰以及压缩面有关。分析了隔离段中的激波的反射状况，提出了隔离段二次流的形成和发展规律以及影响因素。在隔离段进口附近存在一分离区。隔离段内二次涡的形成、发展和消失由上下壁面的压差，激波的移动方向(波后气流的横向流动)等决定。研究还表明对于高超声速进气道即使是二元进气道，也要关注二次涡的存在和发展。     

可压缩球涡和平面激波相互作用的数值研究

行进了可压缩球涡和平面激波相互作用的数值研究。一方面，用非定常Ｅｕｌｅｒ方程数值模拟了一种可压缩球涡的自身演化，同时，利用Ｒａｎｋｉｎｅ－Ｈｕｇｏｎｉｏｔ关系，在流场中嵌入运动激波，求解了激波－球涡相干的流动过程。     

脉冲爆轰发动机第一和第二次循环的数值模拟

对脉冲爆轰发动机两次循环工作过程进行轴对称数值模拟,考虑了包含19个基元反应和9种组份的H2-O2-N2详细化学反应动力学机理,通过改进的ISAT方法(动态自适应建表)减少反应化学的计算时间.数值结果显示了脉冲爆轰发动机充气,起爆,卸压和再充气的循环工作过程.根据数值计算结果,详细地讨论了两次循环管内、外流场的结构间的差异.两次循环中管内、外流场结构的差异又导致二者推力和冲量的不同.     

流向涡与斜激波相互作用的试验研究和数值模拟

为研究流向涡与斜激波相互作用在超声速燃烧中的应用,进行了由翼产生的流向涡与楔块产生的二维斜激波相互作用的燃烧室冷流试验研究.在不同激波强度下,纹影仪捕捉到了强、中、弱不同的涡/波作用现象.仿真与试验结果符合得较好.试验结果表明:对马赫数为2.3流场,在翼攻角10°时,能产生强流向涡,此工况下,锲角越大,涡/波作用越强.仿真结果表明:马赫数对涡/波作用影响较大,总压影响不明显,总温可影响亚声速回流区的尺寸.      

考虑顶隙的压气机单转子内三维跨音速粘性流场结构的数值模拟

借助于考虑叶片顶部间隙的压气机转子内跨音速湍流流场的三维数值解,得到了该转子流场中的复杂涡系与波系结构。给出了叶片表面及轮毂面上的壁面流型,分析并揭示了转子叶道内的马蹄涡、通道涡、角涡、尾涡等复杂涡系及其与激波的相互干涉;特别是,分析并揭示了转子叶道内顶隙射流的特点,得到了顶隙附近激波与泄漏流动的干涉形态及泄漏流动的发展,并着重给出了其与激波的干涉形态。     

超声速弱欠膨胀冲击射流流场结构

为了对超声速弱欠膨胀冲击射流的流场结构细节进行研究,使用大涡模拟方法对其进行了数值模拟。利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-kutta法。亚格子尺度模型采用的是修正Sm agorinsky涡粘性模型。通过与经典的冲击射流实验比较,证明了程序的可靠性。数值模拟得到了剪切层以及壁面射流中的涡结构和主射流中的激波结构,并且在此基础上对涡合并和板前激波和涡干扰现象进行了深入研究。发现涡合并现象主要出现在流场的上游,越往下游出现的几率越小;涡和板前激波的相互作用会引起激波位置和强度以及冲击平板上冲击区的压强的显著变化,同时也会导致涡的变形。     

二维喷管的初始流动

基于可压缩Navier-Stokes方程,采用大涡模拟方法与高精度混合WENO/TCD格式,对Ma=1.4的超声速平面射流初始流场进行了数值研究.数值结果清晰地描述了超声速平面射流初始流场的结构特征,包括主涡环与激波结构以及它们演变过程.因主涡环内存在涡导激波对,激波与涡相互作用加速射流剪切层失稳,使剪切层首次卷起形成小涡的位置出现在涡导激波对后,此与亚声速射流情况不同.小涡串卷起成后,相继与涡导激波对相互作用,使激波出现明显的变形并加速主涡环失稳.      

跨声速压气机转子叶尖非定常流场数值研究

采用时间精确求解方法对某高负荷跨声速轴流压气机转子在98%设计转速下的叶尖非定常流场进行了数值研究.结果显示,激波自身振荡不明显,叶尖区域流动的非定常性主要来源于叶尖泄漏涡的破碎及其与激波之间的相互干涉.对比设计状态与失速状态下叶尖泄漏涡的特点发现:在近失速点时,叶尖区域间歇性出现前缘溢流.分析表明:叶尖泄漏涡在激波后破碎是造成堵塞的主要原因,也是造成spike型失速初始扰动的原因.      

跨声速风扇转子间隙流动结构分析

为探索跨声速转子间隙流动结构,归纳间隙泄漏涡（TLV）和激波相互作用机理,以跨声速转子为研究对象,数值模拟不同间隙下流动特性.此外,着重探究TLV与激波的相互作用,斜激波受到间隙泄漏流的干扰,被削弱打断向上游凹曲;建立三维模型,加深对跨声速叶尖区域流动的全面认识.研究表明:随间隙增大,相同流量下,效率越小,压比也越小;TLV强度更大,偏离吸力面程度增大,沿周向和展向影响范围都越大.压力差提供泄漏流迁移动力,间隙提供泄漏流形成通道.选取h/c=1.0%间隙,随着间隙高度的增大,泄漏流周向运动趋势更明显且二次泄漏现象更剧烈;沿泄漏流方向无量纲流向涡量有少量减少,无量纲螺旋度较高,集中涡特性明显.