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

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

爆轰波穿越惰性气团时的透射激波参数分析

爆轰波穿越惰性气团再进入可燃混合气体时,透射激波的参数主要受到三个基本过程的影响:头激波两次穿越惰性气团与可燃气体之间接触面导致的激波强度下降、放热反应中止、爆轰波后稀疏波对激波的衰减作用。放热反应中止是透射激波衰减的主要原因,不同种类的惰性气体对透射激波参数影响不大。     

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

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

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

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

基于侧向膨胀影响爆震波的自持机理

为了研究受侧向膨胀影响时爆震波的传播特性及自持机理,在实验段对比分析了当量比(0.70~2.25)和波前预混气高度(1,2,3cm)对爆震波自持传播能力的影响.实验表明:波前预混气高度越高、预混气活性越强,则爆震波抵御侧向膨胀影响的能力越强,速度亏损越小,自持传播能力越强.运用Fay流体扩张理论,Dabora和Murray速度亏损理论,并结合Zeldovich-von Neumann-D?ring(ZND)模型对受侧向膨胀影响的爆震波激波角、界面角和速度亏损进行理论预测,证明Dabora的理论预测与实验值吻合很好,且发现若要受侧向膨胀影响的爆震波自持传播,则其速度亏损的极限为7.0%~11.0%.      

狭缝内爆轰波传播模式的实验研究

为研究狭缝内爆轰波起爆及各种传播模式的特点,实验得到了不同初始压力下(p0=10~60k Pa)化学当量比的乙烯/氧气预混气体在高度1.0mm,宽度20mm狭缝内的爆轰性能。采用烟膜片记录爆轰波运行轨迹,高速摄影捕捉火焰面。结果表明:随着初始压力的降低,实验依次得到稳定爆轰模式、"结巴"爆轰模式、驰振爆轰模式和低速稳定传播模式,中间两者为不稳定模式。不稳定爆轰传播模式具有强烈的速度震荡,在一个周期内,烟膜板中周期性的出现过驱爆轰的胞格结构,随后"结巴"爆轰以单头形式与侧壁面发生碰撞并向前传播,驰振爆轰短暂熄灭。低速稳定传播速度为45%的C-J速度,烟膜板中无胞格结构出现。狭缝内爆轰波速度亏损值为10%~15%,高于大尺寸通道中的速度亏损。     

聚焦角度对空气中激光维持爆轰波的影响

吸气式激光推进的激光能量沉积过程主要发生在激光支持的爆轰波（LSD波）的形成与演化阶段,分析LSD波的传播规律和影响因素对揭示能量沉积机理起到基础性作用。利用纳秒四分幅高速相机和延长光路的纹影系统,拍摄了不同聚焦角度下CO₂激光经透镜聚焦击穿空气流场的纹影照片,结合一维理论公式分析了LSD波速度。结果表明：在本文的实验条件下,吸气式激光推进击穿空气经历了LSD波和激光支持的燃烧波（LSC波）两个阶段, LSD波阶段的持续时间约7.1 μs;3种聚焦角度下LSD波速度与理论公式预测的整体趋势比较吻合;聚焦角度越小,LSD波的初始速度越大,聚焦角度θ=4.29°时实验数据达到约14 km/s。     

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

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

气相爆轰波马赫反射非自相似性特征的实验

为了研究在CJ(Chapman-Jouguet)状态下气相爆轰波马赫反射的非自相似性特征,搭建了由驱动段、传播段以及观测段组成的矩形管道;利用纹影和烟膜实验方法分别对3种预混气发生马赫反射现象进行了实验研究.实验结果表明:对稳态或非稳态气体,测量了马赫反射三波点位置高度,验证了马赫反射三波点轨迹线为波动的曲线,即自相似性对爆轰波马赫反射失效;非自相似性的重要特征是三波点轨迹线在楔角初始位置时基本遵循无反应冲击波理论计算结果,随后逐渐偏离并向有反应冲击波转变,之后接近平行于反应冲击波理论计算并稳定地处于两个理论值之间.最后给出在初始压力为10kPa和楔角为30°的条件下,对于所研究的3种气体,其马赫反射三波点发生转变的位置高度分别约为0.8,1.05cm以及0.5cm,可见相同初始条件下非稳态气体的马赫反射发生转变比稳态气体提前.      

预混气爆轰马赫反射实验研究

为研究斜爆轰马赫反射的实验现象以及不稳定性对马赫反射的影响,在矩形爆轰管道中,对3种预混气的马赫反射现象进行了实验研究,取得烟膜记录。烟膜结果显示,三波点高度在楔角附近与无反应三激波理论类似,随后逐渐偏向反应三激波理论,三波点实际高度由反应实际放热量决定,三种气体偏转始末距离与结束偏转距离之比平均值分别为0.42,0.49,0.31,与三者活化能大小关系相符合。胞格结构随初始压力增大而渐趋规律,胞格尺寸也逐渐减小,由CJ区域到过驱区域胞格结构存在明显的转变。横波是影响斜爆轰稳定的重要因素,压力增大可抑制分子的无规则运动,使横波出现的位置更加集中,胞格更为规律。预混气爆轰不稳定度大的气体横波出现的时空位置更不易确定,胞格结构更不规律,胞格尺寸较小。     

A self-similar solution of a curved shock wave and its time-dependent force variation for a starting flat plate airfoil in supersonic flow

The problem of aeroelasticity and maneuvering of command surface and gust wing interaction involves a starting flow period which can be seen as the flow of an airfoil attaining suddenly an angle of attack. In the linear or nonlinear case, compressive Mach or shock waves are generated on the windward side and expansive Mach or rarefaction waves are generated on the leeward side. On each side, these waves are composed of an oblique steady state wave, a vertically-moving one-dimensional unsteady wave, and a secondary wave resulting from the interaction between the steady and unsteady ones. An analytical solution in the secondary wave has been obtained by Heaslet and Lomax in the linear case, and this linear solution has been borrowed to give an approximate solution by Bai and Wu for the nonlinear case. The structure of the secondary shock wave and the appearance of various force stages are two issues not yet considered in previous studies and has been studied in the present paper. A self-similar solution is obtained for the secondary shock wave, and the reason to have an initial force plateau as observed numerically is identified. Moreover, six theoretical characteristic time scales for pressure load variation are determined which explain the slope changes of the time-dependent force curve.     

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

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

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

在来流马赫数为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.     

二维凹面腔内暂冲式激波聚焦机理的实验

针对两级脉冲爆震发动机中第二级凹面腔内激波聚焦起爆能量和周期的研究,设计了二维凹面腔内暂冲式激波聚焦实验系统。通过分析粒子图像测速系统捕捉的照片和动态压力数据,获得了凹面腔内激波聚焦高压区压力特性和多循环流场演化特征。并且开展了A和B两种型面凹面腔的对比实验,进一步分析了多循环激波聚焦的形成机理。结果表明:入射激波第一次在凹面腔底部聚焦后形成的聚焦反射激波在喷口射流边界上发生反射,产生的反射激波再次在凹面腔底部聚焦,形成了第二次激波聚焦,如此往复便形成了周期性的激波聚焦,聚焦频率高达7~10kHz;并且,速度场和射流强度随入射激波碰撞而减弱,在排气过程中又增强,呈现周期性变化。      

连续旋转爆震波细致结构及自持机理

基于化学非平衡流解耦方法和详细化学反应机理,在0.1mm的网格尺度下对连续旋转爆震波内流场和爆震波平掠惰性气体界面的传播过程进行了数值模拟。详细分析了连续旋转爆震波的空间结构,由于在出口侧无几何约束,受一系列稀疏波的影响,旋转爆震波强度受到不同程度的削弱,形状也发生了相应的变化。爆震波平掠惰性气体界面时的流场分布与连续旋转爆震波结构极为类似,对比研究了波前可燃气体层高度对爆震波结构和传播过程的影响,当爆震波高度较小时,其抵御惰性气体侧稀疏波衰减的能力弱,波阵面严重变形,爆震波解耦。若爆震波要稳定自持传播,其必须要达到一定的临界高度来抵御稀疏波的影响。     

轴向入射激波反射聚焦的实验和数值模拟

为研究两级脉冲爆震发动机凹腔入口宽度对第2级凹腔中激波聚焦过程的影响,设计轴向入射激波聚焦的实验系统,对不同入口宽度下激波反射聚焦的过程进行了实验和数值模拟。采用纹影系统拍摄凹腔中的流场结构,并测量了凹腔顶点处的动态压力。采用保持强稳定性（SSP）的Runge Kutta格式、weighed essential non oscillation（WENO)格式和块结构网格自适应加密（SAMR）算法对激波反射、衍射和聚焦的过程进行了数值模拟,数值模拟结果和实验吻合较好。通过分析比较相同马赫数激波通过不同宽度凹腔的反射聚焦过程,发现较小宽度凹腔入口下激波反射和聚焦较强,顶点的压力峰值较高,激波聚焦诱导的射流较强。      

AUSM+格式的改进

为了提高AUSM 格式对流场计算中激波和剪切流的分辨率,减小数值振荡和激波后数值过冲问题,本文在原有AUSM 格式的基础上经过数值研究分析,对AUSM 格式进行了声速、特征马赫数的处理方式,对网格界面马赫数等进行了改进,并用改进后的AUSM 格式求解了含有脱体激波、斜激波、剪切流、膨胀扇、及湍流边界层等流动特征的复杂流动.经典型算例检验证实,改进后的AUSM 格式对激波和剪切流的分辨率得到了加强,求解的准确性和收敛性有一定的提高.改进后的AUSM 格式继承了原AUSM 格式优点,计算处理问题能力增强,可以有效地应用于亚声、跨声和高超声流动的工程实际问题中.     

基于前体激波的内转式进气道一体化设计

在腹部进气的乘波前体/内转式进气道的一体化设计中,为使进气道捕获截面和唇口型线的形状与飞行器前体激波较好匹配,提出一种基于前体激波形状的一体化设计方法。首先,计算乘波前体流场并提取前体激波形状;其次,将进气道捕获型线（ICC）投影在前体激波曲面上,得到可全流量捕获的进气道唇口型线（IFCC）;再次,给定进气道基本流场的中心体轴线位置,确定基本流场的入射激波形状;然后,给定基本流场的沿程压缩规律,应用特征线法确定进气道的基本流场;最后,将ICC顺来流方向投影至进气道入射激波曲面上,经流线追踪和黏性修正得到最终的进气道型面。数值模拟结果表明,对于典型飞行器前体,在设计马赫数为7.0的条件下,应用该方法得到的进气道流量捕获系数达0.976,隔离段出口截面的马赫数、压比和总压恢复系数分别为3.17、38.9和0.487。