Search

The Jones-Wilkins-Lee (JWL) equation of state (EOS) is the most popular equation of state used in hydrocodes to model detonation products thermodynamics resulting from high explosives. However, the JWL EOS presents several major difficulties. Namely, its range of validity and convexity is limited as its parameters are adjusted on a given reference curve. When used in conditions where the thermodynamic state varies significantly, computational failure may happen. It occurs frequently in multiphase flow computations when heat and mass transfers with other phases are present, as the resulting thermodynamic state is far from the reference curve used for adjustment. Moreover, the detonation products composition is absent in the JWL formulation. This is another important limitation when additional physics are involved, such as post-combustion, phase transition, or when dealing with non-ideal explosives. In this work, the Becker-Kistiakowsky-Wilson (BKW) EOS is considered instead. The BKW EOS is widely used in thermochemical codes, as parent EOS to provide the reference data used to fit JWL EOS parameters. The BKW EOS present several advantages. Its formulation considers variable gas species in the detonation products. No specific adjustment is needed for a given explosive, and only the detonation products composition requires computational efforts. Its range of validity is large, as the various parameters have been adjusted for each gas species separately. However, the BKW EOS also present difficulties which limit its application in hydrocodes. These difficulties are the aim of the present work. The BKW EOS is currently restricted to thermochemical codes as its mathematical formulation is computationally expensive. It also requires the gas chemical composition to be computed via an appropriate chemical equilibrium solver. In the present work, the computational cost of the BKW EOS is reduced thanks to the thermodynamic relaxation method of Neron et al. (2023). The variable composition of the detonation products is accounted for through specific relaxation terms computed once, with the help of conventional thermochemical codes. The condensed species, such as solid carbon, are modelled with the Cochran-Chan (CC) EOS, having a simpler thermodynamic formulation and being thermodynamically consistent compared to usual EOSs used for condensed species. This multiphase formulation made of BKW gas mixture with variable composition and condensed phases is considered in temperature and pressure equilibrium. It is embedded in a multiphase flow formulation, with the help of the thermodynamic relaxation method mentioned above. It shows comparable computational cost as the JWL formulation, with significantly enhanced physics capabilities and improved robustness., Comment: Journal of Computational Physics, 2024

The Noble-Abel (NA) equation of state (EOS) is widely used in interior ballistics of guns as well as rocket propulsion computations. Its simplicity and accuracy are key points for intensive computations with hyperbolic two-phase flow models considered in interior ballistics codes. An alternative is examined in the present contribution through a first-order virial (VO1) equation of state. Appropriate methods for the determination of related parameters, such as specific gas constant, covolume and condensed material energy for both formulations (NA and VO1) are presented. Combination of closed bomb vessel experiments and thermochemical code computations are needed. An extended VO1 EOS with temperature dependent specific heat is examined. Then extension to multiple reactive materials is addressed. Examples are examined for each formulation (NA and VO1) and comparisons are done with the Becker-Kistiakowsky-Wilson (BKW) EOS as reference. Several conclusions emerged. First, consideration of specific heat temperature dependance in interior ballistics of guns computations appeared insignificant. Second, VO1 appeared more accurate than NA, particularly when gas density comes out of the range used for the EOS parameters determination. Last, regarding mixtures of condensed reactive materials, producing burnt gas mixtures, NA appeared again less accurate than VO1. However, its formulation is explicit, while VO1 requires numerical solving of a non-linear equation, with consequences on computational cost.