均匀自由场中球形TNT装药爆炸数值计算分析

为更加精确地描述宽范围内均匀自由场中球形装药爆炸过程的特性,采用球形非等温爆炸模型,考虑装药爆炸近区高温高压条件下真实空气的特殊性质,引入真实空气量纲一拟合状态方程和JWL爆轰产物状态方程,建立拉格朗日球形坐标系下完备的量纲一流体动力学方程组;采用对点源强爆炸自模拟数值解的结果进行样条插值作为初始条件,将人工黏性法与激波装配法相结合,采用显式中心差分格式进行数值求解。通过算例计算,得到冲击波不同时刻的波形曲线和主激波峰值超压衰减曲线,通过拟合得到宽范围内的主激波峰值超压计算公式,计算值与有关文献以及实验较为吻合,验证了公式的准确性。结果表明,采用真实空气状态方程得到的结果适用于描述装药空中爆炸近区和远区流场特性。