In the present paper, we proposed a two-dimensional finite difference method (FDM) of characteristic lines to address problems of the non-isentropic steady flow of cylindrical explosive underwater explosion. This method describes the non-isentropic effect by adding an entropy-related variable along the flow line to the pressure-related equation along the Mach line, so that both the isentropic flow and the non-isentropic flow can be described in the same equations of the characteristics. Based on the features of the near-field shock wave we firstly modeled the underwater explosion with an infinitely long cylindrical explosive, then discretized those equations using this finite difference method and constructed an appropriate grid to ensure the numerical convergence, and finally calculated the underwater near-field shock wave for several explosives by programming. The numerical examples showed that the results of this method are consistent with those of the commercial finite element software AUTODYN and those of experiments, suggesting that the FDM of characteristics can capture the shock wave front with relatively high accuracy, and confirming that this method is applicable to solving problems in cylindrical explosive underwater explosion.%针对二维定常可压缩超声速非等熵柱状流,提出一种特征线差分解法,通过在沿马赫线的相容方程中添加沿流线的熵变项以描述非等熵效应,得到等熵流和非等熵流均适用的三族特征线方程组.根据水下爆炸近场特点,建立无限长柱状装药的定常模型,将三族特征线方程组用有限差分法离散求解,通过构造合适的网格保证计算格式可以数值上收敛,由此编制程序并计算几种柱状炸药的水下爆炸近场冲击波.对比有限元模拟结果和实验结果发现,特征线差分法可以比较准确地捕捉冲击波形状并计算冲击波后流场,从而验证了所提出的三族特征线差分法的准确性.
展开▼