# One-Dimension Nonisentropic Model for the Flow of Aluminized Explosive Products

Authors: Ji Duan; Cheng Wu; Fengjiang An; Shasha Liao; Hengzhu Liu; Lihui Dai
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

A new analytic model of aluminized explosive products based on the method of characteristics for planar isentropic flow is proposed herein. The contribution of Al oxidation in the explosion products is investigated analytically. The flow behind the detonation front cannot be treated as isentropic due to the Al oxidation in the products. To solve the nonisentropic flow field of aluminized explosives products analytically, the assumption of local isentropic process is proposed. Based on this assumption, the flow field behind the detonation front of aluminized explosive is a function of only the reacted aluminum mass fraction in each time range. The metal plate test was conducted with the metal plate driven by RDX/Al/wax (76/20/4) and RDX/LiF/wax (76/20/4). The reacted aluminum mass can be obtained indirectly from the experiment results. The reacted aluminum mass was then applied to the analytic model, and the velocity of metal plate driven by RDX/Al/wax (76/20/4) and RDX/LiF/wax (76/20/4) was calculated. The final velocity of the metal plate driven by RDX/Al/wax was 7.8% higher than that driven by RDX/LiF/wax.


---

## Body

## 1. Introduction

The addition of aluminum to condensed explosives to increase the total energy release of the explosive is common practice. Aluminum in its powdered form generally does not react quickly enough to contribute to the detonation front itself [1]; however, it can react in the products of the condensed phase explosive or in the surrounding atmosphere, significantly contributing to the strength and acceleration of the blast wave.Often, a multiphase computational fluid dynamics (CFD) model is used (Frost et al. [2, 3], Zhang et al. [4], Ripley et al. [5], Milne et al. [6], Cooper et al. [7], Massoni et al. [8], Kim et al. [9], etc.) to predict the conditions for metal particle reaction in explosives. Mesoscale modeling with hydrocodes is also widely used (Frost et al. [10], Ripley et al. [11], Zhang et al. [12], Milne et al. [6], etc.). Although these models can, in principle, describe a great amount of detail about the phenomena of interest, developing an analytic model that captures the key elements of the problem in a way that makes the dominant features easily discernible would be preferred. Just like the classical Seshadri formulation [13] analytically describes the structure of premixed flames propagating in a uniform cloud of fuel particles. To establish the analytic model, the author simplified the problem and made some assumptions (gravitational effects, diffusion caused by pressure gradient, and heat transport by radiation is negligible). The Seshadri formulation is used to analyze the flame propagating under normal pressure and temperature. It is not suitable to analyze the combustion of Al and detonation products under high pressure and temperature. In this paper, we have developed a new model based on the classical theory for ideal explosive products (Taylor [14]) that incorporates the oxidation of Al in the products, allowing us to analytically investigate the contribution of Al oxidation in the detonation product.The flow behind the detonation front of the aluminized explosive is more complex than the flow behind the detonation front of an ideal explosive. To analytically solve the flow behind the detonation front of the aluminized explosive, it is necessary to make some assumptions. For example, we propose an assumption called local isentropic process, which enables the conclusion that the flow field behind the detonation front of aluminized explosive is only a function of the reacted aluminum mass fraction at each time. A metal plate test was conducted to obtain the velocity of the plate and, indirectly, the mass fraction of reacted aluminum powder. Applying the model and the described assumptions, we calculated the velocity of metal plate driven by the RDX/Al/wax (76/20/4) and RDX/LiF/wax (76/20/4). The final velocity driven by RDX/Al/wax is 7.8% higher than the final velocity driven by RDX/LiF/wax. We then compared the test result with the result calculated by our analytic model and found that the final velocity of metal plate driven by the RDX/Al/wax (76/20/4) calculated by the analytic model is 4.5% higher than the test result.


## 2. The Change of Entropy of the Detonation Products

The flow of the expanding detonation products is in general a highly complex problem. In the classic theory for the flow of detonation products, it is assumed that the flow behind the detonation front of explosive contained in a tube strong enough to confine the motion to one dimension and detonation products expand in a vacuum withγ=3 (Stanyukovich [15]). The flow of detonation products is isentropic approximately. As a result, the flow of the expanding detonation products can be solved analytically. Due to the reaction of Al particles with the products, however, the classic theory for the flow of detonation products is not suitable for the nonisentropic flow of detonation products of aluminized explosive. In this section, the change of entropy of the detonation products caused by chemical reaction will be discussed.In the present study, it is assumed that the flow of detonation products in the tube is adiabatic. In the second law of thermodynamics, entropy is an extensive state function under a reversible process:dS=dQ/T. However, the expansion of the detonation products with chemical reactions is irreversible process; thus(1)dS>dQT>0,where Q is the heat released by chemical reaction and T is the temperature of the products.According to (1), the entropy of detonation products behind the detonation front increases with mass fraction of reacted aluminum powder. Actually, the change of entropy caused by the chemical reaction contains two parts: entropy flow and entropy generation. The entropy flow is caused by the change of heat in the thermodynamic system and the entropy generation is caused by the irreversible factors. In this section, the heat released by chemical reaction is the main reason for the change of entropy in the products. So the entropy generation caused by the irreversible factors will not be discussed. In order to model the flow behind the detonation front of the aluminized explosive, the contribution of the Al oxidation in the products must be considered. However, the flow of detonation products of aluminized explosive is more complex than the flow of detonation products of ideal explosive. Some assumptions are necessary to solve the model analytically and will be discussed in the next section.


## 3. The Assumption of Analytic Model

To solve the flow behind the detonation front and the contribution of Al oxidation in the products analytically, we simplified the problem. Consider the detonation of aluminized explosive contained in a tube so strong that the motion is confined to one dimension. The flow field generated by the high-explosive detonation products is modeled in a highly idealized manner: the products are treated as a perfect gas expanding into vacuum, with planar geometry. This simplification is based on the classical one-dimension model for ideal explosive. The Al particles are uniformly distributed in the products.


### 3.1. The Assumption of Al Particles in Detonation Products

In this paper, it is assumed that none of the aluminum particles reacted during the detonation and the Al oxidation occurs in the products behind the detonation front [16, 17]. For detonation of aluminized explosives, the duration for energy release of the reaction of explosive components is generally less than 0.1 μs, while the energy release for Al oxidation is in the order of microseconds to several milliseconds [18, 19]. Therefore, the reaction rate of explosive components is much quicker than the rate of Al oxidation. So the assumption mentioned above is valid for aluminized explosive.It is assumed that Al particles in the flow field do not affect the flow before Al particles oxidation. This assumption is valid for the small particle loading in the high explosive (e.g., Rudinger [20]). The small particles could be viewed as small perturbation to the flow field of detonation products.The afterburning of the aluminum is modeled following the work of [21], which was originally applied to the combustion of aluminum in a gas phase. We adapted the Noble-Abel equation to determine the EOS of gas phase aluminum (Kim et al. [22]).


### 3.2. The Assumption of Local Isentropic Process

The flow behind the detonation front of aluminized explosive cannot be treated as an isentropic process approximately due to the chemical reaction of Al particles in the detonation products. We propose an assumption of local isentropic process to analyze the nonisentropic flow behind the detonation front of aluminized explosive.The rate of aluminum oxidation is relatively slow, so we divided the products expansion process into many small ranges along the time axis. In each short time range, we assume that the process is approximately isentropic (the mass fraction of reacted aluminum is approximately constant). Using this method simplifies the flow behind the detonation front of aluminized explosive and enables analytical evaluation. We call this assumption local isentropic process.Note that the detonation of aluminized explosive is extremely violent and fast; it is difficult to obtain detailed quantitative data related to this phenomenon directly. In this paper, the comparison of the velocities determined experimentally and obtained through the model is conducted. The model correctly described the contribution of Al oxidation in the detonation products of aluminized explosive. The assumption of local isentropic process is indirectly validated by this result.


## 3.1. The Assumption of Al Particles in Detonation Products

In this paper, it is assumed that none of the aluminum particles reacted during the detonation and the Al oxidation occurs in the products behind the detonation front [16, 17]. For detonation of aluminized explosives, the duration for energy release of the reaction of explosive components is generally less than 0.1 μs, while the energy release for Al oxidation is in the order of microseconds to several milliseconds [18, 19]. Therefore, the reaction rate of explosive components is much quicker than the rate of Al oxidation. So the assumption mentioned above is valid for aluminized explosive.It is assumed that Al particles in the flow field do not affect the flow before Al particles oxidation. This assumption is valid for the small particle loading in the high explosive (e.g., Rudinger [20]). The small particles could be viewed as small perturbation to the flow field of detonation products.The afterburning of the aluminum is modeled following the work of [21], which was originally applied to the combustion of aluminum in a gas phase. We adapted the Noble-Abel equation to determine the EOS of gas phase aluminum (Kim et al. [22]).


## 3.2. The Assumption of Local Isentropic Process

The flow behind the detonation front of aluminized explosive cannot be treated as an isentropic process approximately due to the chemical reaction of Al particles in the detonation products. We propose an assumption of local isentropic process to analyze the nonisentropic flow behind the detonation front of aluminized explosive.The rate of aluminum oxidation is relatively slow, so we divided the products expansion process into many small ranges along the time axis. In each short time range, we assume that the process is approximately isentropic (the mass fraction of reacted aluminum is approximately constant). Using this method simplifies the flow behind the detonation front of aluminized explosive and enables analytical evaluation. We call this assumption local isentropic process.Note that the detonation of aluminized explosive is extremely violent and fast; it is difficult to obtain detailed quantitative data related to this phenomenon directly. In this paper, the comparison of the velocities determined experimentally and obtained through the model is conducted. The model correctly described the contribution of Al oxidation in the detonation products of aluminized explosive. The assumption of local isentropic process is indirectly validated by this result.


## 4. Model of Aluminized Explosive Products

In this study, it is assumed that the detonation of aluminized explosive contained in a tube is strong enough to confine the flow of the detonation products to one dimension. As mentioned in the assumption of aluminized explosive products, the high-explosive detonation products are treated as a perfect gas and the flow field generated by the high-explosive detonation products is not affected by Al particles before their oxidation. The state equation of high-explosive detonation products is(2)Pexplosive=bρRT,where b is the initial high-explosive mass fraction of aluminized explosive.It is also assumed that the Al particles are uniformly distributed in the explosive and the Al particles exist in gas phase in the products. To determine the EOS of gas phase aluminum, we adapted the Noble-Abel equation of the form [22](3)PAl=aλAlρRT1-nAnnAn:～0.75,where R denotes gas constant, λAl is the reacted aluminum mass fraction, ρ is the density of aluminized explosive products, a is the initial aluminum mass fraction of aluminized explosive, n is the number of moles per unit volume, and An is an empirical constant. The value of nAn is obtained from [22], due to the fact that aluminized explosive studied in [22] is similar to our study. (Note that the influence of value of nAn on our model is weak.)The standard mixture rule applies to the aluminized explosive products. The state equation of aluminized explosive products is(4)Ptotal=Pexplosive+PAl=b+aλAl1-nAn·ρRT=AλAl·ρRT.Based on the first law of thermodynamics, the equation can be obtained as follows:(5)Cp=Cv+dPvdT,where Cp is specific heat of products at constant pressure, Cv is specific heat of products at constant volume, P is pressure of products, and v is specific volume.By inserting (4) in (5), it can be found that(6)Cp=Cv+db+aλAl/1-nAn·RTdT.Following the assumption of local isentropic process, the expansion of products in the tube was divided into several small time ranges. The rate of Al oxidation is slow relative to the short time ranges and the mass fraction of reacted aluminum is approximately constant. In each time range, (6) can be simplified as follows:(7)Cp=Cv+b+aλAli1-nAn·R.The adiabatic index of the products is(8)γ=CpCv.By inserting (8) in (7), it can be found that(9)Cp=γγ-1b+aλAli1-nAn·R.Based on the first law of thermodynamics, the following equation can be obtained:(10)dh=de+Pdv+vdP=CpdT.Equation (10) can be simplified to(11)dS=dqT=CpdTT-vTdP.The state equation of aluminized explosive products (4) can be inserted into (11), and the following equation can be obtained:(12)dS=CpdTT-b+aλAli/1-nAn·RPdP.By inserting (9) in (12), it can be found that(13)dS=CpdTT-γ-1γCpdPP.Based on the assumption of local isentropic process, the flow of detonation products is approximately isentropic in each time range. Thus(14)TPγ-1/γ=Ci.By inserting (4) in (14), it can be found that(15)Ptotalργ=b+aλAli1-nAn·CiRγ,where Ci is constant and λAli is the mass fraction of reacted aluminum in time range i.The equations which represent the conservation of mass and momentum behind the detonation front of the aluminized explosive are(16)∂ρ∂t+u∂ρ∂x+ρ∂u∂x=0,∂u∂t+u∂u∂x+1ρ∂p∂x=0.In each time range, using the method of characteristics for planar isentropic flow, (16) can be converted to ordinary differential equation. In time range i, the Riemann invariant along a right-running characteristic is(17)ui+2γ-1ci=constant.Whenγ=3, u+c=constant. Therefore, the right-running wave has a constant slope in time range i:(18)dxdt=ui+ci.The same can be shown for left-running characteristics. However, the slope of the characteristic is different in every time range due to the reaction of Al oxidation in the detonation products. Therefore, the characteristics of aluminized explosive and ideal explosive differ, as shown by the sketch in Figure1.Figure 1
Sketch of characteristics for aluminized explosive and ideal explosive.From the sketch of characteristics for aluminized explosive and ideal explosive (Figure1), the influence of chemical reaction in the products is apparent. Based on the assumption of local isentropic process, the characteristic for the flow of detonation products of aluminized explosive is divided into several segments along time axis. In each time range, the flow behind the detonation front of aluminized explosive is approximately isentropic. However, the total flow behind the detonation of aluminized explosive is nonisentropic.For the flow of the detonation products of aluminized explosive, it can be obtained thatS1≠S2≠S3≠S4≠S5≠S6 which is depicted in Figure 1. Parameters of the flow such as sound speed c and pressure of products p are not unchanged along the characteristic line. The parameters of the flow field at time t1 will be analyzed first. Because mass is conserved and the volume is unchanged at the time t1, the chemical reaction does not influence the density. In this study, it is assumed that the energy released by the chemical reaction mainly influences the pressure and sound speed of the products, while the particle velocity and approximate temperature of the products remain unchanged. According to (4), it can be obtained that(19)p2=AλAl2AλAl1p1,ρ2=ρ1,u2=u1,c2=γp2ρ2,where p, ρ, u, and c are pressure, density, particle velocity, and sound speed in the products, respectively, and the subscripts represent the time range. Thus, the parameters of detonation products along the characteristic in the time range i+1 can be calculated as follows:(20)pi+1=AλAli+1AλAlipi,ρi+1=ρi,ui+1=ui,ci+1=γpi+1ρi+1.Based on the assumption of local isentropic process, the parameters of the flow along the characteristic are only related to the mass fraction of reacted aluminumλAli. Therefore, the flow field behind the detonation front of aluminized explosive can be calculated analytically when the mass fraction of reacted aluminum is known.


### 4.1. The Metal Plate Acceleration Test and Mass Fraction of Reacted Aluminum

The metal plate acceleration test was conducted to obtain the acceleration ability of the explosive and the mass fraction of reacted aluminum. The metal plate was driven by two kinds of explosives in this experiment. The explosives and charge size are shown in Table1. Sample 1 and sample 2 accelerate the same copper plate, which has a charge size of Φ40 mm × 1 mm. The experimental setup is shown in Figure 2.Table 1
Components and charge size of the explosive.
Sample number
Proportion of components (wt.%)
Diameter of Al particle (μm)
Charge size (mm)
1
RDX/Al/wax 76/20/4
5
Φ40 × 15
2
RDX/LiF/wax 76/20/4
5
Φ40 × 15Figure 2
Schematic diagram of the metal plate acceleration test.The reaction time of explosion is on the order of a microsecond, making it difficult to obtain detailed quantitative data of the mass fraction of reacted aluminum directly. Unlike Al powder, LiF always remain inert. Based on the test result, the mass fraction of reacted aluminum of detonation products can be obtained indirectly. The difference in kinetic energy of the metal plates driven by RDX/Al/wax and RDX/LiF/wax is the energy contributed by Al oxidation in the detonation products. Thus(21)σQAlωaλt=12mVAl2t-VLiF2t,where σ is the efficiency of Al oxidation reaction heat in detonation products, QAl is the Al oxidation reaction heat, λ is the mass fraction of reacted aluminum, ω is the mass of aluminum explosive, a is the mass fraction of aluminum in the explosive, m is the mass of metal plate, and VAl and VLiF are the velocities of the metal plates driven by RDX/Al/wax (sample 1) and RDX/LiF/wax (sample 2), respectively. The results of metal plate acceleration test are shown in Figure 3.Figure 3
Velocity of the metal plate driven by RDX/Al/wax and RDX/LiF/wax.The reactions of Al with detonation products mainly contain three parts:Al+1.125CO2→0.5Al2O3+0.75CO+0.37C△H1=0.5×(-1675.7)+0.75×(-110.525)-1.125×(-393.509)=478.046kJ/molAl+1.5CO→0.5Al2O3+1.5C△H2=0.5×(-1675.7)-1.5×(-110.525)=672.063kJ/molAl+1.5H2O→0.5Al2O3+1.5H2△H3=0.5×(-1675.7)-1.5×(-241.83)=475.105kJ/molQAl=△H1+△H2+△H3/27/3=20.126kJ/gWhenQAl=20.126kJ/g and σ=0.18 [23], the mass fraction of reacted aluminum in the products can be calculated, as shown in Figure 4.Figure 4
The mass fraction of reacted aluminum in detonation products.


## 4.1. The Metal Plate Acceleration Test and Mass Fraction of Reacted Aluminum

The metal plate acceleration test was conducted to obtain the acceleration ability of the explosive and the mass fraction of reacted aluminum. The metal plate was driven by two kinds of explosives in this experiment. The explosives and charge size are shown in Table1. Sample 1 and sample 2 accelerate the same copper plate, which has a charge size of Φ40 mm × 1 mm. The experimental setup is shown in Figure 2.Table 1
Components and charge size of the explosive.
Sample number
Proportion of components (wt.%)
Diameter of Al particle (μm)
Charge size (mm)
1
RDX/Al/wax 76/20/4
5
Φ40 × 15
2
RDX/LiF/wax 76/20/4
5
Φ40 × 15Figure 2
Schematic diagram of the metal plate acceleration test.The reaction time of explosion is on the order of a microsecond, making it difficult to obtain detailed quantitative data of the mass fraction of reacted aluminum directly. Unlike Al powder, LiF always remain inert. Based on the test result, the mass fraction of reacted aluminum of detonation products can be obtained indirectly. The difference in kinetic energy of the metal plates driven by RDX/Al/wax and RDX/LiF/wax is the energy contributed by Al oxidation in the detonation products. Thus(21)σQAlωaλt=12mVAl2t-VLiF2t,where σ is the efficiency of Al oxidation reaction heat in detonation products, QAl is the Al oxidation reaction heat, λ is the mass fraction of reacted aluminum, ω is the mass of aluminum explosive, a is the mass fraction of aluminum in the explosive, m is the mass of metal plate, and VAl and VLiF are the velocities of the metal plates driven by RDX/Al/wax (sample 1) and RDX/LiF/wax (sample 2), respectively. The results of metal plate acceleration test are shown in Figure 3.Figure 3
Velocity of the metal plate driven by RDX/Al/wax and RDX/LiF/wax.The reactions of Al with detonation products mainly contain three parts:Al+1.125CO2→0.5Al2O3+0.75CO+0.37C△H1=0.5×(-1675.7)+0.75×(-110.525)-1.125×(-393.509)=478.046kJ/molAl+1.5CO→0.5Al2O3+1.5C△H2=0.5×(-1675.7)-1.5×(-110.525)=672.063kJ/molAl+1.5H2O→0.5Al2O3+1.5H2△H3=0.5×(-1675.7)-1.5×(-241.83)=475.105kJ/molQAl=△H1+△H2+△H3/27/3=20.126kJ/gWhenQAl=20.126kJ/g and σ=0.18 [23], the mass fraction of reacted aluminum in the products can be calculated, as shown in Figure 4.Figure 4
The mass fraction of reacted aluminum in detonation products.


## 5. The Flow of Aluminized Explosive Products behind Metal Plate

The motion of explosion products is confined to one dimension; thus when the detonation front reaches the metal plate, there will be a reflected compression wave propagating into the products. The nonsimple wave region behind the reflected compression wave and the nonsimple wave region can be represented by the following equations:(22)x=ui+cit+Fi,(23)x=ui-cit+Fi′,where the subscript i represents the flow in the time range i. Fi is a parameter related to the mass fraction of aluminum λi. The mass fraction of aluminum in time range i is treated as approximately constant so Fi is also a constant in time range i. Fi′ is also related to the motion of metal plate. Therefore, the motion of metal plate must be calculated first in solving the flow of detonation products.The equation of conservation of momentum for the metal plate is(24)MdVdt=Arpm,where M is the mass of the metal plate, V is the velocity of the metal plate, Ar is the surface area of the metal plate, and pm is pressure of detonation products on the surface of metal plate. Considering the equation of the state, the following equation can be obtained:(25)pmpi=cmciγ,where cm is the sound velocity on the surface of the metal plate and pi and ci are the initial condition in the time range i.By inserting (25) in (24), it can be found that(26)dVdt=ArpiMciγcmγ.When the compression waves reach the metal plate, the velocity of detonation products changes fromui to um which is equal to V. The velocity of sound also changes from ci to cm. Thus(27)ui+ci=um+cm.The derivation of (22), then, is(28)dxdt=um+cm+dumdt+dcmdtt.By inserting (26) in (28), it can be found that(29)dcmdt+cmt+ψcmγ=0,where ψ=Arpi/Mciγ.In this study,γ≈3 and solve the differential equation (29) as follows:(30)1+2ψcm2tcm2t2=ϑ1=constant.In time range 1, the initial conditions are(31)p1=pH=14ρ0D2,c1=cH=34D.When the first compression wave reaches the metal plate, the boundary conditions are(32)t=LD,um=0,cm=D,where ρ0 is initial density of explosive, pH and cH are C-J parameters, L is the length of explosive, and D is the detonation velocity.By inserting the boundary conditions in (30), it can be found that(33)ϑ1=1+ηL2=constant,where η=32/27m/M.By inserting (33) in (30), it can be found that(34)cm=Ltξ1,where ξ1=1+η1-L/Dt-0.5.By inserting (34) in (26) and integrating, it can be found that(35)V=D1+2ξ1-1ηξ1-Lξ1Dt.The initial and boundary conditions of the time rangei+1 can be obtained from the properties of the time range i. The velocity of the metal plate driven by aluminized explosive in time range i+1 can be calculated:(36)V-Vi=AλAliAλAl01/4Lξi+1titi+2Dη1ξi+1ti-Lξi+1t-2Dηξi+1,ξi+1=AλAliAλAl0L2ϑi+1-ηLDt-0.5,ϑi+1=1+AλAl0/AλAliηcm2titi/LDcm2titi2=constant,where ti is the initial time and Vi is the initial velocity of the metal plate in the time range i+1.In time range 1 (t0～t1), the nonsimple wave region behind the metal plate can be represented by the following equations:(37)x=u1+c1t+F1,(38)x=u1-c1t+F1′,where F1=0.Based on the characteristic method, it can be obtained that(39)u1+c1=um+cm.By inserting (34) and (39) in (37), it can be found that(40)um=x-F1t-Lξ1t.By subtracting (34) from (40), it can be obtained that(41)um-cm=u1-c1=x-F1t-2Lξ1t.By inserting (41) in (38), it can be found that(42)F1′=F1+2Lξ1.In time range 2 (t1～t2), the nonsimple wave region behind the metal plate can be represented by the following equations:(43)x=u2+c2t+F2,x=u2-c2t+F2′.As mentioned in the section of the effect of Al oxidation, the density and velocity of particles along a characteristic remain unchanged. At the timet1, the following relationship can be obtained:(44)x=u1+c1t1=u2+c2t1+F2.ThereforeF2 can be obtained:(45)F2=1-AλAl1AλAl0c1t1,F2′=F2+2Lξ2.In time rangei (ti-1～ti), the nonsimple wave region behind the metal plate can be represented by the following equations:(46)x=ui+cit+Fi,x=ui-cit+Fi′,Fi=1-AλAli-1AλAli-2ci-1ti-1,Fi′=Fi+2Lξi.Therefore, the flow of the aluminized explosive products behind the metal plate can be solved by the described one-dimension nonisentropic model.


## 6. The Application of the Model to the Aluminized Explosive

The one-dimension nonisentropic model is applied to calculate the velocity of the metal plate driven by RDX/Al/wax (76/20/4) and RDX/LiF/wax (76/20/4). The characteristics of the explosives are shown in Table2.Table 2
Characteristics of the explosive and the charge size.
Explosive
Aluminum content (wt.%)
Density (g⋅cm−3)
Detonation velocity (m⋅s−1)
Charge size (mm)
RDX/Al/wax 76/20/4
20
1.82
8000
Φ40 × 15
RDX/LiF/wax 76/20/4
—
1.82
8000
Φ40 × 15The initial and boundary conditions are shown in Table3.Table 3
The initial condition and boundary condition in time range 1.
Initial condition
p
H (GPa)
28.8
ρ
H (g⋅cm−3)
2.4
u
H (m/s)
2000
c
H (m/s)
6000
Boundary condition
u
m (m/s)
0
c
m (m/s)
8000
(a)pH, ρH, uH, and cH represent the pressure, density, velocity of particle in the products, and sound speed on the detonation front.
(b)um and cm represent the velocity of particle in the products and sound speed behind the metal plate.Based on the assumption of local isentropic process, the process of products expansion was divided into 35 time ranges on average. The velocity of the metal plate and sound speed of the detonation products behind the metal plate were calculated and are depicted in Figure5.Figure 5
Velocity of metal plate and sound speed of the detonation products behind the metal plate, as calculated by the one-dimension nonisentropic model.When calculated by the one-dimension nonisentropic model, the final velocity of the metal plate driven by RDX/Al/wax was 3570 m/s and the velocity when driven by RDX/LiF/wax was 3310 m/s. The final velocity of the metal plate driven by RDX/Al/wax was 7.8% higher than that of the RDX/LiF/wax. When experimentally tested, the final velocity of metal plate driven by RDX/Al/wax was 3418 m/s, 9.5% higher than that of the RDX/LiF/wax. The final velocity calculated by the model is higher than the experimentally determined velocity (Table4).Table 4
Comparison of the final velocities determined experimentally and through the model.
Explosive
Final velocity (m/s)
Experimental
Model
RDX/Al/wax
3418
3570
RDX/LiF/wax
3121
3310Al oxidation increases the acceleration ability of the explosive, which can be seen when compared to RDX/LiF/wax (Figure5). The one-dimension nonisentropic model analytically describes the contribution of Al oxidation in the products. The problem was simplified and the assumption of local isentropic process was employed to analyze the contribution of Al oxidation in the products. The result calculated by the model cannot be as precise as the result which is calculated by CFD; however the objective of the present study is to develop a simplified analytical model to describe the contribution of Al oxidation, in the same fashion of the classical isentropic flow model for detonation products of ideal explosive.Note that the model will produce large error when the interval of time range is divided too largely. We divided the process of detonation products expanding into 10, 15, 20, and 25 time ranges averagely and calculated the metal plate velocity, respectively. When the process is divided into 20, 25, and 35 time ranges averagely, the velocities of metal plate are almost the same. But the results on 10 and 15 time ranges are not good. If each time range is too large, the mass fraction of reacted aluminum cannot be treated as constant approximately and the result calculated by our model will produce large error.


## 7. Summary

In this study, a simple model for aluminized explosive products was constructed and used to calculate the motion of a metal plate driven by aluminized explosive. This model operates on the assumption that the Al particles have no effect on the flow field before Al particles oxidation. This one-way coupling assumption is valid for small particle loading in the high explosive. The detonation of aluminized explosive contained in a tube is strong enough to conifer the motion to one dimension. It also assumes that the process of products expanding in the tube is adiabatic. In order to solve the complete flow field analytically, the expansion process was divided into several time ranges. In the short time ranges, we assume that the process is approximately isentropic. Based on the assumptions, the model can be solved analytically. The metal plate acceleration test was conducted to indirectly obtain the velocity of metal plate and the mass fraction of reacted aluminum of detonation products. Finally, the model was applied to calculate the motion of the metal plate, and the model and experimental results were compared. The analytic model result was slightly higher than the test result. The model was a simple method that captured the contribution of the Al oxidation in the products to acceleration ability. Moreover, the flow of aluminized explosive products behind the metal plate can be solved based on the model. The application of the analytic model for the detonation products of aluminized explosive is analogous to the classical isentropic flow model for detonation products of ideal explosives.



---

*Source: 1717095-2017-04-04.xml*