Understanding the differences in the shock compression and detonation response of insensitive high explosives (IHEs) and conventional HEs (CHEs) is a long-standing need in HE science and technology. Having previously examined 1,1-diamino-2,2-dinitroethene (FOX-7) IHE single crystals [Winey et al., J. Appl. Phys. 130, 015902 (2021)], the shock and detonation response of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX)—a widely used CHE—was determined using wave profile measurements in ∼250 μm thick single crystals shock compressed to 63 GPa. In marked contrast to FOX-7, RDX single crystals shocked along the [100] and [111] orientations showed wave profile features consistent with chemical decomposition onset at 15 GPa. These features were more pronounced for [100] RDX, suggesting a higher decomposition rate compared to [111] RDX. At 51 GPa and above, flat-topped single waves were observed for both orientations, showing the classic Chapman–Jouget (C–J) detonation response in which the decomposition is completed within the detonation front. The Hugoniot states and sound speeds determined for the detonation products were similar for both [100] and [111] orientations, showing that the overdriven detonation response for RDX does not depend on crystal orientation. The C–J pressure for RDX single crystals (35 GPa)—determined experimentally—is comparable to that of FOX-7. However, compared to FOX-7, chemical decomposition onset for RDX occurs at much lower pressures and the overdriven C–J detonation response occurs at higher pressures. The present findings constitute the first experimental comparison of the shock and detonation response of conventional and insensitive HE single crystals over a broad pressure range below and above the C–J pressure.

Achieving high performance and safety have been key—and sometimes competing—objectives throughout the long history of high explosive (HE) development and their use.1–5 The desire for high performance (detonation wave velocity, pressure, and reaction products) has motivated the development of HE crystals that are denser and more energetic, while enhanced safety (avoiding unwanted detonation) has been the motivation for developing insensitive high explosives (IHEs) that are less sensitive to shock wave initiation than conventional HEs (CHEs). A continuing goal of HE scientific activities is attaining the optimal combination of high performance and insensitivity.

Recently, wave profile measurements were reported6 for 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals shock compressed to 64 GPa—almost twice the calculated Chapman–Jouget (C–J) pressure.7 These wave profiles—the first reported for HE single crystals shocked above the C–J pressure—in thin (∼250 μm) crystals revealed steady overdriven detonations that exhibited classic C–J detonation behavior:8–10 completion of chemical decomposition within the detonation front—an unexpected finding for an IHE. Determination of the detonation products equation of state (EOS) and the C–J pressure (35 GPa) provided6 an experimental demonstration of performance similar to or better than the theoretically calculated response of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX)7—a commonly used CHE. The high performance of FOX-7 relative to other IHEs, such as 2,4,6-triamino-1,3,5-trinitrobenzene (TATB), together with the observed shock insensitivity of FOX-7 to ∼30 GPa, revealed a near-optimal combination of high performance and insensitivity not previously reported for HE crystals.

Gaining insight into the differences in the shock compression and detonation response of IHEs and CHEs is important for guiding research efforts to further enhance the desired combination of high performance and insensitivity. In particular, wave profile measurements for CHE single crystals at high pressures—to complement the previous FOX-7 study6—are needed to better understand the intrinsic response of the HEs, avoiding the complexities inherent to the response of composite (or polymer-bonded) HE formulations. However, the peak pressures reached in previous wave profile measurements for shock compressed CHE single crystals—pentaerythritol tetranitrate (PETN),11 octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX),12 and RDX13–15—have not exceeded 5 GPa (or ∼13 GPa for unsteady waves in PETN).16 Hence, wave profiles in CHE single crystals are needed at pressures comparable to the pressures examined in FOX-7.6 

To address this need, we conducted well-characterized plane shock wave experiments to measure wave profiles in RDX single crystals to 63 GPa. Because the shock response of RDX single crystals can be examined along more than one orientation (unlike FOX-7), wave profiles were measured for [100]- and [111]-oriented RDX. The experiments were undertaken to address the following questions:

  1. How does the shock sensitivity of RDX single crystals compare to that of FOX-7? In particular, at what pressures do shock-induced chemical decomposition features appear in the measured RDX wave profiles?

  2. At what pressures do steady detonations develop within thin (∼250 mm) RDX crystals? Will such detonations exhibit an idealized C–J behavior at high pressures, similar to FOX-7?

  3. What is the response of RDX detonation products (products Hugoniot curve and EOS) at pressures above the C–J state and how does it compare to that of FOX-7?

  4. Does crystal orientation play a significant role in the shock compression and detonation response of RDX single crystals?

The experimental methods used to obtain wave profiles are presented in Sec. II. The experimental results, including the measured wave profiles, are presented in Sec. III. In Sec. IV, the results are analyzed and discussed to address the scientific questions noted above. The main findings of this work are summarized in Sec. V.

The experimental configuration, shown schematically in Fig. 1, and the experimental methods used in this work are similar to those described in Refs. 6 and 17. Hence, only a brief summary of the methods is presented here.

FIG. 1.

Configuration for plate impact experiments on RDX single crystals. Approximate thicknesses of the Al buffers and the RDX samples are shown here; the actual thicknesses are listed in Table I. TPX is a widely used trade name (Mitsui Chemicals, Inc.) for the polymer polymethylpentene.

FIG. 1.

Configuration for plate impact experiments on RDX single crystals. Approximate thicknesses of the Al buffers and the RDX samples are shown here; the actual thicknesses are listed in Table I. TPX is a widely used trade name (Mitsui Chemicals, Inc.) for the polymer polymethylpentene.

Close modal
TABLE I.

Parameters for plate impact experiments on RDX single crystals.

Experiment numberRDX orientationRDX thickness (μm)Impactor/buffer materialsBuffer thickness (μm)Impact velocity (mm/μs)
1 (22-605) [111] 256 ± 2 Cu/Al 823 ± 2 1.690 ± 0.004 
2 (22-601) [100] 259 ± 2 Cu/Al 820 ± 2 1.694 ± 0.001 
3 (22-603) [100] 255 ± 2 Cu/Al 837 ± 2 2.140 ± 0.001 
4 (22-604) [111] 256 ± 2 Cu/Al 862 ± 2 2.148 ± 0.002 
5 (22-2SH13) [111] 258 ± 2 Al/Al 848 ± 2 4.016 ± 0.005 
6 (22-2SH01) [100] 259 ± 2 Al/Al 888 ± 2 4.083 ± 0.003 
7 (22-2SH12) [111] 258 ± 2 Al/Al 852 ± 2 4.560 ± 0.005 
8 (20-2SH36) [100] 254 ± 2 Al/Al 892 ± 2 4.574 ± 0.004 
9 (22-2SH11) [111] 255 ± 2 Al/Al 843 ± 2 5.365 ± 0.006 
10 (20-2SH41) [100] 248 ± 2 Al/Al 892 ± 2 5.422 ± 0.004 
11 (20-2SH54) [100] 257 ± 2 Al/Al 818 ± 2 6.244 ± 0.004 
12 (22-2SH02) [111] 255 ± 2 Al/Al 818 ± 2 6.292 ± 0.008 
Experiment numberRDX orientationRDX thickness (μm)Impactor/buffer materialsBuffer thickness (μm)Impact velocity (mm/μs)
1 (22-605) [111] 256 ± 2 Cu/Al 823 ± 2 1.690 ± 0.004 
2 (22-601) [100] 259 ± 2 Cu/Al 820 ± 2 1.694 ± 0.001 
3 (22-603) [100] 255 ± 2 Cu/Al 837 ± 2 2.140 ± 0.001 
4 (22-604) [111] 256 ± 2 Cu/Al 862 ± 2 2.148 ± 0.002 
5 (22-2SH13) [111] 258 ± 2 Al/Al 848 ± 2 4.016 ± 0.005 
6 (22-2SH01) [100] 259 ± 2 Al/Al 888 ± 2 4.083 ± 0.003 
7 (22-2SH12) [111] 258 ± 2 Al/Al 852 ± 2 4.560 ± 0.005 
8 (20-2SH36) [100] 254 ± 2 Al/Al 892 ± 2 4.574 ± 0.004 
9 (22-2SH11) [111] 255 ± 2 Al/Al 843 ± 2 5.365 ± 0.006 
10 (20-2SH41) [100] 248 ± 2 Al/Al 892 ± 2 5.422 ± 0.004 
11 (20-2SH54) [100] 257 ± 2 Al/Al 818 ± 2 6.244 ± 0.004 
12 (22-2SH02) [111] 255 ± 2 Al/Al 818 ± 2 6.292 ± 0.008 

The RDX single crystals were grown at the Los Alamos National Laboratory (LANL) from stock produced by the Woolwich (HMX-free) process with additional purification before crystallization. The α-phase crystals (Pbca space group, ρ0 = 1.806 g/cm3)18 were oriented along the [100] or [111] crystallographic direction, cut, and provided to us by Dr. Daniel E. Hooks of LANL. Typical lateral dimensions of the prepared samples were 7–8 mm. Each crystal was ground and polished to a thickness of ∼250 μm, comparable to the sample thickness used in the FOX-7 experiments.6 The polished RDX single crystals were bonded to a 1050 Al buffer19 and a [100]-oriented LiF window20 using epoxy. Prior to bonding, an Al mirror was deposited on the sample side of the LiF window.

Plate impact experiments were conducted using the configuration shown in Fig. 1. Using a two-stage gun,21 1050 Al flyers were impacted onto the Al/RDX/LiF target assemblies (a powder gun was used to launch C101 Cu flyers20 for experiments 1–4); Table I shows the measured impact velocities, together with other experimental parameters. A velocity interferometer system (VISAR)22 having 0.7 ns time resolution was used to measure particle velocity histories at the RDX/LiF interface. In addition, velocity histories were recorded at three locations at the Al buffer/RDX interface to measure the impact tilt and to determine the shock wave arrival time.

Twelve plate impact experiments were conducted on RDX single crystals shock compressed to pressures ranging from 15 to 63 GPa. The measured particle velocity histories at the RDX/LiF interface for all 12 experiments are shown in Fig. 2. The measured wave profiles for many of the experiments show a large jump in particle velocity to a near-constant state, followed by a small jump at a later time. The small jump is due to a reflected wave that originates at the RDX/LiF interface and then reverberates between the Al buffer and the LiF window.

FIG. 2.

Wave profiles for shock compressed RDX single crystals (∼250 μm thick) measured at the RDX/LiF interface. The shock wave entered the samples at time zero. To facilitate comparison of the wave profiles, including transit times through the samples, the time axis has been normalized by the sample thickness. The Hugoniot pressure for experiments 7 and 8 was not determined due to the absence of a constant state following the initial jump.

FIG. 2.

Wave profiles for shock compressed RDX single crystals (∼250 μm thick) measured at the RDX/LiF interface. The shock wave entered the samples at time zero. To facilitate comparison of the wave profiles, including transit times through the samples, the time axis has been normalized by the sample thickness. The Hugoniot pressure for experiments 7 and 8 was not determined due to the absence of a constant state following the initial jump.

Close modal

At 15 GPa (experiments 1 and 2), the measured velocity profiles for both RDX orientations show an initial jump followed by a slow increase in particle velocity; this velocity increase is more rapid for [100] RDX than for [111] RDX (Fig. 3). The observed slow velocity increase indicates the onset of shock wave-induced chemical decomposition in RDX single crystals, with the reaction progressing more rapidly for [100] RDX.

FIG. 3.

Measured RDX wave profiles at 15 GPa, together with the 15 GPa FOX-7 wave profile reported previously (Ref. 17). The time axis is normalized by the sample thickness. The RDX wave profiles continue to rise after the initial jump, but the FOX-7 profile does not; the [100] RDX profile rises faster than the [111] profile.

FIG. 3.

Measured RDX wave profiles at 15 GPa, together with the 15 GPa FOX-7 wave profile reported previously (Ref. 17). The time axis is normalized by the sample thickness. The RDX wave profiles continue to rise after the initial jump, but the FOX-7 profile does not; the [100] RDX profile rises faster than the [111] profile.

Close modal

At 20 GPa (experiments 3 and 4), the large initial jump is again followed by a slow increase in velocity, with the velocity increasing more rapidly for [100] RDX than for [111] RDX (Fig. 4). However, for both orientations, the increase in velocity at 20 GPa is somewhat more rapid than at 15 GPa—especially for the [111] orientation—suggesting that the decomposition rate increases with shock pressure.

FIG. 4.

Measured RDX wave profiles at 20 GPa, together with the 21 and 25 GPa FOX-7 wave profiles reported previously (Refs. 6 and 17). The time axis is normalized by the sample thickness. The RDX wave profiles continue to rise after the initial jump, but the FOX-7 profiles do not; the [100] RDX profile rises somewhat faster than the [111] RDX profile.

FIG. 4.

Measured RDX wave profiles at 20 GPa, together with the 21 and 25 GPa FOX-7 wave profiles reported previously (Refs. 6 and 17). The time axis is normalized by the sample thickness. The RDX wave profiles continue to rise after the initial jump, but the FOX-7 profiles do not; the [100] RDX profile rises somewhat faster than the [111] RDX profile.

Close modal

At 31–32 GPa (experiments 5 and 6), approaching the calculated C–J pressure (∼35 GPa),7 the initial jump and slow increase in velocity are followed by a second velocity jump and significant velocity fluctuations (Fig. 5). The second velocity jump occurs at a significantly earlier time than that expected for the reflected reverberation wave noted earlier. Therefore, the second velocity jump for experiments 5 and 6 suggests the onset of a rapid shock-induced chemical reaction that arises subsequent to the initial slow reaction onset. The second jump occurs at a later time for [111] RDX, compared to [100] RDX. The onset of a rapid chemical reaction is also indicated in experiments 7 and 8, where the measured wave profiles show an initial spike followed by a significant reduction in velocity (Fig. 5); however, the measured profiles for [100] RDX and [111] RDX are quite similar. Further discussion of experiments 5–8 is presented in Sec. IV.

FIG. 5.

Measured RDX wave profiles near the calculated Chapman–Jouget pressure (35 GPa, Ref. 7), together with reported FOX-7 wave profiles at similar pressures (Ref. 6). The time axis is normalized by the sample thickness. The Hugoniot pressure for experiments 7 and 8 was not determined due to the absence of a constant state following the initial jump.

FIG. 5.

Measured RDX wave profiles near the calculated Chapman–Jouget pressure (35 GPa, Ref. 7), together with reported FOX-7 wave profiles at similar pressures (Ref. 6). The time axis is normalized by the sample thickness. The Hugoniot pressure for experiments 7 and 8 was not determined due to the absence of a constant state following the initial jump.

Close modal

At higher pressures (51 GPa and above), the measured profiles are flat-topped single waves. These flat-topped profiles show that the shock wave-induced chemical reaction is completed within the shock front, implying reaction times less than the 0.7 ns time resolution of our measurements. Completion of the reaction within the shock front is consistent with the idealized C–J detonation theory.

The analyses conducted on our measured RDX results—and presented in this section—are similar to those carried out previously for FOX-7 single crystals.6,17 However, unlike the previous FOX-7 work, the RDX results presented here encompass two different crystal orientations ([100] and [111]). Therefore, in this section, the results for [100] RDX and [111] RDX will be compared with each other and with the previous results for FOX-7 single crystals.6,17

For each experiment, the shock compression end states (or Hugoniot states) were determined from the measured shock wave velocities and the measured impact velocities using standard impedance matching methods.23 The results are tabulated in Table II and shown in Fig. 6.

FIG. 6.

Measured Hugoniot states for RDX single crystals in the shock velocity–particle velocity plane (a) and in the pressure–volume plane (b). The red and blue solid circles are the measured Hugoniot states for [100] RDX and [111] RDX, respectively. The red solid triangle is the calculated Chapman–Jouguet state for RDX from Ref. 7. The red solid and red dashed curves are the JWL Hugoniot curve and isentrope, respectively, determined from a fit to the measured [100] and [111] RDX products Hugoniot states, together with calculated results at the C–J state and on the expansion isentrope (Ref. 7). The red dotted-dashed curve was determined from a nonlinear shock velocity–particle velocity fit to the measured RDX Hugoniot states below 2.6 mm/μs particle velocity. The black dotted-dashed curve is the calculated Hugoniot curve for unreacted RDX from Ref. 28. The green solid curve is the JWL Hugoniot curve for FOX-7 detonation products (Ref. 6) and the green dotted-dashed curve is the linear shock velocity–particle velocity fit to the measured unreacted FOX-7 Hugoniot states (Ref. 17).

FIG. 6.

Measured Hugoniot states for RDX single crystals in the shock velocity–particle velocity plane (a) and in the pressure–volume plane (b). The red and blue solid circles are the measured Hugoniot states for [100] RDX and [111] RDX, respectively. The red solid triangle is the calculated Chapman–Jouguet state for RDX from Ref. 7. The red solid and red dashed curves are the JWL Hugoniot curve and isentrope, respectively, determined from a fit to the measured [100] and [111] RDX products Hugoniot states, together with calculated results at the C–J state and on the expansion isentrope (Ref. 7). The red dotted-dashed curve was determined from a nonlinear shock velocity–particle velocity fit to the measured RDX Hugoniot states below 2.6 mm/μs particle velocity. The black dotted-dashed curve is the calculated Hugoniot curve for unreacted RDX from Ref. 28. The green solid curve is the JWL Hugoniot curve for FOX-7 detonation products (Ref. 6) and the green dotted-dashed curve is the linear shock velocity–particle velocity fit to the measured unreacted FOX-7 Hugoniot states (Ref. 17).

Close modal
TABLE II.

Experimental results for shock compressed RDX single crystals.

Experiment numberRDX orientationShock velocity (mm/μs)Particle velocity (mm/μs)aPressure (GPa)aV/V0a
[111] 5.46 ± 0.07 1.500 ± 0.008 14.8 ± 0.1 0.725 ± 0.005 
[100] 5.45 ± 0.07 1.504 ± 0.007 14.8 ± 0.1 0.724 ± 0.005 
[100] 6.03 ± 0.09 1.858 ± 0.009 20.2 ± 0.2 0.692 ± 0.006 
[111] 6.04 ± 0.09 1.864 ± 0.009 20.3 ± 0.2 0.691 ± 0.006 
[111] 6.96 ± 0.11 2.464 ± 0.013 31.0 ± 0.3 0.646 ± 0.007 
[100] 7.11 ± 0.11 2.503 ± 0.013 31.7 ± 0.3 0.643 ± 0.007 
[111] 8.67b c c c 
[100] 8.67b c c c 
[111] 9.13 ± 0.18 3.097 ± 0.022 51.1 ± 0.7 0.661 ± 0.009 
10 [100] 9.12 ± 0.18 3.136 ± 0.022 51.6 ± 0.7 0.656 ± 0.009 
11 [100] 9.52 ± 0.19 3.624 ± 0.025 62.3 ± 0.8 0.619 ± 0.010 
12 [111] 9.55 ± 0.19 3.652 ± 0.026 63.0 ± 0.9 0.618 ± 0.010 
Experiment numberRDX orientationShock velocity (mm/μs)Particle velocity (mm/μs)aPressure (GPa)aV/V0a
[111] 5.46 ± 0.07 1.500 ± 0.008 14.8 ± 0.1 0.725 ± 0.005 
[100] 5.45 ± 0.07 1.504 ± 0.007 14.8 ± 0.1 0.724 ± 0.005 
[100] 6.03 ± 0.09 1.858 ± 0.009 20.2 ± 0.2 0.692 ± 0.006 
[111] 6.04 ± 0.09 1.864 ± 0.009 20.3 ± 0.2 0.691 ± 0.006 
[111] 6.96 ± 0.11 2.464 ± 0.013 31.0 ± 0.3 0.646 ± 0.007 
[100] 7.11 ± 0.11 2.503 ± 0.013 31.7 ± 0.3 0.643 ± 0.007 
[111] 8.67b c c c 
[100] 8.67b c c c 
[111] 9.13 ± 0.18 3.097 ± 0.022 51.1 ± 0.7 0.661 ± 0.009 
10 [100] 9.12 ± 0.18 3.136 ± 0.022 51.6 ± 0.7 0.656 ± 0.009 
11 [100] 9.52 ± 0.19 3.624 ± 0.025 62.3 ± 0.8 0.619 ± 0.010 
12 [111] 9.55 ± 0.19 3.652 ± 0.026 63.0 ± 0.9 0.618 ± 0.010 
a

Hugoniot states were determined from the measured shock wave velocities and impact velocities using impedance matching (Ref. 23).

b

Because the wave propagation in experiments 7 and 8 was likely unsteady, the measured shock velocity represents an average value.

c

Hugoniot states were not determined for experiments 7 and 8 because of the absence of a near-constant state following the initial jump.

For particle velocities (up) greater than 2.6 mm/μs, the measured shock wave velocities (Us) shown in Fig. 6(a) are greater than the calculated C–J detonation velocity for RDX (8.76 mm/μs).7 These results, together with the wave profiles presented in Fig. 2, show that the development of steady overdriven detonation waves in RDX single crystals requires less than 250 μm of run distance from the impact surface at sufficiently high pressures, similar to the response reported previously for FOX-7.6 These findings for RDX and FOX-7 are in contrast to the multi-millimeter run distances reported for establishing steady detonations in TATB-based composite IHEs.24,25

For particle velocities less than 2.6 mm/μs, the measured shock velocities are significantly lower than the calculated C–J detonation velocity.7 As shown in Fig. 6(a) (red dotted-dashed curve), the Us–up results for up < 2.6 mm/μs are fitted well using the nonlinear relationship26 

Usc0=1.290.29exp(4.0upc0)+1.51upc0,
(1)

where c0 = 2.525 mm/μs is the ambient bulk sound speed determined using the RDX elastic constants;27 this finding is in contrast to the linear Us–up relationship calculated for unreacted RDX (black dotted-dashed curve).28 In addition, a linear Us–up relationship was also determined previously for unreacted FOX-7 (green dotted-dashed curve).6,17 Although the reason for the nonlinear Us–up relationship measured for RDX is not clear, the measured wave profiles shown in Figs. 3–5 suggest that the Hugoniot states in Fig. 6 likely do not correspond to completely unreacted RDX. For experiments 5 and 6, the second jump in the RDX wave profiles—which contrasts with the large spikes observed for FOX-7 (Fig. 5)— occurs subsequent to the arrival of the shock wave at the LiF window. Therefore, the onset of a rapid chemical reaction—indicated by the second jump—is localized near the RDX/LiF interface and results from the increased pressure and temperature due to the reflection of the shock wave from the LiF window.29 For experiments 7 and 8, the measured shock velocities (not shown in Fig. 6; see Table II) are significantly lower than the calculated C–J detonation velocity. However, unlike the other experiments, these velocities cannot be associated with a nearly flat-top wave.

Figure 6(b) shows the pressure (P)–volume (V/V0) states corresponding to the Us–up results in Fig. 6(a). Due to chemical energy release in the detonation wave front, the PV/V0 Hugoniot states at 51 GPa and above have significantly larger volumes compared to the extrapolated Hugoniot curve for RDX (red dot-dashed curve). The results show that the PV/V0 Hugoniot responses for [100] RDX and [111] RDX are similar over the range of pressures examined here. We note that Hugoniot states were not determined for experiments 7 and 8 because a near-constant state following the initial jump was not observed.

The measured Hugoniot states at 51 GPa and above arise from steady detonations. Therefore, they provide important insight into the detonation products response, including the products equation of state (EOS). We used the well-established Jones–Wilkins–Lee (JWL) EOS9,30–32 to fit the measured RDX Hugoniot states, together with the calculated results for RDX at the C–J state and along the expansion isentrope.7 The JWL EOS is defined by9,30–32

P=A(1ωV0R1V)eR1V/V0+B(1ωV0R2V)eR2V/V0+ωV0VE,
(2)

where P is the pressure; V is the specific volume; E is the internal energy per unit volume; and A, B, R1, and R2 are constants. In Eq. (2), ω is traditionally used to denote the Grüneisen parameter—usually denoted elsewhere by Γ and defined as Γ=V(P/E)V23,33—which is assumed to be constant in the JWL formulation.9,30–32

The Hugoniot curve and the isentrope determined using the JWL EOS are shown in Fig. 6 (red solid and red dashed curves, respectively). The JWL fit presented here results in Γ = 0.51, which is consistent with the values obtained from JWL fits for overdriven detonations in FOX-7 single crystals6 and in composite HEs.34–36 The parameter values obtained from our JWL fit are listed in Table III.

TABLE III.

RDX parameters for the Jones–Wilkins–Lee equation of state (Refs. 9 and 30–32).

A (GPa)B (GPa)R1R2ωE0 (kJ/cm3)
6.918 × 105 671.4 19.52 4.20 0.51 −10.88a 
A (GPa)B (GPa)R1R2ωE0 (kJ/cm3)
6.918 × 105 671.4 19.52 4.20 0.51 −10.88a 
a

Reference 7.

Also shown in Fig. 6 is the Hugoniot curve determined previously for the FOX-7 detonation products (green solid curve).6 The FOX-7 products Hugoniot curve lies above that of RDX in the PV/V0 plane, suggesting higher performance for overdriven detonations.

As discussed in Sec. III, the small jump in particle velocity in the later portions of the measured wave profiles (Fig. 2) is due to a reflected wave that originates at the RDX/LiF interface and then reverberates between the Al buffer and the LiF window. Therefore, for experiments at 51 GPa and above, the measured arrival time of the reverberating wave provides information about the sound speed in the RDX detonation products, which provides an important constraint on the detonation products EOS.

To examine this constraint, we carried out numerical simulations to match the wave profiles—including the small particle velocity jump at later times—measured at 51 GPa and above. The wave profile calculations utilized the well-established finite-difference, artificial viscosity approach37 in a wave propagation code.38 Because the chemical decomposition in RDX shock-compressed to 51 GPa and above is completed within the shock front, no further reaction was considered and the detonation products response was described using the Mie–Grüneisen EOS,23,33

P=PH+ΓV(EEH),
(3)

where PH and EH are the pressure and specific energy, respectively, on the Hugoniot curve. The Grüneisen parameter Γ was assumed constant, consistent with the JWL EOS formulation. To calculate wave profiles, the Hugoniot curve used to determine PH and EH in Eq. (3) incorporated a smooth transition from the RDX Hugoniot curve [Eq. (1)] to the RDX detonation products Hugoniot curve, as shown in Fig. 7. The smooth transition was incorporated for numerical convenience and has no bearing on the findings.

FIG. 7.

Hugoniot curves for RDX. The red solid curve is the JWL Hugoniot curve for the RDX detonation products (Fig. 6) and the red dotted-dashed curve is the RDX Hugoniot curve (Fig. 6). The black dashed curve was used in numerical simulations to calculate RDX wave profiles at 51 GPa and above; the corresponding RDX products Hugoniot states are shown as black circles.

FIG. 7.

Hugoniot curves for RDX. The red solid curve is the JWL Hugoniot curve for the RDX detonation products (Fig. 6) and the red dotted-dashed curve is the RDX Hugoniot curve (Fig. 6). The black dashed curve was used in numerical simulations to calculate RDX wave profiles at 51 GPa and above; the corresponding RDX products Hugoniot states are shown as black circles.

Close modal

The wave profiles calculated using Γ = 0.51 (red curves) are shown in Fig. 8, together with the measured wave profiles (black curves). The calculated profiles provide a good match to the measured shock wave arrival times and the measured particle velocities. However, the calculated arrival time of the small reverberating wave is early compared to the measured profile, indicating that the calculated detonation products sound speed is too high. As shown in Fig. 8, the Mie–Grüneisen EOS with Γ = 1.05 (green curves) provides a somewhat better match to the measured wave profiles. However, the calculated reverberating wave arrival time is too late for experiments 9 and 10 and too early for experiments 11 and 12.

FIG. 8.

Calculated and measured wave profiles for shock compressed RDX. The shock wave entered the samples at time zero. The time axis has been normalized by the sample thickness. The black curves are the measured wave profiles from experiments 9–12. The red curves are wave profiles calculated using the Mie–Grüneisen equation of state with Γ = 0.51; the green curves were calculated similarly, but using Γ = 1.05. The blue dashed curves and blue solid curves were calculated using Γ = 0.95 and Γ = 1.25, respectively. The curves for experiments 10 and 12 have been offset vertically by 0.15 mm/μs for visual clarity.

FIG. 8.

Calculated and measured wave profiles for shock compressed RDX. The shock wave entered the samples at time zero. The time axis has been normalized by the sample thickness. The black curves are the measured wave profiles from experiments 9–12. The red curves are wave profiles calculated using the Mie–Grüneisen equation of state with Γ = 0.51; the green curves were calculated similarly, but using Γ = 1.05. The blue dashed curves and blue solid curves were calculated using Γ = 0.95 and Γ = 1.25, respectively. The curves for experiments 10 and 12 have been offset vertically by 0.15 mm/μs for visual clarity.

Close modal

These results show that the wave profiles at 51 GPa and above cannot all be matched by calculations using a single Γ value, in contrast to previous calculations for FOX-7.6 As shown in Fig. 8, a good match to all of the wave profiles at 51 GPa and above is obtained from calculations using Γ = 0.95 (blue dashed curves) for experiments 9 and 10 and Γ = 1.25 (blue solid curves) for experiments 11 and 12. Both of these Γ values are significantly larger than that obtained from the JWL fit for RDX (Γ = 0.51) and are somewhat smaller than the Γ value obtained previously for FOX-7 (Γ = 1.30).6 

As noted previously,35 the sonic condition8,9 at the C–J point of a steady detonation implies that the shock wave velocity, the Lagrangian sound speed of the detonation products, and the sum of the Eulerian sound speed and the particle velocity are all equal at the C–J state. Similar to the previous FOX-7 study,6 this condition was used, together with sound speed determinations for the detonation products, to determine the C–J state for RDX. The sound speeds in the RDX detonation products were determined using33,35

cE2=V2dPHdV[1Γ2(V0V1)]+VΓ2PH,
(4)

where PH is the Hugoniot pressure, V is the specific volume, dPH/dV is the slope of the Hugoniot curve at PH, and Γ is the Grüneisen parameter. In Eq. (4), the JWL Hugoniot curve for RDX in Fig. 6 was used, together with Γ = 0.95 for experiments 9 and 10 and Γ = 1.25 for experiments 11 and 12. We note that cE in Eq. (4) is the Eulerian sound speed, the acoustic wave velocity in the compressed fluid. The Lagrangian sound speed is defined by

cL=cEV0V.
(5)

As shown in Fig. 9, the shock wave velocity, the Lagrangian sound speed, and the sum of the Eulerian sound speed and the particle velocity all converge at the Us–up state corresponding to a C–J pressure of 35 GPa—providing experimental confirmation of the calculated C–J pressure.7 Thus, the calculated C–J state is shown to be consistent with the measured wave profiles at 51 GPa and above, similar to the previous finding for FOX-7.6 

FIG. 9.

Wave velocities for shock compressed RDX. The black solid circles are the RDX products Hugoniot states at 51 GPa and above. The red solid curve is the JWL Hugoniot curve for the reaction products (Fig. 6) and the red dotted-dashed curve is the RDX Hugoniot curve (Fig. 6). The blue solid circles and blue open circles are the Lagrangian and Eulerian sound speeds, respectively, determined from the RDX reaction products EOS with Γ = 0.95 at 51–52 GPa and Γ = 1.25 at 62–63 GP. The blue solid square is the Chapman–Jouguet state determined using the sound speeds and shock velocities.

FIG. 9.

Wave velocities for shock compressed RDX. The black solid circles are the RDX products Hugoniot states at 51 GPa and above. The red solid curve is the JWL Hugoniot curve for the reaction products (Fig. 6) and the red dotted-dashed curve is the RDX Hugoniot curve (Fig. 6). The blue solid circles and blue open circles are the Lagrangian and Eulerian sound speeds, respectively, determined from the RDX reaction products EOS with Γ = 0.95 at 51–52 GPa and Γ = 1.25 at 62–63 GP. The blue solid square is the Chapman–Jouguet state determined using the sound speeds and shock velocities.

Close modal

The results and analyses presented here for RDX single crystals—a conventional HE—complement those reported previously for FOX-7 IHE single crystals,6,17 providing important insight into the differences in the shock compression and detonation response of IHEs and conventional HEs. The comparisons presented here between the shock compression responses of RDX and FOX-7 are particularly useful because they involve the intrinsic response of conventional HE and IHE single crystals, rather than the significantly more complex response of composite (or polymer-bonded) HE formulations.

The wave profile features indicating the onset of a chemical reaction, observed for RDX at 15 and 20 GPa (Figs. 2–4), demonstrate the considerably higher shock sensitivity of RDX single crystals relative to FOX-7 single crystals, which showed no sign of chemical decomposition to ∼30 GPa (Figs. 3 and 4).6,17 Approaching the C–J pressure, the onset of relatively slow reaction rates in the shock compressed RDX crystals gives way to much faster reactions (Fig. 5). Compared to FOX-7,6 the energy release in RDX is slower and occurs over a considerably wider pressure range.

The flat-topped single waves observed at 51 GPa and above (Fig. 2) show that RDX, like FOX-7, can exhibit steady overdriven detonation waves within 250 μm of run distance from the impact surface at sufficiently high pressures, and these detonation waves are well-described by the classic Chapman–Jouguet (C–J) detonation theory.8–10 Interestingly, the shock pressures required to observe this behavior are higher for RDX than for FOX-7, despite the significantly greater sensitivity of RDX to the onset of chemical decomposition at lower pressures.

Because of the C–J detonation exhibited by the RDX single crystals, the measured wave profiles provided significant insight into the RDX detonation products response. In particular, the detonation products EOS for RDX single crystals—a key determinant of HE performance9,39—was determined to pressures almost twice that of the C–J state, similar to the analysis carried out previously for FOX-7.6 The FOX-7 products Hugoniot curve lies above that of RDX in the PV/V0 plane (Fig. 6), suggesting higher performance for FOX-7 at overdriven detonation conditions.

Our arrival time measurements for the small reverberation wave observed in the measured profiles (Fig. 2) provided additional constraints on the detonation products EOS. Numerical simulations using a single constant Grüneisen parameter were unable to match the measured wave profiles at 51 GPa and above (Fig. 8), in contrast to previous results for FOX-7.6 The Grüneisen parameters used to match the overdriven detonation wave profiles for RDX at 51 GPa (Γ = 0.95) and at 63 GPa (Γ = 1.25), while somewhat smaller than that used in similar calculations for FOX-7,6 were significantly larger than those required to match the response of detonation products for RDX, FOX-7,6 and other HEs31,34–36 at relatively large expansions (V/V0 > 3). Therefore, our results for RDX—together with the previous results for FOX-76—suggest that the constant Γ assumption incorporated in the JWL EOS9,30–32 is not a good assumption for modeling the detonation products response of HEs over volume ranges that encompass both highly compressed and highly expanded states.

Using the detonation products EOS, sound speeds were determined for shock compressed RDX and were used, together with the detonation products Hugoniot curve, to determine the C–J state for RDX. The resulting C–J pressure (35 GPa) for RDX single crystals is the same as that reported previously for FOX-7 single crystals,6 providing experimental confirmation for previous calculations.7 

Our results have also provided insight into the role of crystal orientation on the shock compression response of RDX crystals. Figures 2–5 show that the chemical decomposition develops more rapidly for [100] RDX, compared to [111] RDX, indicating that [100] RDX is somewhat more sensitive under shock compression. However, the orientation dependence reported here for the RDX response is significantly less pronounced than that reported previously for PETN.16 At pressures above the C–J state, the Hugoniot response and detonation products EOS are similar for [100] RDX and [111] RDX. This finding is consistent with the hypothesis that the chemical reactions involved in the overdriven detonation response depend primarily on the structure of the RDX molecules and not on the crystal orientation. Although a corresponding examination of the role of crystal orientation on the shock compression and detonation response of FOX-7 single crystals is very desirable, crystal orientations other than {101} are not available—likely due to the layered crystal structure of FOX-7.40 

The above results and findings demonstrate the usefulness of comparing the response of different HE single crystals to shock compression below and above the C–J pressure. Extension of the same approach to other HE single crystals is expected to provide fruitful insights. Also, the experimental results and findings presented here are expected to motivate theoretical studies to further understand the observed differences in the response of RDX and FOX-7 single crystals over a broad range of pressures below and above the C–J pressure. Such efforts are needed to provide a more detailed understanding of detonation performance and shock insensitivity for IHE and CHE single crystals and to guide the development of HE crystals that further enhance the high performance and insensitivity observed for shock compressed FOX-7.6 

N. Arganbright and K. Zimmerman are acknowledged for their expert assistance with the plate impact experiments. The RDX single crystals used in this work were provided to us by Dr. Daniel E. Hooks of Los Alamos National Laboratory. This work was supported by the Office of Naval Research (Grant No. N00014-19-2047) and the Department of Energy/National Nuclear Security Administration (Cooperative Agreement No. DE-NA0003957).

The authors have no conflicts to disclose.

J. M. Winey: Conceptualization (equal); Formal analysis (lead); Supervision (equal); Writing – original draft (lead); Writing – review & editing (equal). Y. Toyoda: Investigation (lead). Y. M. Gupta: Conceptualization (equal); Supervision (equal); Writing – review & editing (equal).

The data that support the findings of this study are available within the article.

1.
G. W.
MacDonald
,
Nature
95
,
508
(
1915
).
2.
3.
T. R.
Gibbs
and
A.
Popolato
, in
LASL Explosive Property Data
(
University of California
,
Los Angeles
,
1980
).
4.
B. M.
Dobratz
and
P. C.
Crawford
,
LLNL Explosives Handbook: Properties of Chemical Explosives and Explosive Simulants
(
Lawrence Livermore National Laboratory Report UCRL-52997
,
1985
).
5.
T. M.
Klapötke
,
Energetic Materials Encyclopedia
(
de Gruyter
,
Berlin
,
2018
).
6.
J. M.
Winey
,
Y.
Toyoda
, and
Y. M.
Gupta
,
J. Appl. Phys.
130
,
015902
(
2021
).
7.
P.
Samuels
,
K.
Spangler
,
D.
Iwaniuk
,
R.
Cornell
,
E. L.
Baker
, and
L. I.
Stiel
,
AIP Conf. Proc.
1979
,
150033
(
2018
).
8.
9.
W.
Fickett
and
W. C.
Davis
,
Detonation
(
University of California
,
Los Angeles
,
1979
).
10.
A. N.
Dremin
,
Toward Detonation Theory
(
Springer
,
New York
,
1999
).
11.
J. J.
Dick
and
J. P.
Ritchie
,
J. Appl. Phys.
76
,
2726
(
1994
).
12.
J. J.
Dick
,
D. E.
Hooks
,
R.
Menikoff
, and
A. R.
Martinez
,
J. Appl. Phys.
96
,
374
(
2004
).
13.
D. E.
Hooks
,
K. J.
Ramos
, and
A. R.
Martinez
,
J. Appl. Phys.
100
,
024908
(
2006
);
see also the erratum:
D. E.
Hooks
,
K. J.
Ramos
, and
A. R.
Martinez
,
J. Appl. Phys.
109
,
089901
(
2011
).
14.
M. J.
Cawkwell
,
K. J.
Ramos
,
D. E.
Hooks
, and
T. D.
Sewell
,
J. Appl. Phys.
107
,
063512
(
2010
).
15.
K. J.
Ramos
,
D. E.
Hooks
,
T. D.
Sewell
, and
M. J.
Cawkwell
,
J. Appl. Phys.
108
,
066105
(
2010
).
16.
J. J.
Dick
,
J. Appl. Phys.
81
,
601
(
1997
).
17.
J. M.
Winey
,
Y.
Toyoda
, and
Y. M.
Gupta
,
J. Appl. Phys.
127
,
155901
(
2020
).
18.
C. S.
Choi
and
E.
Prince
,
Acta Cryst. B
28
,
2857
(
1972
).
19.
D.
Choudhuri
and
Y. M.
Gupta
,
AIP Conf. Proc.
1426
,
755
(
2012
).
20.
P. A.
Rigg
,
M. D.
Knudson
,
R. J.
Scharff
, and
R. S.
Hixson
,
J. Appl. Phys.
116
,
033515
(
2014
).
21.
A. H.
Jones
,
W. M.
Isbell
, and
C. J.
Maiden
,
J. Appl. Phys.
37
,
3493
(
1966
).
22.
L. M.
Barker
and
R. E.
Hollenbach
,
J. Appl. Phys.
43
,
4669
(
1972
).
23.
J. M.
Walsh
,
M. H.
Rice
,
R. G.
McQueen
, and
F. L.
Yarger
,
Phys. Rev.
108
,
196
(
1957
).
24.
C. M.
Tarver
,
J. W.
Kury
, and
R. D.
Breithaupt
,
J. Appl. Phys.
82
,
3771
(
1997
).
25.
W. L.
Seitz
,
H. L.
Stacy
,
R.
Engelke
,
P. K.
Tang
, and
J.
Wackerle
, in
Proceedings of the Ninth Symposium (International) on Detonation
, edited by
E. L.
Lee
and
J. M.
Short
(
Office of Naval Research
,
Arlington
,
1989
), pp.
657
669
.
26.
R. W.
Woolfolk
,
M.
Cowperthwaite
, and
R.
Shaw
,
Thermochim. Acta
5
,
409
(
1973
).
27.
B.
Sun
,
J. M.
Winey
,
N.
Hemmi
,
Z. A.
Dreger
,
K. A.
Zimmerman
,
Y. M.
Gupta
,
D. H.
Torchinsky
, and
K. A.
Nelson
,
J. Appl. Phys.
104
,
073517
(
2008
).
28.
D. C.
Sorescu
and
B. M.
Rice
,
J. Phys. Chem. C
120
,
19547
(
2016
).
29.
S. A.
Sheffield
,
D. D.
Bloomquist
, and
C. M.
Tarver
,
J. Chem. Phys.
80
,
3831
(
1984
).
30.
E. L.
Lee
,
H. C.
Hornig
, and
J. W.
Kury
, Lawrence Radiation Laboratory Report UCRL-50422 (
1968
).
31.
E.
Lee
,
M.
Finger
, and
W.
Collins
, Lawrence Livermore National Laboratory Report UCID-16189 (
1973
).
32.
E. L.
Lee
and
C. M.
Tarver
,
Phys. Fluids
23
,
2362
(
1980
).
33.
R. G.
McQueen
,
S. P.
Marsh
,
J. W.
Taylor
,
J. N.
Fritz
, and
W. J.
Carter
, in
High-Velocity Impact Phenomena
, edited by
R.
Kinslow
(
Academic
,
New York
,
1970
), pp.
293
417
.
34.
L. G.
Green
,
C. M.
Tarver
, and
D. J.
Erskine
, in
Proceedings of the Ninth Symposium (International) on Detonation
, edited by
E. L.
Lee
and
J. M.
Short
(
Office of Naval Research
,
Arlington
,
1989
), pp.
120
682
.
35.
J. N.
Fritz
,
R. S.
Hixson
,
M. S.
Shaw
,
C. E.
Morris
, and
R. G.
McQueen
,
J. Appl. Phys.
80
,
6129
(
1996
).
36.
C. M.
Tarver
,
J. Phys. Chem. A
124
,
1399
(
2020
).
37.
M. L.
Wilkins
, in
Methods in Computational Physics, Volume 3
, edited by
B.
Alder
,
S.
Fernbach
, and
M.
Rotenberg
(
Academic
,
New York
,
1964
), pp.
211
263
.
38.
Y. M.
Gupta
,
COPS Wave Propagation Code
(
SRI International
,
Menlo Park, CA
,
1976
).
39.
M.
Short
,
J. J.
Quirk
,
C. B.
Kiyanda
,
S. I.
Jackson
,
M. E.
Briggs
, and
M. A.
Shinas
, in
Proceedings of the Fourteenth International Detonation Symposium
, edited by
S.
Peiris
,
C.
Boswell
, and
B.
Asay
(
Office of Naval Research
,
Arlington
,
2010
), pp.
769
778
.
40.
U.
Bemm
and
H.
Östmark
,
Acta Crystallogr. C
54
,
1997
(
1998
).