基于详细氢氧化学动力学模型,建立了描述氢氧爆轰的多组分反应欧拉方程组.针对建立的反应欧拉方程组,数值方法上采用3阶Additive Runge-Kutta方法对时间项进行积分,采用5阶精度的加权本质无振荡(WENO)格式对空间对流项进行离散,自主研发了大规模高精度计算程序.该程序能够处理化学反应源项引起的刚性问题,且能节省计算时间和计算内存.对半球型、半椭球型、圆锥型3种结构形式凹面腔内的激波聚焦起爆过程进行了数值模拟,数值模拟研究得到了不同结构形式凹面腔内的激波聚焦起爆过程.
Based on detailed hydrogen and oxygen chemical kinetics model,the multi-species reactive Euler equations were established to describe the hydrogen and oxygen detonation.Using third-order Additive Runge-Kutta methods in time discretization,using fifth-order weighted essentially non-oscillatory(WENO)scheme in spatial discretization,a high-solution large-scale program was developed.It was verified that the program could solve stiff source term caused by multispecies chemical reaction, which could also save computation time and computation memory.With this program,the process of detonation initiation with convergence of shock waves was numerically simulated for 3types of concave cavity,hemispherical concave cavity,semi-ellipsoidal concave cavity and conical concave cavity.The simulation results show the characteristic of detonation initiation in different types of concave cavity.