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.     

正向爆轰驱动高焓激波风洞的数值模拟

对充满氢氧可燃气体、带扩容腔的正向爆轰驱动的激波风洞进行了数值模拟。计算采用了欧拉方程,频散可控耗散差分格式（DCD）和改进的二阶段化学反应模型。在扩容腔附近采用二维轴对称计算模型,而在驱动段和被驱动段的直管道部分则采用一维计算模型。本文分析了爆轰波在管道中的传播、反射和绕射过程。计算结果表明扩容腔的尺寸对爆轰波的传播、反射、汇聚等起着决定性的作用;带扩容腔的正向爆轰驱动的激波风洞能够得到平稳的持续时间较长的气流,提高了实验的精确度和可重复性。     

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

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

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

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

飞行Ma12条件超燃发动机流场及燃烧特征分析

为提升针对高马赫数发动机的模拟能力,对计算方法进行了可压缩性修正,并针对飞行Ma12条件下超燃冲压发动机进行了多状态三维数值模拟,分析了发动机内波系、参数以及燃烧性能特征.研究结果表明:(1)修正后的方法计算所得激波位置及强度与试验值吻合,在激波串模拟、高马赫数发动机模拟上均展现了更优的能力.(2)发动机内形成激波与反...     

Micromechanical modeling of shock-induced chemical reactions in heterogeneous multi-material powder mixtures

Shock waves create a unique environment of high pressure, temperature, strain rate, and gradients thereof. Chemical reactions that occur in this regime can lead to the synthesis of new materials that are not possible under conventional conditions. Shock-induced chemical reactions (also known as shock synthesis) is difficult to study experimentally due to the small time and space scales over which it occurs. To aid in our understanding, numerical simulations of shock-induced chemical reactions (SICR) are performed on the mesoscopic level. The model reaction, Nb+2Si⇌NbSi₂, is chosen for this research because it has been studied extensively. Comparisons of the calculated results with the available experimental data are presented to validate the modeling process, and additional calculations are used to determine the reaction threshold.     

Coarse-graining atomistic dynamics of brittle fracture by finite element method

In this work we present the finite element (FE) implementation of an atomistic formulation of balance equations and its application to coarse-grained (CG) simulation of dynamic fracture. First, we simulate a notched specimen that contains about 1.8 million atoms by the CG-FE method, and we compare the CG-FE results with that by all-atom molecular dynamics (MD) simulations. We find that CG-FE simulations with about 5% degrees of freedom of the MD simulation can capture the essential dynamic features, not in exact correspondence, but qualitatively and quantitatively similar to that obtained by MD simulations. We then proceed to simulate a series of micron-sized specimens by the CF-FE method. We find that it is the interaction of the forward propagating crack with the stress waves being reflected back by the boundaries of the specimen that triggers the dynamic instability and hence the branching of cracks in micron-sized specimens. The potential application of the method and future work for improvements are discussed.     

FE Formulation of Detonation Important by the Stabilized Galerkin Method

We introduce a finite element formulation of Detonation problems. The detonation involves supersonic combustion waves which emanate from an exothermic chemical reaction occurring behind propagating shock waves. The development of geometrically versatile schemes is required to meet practical needs. The Galerkin/Least-squares method with entropy variables possesses good accuracy and stability properties. As a first attempt to handle chemical reaction model, we use a simple two front model which has a good performance and a flexibility in adapting to the experimental data. It is emphasized that the process of generation and propagation of the shock triple points has been demonstrated.     

固体塑性实验室时间尺度的分子动力学模拟概述

随着超级计算机软硬件的飞速提升,基于经验势函数的分子动力学模拟在解析固体塑性的微观机制方面发挥着关键作用.但是,由于传统分子动力学基于牛顿运动方程数值积分,积分时间步长通常为飞秒量级,其模拟的时间尺度通常限于纳秒量级,从而为模拟长时间尺度固体塑性机制带来了巨大的挑战.本文从分子动力学模拟的时间尺度限制切入,介绍目前国际...     

爆轰波平掠惰性气体界面及其解耦现象的数值研究

采用基元化学反应模型和迎风TVD格式,数值研究了爆轰波平掠惰性气体界面时的物理现象及其作用机制,并用点隐方法克服化学反应源项引起的计算刚性。数值结果显示,当爆轰波平掠过惰性气体界面时,形成了爆轰波、界面、透射激波以及稀疏波相互作用的现象。在高N2比例稀释的可燃混合气体情况下,当爆轰波平掠过特定惰性气体界面时,它与惰性气体界面相互作用产生的稀疏波可以导致爆轰波的解耦。     

激波管双爆轰驱动段性能的数值模拟研究

采用频散可控的耗散格式(DCD),求解Euler方程和一种改进的二阶段化学反应模型, 对氢氧反向-正向双爆轰驱动段激波管进行了数值模拟. 计算结果表明：当辅驱动段与主驱动 段初始压力比小于临界值时,Taylor波仍会出现,但波扇夹角较单一前向爆轰驱动段小,入 射激波马赫数衰减率变小;当初始压力比等于临界值时,主驱动段中的Taylor波完全被消除, 入射激波马赫数不再衰减. 当初始压力比大于临界值时,在主驱动段中能产生过驱动爆轰波, 不仅Taylor波被完全消除,而且驱动能力较单一前向爆轰驱动段强.     

环形激波绕射, 反射和聚焦的数值模拟研究

应用频散可控耗散格式对环形激波在圆柱形激波管内绕射、反射和聚焦的问题进行了数值模拟研究，研究结果表明环形激波形成强烈聚焦的关键因素是环形激波在圆柱形管道中向对称轴运动时，绕射激波就不断加速而不作通常情况下的衰减；不同马赫数的环形激波绕射也产生不同马赫数及形状的准柱形激波，导致聚焦效果和位置的差异；另外，环形激波聚焦于一个点而圆柱形激波聚焦于一条线，两者有本质不同。     

Numerical simulations of the process of multiple shock–flame interactions

Based on a weighted essentially nonoscillatory scheme, the multiple interactions of a flame interface with an incident shock wave and its reshock waves are numerically simulated by solving the compressible reactive Navier–Stokes equations with a single-step Arrhenius chemical reaction. The two-dimensional sinusoidally perturbed premixed flames with different initial perturbed amplitudes are used to investigate the effect of the initial perturbation on the flame evolutions. The results show that the development of the flame interface is directly affected by the initial perturbed amplitudes before the passages of reshock waves, and the perturbation development is mainly controlled by the Richtmyer–Meshkov instability(RMI). After the successive impacts of multiple reshock waves, the chemical reaction accelerates the consumption of reactants and leads to a gradual disappearance of the initial perturbed information. The perturbation developments in frozen flows with the same initial interface as those in reactive flows are also demonstrated.Comparisons of results between the reactive and frozen flows show that a chemical reaction changes the perturbation pattern of the flame interface by decreasing the density gradient,thereby weakening the baroclinic torque in the flame mixing region, and therefore plays a dominant role after the passage of reshock waves.     

激波绕射碰撞加速诱导爆轰的数值模拟

基于单步化学反应的Euler方程和对激波(爆轰波)、接触间断具有良好捕捉效果的Roe/HLL混合格式以及自适应网格技术,模拟了激波在方形管中与方块障碍物相互作用,并发生绕射碰撞来诱导爆轰的过程.结果表明,弱激波在绕经方块时,形成上、下绕射激波并在方块尾部发生碰撞,生成局部高温高压点,可加快爆轰的形成;而当管内阻塞比超过...     

爆炸驱动典型活性材料能量释放特性研究

为了研究活性材料爆炸驱动反应特性,基于粉末压制成型工艺,制备了Al/PTFE、Al/Ni两种典型的活性材料及Al₂O₃/PTFE、Al₂O₃/PTFE/W惰性材料。通过爆炸驱动试验,并结合高速摄影、远红外热像仪以及峰值超压测试技术,分析了不同活性材料壳体装药爆炸火球、温度场分布及空气冲击波峰值超压等特性。同时,在炸药爆炸空气冲击波峰值超压经验计算模型中考虑了活性材料释放的化学能,分析了反应释放能量对空气冲击波的影响规律。结果表明：活性材料在爆炸驱动过程中经历了强加载条件下反应、产生碎片并向四周飞散、撞击钢板及后续反应等阶段。活性材料对炸药爆炸产生的空气冲击波具有强化作用,爆炸加载瞬间材料仅发生了部分化学反应。

     

Hot spot formation from shock reflections

Heterogeneities sensitize an explosive to shock initiation. This is due to hot-spot formation and the sensitivity of chemical reaction rates to temperature. Here, we describe a numerical experiment aimed at elucidating a mechanism for hot-spot formation that occurs when a shock wave passes over a high-density impurity. The simulation performed is motivated by a physical experiment in which glass beads are added to liquid nitromethane. The impedance mismatch between the beads and the nitromethane results in shock reflections. These, in turn, give rise to transverse waves along the lead shock front. Hot spots arise on local portions of the lead front with a higher shock strength, rather than on the reflected shocks behind the beads. Moreover, the interactions generated by reflected waves from neighboring beads can significantly increase the peak hot-spot temperature when the beads are suitably spaced.     

A chemical shock tube driven by detonation

A chemical shock tube driven by a detonation driver is described in the present paper. This shock tube can produce a single controlled high-temperature pulse for studies of gas-phase reaction kinetics, but the difficulty associated with the timing for the rupture of diaphragms in the conventional chemical shock tube is overcome, because the detonation wave in the driver section can be predicted correctly and shows a good repeatability. In addition, this shock tube is capable of providing higher temperature conditions for the test gas than the conventional high-pressure shock tube, owing to the inherently high-driving capability of the detonation driver. The feasibility of this shock tube is examined by numerical simulations and preliminary experiments.     

On the interference of a weak and a strong shock wave in a dissociating nitrogen flow

The influence of the nitrogen dissociation on the interactions due to the interference of two planar shock waves in a hypersonic high enthalpy flow is theoretically investigated for infinite reaction rates. The two limiting cases of infinitely slow and infinitely fast reactions are modelled as a perfect gas and an ideal dissociating gas in chemical equilibrium.To investigate the influence of finite reaction rates on the interactions of shock waves, experiments are performed in the high enthalpy shock tunnel Göttingen (HEG) with a wind tunnel model consisting of a wedge type shock generator and a transversally mounted cylinder. The pressure and heat transfer loads resulting from the shock wave interferences are measured and the flow field is visualized by means of interferograms. The experimental results are compared with the results of a numerical simulation for a dissociating nitrogen flow and with the experimental results for a perfect gas flow.     

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

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