钝感炸药爆轰波传播和驱动的数值模拟研究  

Numerical simulation research on propagation and driven of non-ideal detonation wave of insensitive high explosive

作　　者：姜洋[1] 钟敏[1] 谭多望[1] 李平[1] 柏劲松[1] 

机构地区：[1]中国工程物理研究院流体物理研究所,绵阳621900

出　　处：《计算力学学报》2011年第B04期171-175,共5页Chinese Journal of Computational Mechanics

基　　金：国家自然科学基金资助项目(10972045)

摘　　要：为实现对非理想爆轰波传播和爆轰驱动问题的精确计算,研究了与贴体坐标下Level Set方法相结合的爆轰冲击波动力学(DSD)计算方法;根据Hamilton-Jacobi方程的Godunov差分格式,提出了非正交的贴体坐标系中Level Set函数方程的差分格式及其相应的数值方法;并将该DSD方法与二维流体动力学程序TDY2D耦合后得到了编码TDY DSD。分别运用DSD和TDYDSD程序对钝感炸药JB-9104爆轰波传播和驱动飞片模型进行了数值模拟计算。计算结果表明了方法和程序的正确性与有效性,能够对复杂二维构型中非理想爆轰问题进行高精度数值模拟。In order to compute the problem accurately about the propagation of non-ideal detonation wave and the acceleration by detonation products, the Detonation Shock Dynamics （DSD） method which was combined with the Level Set method in body-fitted coordinate was studied in this paper. According to the Hamilton-Jacobi formulation of the Godunov＇s scheme, the finite difference method and algorithm for the Level Set equation in non-orthogonal body-fitted coordinate were researched. Then the corresponding code DSD was incorporated into the 2D hydrodynamic code TDY2D to form a code TDY DSD. The code DSD and TDY DSD were used to compute these models about the propagation of detonation wave and the acceleration by detonation products of insensitive high explosive JB-9014. The computation results show that these codes are correct and effective. The precision is high in calculating the problem about non-ideal detonation in complex configuration.

关 键 词：非理想爆轰 DSD方法 贴体坐标系 数值模拟 

分 类 号：O241[理学—计算数学] O382[理学—数学]