基于多核计算机集群系统的网格分区策略

网格分区技术是提高并行计算效率的重要手段。基于图的分区技术(如Metis等)已经可以完全做到负载平衡,但是随着计算机集群系统的发展,节点间的通讯成为并行计算中的重要时间消耗部分。为此,提出一种适用于非结构网格的新分区策略(MC Partition method),减少了节点间的通讯规模,结果表明,这种分区方法可以较大的提高并行计算效率。     

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

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

连续旋转爆轰燃烧室增压特性的数值研究

为了分析爆轰燃烧室几何尺寸对增压特性的影响,利用二维欧拉方程数值研究了当量H2/Air在连续旋转爆轰燃烧室中的燃烧流场情况。研究表明:在准稳定状态下爆轰燃烧流场特征参数随时间呈周期性振荡,且其振幅保持不变;横波及局部爆炸波的存在是爆轰波后侧附近总压沿x轴呈现波动变化的主要原因;影响连续旋转爆轰燃烧室增压特性的直接因素为斜激波的高度及爆轰波强度,减少轴向尺寸或增加周向尺寸可增强其增压特性。连续旋转爆轰燃烧室的增压比高达2.52,与传统的燃烧室相比,大大提高了燃气的做功能力。     

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

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

斜爆轰波的波角和法向速度-曲率关系初探

为研究斜劈诱导斜爆轰波的波阵面弯曲效应，以期为斜爆轰的不稳定性及其演化规律提供新的见解，基于加权本质无振荡（WENO）格式空间离散和附加Runge-Kutta方法时间离散的求解器，针对不同的化学反应参数（释热量、放热速率和化学反应区参考长度）条件，开展斜爆轰波的数值计算研究。结果表明斜爆轰波沿波阵面的波角变化可分为3个区域：区域I，波角平滑减小；区域II，波角跃升后衰减；区域III，波角有规律振荡。波阵面法向速度-曲率关系在区域I呈现准垂直直线变化趋势，并伴随着爆轰波强度的不断衰减；在区域III则呈现出"D"形曲线，即由极曲线段、光滑水平变化段和拟线性变化段组成，为类胞格结构的周期性演变；区域II可认为是以上两个区域特征的耦合。不同的化学反应参数对斜爆轰波波阵面的弯曲效应影响存在较大差别。     

双楔面诱导的斜爆轰波阵面的磁流体控制

基于高超声速飞行器爆轰推进发展需要,数值研究了由双楔面诱导的斜爆轰波类型及其磁流体(MHD)控制的可行性,探讨了后楔倾角变化对双楔斜爆轰流场与MHD控制的影响.结果表明:对于前、后楔倾角分别为15°与20°的双楔面诱导的斜爆轰,Lorentz力可以将不同来流马赫数条件下的稳定斜爆轰波阵面恢复到设计位置,并可以使不稳定斜爆轰波趋于稳定,但是,无法将不稳定斜爆轰波恢复到设计位置.而当后楔倾角稍微增加时,若爆轰流场稳定,同样可对斜爆轰波阵面进行有效控制.但当后楔倾角大于25°时,稳定斜爆轰波流场失稳,此时MHD无法控制流场稳定性.但若增加气体反应速率,虽然斜爆轰波不稳定,MHD控制却可使其趋于稳定,但无法将爆轰波阵面恢复到设计位置.      

超声速弹丸诱导爆轰波的无网格数值模拟

为研究锥角对超声速弹丸诱导爆轰波形态的影响,发展了耦合有限速率化学反应模型的最小二乘显式无网格算法,其流体动力学采用含化学反应源项的多组分Euler方程建模,对流项和时间项分别采用多组分HLLC(Harten-Lax-van Leer-Contact)格式和四阶Runge-Kutta法计算。对尖劈诱导斜爆轰以及激波诱导燃烧流场进行了模拟,验证了算法的有效性。最后对等当量比甲烷/空气预混气体中,不同锥角弹丸诱导爆轰波流场进行了模拟,云图同实验阴影照片吻合较好,结果表明当锥角处于70°~100°时,易形成驻定斜爆轰波;锥角较小不利于可燃混合气体的点燃,仅能形成驻定斜激波;过大的锥角将导致爆轰波的脱体。     

旋转爆轰波中多波流动模式的数值研究

为了研究旋转爆轰燃烧室内复杂的波系结构,特别是爆轰波多波传播模式的影响因素,采用两步诱导-放热总包反应模型,对简化的二维旋转爆轰波进行数值模拟。在不考虑环形爆轰燃烧室曲率的情况下,研究了进口总温和周向尺寸对单波及多波流场结构的影响。结果表明,一定范围内进口总温的增加会使得爆轰波的波头数目增加,且双波流场结构中存在双波对撞和双波同向两种传播形式,双波同向传播和三波同向传播之间则稳定存在着两组双波对撞流场,且有对应的温度范围;对于双波对撞模式,持续减小周向尺寸能够使得流场结构转变成单波模式;周向尺寸的增加则会使得双波对撞向双波同向传播模式转变。     

高速飞行弹丸诱导的三维爆轰波数值模拟

采用三维非定常化学反应控制方程组对驻定在带攻角高速飞行弹丸上的斜爆轰波流场进行了数值模拟。对流项和化学反应项采用时间分裂法,分别用TVD格式和基元反应模型进行处理。计算结果表明三维数值模拟可以更准确的描述弹丸飞行攻角对爆轰波的形成及结构的影响,是进一步研究超驱爆轰推进技术的重要手段。     

曲面激波诱导斜爆轰的数值模拟

为了研究恰当比预混氢气-空气斜爆轰流场的波系结构和流动特征,基于带化学反应的Navier-Stokes方程,对弹头及楔穿越预混气体时诱导的斜爆轰进行了数值模拟。对流项的离散采用Steger-Warming格式,时间项采用二阶Runge-kutta方法。结果表明,对于弹头：（1）在亚爆轰条件下,能够模拟氢气-空气预混诱导爆轰流场的精细结构;（2）在超爆轰条件下,通过精细调整网格,能够很好地分辨强烈耦合的激波和燃烧波,且与Lehr实验吻合良好。对于楔结构：捕获到了清晰的三波点及其复杂精细的斜爆轰流场结构,预测的诱导燃烧距离、激波角和斜爆轰角与实验吻合良好;通过对流场波系结构变化过程的研究,获得了流场三波点随时间的演化过程。     

Investigation of hot jet on active control of oblique detonation waves

The hot jet injection is utilized to actively control the oblique detonation wave, such as initiating and stabilizing an oblique detonation wave at a desired position that is shorter than the length of induction zone, and adjust the height of the oblique detonation wave at the exit of combustor when the oblique detonation wave engine is working on off-design flight conditions. The fifth order Weighted Essentially Non-Oscillatory (WENO) scheme and a two-step reversible reaction mechanism of the stoichiometric H₂/Air are adopted in the simulations. With the help of hot jet injection, the transition from inert oblique shock wave to the oblique detonation wave immediately occurs near the position of hot jet injection, and consequently the length of combustor can be reduced. The angle of oblique detonation wave also decreases as the hot jet injection approaches the nose of wedge. Additionally, the height of the oblique detonation wave at the exit of combustor can be flexibly adjusted, and also depends on the injection position and the strength of the hot jet. If the velocity of the hot jet is too weak to directly trigger the overall oblique detonation wave at the position of injection, increasing the injection pressure will improve the strength of the hot jet and results in a successful transition.     

尖劈诱导的斜爆轰波的精细结构及其影响因素

为了揭示斜爆轰的精细结构及其影响因素,基于带化学反应的Euler方程,利用五阶WENO格式,对尖劈诱导的斜爆轰进行了数值模拟。结果表明,在劈面的压缩下,斜爆轰波阵面被分为三种结构,依次是类ZND结构、单三波点阵面结构和双三波点阵面结构。同时,根据三波点的轨迹绘制了胞格结构,其中单三波点的轨迹为平行直线,而双三波点的轨迹为倾斜的蜂窝状结构。讨论了来流马赫数和尖劈角度对劈面压缩性的影响,随着马赫数和劈角的增加,尖劈的压缩更厉害,爆轰波阵面也就越光滑。     

突跃型与平滑型斜爆震波起爆机制数值模拟

采用带有化学反应的Euler方程,对突跃型与平滑型2种形态的斜爆震波(ODW)在起爆机制和结构特点等方面的差异进行了研究.结果表明:斜激波(OSW)后是否存在亚声速区是判定斜爆震形态的依据.当波后火焰抬升斜激波面使得亚声速区存在时,横向激波与三波点结构就会出现并形成突跃型斜爆震.当波后所有区域皆为超声速区时,波后火焰燃烧无法影响到上游激波面,则会形成平滑型斜爆震.通过斜激波关系式给出了横向激波形成所需的临界斜激波角度,只有当斜激波角度大于此临界角度时才能形成突跃型斜爆震波,反之则形成平滑型斜爆震波.      

非预混条件下的旋转爆轰燃烧室双波头演化过程数值模拟

针对旋转爆轰燃烧室双波头演化过程中流场结构变化的问题,对非预混条件下的旋转爆轰燃烧室从起爆到形成稳定的双波头过程进行了数值模拟研究。研究结果表明,从起爆到形成稳定爆轰过程,燃烧室主要经历了起爆、爆轰波对撞和稳定爆轰三个阶段;在爆轰波对撞阶段,首次对撞是两个爆轰波间的对撞,由于对撞点处缺少新鲜混合气,从而在对撞结束后衰减为两个压力波。第二次对撞是两个压力波间的对撞,因为在第二次对撞点附近存在新鲜混合气来支撑爆轰波的持续传播,故对撞结束后产生了一个爆轰波和一个较弱的压力波;第二次对撞发生后,燃烧室内的压力波反射叠加并形成局部高压区,此高压区压缩气体使气体温度升高,高温气体引燃混合气后,最终发展成为第二个爆轰波;稳定阶段,两个爆轰波均能稳定自持传播,爆轰波峰面压力可达1.45MPa,波后温度为2500K,爆轰波速度稳定在1738m/s,产生的推力与比冲分别为79.76N和2312.15s;斜激波的存在使燃烧室出口平面流场产生了较大波动。      

斜爆震波起爆特性及其波系结构研究

本文通过对斜爆震波起爆特性及其波系结构的研究进展进行综合评述,分析其研究现状与存在的不足,指导未来斜爆震波的基础研究。首先总结了斜爆震波起爆和驻定特性研究进展,重点回顾了斜爆震波起爆准则和驻定条件,接着对斜爆震波的诱导区结构和波面胞格结构研究进行了总结,随后对斜劈构型对斜爆震波的影响研究现状进行了概述。综述表明,斜爆震波的起爆、驻定和波系结构还需持续开展研究,特别需要关注的是平面斜爆震波的起爆准则和驻定条件,在真实流动条件下更深入研究斜爆震波宏观结构和波面结构的形成机理,探索受限空间内斜爆震波起爆、反射和驻定特性。     

冲压加速器燃烧流场的数值模拟

本文利用ＥＮＮ格式,通过求解包含非平衡效应的完全Ｎ－Ｓ方程,数值模拟了超声速燃烧流动中很有兴趣的Ｈ２／Ａｉｒ斜爆轰波冲压加速器的绕流。     

A theoretical method for solving shock relations coupled with chemical equilibrium and its applications

In this study, a theoretical method is proposed to solve shock relations coupled with chemical equilibrium. Not only shock waves in dissociated flows but also detonation waves in combustive mixtures can be solved. The global iterative solving process is specially designed to mimic the physical and chemical process in reactive shock waves to ensure good stability and fast convergence in the proposed method. Within each global step, the single-variable equations of normal and oblique shock relatio...