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

针对斜爆轰发动机简化模型,采用多组分基元化学反应流动数值模拟技术,数值分析了斜爆轰波系在受限空间内的宏观结构特征及其演变,并进一步分析了爆轰波系结构转变的临界条件.研究结果发现,随着楔面角度的增加,依次出现四种结构:激波诱导燃烧、斜爆轰波双规则反射、回流区马赫反射、楔面燃烧.对于稳定的波系结构,楔面压缩角同时存在下临界...     

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

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

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

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

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

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

激波诱导燃烧模拟中轴对称边界数值异常研究

对钝体激波诱导燃烧模拟中引起对称轴附近出现数值异常现象的几个影响因素进行了细致讨论.在飞行速度低于爆轰速度情形的模拟中,在一定条件下靠近对称轴的位置会出现明显的数值异常——燃烧阵面扭曲.分别从边界处理方法、计算格式、化学反应机制和计算网格等4个方面对该异常产生的原因进行了详细地讨论.研究表明4个因素都会导致数值异常产生,深入的分析表明计算格式是其中最根本的因素,而这种燃烧阵面扭曲的数值异常应该与粉刺问题(Carbuncle)类似.针对如何克服轴对称边界数值异常提出了一些合理的措施,并给出了验证算例的对比与分析.      

受限空间内激波与火焰作用的三维计算

激波与火焰的相互作用常发生在超声速燃烧和燃烧转爆轰过程中,为深入了解这一现象,采用带化学反应的三维Navier-Stokes方程,对平面入射激波及其反射激波与火焰的作用进行了计算研究,其中燃烧过程使用单步反应模型描述。研究结果显示:在受限空间内模拟激波与火焰作用,能更好地符合实验结果,从而体现出受限空间的三维效应;火焰在入射激波的作用下主要经Richtmyer-Meshkov不稳定而发生变形,此时火焰变形以物理作用为主,燃烧膨胀效果相对不明显;当反射激波与变形火焰再次作用后,火焰迅速膨胀变形,放热率维持在较高水平,此时化学反应过程起主要作用;在反射激波的作用下,变形火焰复杂三维涡结构的形成能强化热量与质量的输送,提高燃烧速率。     

激波绕射双方块加速诱导爆轰的数值研究

基于二维单步化学反应Euler方程,应用Roe/HLL混合格式以及自适应网格加密(AMR)技术,对方管中弱激波绕射双方块障碍物加速诱导爆轰现象进行了数值研究。计算结果表明,弱激波与双方块相互作用并发生绕射、碰撞可有效加速诱导爆轰,且在管道阻塞比为45%左右时可达到最佳效果,与相关实验结果相符。另外,研究了方块个数及排列方式对激波绕射加速诱导爆轰的影响,得到了最佳方块数量与方块排列方式。     

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

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

射流涡发生器对激波边界层作用诱导的流体分离控制大涡模拟研究

为了研究射流涡发生器对激波边界层作用所诱导的流动分离控制机理及其流场特性,基于大涡模拟(Large eddy simulation)方法和高阶TCD/WENO混合格式,对来流马赫数Ma=2.5情况下,平板上射流涡发生器对激波与边界层相互作用所诱导流场进行了数值模拟。结果表明,射流涡发生器对激波边界层的流体分离有一定的抑制作用,与无控制情况相比,射流作用下进出口总压恢复系数由85.9%提高到94.6%。射流尾涡主要集中于一环状区域内,在该区域内,入射激波与马蹄涡、桶形激波上方的涡管以及剪切涡相互作用,导致整体尾流被激波往下压缩。同时在激波的压缩下,各涡之间相互缠绕、挤压合并,形成多个流向小涡结构,将边界层外的高速流体卷入边界层内,从而增加边界层底层能量,达到抑制流动分离的目的。     

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

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

气体动理论BGK格式的网格自适应方法

为了提高气体动理论BGK（Bhatnagr-Gross-Krook）格式在超声速流动问题计算时激波捕捉的准确性与计算效率,提出了一种适用于气体动理论BGK格式的网格自适应加密方法。该方法采用基于四边形的链表技术来描述网格的拓扑结构,在物理量重构过程中,使用了在四边形网格中表现优异的van Leer限制器,以保证粗细网格过渡处物理量重构的精度。用跨声速翼型绕流（马赫数Ma=0.85）、超声速前台阶流（Ma=3）和高超声速圆柱绕流（Ma=8.03）等多个典型算例验证了BGK自适应网格方法。计算结果表明,自适应网格BGK方法在保证数值精度的前提下,可大幅度提高计算效率。这为该方法用于高效地解决复杂问题提供了一种选择。     

基于自适应网格加密的超声速可燃气热射流起爆详细反应数值模拟

采用块结构自适应网格加密开源程序AMROC,在高性能计算集群中,进行氢气、氧气、氮气详细反应机理的二维超声速热射流爆震起爆自适应网格精细数值模拟,研究在热射流的持续喷射作用下超声速可燃气热射流爆震起爆,以及形成的精细的爆震胞格结构.结果表明:超声速可燃气流中热射流类似气动斜劈,形成局部激波诱导燃烧.超声速可燃气中持续的热射流喷射会形成过驱爆震,并导致不规则爆震胞格的生成.热射流的扰动压缩作用对过驱爆震的形成以及不规则爆震胞格的产起着关键的作用.热射流的扰动压缩波以当地声速通过爆震波后的亚声速区域作用于爆震波,使爆震波处于持续过驱状态,并形成不规则的爆震胞格结构.     

化学动力学模型对H2/AIR超燃模拟的影响

为了研究超燃流动中不同化学反应模型的影响，提出流动项与化学反应生成源项解耦处理的化学非平衡流动计算方法，从薄层近似N－S方程出发，采用ENO差分格式，数值模拟了超声速H2／Air 预混合气体中钝体激波诱导的振荡燃烧流场，针对飞行速度大于和接近爆轰速度两种不同流动状态，比较Evans-Schexnayder提出的7组分8反应机理（ES模型），Jachimowski提出的13组分32反应机理（J模型）及其简化为9组分19反应机理（J－1模型）对流场参数的影响，计算表明反应机理引起压力场变化较小，对温度等参数影响较大，对于飞行马赫数低于7的问题，取不包括N原子反应的Jachimowski的9组分19反应机理较合适。     

压气机线性振荡叶栅气弹试验研究（一）：非定常气动响应

 为了加强对叶轮机三维气动弹性机理的理解,为数值模拟方法提供试验验证数据,开展了压气机线性叶栅气动弹性试验研究。本研究分为两部分,本文为第1部分,研究压气机线性振荡叶栅的气弹稳定性,第2部分重点辨别叶尖间隙对此振荡叶栅气弹稳定性的影响。组建了线性叶栅试验台,其试验段由7个真实可控扩散叶型组成,中间叶片由凸轮机构驱动模拟三维一阶弯曲振动。使用外置压力传感器在6个叶高处测量叶片表面的定常压力分布和3个折合频率下叶片表面非定常压力分布。由管传递函数修正连接管造成的非定常信号失真。运用影响系数法构造协调叶栅任意叶片间相位角下的非定常气动响应结果。试验结果表明：非定常气动响应具有完全的三维特征;折合频率和叶片间相位角对叶栅的气动弹性稳定性具有决定性影响。     

High-order compact finite volume methods on unstructured grids with adaptive mesh refinement for solving inviscid and viscous flows

In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The high-order variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4est based on the forest of octrees is adopted. Several two- and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.     

NUMERICAL INVESTIGATION OF UNSTEADY DROPLET COMBUSTION PROCESS IN HOT QUIESCENT AIR ENVIRONMENTS (HQAE)

ＮＵＭＥＲＩＣＡＬＩＮＶＥＳＴＩＧＡＴＩＯＮＯＦＵＮＳＴＥＡＤＹＤＲＯＰＬＥＴＣＯＭＢＵＳＴＩＯＮＰＲＯＣＥＳＳＩＮＨＯＴＱＵＩＥＳＣＥＮＴＡＩＲＥＮＶＩＲＯＮＭＥＮＴＳ（ＨＱＡＥ）ＺｈｕＧｕａｎｇｓｈｅｎｇ，ＬｉｎＪｕｎｙｉ，ＺｈａｎｇＹｕｅｌｉｎ...     

不同直径圆球诱导燃烧的振荡机制与频率特性

为了分析不同直径圆球诱导振荡燃烧的规律,并揭示圆球大小在振荡燃烧现象中所发挥的深层次作用,采用二维轴对称Euler方程和基元反应模型,对不同直径的圆球在H2/air预混气体中诱导振荡燃烧的现象开展数值模拟研究.研究发现,振荡频率并不是简单地随直径增大而逐渐从高频向低频连续过渡,而是存在两次突变,形成了超高频、高频以及低...     

模拟存在运动激波流场的无网格自适应算法

将无网格自适应算法发展用于求解存在运动激波的非定常流动问题。为了降低初始布点要求,同时提高对运动激波的分辨率,提出了一种结合点云加密和点云粗化的动态点云技术。采用该技术,在非定常计算的每一个时间步,首先通过基于压强梯度的流场探针识别出激波所在位置,然后对激波附近点云进行加密,从而提高对激波的分辨率,同时对远离激波的点云进行粗化,尽可能减少冗余节点,提高计算效率。采用所提出的无网格自适应算法求解了N-S方程,首先对一维运动激波问题进行了数值模拟,通过将数值模拟结果与精确解比较,验证了所提算法求解存在运动激波的非定常流动的可行性,在此基础上,将方法推广用于二维非定常流动问题的求解,分别对瞬时激波碰撞圆柱以及双弧翼型激波振荡流动进行了数值模拟,计算结果表明,方法可以显著提高对激波的分辨率,同时相对于直接采用较密布点的非自适应算法,自适应算法在计算效率上具有优势。     

Large eddy simulation of hydrogen/air scramjet combustion using tabulated thermo-chemistry approach

Large eddy simulations(LES) have been performed to investigate the flow and combustion fields in the scramjet of the German Aerospace Center(DLR). Turbulent combustion is modeled by the tabulated thermo-chemistry approach in combination with the presumed probability density function(PDF). A b-function is used to model the distribution of the mixture fraction, while two different PDFs, d-function(Model I) and b-function(Model II), are applied to model the reaction progress. Temperature is obtained by solving filtered energy transport equation and the reaction rate of the progress variable is rescaled by pressure to consider the effects of compressibility. The adaptive mesh refinement(AMR) technique is used to properly capture shock waves,boundary layers, shear layers and flame structures. Statistical results of temperature and velocity predicted by Model II show better accuracy than that predicted by Model I. The results of scatter points and mixture fraction-conditional variables indicate the significant differences between Model I and Model II. It is concluded that second moment information in the presumed PDF of the reaction progress is very important in the simulation of supersonic combustion. It is also found that an unstable flame with extinction and ignition develops in the shear layers of bluff body and a fuelrich partially premixed flame stabilizes in the central recirculation bubble.