# Simulation of combustion-to-explosion transition processes in Octogen-based explosive

### Abstract

*LX-10*, enclosed in a solid assembly. The calculations were carried out in the

*ANSYS Autodyn*finite element analysis program using the Powder Burn combustion model. The satisfactory agreement between the calculated and experimental data was obtained.

#### For citation:

Zuev Yu.S., Karmanov N.I. Simulation of combustion-to-explosion transition processes in Octogen-based explosive. *Journal of «Almaz – Antey» Air and Space Defence Corporation*. 2017;(4):73-79.
https://doi.org/10.38013/2542-0542-2017-4-73-79

## Introduction

During their lifetime cycle, products containing explosives may be exposed to various unanticipated emergency impacts (falling, fire action, bullet impact, etc.). The most frequently occurring emergency is fire. It is known that at explosive ignition, various outcomes can be possible: product containment failure and explosive burning out, convective combustion, low-velocity explosive transformation, or detonation [1].

The most severe consequences of accidents are associated with detonation. In this case, high-velocity fragments are produced, which may cause detonation in products located nearby. At the same time, in case of ignition of octogen-based explosives, a low-velocity explosive transformation scenario can be possible. In this case, larger fragments will be formed, but flying at a lower velocity than fragments produced by detonation. In view of this, it is vital to be able to predict product reaction to the impact of thermal field generated by fire.

The product reaction can be predicted using advanced finite-element analysis software applications, such as *ANSYS Autodyn*. The *ANSYS Autodyn* software features a fairly flexible combustion model *Powder Burn*, suitable for describing combustion processes in a wide range of pressures typical for transient conditions of explosive transformation (0.002...0.030 Mbar).

## Powder Burn model structure

The model includes equations for the state of solids and gaseous products of decomposition, as well as for solid burnout rate [2]. A solid specimen is presumed to consist of particles of certain shape and corresponding characteristic size r_{0}.

The equation for the state of gaseous products of decomposition can be set in the form of the *JWL* equation. To describe the behaviour of a solid, the compaction equation can be used, for which it is necessary to set two tabular functions: the first one is the relationship between solid-phase pressure vs. density, and the second – relationship between sound velocity vs. density. The burnout rate is calculated by the formula

where F – burnout magnitude (varying from zero to unity);

t – time;

H (P_{g}) – tabular dependence of particle layerby-layer combustion rate H vs. gas-phase pressure P_{g} .

Here also, constant G characterises particle size. For example, for a spherical particle G = 3/r_{0} (r_{0} – radius), for a cubical one – G = 1 / r_{0} (r_{0} – rib length). Constants c characterise particle shape. For example, c = 2/3 for a spherical or cubical particle, c = 1/2 – for a cylindrical one.

## Experimental data

Using the *Powder Burn* model, we shall try to reproduce the experiment results [3] with regard to the combustion-to-explosion transition in high-density octogen-containing LX-10 explosive samples confined in a solid assembly.

The experiment setup schematic is shown in Fig. 1. Five cylindrical samples of the LX-10 explosive (octogen – 95 %, Viton binder – 5 %) were placed in a steel tube with inside diameter of 44.93 mm, length of 203.2 mm, and wall thickness of 2.93 mm. One of the explosive samples contained an igniter.

**Fig. 1**. Experiment setup [3]: 1 – igniter; 2 – LX-10 explosive samples; 3 – AerMet 100 steel tube; 4 – cover

After igniter actuation, an explosion occurred. The recorded velocities V of tube and fragments separation in the plane passing through the igniter and perpendicular to the assembly axis are shown in Fig. 2. The velocity increased for about 250 µs (see Fig. 2), which is a proof of non-detonative character of the explosion. A substantial increase in tube velocity, from 10 m/s to 400...700 m/s (before fragmentation), occurred during 30...41 µs. The values of fragments velocity was changing within a range of 800…2100 m/s.

**Fig. 2**. Steel tube and fragment separation velocities measured using Photonic Doppler Velocimetry (PDV) and radars, respectively [3]:

1 – PDV3; 2 – PDV2; 3 – PDV1; 4 – Radar1; 5 – Radar2; 6 – Radar3

## Initial data for calculations

The purpose of the study was to describe the dynamics of steel tube separation based in the experiment [3], using the *Powder Burn* model. For that, the initial data are required, which were taken from papers [3], [4]. For example, the paper [4] provides parameter values for the JWL equation of state both for solid LX-10 explosive and for the gaseous products of LX-10 detonation (Table 1). The JWL equation appears as follows:

where P – pressure;

ρ – density;

E – internal energy;

A, B, R_{1}, R_{2} – parameters;

ρ_{0} – initial density of the LX-10 explosive.

For description of strength characteristics of the LX-10 explosive, a von Mises strength model with the shear modulus equal to 0.05 Mbar and yield strength Y_{0} = 0.002 Mbar was selected.

The parameters for detonation products taken from Table 1 can be used right away, because in the Powder Burn model, the explosive decomposition products are also described by the JWL equation.

Table 1

Parameters for LX-10 explosive (g/cm^{3})

Parameters |
Solid explosive LX-10 |
Detonation products |
---|---|---|

А, Mbar |
7320 |
16,689 |

В, Mbar |
-0,052654 |
0,5969 |

R |
14,10 |
5,900 |

R |
1,410 |
2,100 |

w |
0,8867 |
0,450 |

E |
4,40110 |
0,104 |

For solid behaviour description, the *Powder Burn* model employs the compaction equation. This equation was originally intended for modelling the behaviour of porous substances, such as loose sand or gunpowder of bulk density. It is known that substances in the bulk state become considerably compacted under load, and with load removed, their density drops, although not back to the initial state but to a certain value which is greater than the initial density. The compaction equation shall be adapted to address the problem under consideration, in which the initially high-density low-porosity explosive samples are used. In order to determine the compaction equation parameters, we shall use the JWL equation parameters for a solid explosive LX-10 (see Table 1).

For the compaction equation, it is necessary to set two tabular functions: function P(ρ), which defines the dependence of solid-phase pressure vs. density in the process of loading; function c (ρ^{(A)}), which describes unloading of earlier loaded substance down to a certain pressure and defines the dependence of sound velocity vs. density of the unloaded substance. To calculate the values of tabular functions included in the compaction equation, a curve for the JWL equation should be plotted for the case of blast-shock wave loading with the parameters taken from Table 1 (Fig. 3).

The values of tabular function P(ρ) can be select right away. The curve in Fig. 3 is obtained from the JWL equation for solid explosive, with the points corresponding to tabular function P(ρ). Let us take ten points belonging to the curve that corresponds to the JWL equation. The values of this tabular functions that defines loading of the explosive are given in Table 2.

Table 2

Tabular function values

№ |
P, Mbar |
ρ, g/cm |
№ |
P, Mbar |
ρ, g/cm |
---|---|---|---|---|---|

1 |
0 |
1,865 |
6 |
0,043 |
2,15 |

2 |
2,848·10 |
1,900 |
7 |
0,094 |
2,29 |

3 |
6,888·10 |
1,940 |
8 |
0,169 |
2,40 |

4 |
0,012 |
1,980 |
9 |
0,248 |
2,48 |

5 |
0,019 |
2,030 |
10 |
0,539 |
2,66 |

To fill in tabular function c(ρ^{(A)}), we can proceed as shown in Fig. 3, where the curve corresponds to the JWL equation for a solid explosive, and the points correspond to tabular function P(ρ) (see Table 2). The line passing, for example, through points (ρ_{4} ; P_{4}) and (ρ_{5} ; P_{5}), will intersect the abscissa axis in point (ρ_{4}^{(A)}; P = 0). The tangent of the slope angle of this straight line to the abscissa axis is equal to sound velocity square c_{4} . Unloading of the explosive, initially loaded up to the state (ρ_{5} ; P_{5}), will be occurring along this line. As a result, a pair of values (ρ_{4}^{(A)};c_{4}) is obtained. The rest of the pairs of values can be obtained in a similar way. Table 3 contains the values of tabular function с(ρ^{(A)}), defining explosive unloading.

Table 3

Values of tabular function с(ρ^{(с)})

№ |
c, cm/µs |
р |
№ |
c, cm/µs |
р |
---|---|---|---|---|---|

1 |
0,2870 |
1,8650 |
6 |
0,6036 |
2,0313 |

2 |
0,3178 |
1,8711 |
7 |
0,8257 |
2,1514 |

3 |
0,3575 |
1,8854 |
8 |
0,9937 |
2,2282 |

4 |
0,3742 |
1,8936 |
9 |
1,2715 |
2,3259 |

5 |
0,4472 |
1,9343 |
10 |
1,4718 |
2,4105 |

For steel AerMet 100, the Grüneisen equation of state was used:

where Γ_{0} = 1.84;

ρ_{0} = 7.923 g/cm^{3};

c_{0} = 0.4534 cm/µs;

S_{1} = 1.5 [3].

We introduce the Steinberg – Guinan strength model as well:

where Y – dynamic yield strength;

– equivalent plastic strain;

G_{сд} – dynamic shear modulus;

T – temperature.

The parameters for steel AerMet 100 [3] are as follows. At temperature T_{0} = 300 K and pressure P = 0, the initial shear modulus G_{0} = 0.748 Mbar. The initial yield strength Y_{0} = 0.01016 Mbar at temperature T_{0} = 300 K and pressure P = 0. The maximum yield strength Y_{max} = 0.03 Mbar, β = 2, b = 0, n = 0.5, h = 0.

It is known that for quantitative description of the process of a relatively strong shock wave transformation into a detonation wave in high-density explosives, the *Ignition and Growth* model [4] is used, in which the burnout equation (1) can be written as follows:

Here, ρ – current density;

ρ_{0} – initial explosive density;

I, G_{1}, G_{2}, a, b, c, d, e, g, x, y, z, F_{ig max}, F_{G1max}, F_{G2min} – parameters.

This model is represented in the ANSYS Autodyn and LS-DYNA software applications. For the LX-10 explosive, it is reasonable to use it in pressure ranges of 0.015…0.5 Mbar. To apply the *Powder Burn* model in this pressure range, it is necessary to adjust its parameters to comply with the *Ignition and Growth* model parameters, whose values are given in Table 4 [4].

Table 4

Kinetic parameters of LX-10 for the Ignition and Growth model [4]

a = 0 |
G |
G |
I = 7,43 · 10 |
---|---|---|---|

b = 0,667 |
y = 2 |
z = 2 |
F |

x = 20 |
c = 0,667 |
e = 0,333 |
F |

d = 0,333 |
g = 1 |
F |

In the *Powder Burn* model, combustion rate H(P_{g}) is set by the tabular method, therefore, for each value of pressure P_{g} magnitude it is possible to set such a value of combustion rate H, which will allow to describe on the average the behaviour of curve *dF/dt*, specified by the equation (2). For that, given a specific pressure value, we shall find the mean value of function *dF/dt* in the interval of burnout values F = 0–1. The mean value of function *dF/dt* is computed by the formula

Let us assume that in equation (1) G = 1, and c = 0.001. Then,

The value of parameter c = 0.001 corresponds, for example, to cubical particles burning on one of the six surfaces only. Coefficient G is taken equal to unity for the convenience of determining H(P_{g}). Its actual values for the LX-10 explosive lie within a range of 100…200 1/cm (explosive particle size r_{0} ~ 1/G = 0.005...0.01 cm). Thus, values of table function will contain information on parameter G . Now we can select the desired number of typical pressure values and calculate parameter for each value selected (Table 5).

Table 5

Values of tabular function H(P_{g}) for high pressures (P_{g} ≥ 0,015 Mbar)

№ |
1 |
2 |
3 |
4 |
5 |
6 |
---|---|---|---|---|---|---|

P |
0,015 |
0,030 |
0,050 |
0,080 |
0,120 |
0,450 |

H, µs |
0,019 |
0,077 |
0,216 |
0,613 |
2,792 |
1301 |

Then we select the values of function H(P_{g}) for pressures P_{g} < 0.01 Mbar. The layer-by-layer combustion rates, measured at P_{g} = 0.03 Mbar and P_{g} = 0.12 Mbar, are equal to V = 7∙10^{-4} cm/µs and V = 0.03 cm/µs [5]. Multiplying the combustion rate values by 100, we have 0.07 and 3, respectively, which is close to the values of functions H (0.03 Mbar) = 0.077 µs^{–1} and H (0.12 Mbar) = 2.792 µs^{–1} (see Table 5).

In a formal way, it turns out that the equation (1) must have G ≈ 100. Such a value corresponds to explosive particle size r_{0} ~ 1/G = 0.01 cm and to real sizes of octogen particles. Let us assume that G = 100 for pressures P_{g} < 0.01 Mbar as well. Then we multiply layerby-layer combustion rate values 2∙10^{–5} cm/µs, 5.5∙10^{5} cm/µs, and 1∙10^{–4} cm/µs [5] by 100 at the values of P_{g} equal to 0.0016 Mbar, 0.004 Mbar, and 0.007 Mbar. In so doing, we take the constant G = 1. The values of function H(P_{g}) for pressures P_{g} < 0.01 Mbar (Table 6) will be obtained.

Table 6

Values of tabular function H(P_{g}) for low pressures (P_{g} < 0,01 Mbar)

№ |
1 |
2 |
3 |
4 |
5 |
---|---|---|---|---|---|

P |
10 |
0,0001 |
0,0016 |
0,0040 |
0,007 |

H, µs |
10 |
10 |
0,0020 |
0,0055 |
0,010 |

The first two values of the tabular function are selected such that in the calculation, at the initial moments of time when the pressure is low, there is a weak reaction within the entire explosive. Then, as compression waves or shock waves pass through the explosive, the reaction rate will increase by orders of magnitude.

## Computation results

The computational geometry is shown in Fig. 4. The grid size is 0.1 cm. Computations were conducted in the axisymmetric ALE formulation using the ANSYS Autodyn software (licence No. 774-2013-ША dated 01.10.2013). Ignition would start in several elements in which the initial pressure of 0.002 Mbar was applied for 30 µs. The impact of the igniter, as well as of the channels in the explosive samples in which wires for the igniter were laid, was not considered in the computations.

**Fig. 4**. Computation setup: 1 – cover of *AerMet* 100; 2 – tube of *AerMet* 100; 3 – *LX*-10 explosive; 4 – ignition area

As a result of computation using function H(P_{g}) (see Tables 5, 6), and at G = 1, it was determined that tube separation velocity was increasing from 10 to 100 m/s during 16 µs and from 100 to 700 m/s during 15 µs, with transition to detonation occurring quite early. In the experiment, velocity was increasing from 10 to 100 m/s and from 100 m/s to 700 ms during 19 µs. If G = 0.7, then tube separation velocity increases from 10 to 100 m/s during 26 µs and from 100 to 700 m/s during 21 µs, with detonation occurring closer to the cover. If G = 0.5, then tube separation velocity increases from 10 to 100 m/s during 30 µs and from 100 to 700 m/s during 36 µs with no transition to detonation occurring.

Computations at G = 0.7 are the most consistent with the experiment. The results of this computation are shown in Fig. 5. A comparison between the computational and experimental data is given in Fig. 6.

**Fig. 5**. Computation result at the moment of velocity reaching V = 700 m/s

**Fig. 6**. Comparison between computational and experimental tube motion velocities:

1 – PDV3; 2 – PDV2; 3 – PDV1; 4 – computation data; 5 – Radar1; 6 – Radar2; 7 – Radar3

The computational and experimental displacement of the tube wall is ~1 cm. The maximum computed tube velocity is ~2100 m/s, which corresponds to the maximum velocities of fragments recorded in the experiment. According to [3], destruction of the tube starts at 15 % radial strain. Under such deformation, tube velocity in the computation reached 500 m/s, which corresponds to experimental velocities of 400...700 m/s under which the fragmentation starts.

## Conclusion

Using the *ANSYS Autodyn* finite-element analysis software and its integrated *Powder Burn* combustion model, we managed to satisfactorily describe the results of an experiment for simulating the combustion-to-exposition transition in *LX*-10 explosive samples enclosed in a solid assembly [3]. In the *Powder Burn* model, the burnout rate is selected by the tabular method, therefore a possibility appears to describe energy release for a wide range of pressures (P_{g} = 0.002...0.5 Mbar) where transition from combustion to detonation may occur.

## References

1. Переход горения конденсированных систем во взрыв / А. Ф. Беляев, В. К. Боболев, А. И. Коротков и др.; под ред. П. Ф. Похила. М.: Наука, 1973. 293 с.

2. Atwood А., Friis Е. K., Moxnes J. F. А Mathematical model for combustion of energetic powder materials // 34th International Annual Conference of ICT. Karlsruhe, Federal Republic of Germany, 24–27 June, 2003. Pp. 1–14.

3. McClelland K. М. A., Maienschein J. L., Howard W. M., Nichols A. L., de Haven M. R., Strand O. T. Measurement and ALE3D simulation of violence in a deflagration experiment with LX10 and Aermet-100 alloy // 13th International Detonation Symposium. Norfolk, VA, United States 23–28 July, 2006. Pp. 1–12.

4. Tarver C. M., Urtiew P. A., Chidester S. K., Green L. G. Shock compression and initiation of LX-10 // Propellants, Explosive, Pyrotechnics. 1993. Vol. 18. Pp. 117–127.

5. Nichols A. L., Tarver С. M., McGuire Е. М. ALE3D statistical hot spot model results for LX-17 // Proceedings of Conference of the American Physical Society Topical Conference on Shock Compression of Condensed Matter. Portland, Oregon, July 20–25, 2003. Pp. 1–6.

### About the Authors

**Yu. S. Zuev**Russian Federation

Zuev Yuriy Semenovich – Candidate of Engineering Sciences, Senior Research Fellow, Head of the Department of Technical Mechanics. Science research interests: aircraft and machines dynamics and strength.

Snezhinsk

**N. I. Karmanov**Russian Federation

Karmanov Nikolay Ivanovich – post-graduate student. Science research interests: simulation of combustion-to-explosion transition processes in an explosive.

Snezhinsk

### Review

#### For citation:

Zuev Yu.S., Karmanov N.I. Simulation of combustion-to-explosion transition processes in Octogen-based explosive. *Journal of «Almaz – Antey» Air and Space Defence Corporation*. 2017;(4):73-79.
https://doi.org/10.38013/2542-0542-2017-4-73-79