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

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

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

预混火焰界面的RM (Richtmyer-Meshkov)不稳定导致的界面混合区增长过程在自然界和工程实践中十分常见,但化学反应对其增长的影响机理仍不明确,反应性界面混合区增长速率的预测也未见报道, 因此,开展RM不稳定过程中火焰界面演化和混合区预测的研究十分必要.本文采用带单步化学反应的Navier-Stokes方程和高精度数值格式,研究了正弦形预混火焰界面在平面入射激波及其反射激波作用下的RM不稳定过程.结果表明, 在入射激波作用后的阶段,除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不稳定过程的理论预测方法提供了有益参考.     

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

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

聚心火焰与激波相互作用的数值研究

基于带化学反应的2维轴对称Euler方程,利用带有monotonized centered(MC)限制器的波传播算法,在两端敞开的圆桶中对惰性介质的聚心激波和氢气-空气混合物的聚心火焰与激波的相互作用进行了数值模拟。数值结果表明,在惰性介质中激波在轴心的每次汇聚均可成长为马赫干,马赫干的追赶使激波得到一定程度的增强,但整体还呈下降趋势。在氢气-空气混合物中,燃烧诱导的激波,由于与火焰的反复作用,使激波在轴心处产生马赫干的频率和强度皆高于惰性介质中的情形。同时,火焰在与激波的相互作用过程中发生失稳变形,使其形状呈扁平头部的蘑菇云。     

反射激波冲击重气柱的RM不稳定性数值研究

数值研究了二维气柱在入射激波以及反射激波作用下的Richtmyer-Meshkov(RM)不稳定性发展规律, 采用有限体积法结合网格自 适应技术的VAS2D程序, 精确刻画激波和界面的演化. 入射平面激波的马赫数为1.2, 气柱界面内气体为六氟化硫(SF₆), 环境气体为空气, 激波管的尾端为固壁. 通过改变气柱与尾端之间的距离调节反射激波再次作用已经变形的气柱的时间, 获得不同时刻下已经变形的气柱形态、界面尺寸以及环量演化受到反射激波的影响. 结果表明, 反射激波再次作用气柱时, 气柱所处发展阶段不同, 界面演化规律以及环量随时间的变化也不相同, 反射激波与气柱相互作用过程中的涡量产生和分布与无反射情况差异较大, 揭示了不同情况下界面演化的物理机理.     

Air／SF6界面RM不稳定性诱导的物质混合诊断

不同流体伴随激波的作用,在不同密度介质的界面处往往存在激波诱导的界面不稳定性,即RM(Richtmyer-Meshkov)不稳定性,由此不稳定诱导了物质间的相互混合。文中研究了低马赫(1.23)激波作用Air/SF6界面RM不稳定性问题。Air/SF6初始界面由厚度为1～2μm的硝化纤维薄膜相隔得到,利用阴影测试法研究了Air/SF6界面在1.23马赫数激波冲击下,界面混合宽度随时间的发展过程,以及反射激波作用后混合宽度的再增长。实验结果表明混合宽度线性发展前期与经验公式吻合较好。     

小扰动界面形态对RM不稳定性影响的数值分析

采用Navier-Stokes 方程对入射激波及其反射激波连续诱导小振幅扰动界面的Richtmyer-Meshkov 不稳定性增长过程进行了二维数值模拟,分析了单模和随机多模两种不同初始形态的界面上钉结构和泡结构在反射激波作用前后的发展特性. 研究结果发现:单模扰动的初始界面形态对反射激波前、后界面的扰动增长都有影响,反射激波前的界面形态信息可以通过钉和泡结构之间的反转传递到反射激波过后. 扰动界面上钉结构的发展速度控制了界面混合区总体的发展速度,反射激波前界面上发展成具有完整冠部形态的钉,在反射激波后会反转成复杂的泡结构,此泡结构对反射激波后钉的发展不利. 随机多模界面显示了与单模界面类似的发展规律,但随机多模界面上的复杂泡结构分布的不对称性使得其对钉结构增长的拖曳效应相对要弱,这导致了相似扰动波长下多模随机界面的扰动发展相对单模界面扰动发展要快.      

柱形汇聚激波冲击球形重气体界面的实验研究

采用高速纹影法实验研究了柱形汇聚激波与球形重气体界面相互作用的 Richtmyer-Meshkov不稳定性问题. 激波管实验段基于激波动力学理论设计, 将马赫数为1.2 的平面激波转化为柱形汇聚激波, 气体界面由肥皂膜分隔六氟化硫(内)和空气(外)得到. 采用高速摄影机在单次实验中拍摄激波运动的全过程, 对柱形激波的形成进行了实验验证, 并进一步观测了汇聚激波与球形气体界面相互作用过程中的波系发展和气体界面变形以及反射激波同已变形界面二次作用的流场演化. 结果表明: 当柱形汇聚激波穿过气泡界面以后, 气泡左侧界面极点沿激波传播方向保持匀速运动, 气泡右侧界面发展成为射流结构, 气泡主体发展成为涡环结构; 在反射激波的二次作用下, 流场中无序运动显著增强并很快进入湍流混合阶段.     

NUMERICAL STUDY ON THE EFFECTS OF SMALL-AMPLITUDE INITIAL PERTURBATIONS ON RM INSTABILITY

采用Navier-Stokes 方程对入射激波及其反射激波连续诱导小振幅扰动界面的Richtmyer-Meshkov 不稳定性增长过程进行了二维数值模拟,分析了单模和随机多模两种不同初始形态的界面上钉结构和泡结构在反射激波作用前后的发展特性. 研究结果发现：单模扰动的初始界面形态对反射激波前、后界面的扰动增长都有影响,反射激波前的界面形态信息可以通过钉和泡结构之间的反转传递到反射激波过后. 扰动界面上钉结构的发展速度控制了界面混合区总体的发展速度,反射激波前界面上发展成具有完整冠部形态的钉,在反射激波后会反转成复杂的泡结构,此泡结构对反射激波后钉的发展不利. 随机多模界面显示了与单模界面类似的发展规律,但随机多模界面上的复杂泡结构分布的不对称性使得其对钉结构增长的拖曳效应相对要弱,这导致了相似扰动波长下多模随机界面的扰动发展相对单模界面扰动发展要快.     

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.     

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.     

Interaction of a high-speed combustion front with a closely packed bed of spheres

The head-on collision of a combustion front with a closely packed bed of ceramic-oxide spheres was investigated in a vertical 76.2 mm diameter tube containing a nitrogen diluted stoichiometric ethylene–oxygen mixture. A layer of spherical beads in the diameter range of 3–12.7 mm was placed at the bottom of the tube and a flame was ignited at the top endplate. Four orifice plates spaced at one tube diameter were placed at the ignition end of the tube in order to accelerate the flame to either a “fast-flame” or a detonation wave before the bead layer face. The mixture reactivity was adjusted by varying the initial mixture pressure between 10 and 100 kPa absolute. The pressure before and within the bead layer was measured by flush wall-mounted pressure transducers. For initial pressures where a fast-flame interacts with the bead layer peak pressures recorded at the bead layer face were as high as five times the reflected Chapman–Jouget detonation pressure. The explosion resulting from the interaction developed by two distinct mechanisms; one due to the shock reflection off the bead layer face, and the other due to shock transmission and mixing of burned and unburned gas inside the bead layer. The measured explosion delay time (time after shock reflection from the bead layer face) was found to be independent of the incident shock velocity. As a result, the explosion initiation is not the direct result of the shock reflection process but instead is more likely due to the interaction of the reflected shock wave and the trailing flame. The bead layer was found to be very effective in attenuating the explosion front transmitted through the bead layer and thus isolating the tube endplate. This paper is based on work that was presented at the 21th International Colloquium on the Dynamics of Explosions and Reactive Systems, Poitiers, France, July 23–27, 2007.     

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.     

A high-speed photographic study of the transition from deflagration to detonation wave

Experiments were conducted to investigate the DDT process of the oxyhydrogen gas in the rectangular detonation tube of 3 m long. The repeated obstacle was installed near the ignition plug and the effects of the obstacle on the DDT process were investigated. The behaviour of the combustion and detonation wave were visualized utilizing Imacon high-speed camera with the aid of Schlieren optics. As a result, DDT process was visualized, i.e. (i) multiple shock waves were induced by the expanding combustion wave, because the combustion flame played a role as a piston and compressed the unburned gases. (ii) The acceleration of the combustion wave was occurred and the distance between the shock wave and the combustion flame became shorter. (iii) Eventually, the local explosion was occurred and cause overdriven detonation wave to propagate at the velocity of about 3 kms⁻¹. An abridged version of this paper was presented at the 15th Int. Colloquium on the Dynamics of Explosions and Reactive Systems at Boulder, Colorado, from July 30 to August 4, 1995     

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.     

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

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