3.1. Comparison between Experimental and Simulated Results
The comparison between the experimental and simulated results is presented in
Figure 7. The exposure times of the two X-ray pictures are 16
and 21
, corresponding to the situation when the jet tip reaches 1 CD and 2 CD standoff, respectively. It can be seen from the X-ray that the jet is constantly elongated in the forming process, and the morphology shape of the jet is clear when it reaches 1 CD standoff. At the 2 CD standoff situation, the tip and neck of the jet are not very clear, though the shape and contour of the jet still exist. That may indicate that some materials react at the tip of the jet though the overall shape of the jet is relatively stable at a 2 CD standoff.
At the same time, the morphologies of the reactive material jet formation process shown in the X-ray photographs agree well with the numerical simulation results. In the simulation image, the red part represents Al and the gray part represents the PTFE matrix, respectively. By comparison, it can be found that when the standoff is 1 CD, the head and neck of the jet in the X-ray photo are clearer than the slug, while the Al particles in the jet head and neck are less than the part of the slug in the simulation. It indicates that low content of Al particles in the jet (low material density) would reduce the clarity of the X-ray image. When the standoff is 2 CD, the part with more Al particles in the simulation also has higher definition in the X-ray film. At the same time, compare the two X-ray images, the 1 CD standoff film is brighter than the 2 CD one. That is most likely because the materials in the jet are more diffuse in the 2 CD situation, which leads to lower density than the 1 CD jet, thus the brightness on the 2 CD picture is darker. On the other hand, in the 2 CD standoff case, some materials may react, thus reducing the brightness of the picture. However, according to reference [
6], the average initiation delay time is about 120
, thus the reaction of the reactive materials at 21
is not the main reason for the darker X-ray image. Therefore, it is necessary to analyze the dispersion characteristics of Al particles in the jet forming process for further analysis of the penetration and reaction performance of the reactive material jet.
3.2. Formation of Reactive Material Jet
Figure 8 shows the formation process of a typical reactive material jet. After explosion of the main charge, the reactive material liner is integrally accelerated by the shock wave and detonation products. At 5
and 10
, the inner wall of the liner collapses and gathers into the symmetry axis to form the jet due to its relatively higher speed. Because of the lower density, the PTFE in the liner accelerates faster and becomes the main component of the jet. With the increase of time, the jet containing both PTFE and Al gradually elongates, but the content of PTFE in the jet is significantly higher than that of Al.
To investigate the velocity difference between PTFE and Al during the forming process, the liner and jet in
Figure 8 are divided into six equal-length parts. The curves of the velocity difference between the two components from part 2 to part 6 over time are shown in
Figure 9a. Since the jet head is composed of PTFE with a higher relative velocity, from 5
to 10
, there is a period of time during the jet forming process when the jet head (part 1 when
in
Figure 8) has no Al particles. Thus, the velocity difference in part 1 corresponding to the tip area of the jet is not steady and the data is not included in
Figure 9a. The average velocities of Al and PTFE in the reactive material jet at different times are shown in
Figure 9b.
Overall, the velocity of Al and PTFE increases simultaneously at the initial acceleration process. However, after 5 , the PTFE accelerates faster, up to about 2000 m/s, while the velocity of Al increases more slowly to 1500 m/s. This may be the difference between the density of the two components that contributes to the different particle velocities under the same impact load. However, the two components cannot transfer their stress stably due to the extremely short period of the jet formation process. Therefore, the relative displacement of the two components appears and finally forms the density gradient along the axis of the jet.
For the different parts, the velocity difference curves reach the peak in sequence from part 6 to part 2. It suggests that as the liner collapses, the two components are accelerated and form relative displacement. For each part, the PTFE component moves faster than Al particles, and thus flows forward relative to the Al particles. At the same time, the velocity of the part closer to the jet head is much higher than the backward parts, leading to a larger velocity difference and relative displacement. Finally, fewer aluminum particles remain in the part with higher velocity.
In
Figure 10, the reactive material jet at a standoff of 2 CD is equidistantly divided into 10 parts and the Al particle mass fraction of each part is shown in the calculated Al content curve. The result shows that the Al particle content gradually decreases along the symmetry axis from the slug to the jet tip, which corresponds with the previous statement about the gradient density.
3.3. Granule Size Difference Induced Particle Dispersion
In order to study the influence of the granule size on the particle distribution characteristics of the reactive material jet, four granule sizes (40
, 60
, 80
, and 100
) are chosen to build the discrete models. At the same time, to distinguish the different parts of the jet, the threshold velocity is defined as 2000 m/s in this article to divide the jet into a high-speed section (HSSJ) and a low-speed section (slug). As shown in
Figure 11a, the velocities of all the materials (include Al and PTFE) in HSSJ are over 2000 m/s. Besides,
Figure 11b shows the parts whose velocities are over 2000 m/s and below 2000 m/s in the jets of different standoffs, and the region shape of HSSJ in all standoff cases are basically the same except for stretching over time.
The aluminum content in the HSSJ is analyzed for different granule sizes, at different standoff values. The results are shown in
Table 3. The data in
Table 3 illustrates that the particle content of the HSSJ is rarely affected by the standoffs when the value is more than 1.25 CD. Therefore, the focus of the following analysis will be mainly on the jet under the 2.0 CD standoff condition.
Geometrical morphology of jets formed by granules of different size at 2 CD standoff are shown in
Figure 12. By comparison, smaller granules are more easily dispersed in the jet formed by the reactive material liner. When the granule size is 40
, the aluminum particle can enter the tip area of the jet, while the 100
situation has no particles entering the jet tip. Thus, the dispersion of the particles in the jets increases as the granule size decreases. This effect is related to the different velocities of the granules, as explained in the next sections.
In order to describe the change of aluminum content in each part of the jet, the following Equation (9) is proposed to calculate the relative aluminum content
[
13]:
where
is the aluminum mass fraction in the corresponding part of jet,
is the initial aluminum mass fraction in the reactive material liner.
In
Figure 13, the reactive material jets of different granule diameters are equidistantly divided into 10 parts (the same division method as
Figure 10) and the relative aluminum content
distribution in each part of the cases is shown for each granule size. The plot shows that, as a whole, the aluminum content in the reactive material jet increases at first and then decreases from the head to the tail independently of the granule size. For the four defined granule sizes, there is a maximum Al content in the middle region (parts 4–6) of the jet, increasing with the granule size. On the contrary, for the upper region of the jet (parts 7–10), the Al content decreases and it is lower for the larger granule size condition.
The above results indicate that, in the jet formation process, smaller particles are more likely to enter the front part of the jet (the part with higher velocity), leading to a higher aluminum content in the head of the jet, which is consistent with the particle distribution characteristic in HSSJ in the table above. This is due to the fact that under the same loading condition, larger particles gain lower velocity compared with smaller particles, making it more easily to generate relative displacement. Assuming that the load on aluminum particles in the jet formation process is in a single direction and is proportional to the projection area of the particles, the acceleration of aluminum particles can be calculated by:
where
is the stress per unit area,
is the granule density, and
is the granule diameter. The product (
) is related to mass, therefore the size and mass of granule have the same dependence. Meanwhile, the mass and acceleration are related inversely. Accordingly, the acceleration of the particle is inversely proportional to the granule diameter, and finally results in a higher velocity of the smaller granule.
Figure 14 shows the aluminum particle velocity and the velocity difference (
) profiles over time. The Al particle velocity curves show that the aluminum particles in all cases first accelerate integrally, then decelerate for a period of time, and finally accelerate again. At the same time, larger particles decelerate more than the smaller ones, and finally lead to the increase of the average velocity difference among the four kinds of particles. It should be noted that, in order to show the real speed difference, all the average velocities at each moment in
Figure 14 only calculate material with velocities higher than 0 m/s. Additionally, after 8
, the material in the middle and bottom of the liner is accelerated and calculated into the average velocity of aluminum particles, resulting in a decline section in the velocity curve. Therefore, the deceleration corresponding to the granule size actually reflects its acceleration ability. The velocity difference curves in the
Figure 14 show that, velocity difference increases with the particle diameter, leading to larger relative displacement between the two components. Finally, the content of aluminum granules in the HSSJ decreases while the content of aluminum granules in the slug increases.
The result suggests that because of its weaker speed-up-ability, larger granules would have a higher proportion of particles which concentrate in the middle and bottom area of the jet. On the contrary, smaller particles accelerate faster and can disperse more into the HSSJ. Therefore, the smaller Al particles could be used in the engineering design to increase the density of the jet head which leads to a higher penetration performance. However, a previous study [
10] showed that smaller Al particles speed up the reaction of the reactive material resulting in a weaker penetration performance. Thus, the diameter of the Al particles should be controlled within a relatively appropriate range to obtain a more balanced penetration performance and reaction capacity.
3.4. Evolution of Granules from Liner to Jet
The influence of the initial particle distribution in the liner on the particle distribution characteristics of the jet is studied in this section. A uniform distribution of discrete model is conducted and the reactive material liner is divided into three parts from top to bottom.
Figure 15 shows the movement of the particles from different parts of the liner. Firstly, the particles from the top parts of the liner move to the symmetry axis, and the PTFE matrix accelerates faster than the Al particles due to its lower density. The PTFE then gathers into the jet ahead of the Al particles at 6
. With the collapse of the liner and the elongation of the jet, the particles from the top part of the liner are distributed to the head and tail of the jet, and the particles from the middle and bottom section move to the head and neck of the jet.
The mass fraction and velocity charts of the aluminum particles from different parts of the liner are shown in
Figure 16. In
Figure 16a, the jet formed by the liner with uniform distributed particles is divided into 10 parts, and the mass fractions of the particles from each part are shown as three curves, respectively. As a whole, the particles from three parts of the liner are distributed in sequence from the head to the slug. The bottom particles are mainly distributed in the top and middle area of the jet, and the top particles are mainly distributed in the slug. In HSSJ, most of the aluminum particles are from the bottom and the middle area of the liner, while the top particles can hardly enter the HSSJ.
Figure 16b shows the velocity–time profile of the Al particles from three parts of the liner. The figure shows that, in the formation process, the top particles are first accelerated, and then quickly drop to a relatively low speed. The middle and bottom particles are subsequently accelerated to a higher speed. At the same time, both the curves of the middle and the top particles contain a deceleration period. That is because when the wall of the liner collapses and closes, pressure is transferred from the back forward to the particles closer to the axis, and forcing them into the area with higher velocity. As a result, the particles from other areas slow down. As the jet continues elongating, this energy transfer gradually decreases and the velocity of the middle and top particles tends to be stable.
3.5. Granule Distribution Induced Jet Particle Distribution
According to the above result, the particles from the middle and the bottom part of the liner have a significant impact on the density of the HSSJ. Besides, the powder-pressed manufacturing process also causes uneven density of each part of the reactive material liner, which is due to the concentration effect of the particles [
22,
23].
Figure 17 shows the liner models of different particle distributions. Four liner models with uneven particle distribution (
Figure 17a,b,d,e) are studied and compared with the uniform density one (
Figure 17c). All the liner models are equally divided by area into three parts, and each part is filled with different quantities of Al particles until the overall aluminum particle mass in the liner is
, so that all the cases have the same Al particle mass fraction (
) in the liner. Besides, the Al particle diameter in all samples is set as 40
. At the same time, the particle distribution characteristic of each sample is represented by
, where the middle region of the liners maintains the same concentration in all cases. The five models correspond to configurations with maximum granulate concentration in the bottom region (
), up to the top region (
).
The aluminum mass ratio in HSSJ for different liner granule distributions are shown in
Table 4.
Figure 18 shows the typical images of 2 CD-standoff jet with different initial particle distribution. In all the cases in
Figure 18, the top particles are mainly distributed in the slug, the middle particles are mainly distributed in the back and middle of the jet, and the bottom particles are distributed in the front of the jet. The simulations show clearly that the middle and bottom particles are the main components of aluminum particles in the HSSJ.
Table 4 shows the mass ratio of the particles from different parts of liner to all aluminum material in HSSJ. The results show that with the increase of
(top particles content in liner increases, bottom particles content in liner decreases), the content of bottom particles in HSSJ decreases while the content of the middle particles increases. However, there is a small increase of the content of top particles in HSSJ. As shown in
Figure 18e, the top particles are only distributed in the tip of the jet (the red material) in HSSJ. That is because only a small portion of the top particles which is close to the jet axis can flow into the tip in the formation process. The increase of the portion of the middle particles is due to the decrease of the bottom particle content since the top particles can rarely enter the HSSJ. Besides, the mass ratio of the middle and bottom particles in HSSJ is still above 95%, which means the middle and bottom particles are the key components of the Al in HSSJ. However, in general, the bottom particle density has the greatest influence on the Al particle content in HSSJ.
Figure 19 shows the relative Al content in each part of the jet of different samples. The illustration shows that, with the increase of the
, the position where the maximum relative Al content appears in each sample moves from the middle area to the tail. In the tip of the jet (parts 9 and 10), the sample with more initial particles in the top area of the liner gains more relative aluminum content. That is because the content of the top particles in HSSJ rises with the quantity of the initial top particles in the liner. In the middle parts (parts 4–8), the relative aluminum content in the jet of each sample decreases successively with the increase of
. This is because the bottom particles are mainly distributed in the middle part of the jet, therefore more particles are distributed in the middle part of the jet with the increase of the quantity of the initial bottom particles. On the contrary, the relative aluminum content in the jet of the samples decreases successively with the increase of
in the slug (parts 9 and 10). This is because the top particles are mainly distributed in the slug. As the quantity of particles in the top area of the liner increases, more particles are distributed in the slug.
The change of the Al particle distribution characteristics in the jet with the initial particle distribution in the liner is qualitatively discussed above. It should be noted that, on the premise that the overall Al quantity in the liner remains unchanged, the aluminum content in the HSSJ should be related to the
, which is the ratio of the particle quantity in the top area to that in the bottom area of the liner. Assuming that the relationship between the two can be expressed by
. The fitted curve according to the simulation data is shown in
Figure 20, and the constant
k and
b are fitted as −10.01 and 28.57, respectively. For the fitted curve, it shows a downward trend with the increasing of
.