聚心火焰与诱导激波相互作用及爆燃转爆轰过程

基于带化学反应的二维轴对称Eu ler方程,利用带限制函数的波传播算法,在两端均敞开和一端封闭一端敞开的两种边界情况下,对圆柱型燃烧室中的氢气-空气预混气的聚心燃烧进行了数值模拟,探讨了聚心火焰引发爆轰过程中的波系结构以及火焰的变形。数值结果表明,在两端均敞开的情况下,燃烧诱导的激波在轴心、火焰和壁面之间来回反射,使火焰失稳,加速了火焰的燃烧,但没有转变成爆轰。而在一端封闭的情况下,封闭端反射产生的激波,不断穿越火焰,使火焰严重失稳,加剧了燃烧速度,最终导致爆轰的形成。稳定爆轰阵面的参数与CJ理论计算的结果进行了对比,两者符合的较好。同时,火焰在与激波的作用过程中,形状扭曲变形,对于两端敞开的情况呈对称的扁平头部蘑菇云,而对于一端封闭的情况,则呈封闭端小敞口端大的扁平头部蘑菇云。     

基于成像的RP-3煤油在激波管内自点火过程研究

为深入直观了解RP-3煤油在激波管内的自点火及燃烧过程,利用高分辨率纹影系统和高速摄像机对反射激波后点火燃烧过程进行拍摄,获取了丰富的光学图像。纹影图像显示得益于低压段内高比例氩气作为稀释气体,入射及反射激波均无分叉,为正激波结构,点火燃烧后对密度造成较大影响,流场中出现明显褶皱。针对燃烧过程的高速摄影图像显示激波后燃烧区域基本分为两部分,反射激波后小粒径雾化煤油完全燃烧呈现明亮的蓝色火焰,由于氧气的大量消耗随后将造成煤油液滴间的燃烧出现竞争关系,无法充分燃烧,出现分散的黄色火焰区,后期部分液滴由于贫氧及温度下降造成燃烧不充分而析出碳粒。     

进口斜激波、膨胀波干扰下等直隔离段内的激波串特性

 超燃冲压发动机的隔离段在实际工作中会受到进气道唇罩激波及肩部膨胀扇的显著干扰,本文针对这一特定问题进行了专门研究。提出了唇罩入射激波及肩部膨胀扇的模拟方法,并利用德国Achen的风洞试验对其进行了检验,而后以此研究了入射激波及肩部膨胀扇干扰下隔离段内激波串的基本形态,并分析了出口反压和激波入射位置的影响。仿真结果表明：当激波串在隔离段内不断前移时,受唇罩入射激波及其反射激波的干扰,其高速核心区交替地偏向上下壁面;与无激波入射的情况相比,此时激波串的耐反压能力显著降低,且入射点位置越高,降低幅度越大,管道内的沿程静压分布规律与Waltrup经验公式偏离程度也越来越大。该文结果可为进气道/隔离段的一体化设计提供依据。     

准定常强激波马赫反射波形的数值研究

应用DCD频散控制激波捕捉格式,求解二维、多组分、带有化学反应的Euler方程组,数值模拟了准定常强激波的马赫反射问题.研究结果表明:与经典马赫反射理论相比,在强激波条件下,激波诱导的气体分子振动激发和化学反应使马赫反射的三波点轨迹角变小、马赫杆高度变低、楔顶附体激波倾角变小;马赫杆的相对突出量随入射激波马赫数和楔角的增大而增大,而气体分子的振动、离解等真实气体效应能进一步加剧马赫杆的向前突出.     

径向入射激波聚焦实验和数值模拟研究

为得到不同入射激波马赫数对凹腔内径向入射激波聚焦过程的影响,对马赫数为1.25和2.18的径向入射激波在实验段内反射聚焦的过程进行了实验和数值模拟研究。用高速CCD拍摄了凹腔中气流流场的纹影图片,并采用保持强稳定性的龙格—库塔格式、加权基本无振荡算法和块结构动态自适应网格加密对激波反射聚焦的过程进行了数值模拟,实验和数值模拟结果吻合较好。复杂的激波形状及其之间的相互作用(如激波反射,激波聚焦和不同激波的相互作用)在实验和数值模拟中都能很好地观测到。通过比较不同马赫数入射激波的反射聚焦过程,发现较高的马赫数下聚焦后能产生射流,而较低的马赫数下却没有。同时发现,入射激波马赫数能影响二次激波出现的位置和强度,进而影响聚焦过程流场的紊乱程度。     

斜爆轰波系在受限空间内的演变及其临界条件的数值研究

作为一种面向未来高马赫数吸气推进技术的新概念发动机,斜爆轰发动机的运行依赖其稳定的宏观波系结构。本文针对斜爆轰发动机简化模型,采用多组分基元化学反应流动数值模拟技术,数值分析斜爆轰波系在受限空间内的宏观结构特征及其演变,并进一步分析爆轰波系结构转变的临界条件。研究结果发现,随着楔面角度的增加,依次出现四种结构：激波诱导燃烧,斜爆轰波双规则反射,回流区马赫反射,楔面燃烧。对于稳定的波系结构,楔面压缩角同时存在下临界、亚临界以及超临界三种极限条件,在波系演变过程还伴随激波规则反射和马赫反射的转变。     

一种形状可控的激波增强管道型线设计新方法

激波管所产生的非定常运动激波,若强度和形状能够按照一定的设计要求进行可控条件下的调节,将可望为燃料点火燃烧试验等提供具有独到优势的研究手段。基于激波动力学理论,针对激波管中所产生的平面运动激波,通过设计特定的上下壁面收缩型线,使初始平面运动激波,经收缩段（包括光滑凹形曲线段、斜直线段和光滑凸形曲线段）的变形和强度增加,再以平面波面形状进入较小截面直管段的连续转变过渡,得到了强度增加的平面激波。进一步对所设计的典型型线分别采用数值计算和试验的方法,考核分析激波运动过程中的形状变化,验证了理论方法的可靠性。在此基础上,分析了型线设计的关键参数对激波增强幅度的影响,结果表明,相对于传统激波管方法,本文中所提出的收缩截面方法能更显著地增加平面激波强度;另外,还考察了初始入射激波马赫数对壁面型线和运动激波波面形状的影响,结果表明,对于较强的初始入射激波来说,壁面型线对入射激波强度依赖较小,也就是说,当实际入射激波马赫数即使稍偏离设计状态时,仍然能得到近乎完美的平面形状增强激波。     

定常流场中二维激波反射迟滞现象研究进展

为深入理解超声速飞行器和动力系统流场中存在的激波反射迟滞现象,对近年来有关定常激波反射迟滞文献资料进行了调研和总结。回顾了定常激波反射迟滞现象的提出与证实过程,介绍了改变入射激波形状、改变下游波系结构等可引起激波反射出现迟滞的实现方式,并对粘性效应、三维效应和动态效应等影响激波反射迟滞现象的因素进行了探讨,对激波反射迟滞现象的研究热点和难点问题进行了分析与总结。     

燃料射流隔板对受限超声速反应混合层流场和燃烧特性影响

为了研究燃料射流隔板壁面扰动对于受限超声速反应混合层流场和燃烧特性的影响,基于抽象物理模型激波和膨胀波理论推演,得到了隔板诱发壁面扰动下的受限超声速反应混合层流场结构,并采用数值模拟方法进行了验证。数值结果表明,隔板诱发壁面扰动下的受限超声速反应混合层流场主要由冷态主导,且相对无隔板扰动多了回流区、激波、膨胀波、波与反应混合层相互作用等复杂现象。在此基础上,对于不同隔板厚度计算的结果表明,随着隔板厚度的增加,隔板下游的回流区增大,上下膨胀产生的压力不平衡加剧,反应混合层会产生偏斜,回流区厚度和偏斜距离与隔板厚度成正比。此外,第一道反射激波存在降低点火延迟的作用,且存在一个隔板厚度阈值,阈值以下随着厚度增加点火延迟随之降低。同时,后续的多道反射激波导致混合层发展的局部起伏和局部燃烧增强。     

脉冲爆震燃烧室内缓燃向爆震转捩数值模拟

为了研究脉冲爆震发动机燃烧室内火焰加速及缓燃向爆震转捩过程，利用7组分8反应的氢气详细化学反应机理进行 了2维数值模拟。结果表明：在火焰传播的初始阶段，障碍物、旋涡与火焰的相互作用是主导火焰加速的主要因素。在燃烧区域 产生的膨胀波向前传播并诱导未燃混气流动；障碍物后的回流区及障碍物顶端脱落的旋涡使火焰面拉伸、卷曲，增加火焰面面积， 同时提高燃烧放热强度，使火焰传播速度加快。在火焰发展的后期阶段，激波与障碍物的相互作用会使脱落在已燃混气中的未燃 混气微团发生局部爆炸，产生新的激波进一步推动火焰加速，缩短火焰锋面与前导激波的距离。由于交错型障碍物增加了火焰传 播的距离，其缓燃向爆震的转捩时间比采用对称型障碍物时的更长。     

一种高阶精度大涡模拟技术用于激波/火焰相互作用

基于Smagorinsky亚格子模型和增厚火焰技术,开发了一种高阶精度的反应流大涡模拟求解器,运用于数值研究边界层对激波与火焰相互作用的影响。该求解器的核心在于加入了基于超黏性的激波捕捉技术和时空三阶精度的两步Taylor-Galerkin紧致（TTGC）有限元格式,并通过对一维Shu-Osher问题和二维激波/气泡相互作用问题的计算,验证了求解器对激波、接触间断和湍流脉动等流动细节的捕捉精度,计算结果与实验数据吻合良好。通过对激波管内激波、火焰与边界层相互作用问题的数值模拟,发现由于激波与边界层的相互作用会产生不稳定的激波分叉现象,激波三分叉点传播速度的发展经历了水平匀速运动、小斜率线性增长和大斜率迅速增长3个阶段,由此揭示了激波分叉促进火焰加速的机理。当火焰面传入激波分叉区后,流场不均匀的回流区起到了稳定火焰的作用,一方面分叉结构内火焰面能够为激波的运动持续地供应热量,另一方面局部超声速区域为火焰的快速传播提供动力,使其能紧跟激波。通过对比相同条件下甲烷和乙烯燃烧的数值结果,发现两者爆震点触发的位置都出现在马赫杆后面,热量释放率的变化趋势也大致相同,但乙烯出现爆震的时间比甲烷早。     

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

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

不同半径圆柱诱导CH4/空气预混燃烧的计算研究

本文研究不同半径圆柱诱导CH4/空气预混燃烧流场。采用保自由流5阶WENO格式求解贴体坐标变换后的多组分Euler方程,用基元反应模型描述CH4/空气燃烧。不同于标准WENO格式通量构造方法,该WENO格式数值通量由方程的解进行WENO插值得到,在曲线坐标系下具有保自由流性质。首先给出了保自由流WENO格式精度和保自由流的数值验证,然后计算圆柱诱导CH4/空气预混气燃烧流场,并考察不同半径圆柱的影响,给出燃烧流场压力与温度等值线和云图、压力和温度沿过驻点线分布。结果表明:在考核计算结果网格无关性基础上,该WENO格式可准确地捕捉激波和火焰阵面形状、激波和火焰面驻点距离,得到的计算结果和文献结果相符。增大圆柱半径,激波和火焰面被推向来流方向,激波和火焰面之间距离也减小。和TVD格式相比,5阶WENO格式采用四分之一的网格数可得到近似相同的计算结果。     

动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

吸气式飞行器进气道唇口三维激波/激波干扰

对吸气式飞行器进气道唇口处三维曲面激波/弓形激波干扰流场进行数值模拟,利用典型三维气动干扰试验对采用的数值计算方法进行验证。利用拼接网格技术及逆距离加权插值方法获得入口处流场的非守恒变量,作为激波干扰研究的入口边界条件。数值模拟表明,唇口处激波干扰流动的三维效应十分显著,曲面激波与弓形激波产生斜交,尽管唇口前缘半径很小,但Edney提出的6类激波干扰类型可能沿唇口展向方向同时存在；第Ⅲ和Ⅳ类激波/激波干扰的诱导使得唇口热流分布异常严酷；激波相交处形成斜向“伤疤”状局部高热流条带,峰值热流达到参考热流的4~6倍,可能引起唇口结构的局部烧蚀或破坏。      

超声速气流中的斜激波诱导自点火

在来流马赫数为2.5的条件下对氢气和乙烯燃料的斜激波诱导自点火过程进行了试验研究.采用纹影和高速摄影获得了自点火初始火核的形成和发展过程.试验考察了不同的来流总温、喷注条件、凹腔长度、斜激波强度等因素对自点火过程和点火边界的影响.研究结果表明:入射斜激波能够提高气流的局部参数起到强化自点火的作用.凹腔内形成的低速回流区有助于初始火核形成后的火核逆流传播,并起到防止初始火核因气流振荡而吹脱的稳定源作用.增加喷注压力、增加凹腔长度、提高来流总温有助于自点火的发生.      

斜向环形激波聚焦诱导直接起爆的大涡模拟

基于大涡模拟(LES)数值方法,利用基元反应模型对斜向环形激波聚焦诱爆氢气-空气可燃混合气直接起爆进行了详细的数值模拟.研究结果表明:在进气马赫数为2.5的条件下,斜向环形激波聚焦产生的高能区可以在可燃气体中实现直接起爆形成爆震波,并且燃料能够高效完全反应;而垂直环形激波则没能诱导出爆震燃烧过程.可见,斜向环形入射比垂直环形入射更加有利于燃料混合气聚焦诱导起爆成爆震燃烧.      

Reflection of a moving shock wave over an oblique shock wave

The reflection of a moving shock wave over a wedge immersed in a still gas and the reflection of a wedge induced steady shock wave over symmetrical and asymmetrical reflecting surfaces have received intensive considerations since more than 70 years ago. Here we consider a different shock reflection problem—reflection of a moving shock wave over an initially steady oblique shock wave induced by a wedge immersed in supersonic flow. For the flow condition we considered, five moving triple points, with each connecting an incident shock wave, a reflected shock wave and a Mach stem, are identified. By using the reference frame co-moving with each triple point, the type of each shock wave of this triple point is clarified. The present study is significant in that it treats a new shock reflection problem leading to a new shock reflection configuration and showing potential applications in supersonic flow with unsteady shock interaction.