Rheological Behavior of DNP/HMX Melt-Cast Explosives with Bimodal and Trimodal Particle-Size Distributions

As a matrix for melt-cast explosives, 3,4-dinitropyrazole (DNP) is a promising alternative to 2,4,6-trinitrotoluene (TNT). However, the viscosity of molten DNP is considerably greater compared with that of TNT, thus, requiring the viscosity of DNP-based melt-cast explosive suspensions to be minimized. In this paper, the apparent viscosity of a DNP/HMX (cyclotetramethylenetetranitramine) melt-cast explosive suspension is measured using a Haake Mars III rheometer. Both bimodal and trimodal particle-size distributions are used to minimize the viscosity of this explosive suspension. First, the optimal diameter ratio and mass ratio (two crucial process parameters) between coarse and fine particles are obtained from the bimodal particle-size distribution. Second, based on the optimal diameter ratio and mass ratio, trimodal particle-size distributions are used to further minimize the apparent viscosity of the DNP/HMX melt-cast explosive suspension. Finally, for either the bimodal or trimodal particle-size distribution, if the original data between the apparent viscosity and solid content are normalized, the resultant plot of the relative viscosity versus reduced solid content collapses to a single curve, and the effect of the shear rate on this curve is further investigated.


Introduction
Although additive manufacturing (3D printing) technologies are emerging in the energetic materials area [1], traditional formatting approvals for energetic materials, especially for melt-cast explosives still dominate. Improving the pourability of melt-cast explosives has at least a twofold significance. First, good pourability aids melt-cast explosive suspensions to flow completely into a mold or projectile casing, which helps eliminate voids or other defects in the final explosive charges. Since the voids and other defects can form potential hot spots in the explosive charges upon experiencing external shock or mechanical stimuli [2], good pourability can improve the safety of explosive charges. Second, good pourability (low viscosity) also helps to achieve a high solid content in melt-cast explosives [3], thus, improving the performance of explosive charges. Generally, good pourability can be accomplished by minimizing the viscosity of melt-cast explosive suspensions. Furthermore, the viscosity depends on many factors, such as the temperature, chemical additives, shear rate, particle morphology, and size distribution [4][5][6][7][8][9].
As one of the most important factors affecting viscosity, the particle-size distribution plays a key role in minimizing the suspension viscosity [10][11][12]. Mixing particles of different sizes is a common way to minimize viscosities because it gives a broader particlesize distribution so that finer particles may fit into the spaces between coarser packed particles [13][14][15]. Therefore, a bimodal particle-size distribution is usually better than a unimodal distribution [16,17], and similarly a trimodal or multimodal distribution is better than a bimodal distribution. However, in practical industrial production or even in laboratory investigations, a balance between quality and efficiency is usually considered when minimizing viscosity so that bimodal distributions are more frequently used compared with other distributions [9,18,19]. Furthermore, for bimodal particle-size distributions, the minimization of viscosity may be affected by two crucial process parameters [20][21][22][23]: (i) the ratio λ of the diameter (d 50 in the present study) of coarse particles to that of fine particles and (ii) the ratio ζ of the mass of coarse particles to that of fine particles.
Based on published values [9,16,20,22,[24][25][26][27], λ is usually in the range of 1-10, while ζ is usually in the range of 0.1-4. Generally, if λ is specified, an optimal ζ exists to minimize the viscosity of suspensions under a certain high solid content; vice versa, there exists an optimal λ if ζ is specified. However, such an optimal λ or ζ is generally determined on the basis of limited experimental comparison schemes, so they are usually local optima rather than global optima. Furthermore, if the concepts of λ and ζ are extended to trimodal or multimodal particle-size distributions, it is even harder to experimentally determine their global optimum.
As a new matrix ingredient in melt-cast explosives, DNP is a promising alternative to TNT [28,29]. The melting point and mechanical sensitivity of DNP are comparable to those of TNT [30], while the performance (e.g., the detonation velocity or detonation pressure) of DNP is much better than that of TNT [31,32]. Conversely, the viscosity of DNP is about 80% greater than that of TNT [33], which is disadvantageous for obtaining a suspension with high solid content. Therefore, compared with TNT-based melt-cast explosive suspensions, the minimal-viscosity requirement of DNP-based melt-cast explosive suspensions is more stringent.
Although Zhou has investigated the viscosity minimization of DNP/HMX melt-cast explosive suspensions with bimodal particle-size distributions [33], he considered only a solid content of 60 wt.% and a diameter ratio λ = 35. Since the viscosity of DNP is much greater than that of TNT, it is essential to investigate the rheological behavior of DNP-based melt-cast explosive suspensions at different solid-content and diameter ratios. Moreover, the particle-size distribution should not be limited to bimodal distributions; trimodal distributions should also be investigated and compared.
We, thus, present herein a detailed study of the rheological behavior of DNP/HMX melt-cast explosive suspensions with bimodal and trimodal particle-size distributions. The viscosity of the suspensions is measured using a Haake rheometer. Bimodal particlesize distributions are first investigated during the process of viscosity minimization of the DNP/HMX melt-cast explosive suspensions, then trimodal particle-size distributions are investigated to further minimize the viscosity of the suspensions, and finally the dependence of relative viscosity on reduced solid content is analyzed for both bimodal and trimodal particle-size distributions.

Materials
Both the DNP and HMX used in this study were supplied by the Yinguang Chemical Industry Group Co., Ltd. (Beiyin, Gansu, China) with purities of 99.1% and 99.4%, respectively. As shown in Figure 1, DNP powder is yellow, whereas the crux HMX particles are white. Figure 1b,c show the photographs of HMX particles (sample S 3 ) taken by camera and by scanning electron microscope (SEM), respectively. Moreover, since the HMX particles used in the formulation of explosives generally measure 1-100 µm, six HMX particle-size distributions (labeled S 1 -S 6 in Figure 2) were collected and used to investigate how bimodal and trimodal particle-size distributions affect the viscosity minimization of DNP/HMX explosive suspensions. The mean particle size d 50 of samples S 1 -S 6 was 5.6, 9.2, 38.8, 162.7, 269.6, and 366.6 µm, respectively.

Schemes of Particle Gradation
As mentioned in Section 1, a bimodal particle-size distribution is commonly used to minimize the viscosity of explosive suspensions. Thus, we first investigate how a bimodal particle-size distribution affects the apparent viscosity of DNP/HMX melt-cast explosive suspensions. Based on the optimal bimodal particle-size distribution determined for λ and ζ, we then use a trimodal particle-size distribution to further reduce the apparent viscosity of the DNP/HMX melt-cast explosive suspensions.
Based on the six samples S 1 -S 6 , there are 15 (given by the combination C(6, 2)) possible values for λ. For a practical bimodal particle-size distribution for melt-cast explosives, however, the order of magnitude of the mean particle size is 10 µm ≤ d 50 ≤ 100 µm [9,[33][34][35][36]. Samples with an order of magnitude of the mean particle size of d 50 > 100 µm are seldom used because of serious sedimentation [9]. Furthermore, samples with an order magnitude of the mean particle size d 50 < 10 µm are avoided because they lead to a large viscosity [9].
Moreover, the typical published values for λ are usually in the range of 1-10 [16,22,25,27]. Therefore, only samples S 3 -S 6 with three combinations (S 3 and S 4 , S 3 and S 5 , and S 3 and S 6 ) were used to investigate the viscosity minimization with bimodal particle-size distribution. The values of λ for the three combinations are about 4, 7, and 9 (rounded to integer values), respectively. Compared to λ, the published values for ζ are usually in the range of 0.1-4 [9,20,24,26]. This study compares only five values for ζ: 0.25, 0.5, 1, 2, and 4. Combining the values of λ with those of ζ gives a total of 15 schemes for particle gradation, which are listed in Table 1, where λ 43 is the diameter (d 50 ) ratio of sample S 4 to that of sample S 3 , and ζ 43 is the mass ratio of sample S 4 to that of sample S 3 . Likewise, λ 53 , ζ 53 , λ 63 , and ζ 63 have similar physical meanings.
After investigating the viscosity minimization with bimodal particle-size distribution, viscosity minimization is further investigated with trimodal particle-size distribution. This time, not only are samples S 3 -S 6 used but also samples S 1 and S 2 . This issue is further discussed in Section 3.2.

Viscosity Measurement
The apparent viscosity of the DNP/HMX melt-cast explosive suspensions was measured using an R/S Haake Mars III rheometer (Thermo Fisher Scientific Inc., Waltham, MA, USA). The measuring element has a Couette geometry with a gap of 2.0 mm. Furthermore, the inner and outer radii of the Couette geometry were 8.5 and 10.5 mm, respectively. Therefore, the assumption of a constant shear rate in the gap was respected, the particle migration [37][38][39] was assumed to be negligible, and the material was assumed to be well-mixed.
Moreover, wall slip is commonly observed in rheometry suspensions, and such phenomena are generally associated with many factors, such as low shear rates and smooth walls of the measuring geometries [40,41]. However, the specific analysis of the wall slip phenomenon will not be detailed in this study. An explosion-proof cover was installed outside the original measuring element of the rheometer. As shown in Figure 3, although this cover installs quite simply, it can resist a blast wave if a sample accidentally explodes.
The procedure to measure the apparent viscosity has the following main steps: 1.
In a double-jacketed stainless-steel kettle, DNP powders are first melted at 100 • C (where an oil bath was used due to the risk of explosive samples), following which HMX particles were added incrementally to the molten DNP.

2.
The DNP/HMX mixture was stirred for 15 min at a rate of 500 rev/min to ensure uniform mixing and to eliminate solid agglomerates.

3.
Uniform DNP/HMX suspensions were poured into the measuring element of the rheometer, and the apparent viscosities of the explosive suspensions were measured at shear rates in the range of 1-1000 s −1 .

Apparent Viscosity with Bimodal Size Distribution
To determine the optimal values λ and ζ of the DNP/HMX melt-cast explosive suspensions with a bimodal particle-size distribution, two optional methods were used. One method is to investigate how the shear rate affects the apparent viscosity of these explosive suspensions for a given solid content, and another method is to investigate how the solid content affects the apparent viscosity at a given shear rate.

Shear Thinning Viscous Behavior
Melt-cast explosive suspensions are solid-liquid two-phase mixtures. Generally, the liquid phase exhibits a Newtonian viscous behavior [3,42,43]. However, when the solid phase is dispersed in the liquid phase, the resultant melt-cast explosive suspensions usually exhibit a non-Newtonian (shear thinning) viscous behavior [6][7][8]. Furthermore, the relationship between the apparent viscosity of these suspensions and shear rates approximately follows a power law [9,44]. Figure 4 shows that, for any scheme listed in Table 2, the apparent viscosity η a of the DNP/HMX melt-cast explosive suspension decreases exponentially with the shear rateγ. The slopes between log (η a ) and log (γ) are negative, demonstrating that the DNP/HMX melt-cast explosive suspension presents a non-Newtonian (shear thinning) viscous behavior. Moreover, for any diameter ratio λ (λ 43 ≈ 4, λ 53 ≈ 7, λ 63 ≈ 9), the best parameter for mass ratio ζ is the same for all ζ (ζ 43 = ζ 53 = ζ 63 = 2). However, compared with the cases of λ 43 and λ 63 , the case of λ 53 corresponds to the minimum viscosity of the DNP/HMX melt-cast explosive suspensions. This is further demonstrated in Table 2, when samples S 5 (d 50 = 269.6 µm) and S 3 (d 50 = 38.8 µm, λ 53 ≈ 7) are mixed together with a mass ratio ζ 53 = 2. The viscosity of this DNP/HMX explosive suspension is, thus, minimized.

Mass Apparent Viscosity η a [Pa·s]
Ratio

. Effects of Solid Content
Compared with the shear thinning behavior, the apparent viscosity η a of the DNP/HMX melt-cast explosive suspensions increases with the solid content φ. However, whether based on η a -γ curves or on η a -φ curves, the best parameters for λ and ζ are the same; i.e., when samples S 5 (d 50 = 269.6 µm) and S 3 (d 50 = 38.8 µm; λ 53 ≈ 7) are mixed together with mass ratio ζ 53 = 2, the apparent viscosity of the DNP/HMX melt-cast explosive suspensions is minimized. Further details appear in Figure 5 and Table 3.

Apparent Viscosity with Trimodal Size Distribution
Compared with the bimodal particle-size distribution, trimodal particle-size distribution produces many more combinations, with different diameter ratios combined with different mass ratios. Therefore, it is much more difficult to experimentally determine the optimal values of λ and ζ with trimodal particle-size distributions. One way to reduce this difficulty is to further minimize the viscosity based on the optimal values of λ and ζ obtained from bimodal particle-size distributions.
However, the concepts of λ and ζ must be extended for trimodal particle-size distributions. For bimodal particle-size distributions, λ represents the diameter (d 50 ) ratio of coarse particles to fine particles, and ζ represents the mass ratio of coarse particles to fine particles. In contrast, for trimodal particle-size distributions, λ and ζ represent the diameter and mass ratios for coarse, medium, and fine particles. Furthermore, for trimodal particle-size distributions, particles with an order of magnitude diameter d 50 ∈ 1, 10, and 100 µm are defined as fine, medium, and coarse particles, respectively. Based on this definition, the present samples S 1 and S 2 contain fine particles, S 3 contains medium particles, and S 4 , S 5 , and S 6 contain coarse particles. For bimodal particle-size distributions, sample S 3 contains fine particles (see Section 3.1).
Moreover, similar to the value of λ with bimodal particle-size distributions, the value of λ with trimodal particle-size distributions is also rounded to the nearest integer. For example, λ 63 is rounded to 9, λ 31 is rounded to 7, and λ 631 is rounded to 63 : 7 : 1. In contrast, ζ is not rounded to integer values for either bimodal or trimodal particle-size distributions.

Fixed Diameter Ratio
Given that the optimal diameter ratio for bimodal particle-size distributions has been experimentally determined (λ 53 ), we assume that the optimal diameter ratio remains unchanged for trimodal particle-size distributions (i.e., the diameter ratio between coarse and medium particles is the same as that between medium and fine particles).
Furthermore, for a trimodal particle-size distribution, λ 53 is the diameter ratio between coarse (sample S 5 ) and medium (sample S 3 ) particles. The fine particles must then be sample S 1 , not S 2 , since the diameter ratio λ 31 is rounded to 7 and equal to λ 53 , whereas λ 32 is rounded to 4. Therefore, the diameter ratio for the trimodal particle-size distribution is denoted λ 531 and is rounded to 49 : 7 : 1.
Given the diameter ratio λ 531 , we investigated how different values of the mass ratio ζ 531 affect the apparent viscosity of the DNP/HMX melt-cast explosive suspensions. As shown in Figure 6, the mass ratio ζ 53 = 2 equals the optimal value for bimodal particle-size distributions, whereas the mass ratio ζ 31 varies in the range of 1.25-5, as listed in Table 4. In any case, the apparent viscosity of this explosive suspension is less than that (12.7 Pa·s; see Figure 4 and Table 2) of a bimodal particle-size distribution with the optimal diameter and mass ratios (λ 53 ≈ 7, ζ 53 = 2). However, the viscosity is not minimized with the mass ratio ζ 53 = ζ 31 = 2 (or equivalently ζ 531 = 2 : 1 : 0.5), although, in this case, the apparent viscosity (10.33 Pa·s) is about 19% less than the minimum viscosity (12.7 Pa·s) of bimodal particle-size distributions. Instead, the minimum viscosity (9.19 Pa·s) is obtained with the mass ratio ζ 531 = 2 : 1 : 0.3, which is about 28% less than the minimum viscosity of the bimodal particle-size distributions.

Fixed Mass Ratio
As in the above discussion for fixed diameter ratio, we again assumed that, for trimodal particle-size distributions, the mass ratio between coarse and medium particles is the same as that between medium and fine particles (i.e., ζ cm = ζ mf = 2, or equivalently ζ cmf = 2 : 1 : 0.5, where the subscripts "c," "m," and "f" indicate coarse, medium, and fine particles, respectively).
Following the above definition of coarse, medium, and fine particles gives six combinations of diameter ratios for trimodal particle-size distributions: λ 431 , λ 531 , λ 631 , λ 432 , λ 532 , and λ 632 , as listed in Table 5. As shown in Figure 7, for any of the six diameter ratios, the apparent viscosity of this explosive suspension is less than that of the bimodal particle-size distribution with the optimal diameter and mass ratios (λ 53 ≈ 7, ζ 53 = 2).  However, as with the case in Section 3.2.1, the viscosity is not minimized with the diameter ratio λ 53 = λ 31 ≈ 7 (or equivalently λ 531 ≈ 49 : 7 : 1), instead, the minimum viscosity (9.58 Pa·s) is obtained with the diameter ratio λ 631 , which is about 25% less than the minimum viscosity of bimodal particle-size distributions.

Relative Viscosity Related to Reduced Solid Content
As shown in Section 3.1, the experimental data relating to the apparent viscosity η a and to the solid content φ are characterized by several curves. However, these curves may collapse to a single curve if the solid content is normalized by the maximum solid content φ m . The normalized solid content φ/φ m is called the "reduced solid content." In the literature, the relative viscosity η r is usually plotted as a function of the reduced solid content [10,14]: η r = f (φ/φ m ), where η r ≡ η a /η 0 and η 0 is the Newtonian viscosity of the liquid phase.
For general applications in suspension rheology, the equation proposed by Krieger and Dougherty is a commonly used semi-empirical relation between relative viscosity η r and reduced solid content φ/φ m and is given by [45,46] where α characterizes the divergence when φ approaches φ m . For rigid spheres, α is usually between one and two [46]. However, for the present non-spherical HMX particles, α may well exceed the above range. Moreover, when the parameter is determined by the best fit of Equation (1) to the experimental data, the maximum solid content φ m is obtained simultaneously. Simply plotting (1/η r ) α against φ yields a straight line so that φ m is the x intercept when (1/η r ) α → 0 [47]. This straight line must pass through the point (0, 1) because, when no particles are dispersed in the liquid phase (φ → 0), the apparent viscosity η a of the suspensions equals the Newtonian viscosity η 0 of the liquid phase (i.e., η r → 1).
Combining the apparent viscosity data ( Figure 5) with the Newtonian viscosity of the liquid phase (η 0 for DNP is 16.4 mPa·s at 100 • C [33]), the parameter α is determined to be 3.94, whereas the maximum solid content φ m varies with both λ and ζ (i.e., φ m = φ m (λ, ζ)) [46]. This is physically reasonable. Generally, φ m is attained when neighboring particles are in permanent contact so that the particles are unable to move past one another, in which case the suspension viscosity becomes infinite [48]. Assuming a fixed particle shape, φ m is then approximately characterized by λ and ζ.
As shown by Table 6, the value of φ m varies within the range of 75.1-77.7 wt.%. As part of the background, the apparent viscosities at the solid content of 65 wt.% are also listed in Table 6. Overall, φ m decreases with viscosity, and the largest value of φ m corresponds to the optimal diameter ratio (λ 53 ≈ 7) and mass ratio (ζ 53 = 2) for the minimization of viscosity of the DNP/HMX melt-cast explosive suspensions. These facts demonstrate that the values of φ m obtained are qualitatively correct.
Given the values of λ and ζ, the respective values of φ m are substituted into Equation (1). The experimental relative viscosity η r correlated with reduced solid content φ/φ m is then plotted in Figure 8. Compared with the original η a -φ data scattered on several curves ( Figure 5), the data η r -φ/φ m collapse to a single curve, demonstrating that Equation (1) can be used to characterize the rheological behavior of the DNP/HMX melt-cast explosive suspensions.   Figure 8, it appears that the η r -φ/φ m relation can be uniquely described by Equation (1) provided that the exponent α = 3.94 with the known maximum solid content φ m (λ, ζ). However, the experimental data used to calibrate the parameter correspond to a shear rate of 1 s −1 . Moreover, since the DNP/HMX melt-cast explosive suspensions present a non-Newtonian (shear thinning) viscous behavior (Figure 4), the relative viscosity of this explosive suspension must be a function of the shear rate, i.e., η r = η r (γ).
Therefore, the right-hand side of Equation (1) must also be a function of the shear rate (i.e., either φ m = φ m (λ, ζ,γ) or α = α(γ)). However, as mentioned above, the maximum solid content φ m corresponds to the limiting case and is dominated by the particle shape and particle-size distribution. On this account, the value of φ m is assumed to have physically nothing to do with the shear rate. Consequently, the exponent α must be a function of the shear rate (i.e., α = α(γ)).
Following the same fitting method for parameter α as used for the shear rate of 1 s −1 , the different values of α were determined for three other typical shear rates (10, 100, and 1000 s −1 ), which provided 3.33, 2.73, and 2.01, respectively. Furthermore, the α-γ data can be well fit by whereγ 0 is the reference shear rate (1 s −1 ). Similarly, the α-ln(γ) relation for Inconel feedstocks was also investigated [49]. On the other hand, the value of φ m is the same as listed in Table 6, regardless of whether the shear rate is high or low. The resultant data η r -φ/φ m for the three typical shear rates are plotted in Figures 9-11. As expected, for a given shear rate, these data also collapse to a single curve described by Equation (1). Compared with the two cases with low shear ratesγ = 1 and 10 s −1 (Figures 8 and 9), the experimental data for the two cases with high shear ratesγ = 100 and 1000 s −1 (Figures 10 and 11) deviate more from their respective fits. This may be because the data between log (η a ) and log (γ) deviate from a linear relationship at high shear rates (see Figure 4). Figure 9. Relative viscosity η r as a function of reduced solid content φ/φ m for bimodal particle-size distributions withγ = 10 s −1 . Figure 10. Relative viscosity η r as a function of reduced solid content φ/φ m for bimodal particle-size distributions withγ = 100 s −1 . Figure 11. Relative viscosity η r as a function of reduced solid content φ/φ m for bimodal particle-size distributions withγ = 1000 s −1 .

Effects of Trimodal Distribution on φ m
The experimental data shown in Figures 8-11 were obtained from the rheological behavior of the DNP/HMX melt-cast explosive suspensions with bimodal particle-size distributions. These η r -φ/φ m data are well described by Equation (1). However, for trimodal particle-size distributions, the applicability of Equation (1) needs further investigation. We assumed that the mass ratio ζ cmf between coarse, medium, and fine particles is the same as the optimal value shown in Figure 6 (i.e., ζ cmf = 2 : 1 : 0.3) and considered all six combinations of diameter ratios (ζ 431 , ζ 531 , ζ 631 , ζ 432 , ζ 532 , and ζ 632 ) as listed in Table 7. The corresponding η a -φ rheological data at the four typical shear rates (1, 10, 100, and 1000 s −1 ) were then measured.
Following the same fitting method for the parameter φ m as used for bimodal particlesize distributions, the values of φ m for trimodal particle-size distributions were obtained and are shown in Figure 12. φ m for trimodal particle-size distributions is greater than φ m for bimodal particle-size distributions. The corresponding η r -φ/φ m data are plotted in  Not surprisingly, for a given shear rate, all η r -φ/φ m data for DNP/HMX melt-cast explosive suspensions with trimodal particle-size distributions also collapse to a single curve and are uniquely described by Equation (1). Table 7. Schemes of particle gradation for trimodal size distribution with the optimal mass ratio shown in Figure 6.

Conclusions
The apparent viscosities of the DNP/HMX melt-cast explosive suspensions with both bimodal and trimodal particle-size distributions were measured using an R/S Haake Mars III rheometer. As a matrix ingredient of melt-cast explosive suspensions, the molten DNP presented a Newtonian viscous behavior. However, when HMX particles (solid phase) were dispersed in molten DNP (liquid phase), the resultant DNP/HMX explosive suspensions presented a non-Newtonian (shear thinning) viscous behavior. The viscosity of the DNP/HMX melt-cast explosive suspensions was minimized by both bimodal and trimodal particle gradation for six HMX samples (S 1 -S 6 ) with particle diameters spanning three orders of magnitude (1, 10, and 100 µm).
For bimodal particle-size distributions, only four samples (S 3 -S 6 ) with particle diameters spanning two orders of magnitude (10 and 100 µm) were investigated. The diameter ratio λ and mass ratio ζ (two crucial process parameters) were used to characterize the bimodal particle-size distribution. Three diameter ratios (about 4, 7, and 9) and five mass ratios (0.25, 0.5, 1.0, 2.0, and 4.0) were compared. For the diameter ratios λ 43 , λ 53 , or λ 63 , the best mass ratios ζ 43 , ζ 53 , or ζ 63 were all equal to 2. However, compared with the diameter ratios λ 43 and λ 63 , the diameter ratio λ 53 (about 7) combined with the mass ratio ζ 53 = 2 produced the minimum viscosity of the DNP/HMX melt-cast explosive suspension.
Such an optimal combination of λ and ζ was determined on the basis of η a −γ or η a − φ data. For trimodal particle-size distributions, coarse (S 4 -S 6 ), medium (S 3 ), and fine (S 1 and S 2 ) particles (with order of magnitude diameters of 100, 10, and 1 µm) were all used to further minimize the viscosity of the DNP/HMX melt-cast explosive suspensions. Based on the optimal values of λ and ζ obtained for a bimodal particle-size distribution, two typical particle-gradation schemes for trimodal distribution were compared.
For both bimodal and trimodal particle-size distributions, the relationship between the relative viscosity η r and the reduced solid content φ/φ m can be described by the equation proposed by Krieger and Dougherty (i.e., η r = (1 − φ/φ m ) −α ). Given the suspension composition and particle shape, the maximum solid content φ m depended only on the diameter ratio λ and mass ratio ζ and was independent of the shear rate. In contrast, the exponent α decreased with the shear rate. When the shear rate was 1, 10, 100, or 1000 s −1 , α was 3.94, 3.33, 2.73, or 2.01, respectively. However, at a given shear rate, all the η r -φ/φ m rheological data of the DNP/HMX melt-case explosive suspension collapsed to a single curve, regardless of whether the particle-size distribution was bimodal or trimodal.
However, the best diameter ratio λ and mass ratio ζ obtained from both bimodal and trimodal particle-size distributions were local optima, not global optima, although we compared all possible combinations of particle-size distributions. Moreover, the above equation correlating η r with φ/φ m is essentially semi-empirical, and thus its generalization to DNP/HMX melt-cast explosive suspensions requires further investigation.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.