激波诱导火焰变形、混合和燃烧的数值研究

为了深入研究激波诱导的火焰变形以及由此带来的混合和燃烧变化特性,采用带单步化学反应的Navier-Stokes方程和高网格分辨率,对平面入射激波及其反射激波诱导球形火焰变形的现象进行了二维数值研究,计算结果与实验结果较好吻合。研究结果显示,在反射激波作用前,火焰的变形和皱褶主要受入射激波诱导等物理过程影响;而在反射激波与失稳火焰作用后,燃烧放热率、火焰有效面积和界面长度均迅速增加,控制火焰变形的机制逐渐向化学反应(燃烧)过程过渡;在失稳火焰发展的后期,增强的燃烧过程能够削弱火焰界面的皱褶,进而抑制未燃气和可燃气的混合。由此可以得出结论,激波诱导的火焰界面通过变形可促进界面两侧未燃气与可燃气的混合,进而强化燃烧过程,但燃烧的增强却反过来能抑制混合;认识两者之间的关系有助于利用或控制激波 火焰相互作用过程。     

Explosive-driven shock wave and vortex ring interaction with a propane flame

Experiments were performed to analyze the interaction of an explosively driven shock wave and a propane flame. A 30 g explosive charge was detonated at one end of a 3-m-long, 0.6-m-diameter shock tube to produce a shock wave which propagated into the atmosphere. A propane flame source was positioned at various locations outside of the shock tube to investigate the effect of different strength shock waves. High-speed retroreflective shadowgraph imaging visualized the shock wave motion and flame response, while a synchronized color camera imaged the flame directly. The explosively driven shock tube was shown to produce a repeatable shock wave and vortex ring. Digital streak images show the shock wave and vortex ring propagation and expansion. The shadowgrams show that the shock wave extinguishes the propane flame by pushing it off of the fuel source. Even a weak shock wave was found to be capable of extinguishing the flame.     

Effect of chemical reactivity on the detonation initiation in shock accelerated flow in a confined space

The interactions of a spherical flame with an in- cident shock wave and its reflected shock wave in a confined space were investigated using the three-dimensional reactive Navier-Stokes equations, with emphasis placed on the effect of chemical reactivity of mixture on the flame distortion and detonation initiation after the passage of the reflected shock wave. It is shown that the spatio-temporal characteristics of detonation initiation depend highly on the chemical reactivity of the mixture. When the chemical reactivity enhances, the flame can be severely distorted to form a reactive shock bifurcation structure with detonations initiating at different three-dimensional spatial locations. Moreover, the detonation initiation would occur earlier in a mixture of more enhanced reactivity. The results reveal that the detonations arise from hot spots in the unburned region which are initiated by the shock-detonation-transition mechanism.     

激波冲击火焰的涡量特性研究

激波冲击火焰的现象涉及一系列复杂的物理化学过程,其中涡量的生成与演化对控制火焰发展起重要作用。为系统分析激波冲击火焰过程中的涡量特性,采用二维带化学反应的Navier-Stokes方程对平面入射激波及其反射激波与球形火焰作用的现象进行了数值研究,通过引入并行计算达到高网格分辨率的要求。计算结果表明,斜压项对火焰区内涡量生成起主导作用,压缩项和耗散项在火焰膨胀阶段抑制涡量生成,此外,火焰在激波压缩阶段主要受物理过程而非化学反应过程影响。     

激波诱导火焰失稳与爆轰的条件研究

采用九阶WENO和十阶中心差分格式数值求解激波与火焰作用过程,考察了激波强度、火焰尺寸对激波与球形火焰作用过程的影响。结果表明,增大激波强度或火焰尺寸均可在流场中引发爆轰,但激波强度的影响更大,并且其引发的爆轰可使火焰迅速膨胀,放热率提高,从而影响燃烧特性;此外,爆轰波传播过程中会迅速消耗可燃预混气,合并原有的反射激波,并在流场中形成局部高压区,极大地改变流场结构。     

A three-dimensional numerical study on instability of sinusoidal flame induced by multiple shock waves

The instabilities of a three-dimensional sinusoidally premixed flame induced by an incident shock wave with Mach = 1.7 and its reshock waves were studied by using the Navier–Stokes(NS) equations with a single-step chemical reaction and a high resolution, 9th-order weighted essentially non-oscillatory scheme. The computational results were validated by the grid independence test and the experimental results in the literature. The computational results show that after the passage of incident shock wave the flame interface develops in symmetric structure accompanied by large-scale transverse vortex structures. After the interactions by successive reshock waves, the flame interface is gradually destabilized and broken up, and the large-scale vortex structures are gradually transformed into small-scale vortex structures. The small-scale vortices tend to be isotropic later.The results also reveal that the evolution of the flame interface is affected by both mixing process and chemical reaction. In order to identify the relationship between the mixing and the chemical reaction, a dimensionless parameter, η, that is defined as the ratio of mixing time scale to chemical reaction time scale, is introduced. It is found that at each interaction stage the effect of chemical reaction is enhanced with time.The enhanced effect of chemical reaction at the interaction stage by incident shock wave is greater than that at the interaction stages by reshock waves. The result suggests that the parameter η can reasonably character the features of flame interface development induced by the multiple shock waves.     

铝粉湍流火焰诱导激波现象的实验研究

粉尘湍流火焰诱导激波问题是工业灾害研究中的重要课题．本文在自行设计的大型卧式燃烧管内，对铝粉火焰诱导激波现象进行了实验研究，测试了湍流火焰阵面前压缩波到激波的转捩过程，并将实验结果与数值模拟结果进行了比较．     

Detonation Initiation by Shock Reflection from an Orifice Plate

Results from an experimental investigation of the interaction of a “non-ideal” shock wave and a single obstacle are reported. The shock wave is produced ahead of an accelerated flame in a 14 cm inner-diameter tube partially filled with orifice plates. The shock wave interacts with a single larger blockage orifice plate placed 15–45 cm after the last orifice plate in the flame acceleration section of the tube. Experiments were performed with stoichiometric ethylene–oxygen mixtures with varying amounts of nitrogen dilution at atmospheric pressure and temperature. The critical nitrogen dilution was found for detonation initiation. It is shown that detonation initiation occurs if the chemical induction time based on the reflected shock state is shorter than the time required for an acoustic wave to traverse the orifice plate upstream surface, from the inner to the outer diameter. The similarity between the present results and those obtained from previous investigators looking at detonation initiation by ideal shock reflection produced in a shock tube indicates that the phenomenon is not sensitive to the detailed structure of the shock front but only on the average shock strength.This paper is based on work that was presented at the 20th International Colloquium on the Dynamics of Explosions and Reactive Systems, Montreal, Canada, July 31–August 5, 2005.     

柔性障碍物对甲烷空气爆炸波激励作用的实验研究

为研究柔性障碍物对甲烷空气爆炸波的激励效应,采用双向拉伸聚丙烯（biaxially oriented polypropylene, BOPP）薄膜作为柔性障碍物将管道内甲烷空气预混气体与空气隔开,对比障碍物前后火焰、激波变化,分析膜状柔性障碍物激励效应的机理。实验结果表明:这种具有一定承压能力的柔性障碍物对甲烷爆炸波产生的激励效应不可忽视,在膜片破裂前产生多次激波反射过程,可诱导湍流火焰形成,促使膜前爆炸压力提高,膜片破裂后,火焰在伴流作用下传播速度突增,并加速逐渐逼近前驱冲击波,致使膜后爆炸压力大幅提高;激励效应可使膜片前后最大爆炸压力相差5倍,火焰速度相差7倍;另外在膜片位置2.5 m后增设一道膜片,可增强这种激励效应,而增加膜片的实质是使激波火焰相互作用的次数增加。     

Interaction of a reflected shock from a concave wall with a flame distorted by an incident shock

Observations are presented from calculations where a laminar spherical CH₄/air flame was perturbed successively by incident and reflected shock waves reflected from a planar or concave wall. The two-dimensional axi-symmetric Navier–Stokes equations with detailed chemistry were used. The computational results were qualitatively validated with experiments which were performed in a standard shock tube arrangement. Under the influence of the incident shock wave, a Richtmyer–Meshkov instability is induced in the flame, and the distorted flame finally takes the form of two separated elliptical burning bubbles in the symmetric cross plane. Then, under subsequent interactions with the shock wave reflected from the planar or the concave wall, the flame takes a mushroom-like shape. Transverse waves produced by the shock reflection from the concave wall can compress the flame towards the axis, and the focusing shock generated on the concave wall will lead to a larger mushroom-like flame than that induced by the planar reflection.      

多次激波诱导正弦扰动预混火焰界面失稳的数值研究

激波诱导火焰失稳是实际中常见的现象,为深入研究火焰失稳特性,采用三维单步化学反应的Navier-Stokes方程和9阶weighted essentially non-oscillatory (WENO)的高精度格式,对不同马赫数的入射激波及其反射激波多次诱导正弦型预混火焰界面失稳的现象进行了三维数值模拟,并对计算结果的可靠性进行了验证。研究结果显示,在激波的多次作用下,火焰界面的演变过程主要受Richtmyer-Meshkov (RM)不稳定特性和化学反应特性的双重影响,且随着入射激波强度的增强,上述2种特性均得到进一步强化。为构造体现反应性RM不稳定特性的参数,根据火焰界面混合区平均涡量和化学反应速率,提出了表征界面受不稳定性和化学反应影响的量纲一参数。通过分析发现,在同一入射激波强度下,该参数的对数形式随入射激波和反射激波的多次作用呈基本相同的线性增长趋势;对不同马赫数的入射激波,该参数对数形式的线性增长率也基本一致。这样的变化表明该量纲一参数能够反映反应性RM不稳定过程中火焰界面发展的内在规律。     

INTERFACE EVOLUTIONS AND GROWTH PREDICTIONS OF MIXING ZONE ON PREMIXED FLAME INTERFACE DURING RM INSTABILITY 1)

预混火焰界面的RM (Richtmyer-Meshkov)不稳定导致的界面混合区增长过程在自然界和工程实践中十分常见,但化学反应对其增长的影响机理仍不明确,反应性界面混合区增长速率的预测也未见报道, 因此,开展RM不稳定过程中火焰界面演化和混合区预测的研究十分必要.本文采用带单步化学反应的Navier-Stokes方程和高精度数值格式,研究了正弦形预混火焰界面在平面入射激波及其反射激波作用下的RM不稳定过程.结果表明, 在入射激波作用后的阶段,除RM不稳定本身导致的界面演化为"钉-帽"和"泡"形结构外,化学反应一方面以预混火焰传播的方式促进了界面中"泡"结构的增长,另一方面通过与涡结构的复杂相互作用促进了"钉-帽"结构的增长.化学反应活性越强, 火焰界面的"泡" 结构和"钉-帽"结构的增长越快.在第一次反射激波作用后的阶段,化学反应以相同的火焰传播方式对"泡"和"钉-帽"结构产生影响, 两者效应相抵,因而导致反射激波作用后的阶段中界面混合区增长不受化学反应活性的影响.根据以上分析,分别针对入射激波和第一次反射激波作用后的火焰界面混合区增长速率提出了相应的预测模型,为探索反应性RM不稳定过程的理论预测方法提供了有益参考.     

RM不稳定过程中预混火焰界面演化及混合区增长预测

预混火焰界面的RM (Richtmyer-Meshkov)不稳定导致的界面混合区增长过程在自然界和工程实践中十分常见,但化学反应对其增长的影响机理仍不明确,反应性界面混合区增长速率的预测也未见报道, 因此,开展RM不稳定过程中火焰界面演化和混合区预测的研究十分必要.本文采用带单步化学反应的Navier-Stokes方程和高精度数值格式,研究了正弦形预混火焰界面在平面入射激波及其反射激波作用下的RM不稳定过程.结果表明, 在入射激波作用后的阶段,除RM不稳定本身导致的界面演化为"钉-帽"和"泡"形结构外,化学反应一方面以预混火焰传播的方式促进了界面中"泡"结构的增长,另一方面通过与涡结构的复杂相互作用促进了"钉-帽"结构的增长.化学反应活性越强, 火焰界面的"泡" 结构和"钉-帽"结构的增长越快.在第一次反射激波作用后的阶段,化学反应以相同的火焰传播方式对"泡"和"钉-帽"结构产生影响, 两者效应相抵,因而导致反射激波作用后的阶段中界面混合区增长不受化学反应活性的影响.根据以上分析,分别针对入射激波和第一次反射激波作用后的火焰界面混合区增长速率提出了相应的预测模型,为探索反应性RM不稳定过程的理论预测方法提供了有益参考.      

Numerical investigation of ethylene flame bubble instability induced by shock waves

In this paper, the ethylene/oxygen/nitrogen premixed flame instabilities induced by incident and reflected shock wave were investigated numerically. The effects of grid resolutions and chemical mechanisms on the flame bubble deformation process are evaluated. In the computational frame, the 2D multi-component Navier–Stokes equations with second-order flux-difference splitting scheme were used; the stiff chemical source term was integrated using an implicit ordinary differential equations (ODEs) solver. The two ethylene/oxygen/nitrogen chemical mechanisms, namely 3-step reduced mechanism and 35-step elementary skeletal mechanism, were used to examine the reliability of chemistry. On the other hand, the different grid sizes, Δx × Δy = 0.25 × 0.5mm and Δx × Δy = 0.15 × 0.2mm, were implemented to examine the accuracy of the grid resolution. The computational results were qualitatively validated with experimental results of Thomas et al. (Combust Theory Model 5:573–594, 2001). Two chemical mechanisms and two grid resolutions used in present study can qualitatively reproduce the ethylene spherical flame instability process generated by an incident shock wave of Mach number 1.7. For the case of interaction between the flame and reflected shock waves, the 35-steps mechanism qualitatively predicts the physical process and is somewhat independent on the grid resolutions, while the 3-steps mechanism fails to reproduce the instability of ethylene flame for the two selected grid resolutions. It is concluded that the detailed chemical mechanism, which includes the chain elementary reactions of fuel combustion, describes the flame instability induced by shock wave, in spite of the fact that the flame thickness (reaction zone) is represented by 1–2 grids only.      

Experimental and numerical studies on interactions of a spherical flame with incident and reflected shocks

Observations are presented from experiments and calculations where a laminar spherical CH4/air flame is perturbed successively by incident and reflected shock waves. The experiments are performed in a standard shock tube arrangement, in which a high-speed shadowgraph imaging system is used to record evolutions of the flame. Numerical simulations are conducted by using second-order wave propagation algorithms, based on two-dimensional axisymmetric Navier-Stokes equations with detailed chemical reactions. Qualitative agreements are obtained between the experimental and numerical results. Under actions of incident shock waves, Richtmyer-Meshkov instability responsible for the flame deformation is induced in the flame, and the distoned flame takes a barrel shape. Then, under subsequent actions of the shock wave reflected from a planar wall, the flame takes an inclined non-symmetrical kidney shape in a symmetric cross section, which means a mushroom-like shape of the flame comes finally into being. The vorticity direction in the ring cap has been altered by the reflected shock＇s action, which makes the head of the mushroom-like flame extend quickly to the side wall.     

两相介质中火焰诱导激波现象的研究

基于双流体模型，采用跟踪法处理火焰内边界，ＴＶＤ格式计算气相激波，ＭａｃＣｏ－ｍａｃｋ格式计算颗粒相流场，对两相可燃介质中火焰诱导激波现象进行了理论分析与实验模拟，研究了颗粒相参数对火焰诱导激波现象的影响。     

Regular reflection of a shock wave from a solid wall in a stream of fuel mixture

The two-dimensional stationary problem of regular reflection of a shock wave from a plane solid wall in a fuel gas mixture is examined in the case when the mixture is ignited at the intersection of the incident wave with the wall and a flame front is formed behind the reflected shock wave. The shock waves and the flame front are considered plane surfaces of discontinuity. The fuel mixture and the reaction products are considered perfect, inviscid, and non-heat-conducting gases.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 161–163, July–August, 1978.     

瓦斯/空气预混气体爆燃流场测试系统多参数耦合同步控制实验方法

为了更精确地获得爆炸激波管内瓦斯/空气预混气体爆燃过程中，激波形成过程、压力和火焰传播速度以及火焰与惰性阻燃剂相互作用的流场演化图像。通过分析激波管测试系统中多个目标的时间响应特征及控制方式，利用超高速相机、光电倍增管、时间延时器、固态继电器、电荷放大器和数据采集系统等设备，设计实验方案，分别对激波管中瓦斯/空气预混气体爆燃高压点火系统的响应时间和惰性介质阻燃剂喷射系统的响应时间进行测试。实验结果表明电火花点火的响应时间为微秒量级，而阻燃剂喷射系统的响应时间为毫秒量级，以响应时间为依据，通过设置精确的延迟时间实现多目标同步控制，为完成激波管内瓦斯/空气预混气体爆燃过程的微观流场显示奠定基础。     

Transition of plane overdriven detonation wave to the Chapman-Jouguet regime

The asymptotic laws of behavior for plane, cylindrical, and spherical infinitely thin detonation waves were found in [1, 2] for increasing distance from an igniting source in those cases in which the waves changed into Chapman-Jouguet waves as they decayed. It was shown that the plane overdriven detonation wave approaches the Chapman-Jouguet regime asymptotically, while the transition of the cylindrical or spherical strong detonation wave into the Chapman-Jouguet wave may occur at a finite distance from the initiation source.Similar conclusions are valid for the propagation of stationary steadystate detonation waves which arise with flow of combustible gas mixtures past bodies.However, numerous experiments [3, 4] on firing bodies in a detonating gas show that the overdriven detonation wave which forms ahead of the body decays and decomposes into an ordinary compression shock and a slow combustion front. To establish why the wave does not make the transition to the Chapman-Jouguet regime, in the following we consider the propagation of a plane detonation wave and account for finite chemical reaction rates. We use the very simple two-front model (ordinary shock wave and following flame front). Conditions are found for which transition to the Chapman-Jouguet regime does not occur. We first consider the propagation of an unsteady plane wave and then the steady plane wave. It is found that for all the mixtures used in these experiments transition to the Chapman-Jouguet regime is not possible within the framework of the assumed model.     

Numerical study of three-dimensional developments of premixed flame induced by multiple shock waves

The three-dimensional interactions of a perturbed premixed flame interface with a planar incident shock wave and its reflected shock waves are numerically simulated by solving the compressible, reactive Navier–Stokes equations with the high-resolution scheme and a single-step chemical reaction. The effects of the initial incident shock wave strength (Mach number) and the initial perturbation pattern of interface on the interactions are investigated. The distinct properties of perturbation growth on the flame interface during the interactions are presented. Our results show that perturbation growth is mainly attributed to the flame stretching and propagation. The flame stretching is associated with the larger-scale vortical flow due to Richtmyer–Meshkov instability while the flame propagation is due to the chemical reaction. The mixing properties of unburned/burned gases on both sides of the flame are quantitatively analyzed by using integral and statistical diagnostics. The results show that the large-scale flow due to the vortical motion always plays a dominating role during the reactive interaction process; however, the effect of chemistry becomes more important at the later stage of the interactions, especially for higher Mach number cases. The scalar dissipation due to the molecular diffusion is always small in the present study and can be negligible.