Numerical prediction of steady-state detonation properties of condensed-phase explosives

Within the scope of this study, a computer code named BARUT-X has been developed to calculate the detonation properties of C-H-N-O based condensed-phase explosives using the Chapman-Jouguet (C-J) theory. Determination of the detonation properties is performed in chemical equilibrium and steady-state conditions. Unlike other codes in the literature which use steepest descent optimization method, BARUT-X uses a nonlinear optimization code based on Generalized Reduced Gradient algorithm to compute the equilibrium composition of the detonation products. This optimization code provides a higher level of robustness of the solutions and global optimum determination efficiency. The Becker-Kistiakowsky-Wilson's (BKW) equation of state (EOS) is applied to the high-density gaseous detonation products at high pressures. BARUT-X uses RDX, TNT, BKWR, and BKWN set of constants in the BKW EOS. In addition, the Cowan-Fickett's EOS is applied for the compressible solid carbon in the detonation products. The calculated detonation properties for several condensed-phase explosives by BARUT-X have been compared with those computed by EXPLO5 and FORTRAN BKW codes as well as the experimental data in terms of detonation velocity and detonation pressure. Satisfactory agreement is obtained from these comparisons.