一维粉尘爆轰波结构的数值计算

NUMERICAL CALCULATION OF THE STRUCTURE OF 1-D DUST DETONATION

本文发展了P.Wolanski等人建立的关于一维自持粉尘-空气混合物的爆轰模型,考虑了气、固两相通过击波以后的温度弛豫效应。作者根据文献中所建立的Chapman-Jouguet条件的表述:“在定常爆轰的C-J平面上,我们有M=1和dM=0”,给出了一种确定自持爆速的有效算法。这一方法具有独立于方程组的积分方式,因而优于至今人们给出的其他同类算法。用改进后的模型和算法,作者对粉尘爆轰波结构的基本特征和参数影响作了初步探讨。文中一些理论结果还与国外同类实验结果进行了对比。

The model proposed by wolanski et al.[1] has been developed and a new method for determining the detonation velocity is given here. We considered a model for the reaction zone of dust detonation taking into account the two-phase flow effects and wall losses. The difference from the previous model is that the thermal non-equilibrium between particles and the gas has been considered. For the treatment of C-J condition, using a new result developed in [2] recently(at the C-J plan we have M=l and dM=0), we give an efficient procedure as follous1 ) Assume a value of the detonation velocity D.2 ) Integrate the governing equations starting at the shock.3 ) As M→1 , dM will tend to 0 .4 ) If M is monotony in the neighborhood of M=l then choose a higher value of D and repeat the procedure until M just becomes non-monotony in the neighborhood ofM=l.With the above procedure we have obtained very satisfactory numerical results such as the pressure profiles the temperature distribution and the influence of oxygen concentration etc.