气相爆轰波绕射流场显示研究

采用基于红宝石激光器（波长694.3 nm）的纹影系统,对气相爆轰波绕射进行了初步的流场显示研究。采用单色激光和合适半带宽（15 nm）的滤光片,有效地消除了爆轰波自发光对流场显示的影响。合理设置激光器同步控制系统的触发延时,得到了序列的爆轰波阵面纹影照片。结果表明:图像清晰地显示了爆轰波阵面的诱导激波、横波及化学反应区。当爆轰波在左尖点处绕射,受稀疏波作用,诱导激波与化学反应区明显分离,导致爆轰波衰减为爆燃。分离的诱导激波和折皱的化学反应区在纹影图上清晰可见。诱导激波在垂直支管右壁面反射,诱导二次起爆。畸变爆轰波在水平和垂直支管中均发生马赫反射。提高初压,爆轰波受分叉口几何属性的影响减小,畸变爆轰波在水平和垂直支管下游较易恢复为自持爆轰波。     

可燃气体中激波聚焦的点火特性

数值模拟了二维平面激波从抛物面上反射在可燃气体中聚焦的过程，研究了形 成爆轰波的点火特性. 对理想化学当量比氢气/空气混合气体，在初始压强20kPa的条件下， 马赫数2.6-2.8的激波聚焦能产生两个点火区：第1个点火区是反射激波会聚引起的，第 2个点火区是由入射激波在抛物面上发生马赫反射引起的. 这种条件下流场中会出现爆燃转 爆轰，起爆点分别分布在管道壁面、抛物反射面和第2点火区附近. 起爆机理分别为激波管 道壁面反射、点火诱导激波的抛物面反射和点火诱导的激波与第2点火区产生的爆燃波的相 互作用. 不同的点火和起爆过程导致了不同的流场波系结构，同时影响了爆轰波传播的波动 力学过程.     

气相爆轰波在半圆形弯管中传播现象的实验研究

对气相(2H2/O2/Ar系统)爆轰波在半圆形弯管中的传播现象进行实验研究。用烟迹膜记录了弯管中爆轰波的胞格结构,采用压电传感器测量了沿弯管内外母线指定点的压力时间曲线,得到了爆轰波沿弯管内、外母线的平均速度和胞格尺寸的变化。结果表明:当平面爆轰波进入弯管后,受壁面的几何形状作用,诱导激波阵面发生弯曲。沿诱导激波阵面,自内母线到外母线方向,激波强度逐渐增大。同时,爆轰波后的化学反应区也受到影响,胞格尺寸发生较明显的变化。在本文条件下,当初压p08.00kPa,受扰动的爆轰波在弯管出口下游仍恢复为强度不变的稳定爆轰。胞格记录的三波点迹线表明:受扰动的爆轰波在出口段发生了马赫反射。实验结果还表明:当p0降至5.33kPa,平面稳定爆轰波经过半圆形弯管后,其强度发生衰减并直至出现熄灭。     

气相爆轰在T形管中传播新现象的实验研究

对2H2/O2/Ar系统爆轰波在T形管(截面为40mm×40mm)中传播现象进行了实验研究.用烟迹片记录了T形管中爆轰波的胞格结构,用压电传感器记录了分叉口附近指定点压力时间曲线,得到了爆轰波在分叉口附近的平均速度和胞格图案演变.结果表明:初压P0≥2.67kPa,在水平和垂直支管下游区域(距离分叉口约3.5—6倍方管截面边长),分叉口影响消失,爆轰波恢复稳定,且强度基本保持不变.在分叉口绕射过程中,爆轰波在膨胀区中衰减,诱导激波阵面弯曲.两个支管中发生马赫反射,三波点迹线清晰可见.该传播特性是爆轰波的诱导激波和横波共同作用的结果.分叉口附近的胞格结构先消失再恢复,在无胞格和平衡胞格之间的区域存在细密胞格的过渡区,表征了在诱导激波与化学反应阵面分离后的区域中出现二次点火.P0=2.00kPa,水平支管中稳定自持爆轰能重建,垂直支管中爆轰熄灭.P0<2.00kPa,分叉口上游已不能形成稳定爆轰.还对胞格结构中的几个特征参数进行了测量,并初步分析了P0对这些参数的影响.     

爆轰波与激波对撞的实验研究

对乙炔氧气混合气体中爆轰波与激波的正面对撞现象的实验研究是以高速摄影获取两波对撞的x-t纹影图,以烟迹板记录对撞中的爆轰胞格图案,并基于激波理论和经典CJ爆轰理论求解了两波对撞的稳态解并探寻其规律. 研究发现透射波系包括一道激波和爆轰波,以及紧随爆轰波后的稀疏波区,这一结果对应于一维理论分析中的CJ解. 透射波系基本不受初始压强影响;初始温度也只成比例地改变流场整体速度,温度越高,速度越快;对波系起实质影响作用的是入射激波强度,激波越强,则整个透射流场呈现偏向激波的趋势;理论分析还指出,稀疏波区的出现不可避免,当激波强度趋于声波稀疏波区趋于消失,激波越强则疏波区趋于扩大. 两波对撞存在一个有限的转变阶段,透射爆轰首先减缓,接着迅速迸发为过驱爆轰,然后再逐渐平衡为CJ爆轰. 对于强不稳定的燃气,对撞后爆轰波在空间上的发展极不均衡,一些区域发生火焰面与诱导激波的严重脱离,随后的火焰面失稳发展为诱导激波区内的爆轰波,实验观察到了这种爆轰在烟迹板上留下的极为精细的迹线.      

5.12汶川大地震诱发大型崩滑灾害动力特征初探

对乙炔氧气混合气体中爆轰波与激波的正面对撞现象的实验研究是以高速摄影获 取两波对撞的x-t纹影图,以烟迹板记录对撞中的爆轰胞格图案,并基于激波 理论和经典CJ爆轰理论求解了两波对撞的稳态解并探寻其规律. 研究发现透射波系包括一道激波和爆轰波, 以及紧随爆轰波后的稀疏波区,这一结果对应于一维理论分析中的CJ解. 透射波系基本不受 初始压强影响;初始温度也只成比例地改变流场整体速度,温度越高,速度越快;对波系起 实质影响作用的是入射激波强度,激波越强,则整个透射流场呈现偏向激波的趋势;理论分 析还指出,稀疏波区的出现不可避免,当激波强度趋于声波稀疏波区趋于消失,激波越强则 疏波区趋于扩大. 两波对撞存在一个有限的转变阶段,透射爆轰首先减缓,接着迅速迸发为 过驱爆轰,然后再逐渐平衡为CJ爆轰. 对于强不稳定的燃气,对撞后爆轰波在空间上的发展 极不均衡,一些区域发生火焰面与诱导激波的严重脱离,随后的火焰面失稳发展为诱导激波 区内的爆轰波,实验观察到了这种爆轰在烟迹板上留下的极为精细的迹线.     

激波聚焦诱导点火和爆轰的数值研究

以二维轴对称多组分Euler方程为基础，采用非正交结构化网格和改进的波传播算法，模拟了激波在抛物形反射壁面聚焦反射诱导点火和爆轰的过程，描述了其流场形态。讨论了预混气组成、入射激波强度及反射壁面形状对点火和爆轰的影响。结果表明，激波在抛物形反射壁面顶点处聚焦反射可形成局部高温高压区域，该区域在一定条件下可点燃预混气甚至形成爆轰，其中低稀释剂浓度的预混气、较大的入射激波Mach数和较深的反射壁面有利于可燃预混气形成爆轰。     

激波在不同三角楔面内传播特性的数值研究

激波在收缩管内的反射与聚焦会形成高温高压区,点燃可燃混合气并诱导爆轰,因此对爆轰发动机的点火具有重要意义。本文基于二维N-S方程,结合五阶WENO格式,对马赫数为6的正激波在三角形楔面内的反射与聚焦现象进行了数值研究。结果表明,楔面顶角的变化对激波的反射类型以及聚焦均有明显的影响:随着顶角的增加,激波的反射类型从马赫反射向过渡马赫反射和双马赫反射转变,且壁面上的前向射流更加明显;三波点第一次碰撞产生的高温高压区足够满足可燃混合气体的点火条件,且其温度与压力值随顶角的增加而增大;当激波在楔面上发生临界双马赫反射时,温度与压力达到最大;当顶角增加到一定值时,激波在楔面反射转变为常规反射,不会产生激波对碰,因而没有高温高压区。     

不稳定性对爆轰波楔面马赫反射的影响规律研究

采用一种两步化学反应模型对胞格爆轰波的楔面马赫反射过程进行了数值研究,从而澄清和解释胞格不稳定性对马赫反射发展模式和自相似性的影响。考虑到反应欧拉方程源项的刚性问题,本文采用附加RungeKutta方法耦合非刚性对流项和刚性反应源项,对流项的离散采用五阶精度的WENO格式。计算结果表明,对于稳定胞格爆轰波而言,其马赫反射过程本质上与ZND爆轰波的马赫反射是一致的,整体上不存在自相似性,胞格不稳定性只是造成了三波点轨迹线局部小振幅的波动。在楔面顶点附近,由于马赫杆是强过驱的,爆轰波的马赫反射过程是自相似的。在远场,爆轰波马赫反射的三波点轨迹线渐近的趋向于一条直线,说明重新获得了自相似性。对于不稳定的爆轰波,由于自身的不稳定性可以与马赫反射的强度相匹配,定义其三波点的轨迹是困难的,进行自相似性分析没有意义。     

激波与爆轰波对撞的数值模拟研究

用二阶精度NND差分格式和改进的二阶段化学反应模型模拟了爆轰波与激波的对撞过程,研究了不同强度入射激波对爆轰过渡区域的影响. 当对撞激波较弱时,透射爆轰波演变主要受流动膨胀作用的影响,可划分为对撞影响区、爆轰恢复区和稳定发展区3个阶段. 在爆轰恢复区和稳定发展区,前导激波压力经历一个过冲、然后向稳定爆轰过渡的过程,表现了爆轰波熄爆和再起爆的物理特征. 当对撞激波较强时,可燃混合气体的高热力学参数导致了更高的化学反应活化程度,形成了弱爆轰向稳定爆轰的直接转变.      

Near-wall imaging using toluene-based planar laser-induced fluorescence in shock tube flow

A quantitative thermometry technique, based on planar laser-induced fluorescence (PLIF), was applied to image temperature fields immediately next to walls in shock tube flows. Two types of near-wall flows were considered: the side wall thermal boundary layer behind an incident shock wave, and the end wall thermal layer behind a reflected shock wave. These thin layers are imaged with high spatial resolution (15μm/pixel) in conjunction with fused silica walls and near-UV bandpass filters to accurately measure fluorescence signal levels with minimal interferences from scatter and reflection at the wall surface. Nitrogen, hydrogen or argon gas were premixed with 1–12% toluene, the LIF tracer, and tested under various shock flow conditions. The measured pressures and temperatures ranged between 0.01 and 0.8 bar and 293 and 600 K, respectively. Temperature field measurements were found to be in good agreement with theoretical values calculated using 2-D laminar boundary layer and 1-D heat diffusion equations, respectively. In addition, PLIF images were taken at various time delays behind incident and reflected shock waves to observe the development of the side wall and end wall layers, respectively. The demonstrated diagnostic strategy can be used to accurately measure temperature to about 60 μm from the wall.     

弯管内爆轰波传播的流场显示和数值模拟

采用激光纹影系统拍摄了爆轰波在不同位置的流场照片. 用二阶附加半隐的龙格- 库塔法和五阶WENO格式 分别离散欧拉方程时间和空间导数项，用基元反应来描述爆轰化学反应过程，获得了压力、 温度、典型组元质量分数分布及数值胞格结构和爆轰波平均速度. 结果表明：受壁面稀疏波 和压缩波影响，爆轰波阵面发生畸变. 但由于弯管曲率半径较大，未出现爆轰波熄灭. 靠近 凹壁面的激波强度大于凸壁面侧，且凹壁面侧的反应区宽度较凸壁面侧要窄. 弯管出口处的 三波点数目较入口处减少，爆轰波衰减. 在出口直段，受扰动的爆轰波可恢复为自持爆轰波. 爆轰波流场、胞格结构、平均爆轰波速度的计算和实验结果定性一致.     

Supersonic flow of combustible gas mixture past sphere with account for ignition delay time

Assume an axisymmetric blunt body or a symmetric profile is located in a uniform supersonic combustible gas mixture stream with the parameters M₁, p₁, and T₁. A detached shock is formed ahead of the body and the mixture passing through the, shock is subjected to compression and heating. Various flow regimes behind the shock wave may be realized, depending on the freestream conditions. For low velocities, temperatures, or pressures in the free stream, the mixture heating may not be sufficient for its ignition, and the usual adiabatic flow about the body will take place. In the other limiting case the temperature behind the adiabatic shock and the degree of gas compression in the shock are so great that the mixture ignites instantaneously and burns directly behind the shock wave in an infinitesimally thin zone, i. e., a detonation wave is formed. The intermediate case corresponds to the regime in which the width of the reaction zone is comparable with the characteristic linear dimension of the problem, for example, the radius of curvature of the body at the stagnation point.The problem of supersonic flow of a combustible mixture past a body with the formation of a detonation front has been solved in [1, 2]. The initial mixture and the combustion products were considered perfect gases with various values of the adiabatic exponent .These studies investigated the effect of the magnitude of the reaction thermal effect and flow velocity on the flow pattern and the distribution of the gasdynamic functions behind the detonation wave.In particular, the calculations showed that the strong detonation wave which is formed ahead of the sphere gradually transforms into a Chapman-Jouguet wave at a finite distance from the axis of symmetry. For planar flow in the case of flow about a circular cylinder it is shown that the Chapman-Jouguet regime is established only asymptotically, i. e., at infinity.This result corresponds to the conclusions of [3, 4], in which a theoretical analysis is given of the asymptotic behavior of unsteady flows with planar, spherical, and cylindrical detonation waves.Available experimental data show that in many cases the detonation wave does not degenerate into a Chapman-Jouguet wave as it decays, bur rather at some distance from the body it splits into an adiabatic shock wave and a slow combustion front.The position of the bifurcation point cannot be determined within the framework of the zero thickness detonation front theory [1], and for the determination of the location of this point we must consider the structure of the combustion zone in the detonation wave. Such a study was made with very simple assumptions in [5].The present paper presents a numerical solution of the problem of combustible mixture flow about a sphere with a very simple model for the structure of the combustion zone, in which the entire flow behind the bow shock wave consists of two regions of adiabatic flow-an induction region and a region of equilibrium flow of products of combustion separated by the combustion front in which the mixture burns instantaneously. The solution is presented only for subsonic and transonic flow regions.     

NO PLIF imaging in the CUBRC 48-inch shock tunnel

Nitric oxide planar laser-induced fluorescence (NO PLIF) imaging is demonstrated at a 10-kHz repetition rate in the Calspan University at Buffalo Research Center??s (CUBRC) 48-inch Mach 9 hypervelocity shock tunnel using a pulse burst laser?Cbased high frame rate imaging system. Sequences of up to ten images are obtained internal to a supersonic combustor model, located within the shock tunnel, during a single ~10-millisecond duration run of the ground test facility. Comparison with a CFD simulation shows good overall qualitative agreement in the jet penetration and spreading observed with an average of forty individual PLIF images obtained during several facility runs.     

A combined OH/acetone planar laser-induced fluorescence imaging technique for visualizing combusting flows

A combined OH/acetone planar laser-induced fluorescence (PLIF) imaging technique that provides simultaneous visualizations of regions of unburned fuel and of combustion in a reacting flow is described. OH marks the location of chemical reaction and of combustion products, and acetone vapor, which is seeded into the fuel stream, marks unburned fuel. A single pulse from an ultraviolet laser is used to simultaneously excite both the OH and acetone, and the fluorescence from each is detected on separate cameras. Acetone spectroscopy and chemistry are reviewed to provide a basis for interpreting acetone fluorescence signals in high-temperature combusting environments. The imaging technique is applied to two nonpremixed turbulent reacting flows to assess the utility of the technique for visualizing the instantaneous flow structure and to illustrate the dependence of the interpretation of the acetone PLIF images on the flow conditions.Support was provided for this work by the Air Force Office of Scientific Research, Aerospace Sciences Directorate, with Julian Tishkoff as Technical Monitor, and is gratefully acknowledged. The contributions of Mr. T. C. Island in operating the supersonic flow facility are also greatly appreciated.     

Structure of the detonation wave front in a mixture of nitromethane with acetone

It is shown that the leading front of an inhomogeneous detonation wave is a shock wave in which wave structures of the type of triple shock configurations are moving. It was experimentally found that the reaction in these inhomogeneities occurs in oblique shock waves. The reaction sites at the wave front are ring-shaped. In a 75: 25 mixture of nitromethane with acetone, up to 70% of the front surface is occupied by the reaction at the sites in the wave front. Measurements of the mass velocity profile indicate that afterburning takes place in the unloading area behind the Jouguet plane. Calculations of the heat release in the reaction mixture with a decrease in the mass velocity indicate that the material that have not reacted in the inhomogeneities can be ignited in the induction zone. It is suggested that the adiabatic flashes are a mechanism that generates inhomogeneities in the detonation wave front.     

Experimental observation of the onset of detonation downstream of a perforated plate

In this study, the onset of detonation downstream of a perforated plate subsequent to the reflection of a Chapman–Jouguet detonation upstream is investigated. The experiments were performed with C₃H₈ + 5O₂ and C₂H₂+2.5O₂+70%Ar. The former has a much more irregular transverse wave pattern whereas the latter is known to have a piecewise laminar structure with a regular cellular structure. The onset of detonation phenomenon was found to be significantly different for the two mixtures. For the high argon diluted mixtures, the onset of detonation occurs in the vicinity downstream of the perforated plate. However, if the onset of detonation does not occur close to the plate, the precursor shock decouples from the reaction zone and a deflagration results. For the propane–oxygen mixtures, the onset of detonation is found to occur relatively far from the perforated plate at critical conditions. The major difference between these two mixtures is that the metastable turbulent reaction front can be maintained for relatively long distances for the propane–oxygen mixture. This turbulent metastable regime is also observed to be able to maintain a relatively constant propagation velocity for many channel widths prior to the onset of detonation. For the propane–oxygen mixtures, the onset is caused by a strong local explosion within the turbulent reaction zone.     

Planar laser-induced fluorescence imaging in shock tube flows

The planar laser-induced fluorescence (PLIF) imaging method was used to perform flow visualization and quantitative planar thermometry in shock tube flow fields using toluene as a fluorescence tracer in nitrogen. Fluorescence quantum yield values needed to quantify PLIF images were measured in a static cell at low pressures (<1 bar) for various toluene partial pressures in nitrogen bath gas. Images behind incident and reflected shocks were taken in the core flow away from regions affected by boundary layers. Temperature measurements from these images were successfully compared with predicted values using ideal shock equations. Measured temperatures ranged between 296 and 800 K and pressures between 0.15 and 1.5 atm. The average temperature discrepancies between measurements and the predicted values behind the incident and reflected shocks were 1.6 and 3.6%, respectively. Statistical analyses were also conducted to calculate the temperature measurement uncertainty as a function of image resolution. The technique was also applied to the study of more complex supersonic flows, specifically the interaction of a moving shock with a wedge. Measured temperatures agreed well with the results of numerical simulations in all inviscid regions, and all pertinent features of the single Mach reflection were resolved.     

超音速来流中爆轰波衍射和二次起爆过程研究

预爆管技术被广泛地应用在爆轰波发动机的起爆过程中，但是在超音速来流中基于预爆管技术起始爆轰波的研究并未被广泛地开展。基于此，本文中数值研究了横向超音速来流对半自由空间内爆轰波的衍射和自发二次起爆、及管道内的衍射和壁面反射二次起爆两种现象的影响。数值模拟的控制方程为二维欧拉方程，空间上使用五阶WENO格式进行数值离散，采用带有诱导步的两步链分支化学反应模型。所模拟的爆轰波具有规则的胞格结构，对应于用惰性气体高度稀释过的可爆混合物中形成的爆轰波。结果表明:在半自由空间内，在本文所模拟的几何尺寸下，爆轰波并未成功发生二次起爆现象，但是爆轰波的自持传播距离随着横向超音速来流强度的增强而增加。在核心的三角形流动区域外，波面诱导产生了更多的横波结构；在管道内，横向的超音速来流在逆流侧对出口气流产生了压缩作用，能有效提高波面压力，因此反射后的激波压力也比较高。在同样的几何尺寸下，爆轰波在静止和超音速(Ma=2.0)气流中分别出现了二次起爆失败和成功两种现象，这是由于在超音速来流中化学反应面的褶皱诱导产生了横波结构，横波与管壁以及其他横波之间的碰撞提高了前导激波的强度，并最终促进了爆轰波在超声速流主管道内的成功起始。