Authors: Hanfei Xie, Xiangrong Zhang, Tao Jiang, Yingzhong Zhu, Lin Zhou
Categories: Article, Krieger–Dougherty relation, apparent viscosity, bimodal and trimodal particle-size distribution, melt-cast explosive suspension, shear rate
Source: Polymers
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.
Keywords: melt-cast explosive suspension, apparent viscosity, bimodal and trimodal particle-size distribution, shear rate, Krieger–Dougherty relation
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 particle-size 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 (d50 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 particle-size 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.
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 S3) 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 S1–S6 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 d50 of samples S1–S6 was 5.6, 9.2, 38.8, 162.7, 269.6, and 366.6 μm, respectively.
Figure 1 Photographs of DNP and HMX samples.
Figure 2 Particle-size distributions of six HMX samples.
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 S1–S6, 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≤d50≤100μm [9,33,34,35,36]. Samples with an order of magnitude of the mean particle size of d50>100μm are seldom used because of serious sedimentation [9]. Furthermore, samples with an order magnitude of the mean particle size d50<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 S3–S6 with three combinations (S3 and S4, S3 and S5, and S3 and S6) 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 (d50) ratio of sample S4 to that of sample S3, and ζ43 is the mass ratio of sample S4 to that of sample S3. 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 S3–S6 used but also samples S1 and S2. This issue is further discussed in Section 3.2.
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.
Figure 3 Photograph of the viscosity-measurement setup.
The procedure to measure the apparent viscosity has the following main
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.
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 S5 (d50=269.6μm) and S3 (d50=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.
Figure 4 Effects of the shear rate γ˙ on the apparent viscosity ηa of DNP/HMX melt-cast explosive suspensions.
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 S5 (d50=269.6μm) and S3 (d50=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.
Figure 5 Effects of the solid content ϕ on apparent viscosity ηa of the DNP/HMX melt-cast explosive suspensions.
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 (d50) 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 d50∈1, 10, and 100μm are defined as fine, medium, and coarse particles, respectively. Based on this definition, the present samples S1 and S2 contain fine particles, S3 contains medium particles, and S4, S5, and S6 contain coarse particles. For bimodal particle-size distributions, sample S3 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.
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 S5) and medium (sample S3) particles. The fine particles must then be sample S1, not S2, 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).
Figure 6 Effects of mass ratio ζ531 on apparent viscosity ηa of the DNP/HMX melt-cast explosive suspensions.
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.
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 λ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).
Figure 7 Effects of diameter ratio λcmf on apparent viscosity ηa of DNP/HMX melt-cast explosive suspensions.
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.
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 The relative viscosity ηr as a function of reduced solid content ϕ/ϕm for bimodal particle-size distributions when γ˙=1s−1.
Based on 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 Figure 9, Figure 10 and Figure 11. As expected, for a given shear rate, these data also collapse to a single curve described by Equation (1).
Figure 9 Relative viscosity ηr as a function of reduced solid content ϕ/ϕm for bimodal particle-size distributions with γ˙=10s−1.
Figure 10 Relative viscosity ηr as a function of reduced solid content ϕ/ϕm for bimodal particle-size distributions with γ˙=100s−1.
Figure 11 Relative viscosity ηr as a function of reduced solid content ϕ/ϕm for bimodal particle-size distributions with γ˙=1000s−1.
Compared with the two cases with low shear rates γ˙=1 and 10s−1 (Figure 8 and Figure 9), the experimental data for the two cases with high shear rates γ˙=100 and 1000s−1 (Figure 10 and Figure 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).
The experimental data shown in Figure 8, Figure 9, Figure 10 and Figure 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 particle-size 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 Figure 13, Figure 14, Figure 15 and Figure 16. 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).
Figure 12 Effects of diameter ratio λcmf on maximum solid content ϕm of DNP/HMX melt-cast explosive suspensions.
Figure 13 Relative viscosity ηr as a function of reduced solid content ϕ/ϕm for trimodal particle-size distributions with γ˙=1s−1.
Figure 14 Relative viscosity ηr as a function of reduced solid content ϕ/ϕm for trimodal particle-size distributions with γ˙=10s−1.
Figure 15 Relative viscosity ηr as a function of reduced solid content ϕ/ϕm for trimodal particle-size distributions with γ˙=100s−1.
Figure 16 Relative viscosity ηr as a function of reduced solid content ϕ/ϕm for trimodal particle-size distributions with γ˙=1000s−1.
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 (S1–S6) with particle diameters spanning three orders of magnitude (1, 10, and 100μm).
For bimodal particle-size distributions, only four samples (S3–S6) 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 (S4–S6), medium (S3), and fine (S1 and S2) 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.
One was the fixed-diameter ratio scheme (λcm=λmf), and the other was the fixed-mass ratio scheme (ζcm=ζmf). For either scheme, the apparent viscosity was less than that obtained from the best bimodal particle-size distribution (λ53≈7, ζ53=2). However, the fixed-diameter ratio scheme (λcmf=λcmf≈49:7:1) produced better results compared with the fixed-mass ratio scheme (ζcmf=2:1:0.5). Combining this diameter ratio (λ531) with the optimal mass ratio (ζ531=2:1:0.3), the corresponding apparent viscosity (when γ˙=1s−1 and ϕ=65 wt.%) was 28% less than that obtained from the best bimodal particle-size distribution.
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.
Conceptualization, H.X. and X.Z.; Data Curation, H.X.; Methodology, X.Z., T.J. and L.Z.; Software, Y.Z.; Validation, H.X.; Writing—Original Draft Preparation, H.X. and X.Z.; Supervision, Y.Z.; Project Administration, L.Z. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
The data presented in this study are available on request from the corresponding author.
The authors declare no conflict of interest.
This study was funded by the National Natural Science Foundation of China under Grant No. 11772060 and the Foundation of State Key Laboratory of Explosion Science and Technology under Grant No. YBKT 23-08.
The data presented in this study are available on request from the corresponding author.