矩形管内临界爆轰动力学数值分析

 对矩形管内临界爆轰动力学特征进行了数值分析。采用基元反应描述爆轰化学反应过程,采用二阶附加半隐的龙格-库塔法和5阶WENO格式求解二维反应欧拉方程。对于25%氩稀释化学计量比的氢氧预混气体,当管道宽度为30 mm、初温为300 K时,产生临界爆轰的预混气体初压为3.5 kPa。在此临界条件下,获得了临界爆轰胞格结构、沿壁面的速度和峰值压力曲线及流场波系演变特征。着重对比分析了矩形管内临界爆轰与普通爆轰在爆轰波速度、平均速度、胞格宽长比、横波结构、未反应气囊及旋涡结构之间的差异,深入认识了临界爆轰的不稳定性和化学反应动力学特征。     

气相爆轰波在分叉管中传播现象的数值研究

数值研究气相爆轰波在分叉管中的传播现象.用二阶附加半隐龙格-库塔法和5阶WENO格式求解二维欧拉方程,用基元反应描述爆轰化学反应过程,得到了密度、压力、温度、典型组元质量分数场及数值胞格结构和爆轰波平均速度.结果表明:气相爆轰波在分叉管中传播,分叉口左尖点的稀疏波导致诱导激波后压力、温度急剧下降,诱导激波和化学反应区分离,爆轰波衰减为爆燃波(即爆轰熄灭).分离后的诱导激波在垂直支管右壁面反射,并导致二次起爆.畸变的诱导激波在水平和垂直支管中均发生马赫反射.分叉口上游均匀胞格区和分叉口附近大胞格区的边界不是直线,其起点通常位于分叉口左尖点上游或恰在左尖点.水平支管中马赫反射三波点迹线始于右尖点下游.分叉口左尖点附近的流场中出现了复杂的旋涡结构、未反应区及激波与旋涡作用.旋涡加速了未反应区的化学反应速率.反射激波与旋涡作用并使旋涡破碎.反射激波与未反应区作用,加速其反应消耗,并形成一个内嵌的射流.数值计算得到的波系演变和胞格结构与实验定性一致.     

Energy release effect of mixture on single spinning detonation structure

Unsteady three-dimensional numerical simulation on a single spinning detonation in a circular tube are presented in order to understand the effects of energy release of the mixture on the detonation structure. Overall structures of the spinning detonations such as the shock structure around the spin head, the long pressure trail, and the track angle on the wall are not affected by these effects because they depend on the specific heat ratio of the products which has approximately a constant value. The calculated averaged detonation velocities on the symmetry axis during one cycle decrease inversely with an exponential curve to become the value lower than the CJ detonation velocity. Those for p₀ = 0.1 MPa and p₀ = 0.01 MPa become approximately 0.98 D_CJ and 0.92 D_CJ, respectively, because the energy release in the CJ state for p₀ = 0.01 MPa is 10% lower than that for p₀ = 0.1 MPa. The state of gas behind the head of spinning detonation is also evaluated by the classical oblique shock theory and equilibrium calculation by using the track angle, shock wave angle, and detonation velocity in order to compare with the present and other researcher’s numerical results. The effects of the energy release in the mixture are large on the strength of the transverse detonation.     

并行建表化学加速算法研究及其在气相斜爆轰模拟中的应用

采用详细化学反应机理对气相斜爆轰问题开展数值计算时,由于组分之间的特征时间尺度相差很大,反应源项的直接积分(direct integration,DI)求解通常存在强烈的刚性及非线性现象,导致计算量很大.为了在不损失计算精度的基础上有效减少化学反应过程的计算时间,针对包含2H2+O2详细机理的二维斜爆轰并行计算,提出两...     

Numerical study on spinning detonations

Spinning detonations in both a circular tube and a square tube are presented in order to reveal characteristics of spinning modes by using three-dimensional simulations with a detailed chemical reaction model. The present results show a feature of a single spinning detonation which was discovered in 1926. The shock patterns in both cases are similar except the pressure trail, however, the shock wave angles and the shock wave lengths are shown to be dependent on the cross section configuration of the tube. The pitch angle, the track angle, the Mach stem angle, and the incident shock angle on the tube wall in the numerical results agree well with those in the experimental ones, and they are independent of the compositions of mixture, tube diameters, and initial pressures.     

基于基元反应模型的H2-O2-N2爆轰数值模拟

 采用12组分、23个化学反应的基元化学反应模型,用5阶加权本质无震荡格式（WENO）、3阶TVD Runge-Kutta格式,对H₂-O₂-N₂混合气体胞格爆轰进行了数值模拟。研究了一维ZND爆轰、自维持爆轰的详细结构以及三波点附近的流动结构。计算结果表明：由横波的压力可以显著促进二维爆轰波波阵面的形成;横波的运动生成三波点,三波点造成了爆轰的自维持传播。     

Numerical investigation on propagation mechanism of spinning detonation in a circular tube

Spinning detonations propagating in a circular tube were numerically investigated with a one-step irreversible reaction model governed by Arrhenius kinetics. The time evolution of the simulation results was utilized to reveal the propagation mechanism of single-headed spinning detonation. The track angle of soot record on the tube wall was numerically reproduced with various levels of activation energy, and the simulated unique angle was the same as that of the previous reports. The maximum pressure histories of the shock front on the tube wall showed stable and unstable pitch modes for the lower and higher activation energies, respectively. The shock front shapes and the pressure profiles on the tube wall clarified the mechanisms of two modes. The maximum pressure history in the stable pitch remained nearly constant, and the single Mach leg existing on the shock front rotated at a constant speed. The high and low frequency pressure oscillations appeared in the unstable pitch due to the generation and decay of complex Mach interaction on the shock front shape. The high-frequency oscillation was self-induced because the intensity of the transverse wave was changed during propagation in one cycle. The high-frequency behavior was not always the same for each cycle, and therefore the low frequency oscillation was also induced in the pressure history.     

炸药爆轰的连续介质本构模型和数值计算方法

假定炸药和爆轰产物处于局部热力学平衡状态, 即它们的压力和温度相同, 利用热力学基本关系建立炸药爆轰过程的连续介质本构模型的一般理论框架. 在此框架下, 炸药爆轰本构模型由一组常微分方程构成, 包括炸药和爆轰产物的状态方程、简单混合法则、化学反应速率方程和能量守恒方程, 易于由成熟的计算方法如梯形法等进行求解. 一组广义Maxwell型非线性固体本构形式的微分方程描述了压力和温度随时间的演化速率与应变率和化学反应速率的关系, 借助简单混合物理论, 其中的系数由炸药和爆轰产物的材料参数确定. 未反应的炸药和爆轰产物采用JWL状态方程, 化学反应率方程采用Lee-Tarver点火-燃烧二项式模型, 模拟PBX-9404炸药的一维冲击波起爆过程和爆轰波传播过程. 计算结果表明了本文给出的本构模型和相应计算方法的有效性. 关键词： 炸药爆轰 本构模型 化学反应率方程 数值模拟     

Dynamics of detonation waves in a channel with variable cross section and filled with bubbly fluid

The flow of bubbly fluid comprising a mixture of bubbles filled with explosive and inert gases, which is driven through a converging channel, was studied. Depending on the velocity of the hummer hitting the bubbly fluid boundary, the flow may be accompanied by the development of detonation waves which compress the bubbles with inert gas.     

隔板深度对激波聚焦起爆影响的数值模拟研究

考虑几何结构参数对激波聚焦触发爆轰波的复杂影响,对H2/Air预混气的环形射流激波聚焦起爆现象开展了数值模拟研究,详细分析了不同隔板深度条件下的激波聚焦过程、流场演化特征以及爆轰波参数变化规律.研究结果表明,凹腔内激波聚焦诱导的局部爆炸以及隔板前缘处射流形成"卷吸涡"是引起爆轰波触发的两个重要机制,而隔板深度是影响环形...     

氢气-空气混合物中瞬态爆轰过程的二维数值模拟

 对高温火团引发的氢气-空气混合气的瞬态爆轰过程进行了二维数值模拟,考虑了H₂-O₂-N₂的详细化学反应动力学机理,该机理包含了19个基元反应和9种组分。采用分裂格式处理带化学反应的Euler方程,其中使用全耦合的TVD格式求解流场,使用基于Gear算法的微分方程解法器求解化学反应过程。计算结果表明：在H₂∶O₂∶N₂＝0.4∶0.4∶0.2（摩尔比）的混合气中,高温气团初始温度为T/T₀=5.3时可诱导爆轰,爆轰波以2 300 m/s的速度传播,同时爆轰波阵面在管壁会形成反射波。还对计算的爆轰波后组分的浓度和温度进行了讨论,为理解爆轰波后结构提供信息。     

爆轰波在弯管内传播过程数值分析

应用基元反应模型和频散可控耗散格式(DCD)对氢氧爆轰波在弯管内的传播过程进行了数值模拟.计算中氢氧混合物化学反应采用了8种组分20个反应方程式.在处理化学反应引起的刚性问题时采用了时间算子分裂的方法.计算结果表明,在弯管小曲率半径壁面附近,由于膨胀稀疏作用,爆轰波强度减弱,在局部出现前导激波与放热反应区的解藕以及二次起爆现象;在弯管大曲率半径壁面上爆轰波在马赫反射和正规反射之间相互转变,使爆轰波加强.弯管内的爆轰现象与弯管曲率半径有关.     

变密度随机涡模型在湍流射流扩散火焰中的应用

采用变密度随机涡模型,对H₂/O₂/N₂湍流射流扩散火焰进行数值模拟,湍流过程通过涡的采样、涡的抑制和涡的翻转实现.其中,针对变密度反应流问题,提出一种大涡抑制的新机制,并详细讨论各种参数对模型预测效果的影响.计算结果表明,修改后的模型可以合理预测H₂/O₂/N₂射流火焰结构,能够反映湍流的涡特性;模型中与涡采样和涡抑制有关的参数对预测结果有一定影响.     

Numerical studies on autoignition and detonation development from a hot spot in hydrogen/air mixtures

Detonation development inside spark ignition engines can result in the so called super-knock with extremely high pressure oscillation above 200?atm. In this study, numerical simulations of autoignitive reaction front propagation in hydrogen/air mixtures are conducted and the detonation development regime is investigated. A hot spot with linear temperature distribution is used to induce autoignitive reaction front propagation. With the change of temperature gradient or hot spot size, three typical autoignition reaction front modes are identified: supersonic reaction front; detonation development and subsonic reaction front. The effects of initial pressure, initial temperature, fuel type and equivalence ratio on detonation development regime are examined. It is found that the detonation development regime strongly depends on mixture composition (fuel and equivalence ratio) and thermal conditions (initial pressure and temperature). Therefore, to achieve the quantitative prediction of super-knock in engines, we need use the detonation development regime for specific fuel at specific initial temperature, initial pressure, and equivalence ratio.     

Multiplicity of detonation regimes in systems with a multi-peaked thermicity

The study investigates detonations with multiple quasi-steady velocities that have been observed in the past in systems with multi-peaked thermicity, using Fickett's detonation analogue. A steady-state analysis of the travelling wave predicts multiple states, however, all but the one with the highest velocity develop a singularity after the sonic point. Simulations show singularities are associated with a shock wave which overtakes all sonic points, establishing a detonation travelling at the highest of the predicted velocities. Under a certain parameter range, the steady-state detonation can have multiple sonic points and solutions. Embedded shocks can exist behind sonic points, where they link the weak and strong solutions. Sonic points whose characteristics do not diverge are found to be unstable, and to be the source of the embedded shocks. Numerical simulations show that these shocks are only quasi-stable. This is believed to be due in part to a feature of the model which permits shocks anywhere behind a sonic point.     

Modelling detonation of heterogeneous explosives with embedded inert particles using detonation shock dynamics: Normal and divergent propagation in regular and simplified microstructure

This paper discusses the mathematical formulation of Detonation Shock Dynamics (DSD) regarding a detonation shock wave passing over a series of inert spherical particles embedded in a high-explosive material. DSD provides an efficient method for studying detonation front propagation in such materials without the necessity of simulating the combustion equations for the entire system. We derive a series of partial differential equations in a cylindrical coordinate system and a moving shock-attached coordinate system which describes the propagation of detonation about a single particle, where the detonation obeys a linear shock normal velocity-curvature (D_n–κ) DSD relation. We solve these equations numerically and observe the short-term and long-term behaviour of the detonation shock wave as it passes over the particles. We discuss the shape of the perturbed shock wave and demonstrate the periodic and convergent behaviour obtained when detonation passes over a regular, periodic array of inert spherical particles.     

连续波氟化氘/氟化氢化学激光器TRIP技术混合机理

采用CLNST引入的扰流气注入技术对连续波氟化氘/氟化氢化学激光器进行数值模拟研究。模拟其3维流场有化学反应,但不考虑激射的情况。计算结果显示:加入扰流气后,平行注入的主副气流混合增强明显,气流混合速度明显加快,在喷管出口平面5 mm处即扩散至副喷管中心线位置,且小信号增益系数值明显增强。分析了流场结构得出:扰流气增强气流混合的主要方式在于使反应界面拉伸和扭曲,同时观察到有弱漩涡产生,但沿气流方向较快耗散掉,对流场影响较小。     

Embedded particle size distribution and its effect on detonation in composite explosives

This paper discusses the effects of stochastically varying inert particle parameters on the long-term behaviour of detonation front propagation. The simulation model involves a series of cylindrical high explosive unit cells, each embedded with an inert spherical particle. Detonation shock dynamics theory postulates that the velocity of the shock front in the explosive fluid is related to its curvature. In our previous work, we derived a series of partial differential equations that govern the propagation of the shock front passing over the inert particles and developed a computationally efficient simulation environment to study the model over extremely long timescales. We expand upon that project by randomising several properties of the inert particles to represent experimental designs better. First, we randomise the particle diameters according to the Weibull distribution. Then we discuss stochastic particle spacing methods and their effects on the predictability of the shock wave speed. Finally, we discuss mixtures of plastic and metal particles and material inconsistency among the particles.     

Numerical investigation of laminar cross-flow non-premixed flames in the presence of a bluff-body

A knowledge of flame stability regimes in the presence of cylindrical bluff-bodies of various dimensions is essential to design non-premixed burners. The reacting flow field in such cases is reported to be three-dimensional and unsteady. In the literature, only a few experimental investigations with limited measurements are available. Therefore, in this work, a detailed numerical study of laminar cross-flow non-premixed methane–air flames in the presence of a square cylinder is presented. The flow, temperature, species and reaction fields have been predicted using a comprehensive transient three-dimensional reacting flow model with detailed chemical kinetics and variable thermo-physical properties, in order to get a good insight into the flame stabilisation phenomena. Further, analyses of quantities such as local equivalence ratio, cell Damköhler number, species velocity, net consumption rate of methane, which are not easily obtained through experiments even with detailed diagnostics, have been carried out. The influence of the flow field due to varying inlet velocity of the oxidiser, in the presence of the bluff-body, on flame anchoring location has been analysed in detail. Local equivalence ratio contours obtained from non-reacting flow calculations are seen to be quite useful in analysing the mixing process and in the prediction of flame anchoring locations when the flames are not separated. Cell Damköhler number has been calculated using cell size, species velocity of the fuel, which is a derived quantity, and the net reaction rate of the fuel. The flame zone, which is customarily inferred from the contours of temperature, CO and OH, is also shown to be predicted well by the contour line corresponding to a Damköhler number equal to unity. The net reaction rate of CH₄ and the net rates of two dominant reactions, which consume methane, show clearly the variation in the flame anchoring locations in these three cases. Further, the three-dimensionality of these flames are analysed by plotting the mean temperature contours in y–z planes. Finally, the unsteadiness in the separated flame case is analysed.     

Particle path tracking method in two- and three-dimensional continuously rotating detonation engines

The particle path tracking method is proposed and used in two-dimensional（2D） and three-dimensional（3D） numerical simulations of continuously rotating detonation engines（CRDEs）. This method is used to analyze the combustion and expansion processes of the fresh particles, and the thermodynamic cycle process of CRDE. In a 3D CRDE flow field, as the radius of the annulus increases, the no-injection area proportion increases, the non-detonation proportion decreases, and the detonation height decreases. The flow field parameters on the 3D mid annulus are different from in the 2D flow field under the same chamber size. The non-detonation proportion in the 3D flow field is less than in the 2D flow field. In the 2D and 3D CRDE, the paths of the flow particles have only a small fluctuation in the circumferential direction. The numerical thermodynamic cycle processes are qualitatively consistent with the three ideal cycle models, and they are right in between the ideal F–J cycle and ideal ZND cycle. The net mechanical work and thermal efficiency are slightly smaller in the 2D simulation than in the 3D simulation. In the 3D CRDE, as the radius of the annulus increases, the net mechanical work is almost constant, and the thermal efficiency increases. The numerical thermal efficiencies are larger than F–J cycle, and much smaller than ZND cycle.