引用本文: | 刘世杰,林志勇,孙明波,等.采用不同化学反应源项处理方法的胞格爆轰数值研究[J].国防科技大学学报,2010,32(5):1-6.[点击复制] |
LIU Shijie,LIN Zhiyong,SUN Mingbo,et al.Numerical Simulation of Cellular Detonation Using Different Chemical Reacting Source Term Methods[J].Journal of National University of Defense Technology,2010,32(5):1-6[点击复制] |
|
|
|
本文已被:浏览 7340次 下载 5417次 |
采用不同化学反应源项处理方法的胞格爆轰数值研究 |
刘世杰, 林志勇, 孙明波, 刘卫东 |
(国防科技大学 航天与材料工程学院,湖南 长沙 410073)
|
摘要: |
采用一种改进的化学非平衡流解耦方法对Euler反应流方程进行解耦处理,对流项采用五阶WENO格式进行离散,化学反应源项分别采用梯形公式和α-QSS拟稳态逼近两种方法处理,时间步进采用二阶精度的Runge-Kutta方法,对H2/O2/Ar预混气的胞格爆轰进行了二维数值模拟。两种处理方法都得到了合理的三波点结构,计算结果接近,但α-QSS拟稳态逼近处理方法的化学反应较完全,爆轰波传播速度较高;揭示了爆轰波结构的发展变化过程,横波与横波、横波与壁面碰撞引起的二次起爆对爆轰波的稳定发展至关重要;分析了数值胞格结构,胞格长宽比与参考文献基本一致,验证了本计算方法的有效性。 |
关键词: 胞格爆轰 解耦方法 梯形公式 α-QSS 拟稳态逼近方法 胞格结构 |
DOI: |
投稿日期:2010-03-24 |
基金项目:国家自然科学基金资助项目(90816016) |
|
Numerical Simulation of Cellular Detonation Using Different Chemical Reacting Source Term Methods |
LIU Shijie, LIN Zhiyong, SUN Mingbo, LIU Weidong |
(College of Aerospace and Materials Engineering, National Univ. of Defense Technology, Changsha 410073, China)
|
Abstract: |
An improved de-coupled solver of non-equilibrium flow was used to split the Euler equations for a reacting flow. Inviscid flux was calculated with fifth-order WENO scheme. Trapezoid formula and α-QSS methods were applied to deal with the stiffness generated by the chemical reacting source term of species equations separately. Time integration was performed with two-order TVD Runge-Kutta scheme. Cellular detonation of H2/O2/Ar in a straight tube was numerically studied. Calculated results by different methods were similar, and triple point structures were analyzed. The releasing heat of chemical reaction dealt with by α-QSSmethod is higher than that by trapezoid formula, so the detonation wave propagating speed of α-QSSmethod is higher than that of trapezoid formula too. Development of cellular structure was elaborated. Due to interactions between transverse waves and between transverse wave and tube wall, re-initiation was ignited, which plays a crucial role in the propagating of detonation wave. Ratio of numerical cellular structure length to width is consistent with reference value, which indicates the feasibility of the numerical method. |
Keywords: cellular detonation uncoupled non-equilibrium method trapezoid formula α-QSS cellular structure |
|
|