Abstract

Using the dispersed metal in solid propellants to increase the temperature of combustion products leads to such a problem as the specific impulse loss due to the incomplete combustion of metal particles in the exhaust products. A redistribution of metal loaded into the propellant grain is one of the methods to decrease the specific impulse loss. This paper reports on the ways to obtain the optimum metal particle disposition for the case-bounded propellant grain of tube cross-sectional type. Three different approaches to analyze the metal combustion efficiency are discussed. The influence of the dynamic nonequilibrium of two-phase flow on the optimum metal particles disposition in the propellant grain of tube cross-sectional type is investigated.

1. Introduction

Increasing the energy efficiency of propulsion is an important task of designing solid rocket motors (SRM), one of the component parts of which is the creation of a propellant with the high specific impulse [1]. Composite solid propellants are currently widely used in SRM, which contain aluminium, magnesium, and other metals as fuel additives to increase the specific impulse and combustion stability [2].

It is known the metal additives, on the one hand, lead to improving the energy performance of composite solid propellant due to high temperature of combustion [3], and on the other hand, they are the sources of specific impulse losses due to incomplete combustion of metal particles (especially for small-scale SRM and nozzleless SRM) [4]. The latter can be decreased by the redistribution of metal particles inside a solid propellant in accordance with particularities of a propellant grain and SRM.

The problem of increasing the specific impulse due to the redistribution of metal inside the grain has been dealt with in the last 40 years. For example, in Japan in 1976 it was patented a composition of the solid propellant in which the metal redistribution is carried out by varying the particle size [5]. It was noted that, for the grain of SRM with the channel length to channel diameter ratio equal to 5–30, the redistribution of particle by sizes (44% of channel length is occupied by 50–150 micron particles, 28% by 20–50 micron particles, and 28% by less than 20 micron particles) can give a 4–7% increase in the specific impulse as compared with uniform distribution of the metal powder inside of the propellant. In addition, the metal redistribution accomplished by changing its mass fraction is another way to increase the specific impulse.

A numerical investigation of the influence of the metal particle redistribution in accordance with a metal mass fraction on the combustion efficiency was carried out in paper [6] by assuming that the mass fraction of metal particles varies linearly along the channel. It was shown that there exists an optimal slope of the line providing maximum combustion efficiency for a given mass of metal in the propellant grain. The dependence of the coordinate of the point at which the metal particle leaves the grain surface on the initial size of the particle burning out at the channel outlet was given in the paper [7]. The solution was obtained for the flow in one-dimensional approximation and expressed in terms of modified Bessel functions.

Anyway, the study of the combustion efficiency of metal particles in the combustion chamber of SRM is based on the analysis of the motion of a single particle [6, 8]. The motion of a single particle in the channel of solid propellant grain under the influence of various factors is well studied [913].

It is of interest to reveal the ways of an optimal disposition of metal particles in the propellant grain with a cylindrical channel from the maximum combustion efficiency point of view and to identify main parameters affecting the combustion efficiency.

2. Zero-Dimensional Approach

Combustion products of the metallized composite solid propellant (MCSP) leaving the burning surface are supposed to consist of metal (Al), oxidizing gases (CО2, H2O, and O), and inert gases (I). A chemical reaction takes place between the metal and the oxidizing gas to form a metal oxide and inert gases

Here and are the oxidizing gas and the inert gas, respectively. Equation (1) shows that complete combustion of one mole of metal requires moles of oxidizing gas. This stoichiometric ratio imposes certain restrictions on the permissible content of metal in the MCSP. We define the limit of the metal mass fraction in the MCSP , at which metal can burn completely. Number of gas moles generated at burning of the MCSP equals . Number of metal moles admitted to the flow equals . Therefore, in accordance with (1) for the complete combustion of metal the following equality must hold:

Hence, the limit value of the mass fraction of metal can be found from

Here . Equation (3) shows that depends on and the complex . For values the metal cannot completely burn up, and for oxidizing gases will be sufficient to do it. For different propellants, the molar mass varies from 15 to 30 kg/Kmol, so that the complex varies from 0.8 to 1.6. It follows from (3) that decreasing the molecular weight of gas phase and increasing the mole fraction of oxidizing gas lead to increasing the limit value of .

There are two main time parameters affecting the metal combustion completeness in the SRM, namely, the characteristic burning time of metal particle and the characteristic residence time of metal particle in the combustion chamber of the SRM, .

In the case of zero-dimensional approach, the residence time can be defined as . Assuming the nozzle discharge coefficient and the coefficient of thermal losses equal to unity, we find that and here .

The characteristic burning time of metal particle can be taken from the experimental dependence [14, 15]. Here we use simpler correlation [15]: where 1.062·104 sec/m1.5 for . The condition of the complete combustion of metal one can be written in the form , where . When the ratio is 20 m and is 0.1, then the diameter of particles to be completely burned at the chamber outlet is less than 50 μm, and when equals 0.4, the diameter reaches 100 μm. Decreasing the volume of the combustion chamber leads to decreasing the diameter of particles which can completely burn. The estimate shows that when and , the particle diameter reaches 10 μm and when , the particle diameter is 25 μm.

Thus, the zero-dimensional approach allows us to derive conditions under which the metal particles in the combustion chamber of the SRM will be completely burned: the mass fraction of metal particles in the MCSP has to be less than , determined from (3); the ratio of the burning time of metal particle to the residence time of metal particle in the combustion chamber has to be less than unity, .

3. One- and Two-Dimensional Approaches

The zero-dimensional approach does not allow us to detect the influence of particle distribution along the channel on the combustion completeness of particles. It can be made on the base of one- or two-dimensional approaches.

According to Larson [16], the metal combustion efficiency is defined in terms of mass flow rates as and for the equilibrium flow (equality of gas and particles velocities), the main parameters affecting the efficiency are the ratio of the burning time of metal particle to the mean residence time of metal particle in the combustion chamber of the SRM, , and the ratio of the mole fluxes of metal and oxidizer on the burning surface of propellant, , taking into account the stoichiometric ratio , which can be written as At , we have

The parameter is greater unity in the case of the oxidizer deficiency and is less unity in the case of the oxidizer excess.

As a rule, the MCSP consists of three main components: a binder, an oxidizer, and powdered metal. A mass fraction of metal distributed in the propellant can be expressed through mass densities and volume fractions of the components from where is the volume fraction of metal distributed in the binder, equal to . Since the velocity of combustion products leaving the burning surface of the SRM is expressed in terms of the temperature through the correlation the formula for the parameter is

The necessary conditions for the complete combustion of metal in the propellant grain are and . Following (3), (7) the case corresponds to   .

Analysis of (7)–(10) shows that the increase of the mass fraction of metal leads to the increase in the parameter and the decrease in the parameter and the increase of the mole fraction of oxidizer decreases both and . To evaluate parameters and , letting us take , 3500°K, , , ,  kg/Kmol, and  J/(kg·K), we have , for μm, and for μm.

The residence time of combustion products in the chamber of SRM defined by (4) and used in (10) is an average characteristic, so even if , it does not mean that all the metal from the burning surface of channel is completely burned. Indeed, the metal particle leaving the area located closer to the front end of grain requires more time for the motion before the exiting channel than the particles leaving the area located near the exit of channel. If , there exists a coordinate of a point at the burning propellant surface from which the particle escapes and burns out at the exit of channel. Therefore, at particles burn completely out before the exit, whereas at they do not.

To estimate the coordinate and to define the residence time of metal particle , one can use the approximation of the equilibrium two-phase flow. The rotational flow field in a channel of constant cross-section with a lateral burning surface is given as follows [17]: and the irrotational one is given as follows [18]:

In the case of one-dimensional approach, the continuity equation gives the velocity of product combustion in the form , which coincides with the velocity of the irrotational flow field for two-dimensional case.

One can see from the above correlations that the axial velocity in the case of rotational flow field depends on both and coordinates and in the other cases it depends only on coordinate. Therefore, to find a dependent for the metal particle in the case of the rotational flow field, we have to solve the system of two equations ;  ;  ;   and one equation ; otherwise. Solving the equations gives for the rotational flow field: and for the irrotational flow field

From formulas (13) and (14) it follows that if the particle burning time equals the average residence time of the particle in the channel, , , the value of equals 0.398 for the rotational flow and 0.367 for irrotational flow. It means that particles from more than half of the propellant grain surface do not burn out.

In the limiting case when and , the metal combustion efficiency is yielded by the following.

For one-dimensional flow,

For two-dimensional flow,

So for the uniform distribution of particles along the channel, the value of the metal combustion efficiency is less unity, and for two-dimensional and one-dimensional cases, respectively.

4. Results and Discussions

4.1. Optimum Distribution along the Grain

To provide the equality of the metal combustion efficiency to unity, it is necessary that the residence time of the particle from the point with -coordinate should be greater than its burning time. One can see from the expression (10) for that this may be achieved by variation of or at constant values of and . Let us define the dependent giving the value of equal to unity.

As far as the gas is assumed to be incompressible and flow field is equilibrium, the mole fraction of oxidizer along the particle trajectory can be written as [16] which is valid for both one-dimensional and two-dimensional approaches. The burn rate law for the particle can be represented as Passing to the dimensionless variables , gives Denoting and integrating (19) by the time from 0 to the burning time of particle , we obtain Using previously obtained relations (13), (14) and introducing the variable , we have for two-dimensional flow and for one-dimensional flow Thus, it follows from (21), (22), which can be considered as the optimal distribution of the metal particles, that can be reached by a variation of along at and by a variation of along at . Assuming , we have the following.

For two-dimensional flow,

For one-dimensional flow,

At in both cases, we have —Dirac delta function. Meaning of the latter is that when , the value of parameter (and hence the particle size) should be zero along the entire channel, except for the point .

Figure 1 shows characteristic profiles of at different values of the parameter . The dashed curves correspond to the one-dimensional flow and the solid ones to the two-dimensional flow. Qualitative difference can be clearly seen in the behavior of functions when and for both cases when .

Behavior of along the length of the channel at different values of the parameter is depicted in Figure 2. It is seen that to achieve the equality the metal particles should be absent in the propellant to the right of the intersection point between the abscissa axis and the curve . The increase (decrease) of the value of leads the shift of the intersection point to the left (right) and for one-dimensional flow (dashed curves) the region free of particles are greater than those for two-dimensional flow (solid curves).

4.2. Optimum Distribution across the Grain

During the propellant combustion process, the channel radius decreases that leads to the change of the value of for a given value of . Assuming that the increase of the combustion chamber volume does not affect the temperature in the combustion chamber, the value of is inversely proportional to the channel radius squared in the case of the SRM with the geometry of as follows from (10). The inverse proportionality to the channel radius takes place for SRM with the geometry of . Therefore, when the value of increases, the value of decreases, and dependent and will change. From (21) and (22), one can obtain next correlation where where is the ratio of the current radius of the channel to the initial one. The correlation (25) should be treated either as a relation between the parameters and at constant value of its RHS, which provides , or as an equation for determining the function at given (position of the equal size particles inside the propellant which can completely burn at the channel outlet), or as an equation for determining the function at given (position of propellant portions with the equal mass fraction of metal providing the complete combustion of particles).

Profiles of for 2D flow (solid curves) and for 1D flow (dashed curves) are obtained at , and different values of are depicted in Figure 3. Curves corresponding to different kinds of flow are close to each other near the front end of the grain and are quite different near the rear end of the grain. The distribution of equal size particles () in accordance with their mass content inside the propellant grain is shown in Figure 4.

The increase of (Figure 5) leads to enhancing the region where changes from 0 to 0.9 and increasing a zone free of particles (part of the grain below the curve ). The same behaviour of takes place for and (Figure 6). The fact that the particle residence time in the case of the 1D flow is less than that that in the case of the 2D flow leads to a displacement of profiles (Figure 4) and (Figure 6) downwards at the front end of the grain and upwards at the rear end of the grain with respect to the similar profiles for 2D flow.

Figure 7 shows the dependence of the parameter to complex ratio on the parameter at . It is seen that at (the flux of metal is negligible in comparison with the flux of oxidizing gas) the parameter can be arbitrary due to the arbitrariness of values of and at (stoichiometric ratio of mole fluxes) the value of should be equal to zero.

4.3. Effect of the Velocity Lag

The moving inertial particle has a velocity less than the velocity of the gas whereby the residence time of the particles in the channel increases, which increases the metal combustion efficiency.

Consider the two-dimensional approximation. Following [16], we can write the system of equations describing the motion of a particle in the channel flow field taking into account the fact that the particle velocity lag is negligible at the calculation of the oxidizer concentration along the particle trajectory (so-called King’s model): where

The system of (27) is dimensionless, where the length scale is chosen to be the radius of the channel and the velocity scale is chosen to be .

For rotational flow, and for irrotational flow, ; .

At the initial moment of time, the following conditions are set: ; ; ; .

For irrotational flow, ; for rotational flow, . The condition that the particle resides in the channel is defined by the following inequality: , where is dimensionless length of channel.

To calculate the combustion efficiency of burning particles considering the velocity lag we use the following formula [16]:

Here is a diameter of particle at the exit area which left the burning propellant surface at the point with coordinate . To provide , it is necessary to resolve the system of (27) and to define a value of the coordinate which gives a value of equal to zero.

From the above system of (27)–(30), one can see that the solution depends on, besides and , also three dimensionless parameters , Re, and , where is the relaxation time to average residence time ratio. When , the flow tends to the equilibrium one.

As the propellant burns in the case of a cylindrical channel grain, the outflow velocity of the combustion products from the burning propellant surface changes. It is known [19] the combustion gas velocity at the burning surface is inversely proportional to the channel radius in this case, and the mass flux is proportional to the channel radius raised to the power of , where is an exponent in the burning propellant rate law. Thus, dimensionless parameters , , , Re are also variable. If we denote these parameters at the initial moment of time as , , , , then at another moment of time (in the case of ) they are as follows:

Solving the system of (27)–(30), the dependent which represents the optimum disposition of metal particles in the propellant can be obtained.

The influence of parameters , , on the position of the curve for rotational flow is shown in Figures 8, 9, and 10, respectively. The decrease of values of , and the increase of value of lead to the increase of the deviation of the curve position in the direction of the section from that corresponding to the equilibrium flow. As the shift of the curve in that direction means the increase of the combustion efficiency, then the increase of the value of and the decrease of the values of , are preferable with the energetic point of view in the case of the nonequilibrium flow. It is clear that the closer the curve passes near the curve corresponding to the equilibrium flow, the better the dependent (25) responds to the solution of the system of (27). Therefore, one can say that the dependent (25) describes the solution of the system (27) at large values of and and small .

5. Conclusion

The analytical correlation for the optimum disposition of metal particles in the case-bounded propellant grain of tube cross-sectional type under the assumption of equilibrium two-phase flow is deduced.

Analysis of the dimensionless system of equations describing a motion of burning particle in the channel of circular cross-section with taking into account the velocity lag shows that the optimum disposition of metal particles in the propellant grain depends on five parameters: the relative length of the channel; the particle relaxation time; the molar metal to oxidizer flux ratio at the propellant grain surface; particle combustion to residence time ratio; and Reynolds number of particle.

The increase of the particle relaxation time and the decrease of the particle Reynolds number and the relative channel length are preferable from the combustion efficiency point of view in the case of nonequilibrium flows.

Nomenclature

Nozzle throat area
Particle diameter
Metal combustion efficiency
Functions
Inert gas
Mass rate
Coefficient in the particle burning rate law
Channel length
Dimensionless channel length
Molar mass of species
Exponent in the burning propellant rate law
Oxidizing gas
Exponent in the particle burning rate law
Gas constant
Radius of channel
Reynolds number
Lateral surface area of grain
Time
Temperature
Axial and radial components of velocity
Dimensionless components of velocity
Velocity at the burning surface
Volume
Axial and radial coordinates
Dimensionless coordinates
Mass fraction of metal.
Greek Symbols
Molar fraction of oxidizer
Volume fraction
Dimensionless particle size
Gas dynamic complex
Specific heat capacity ratio of gas phase
Drag function
Bulk density
Dimensionless time
Viscosity
Stoichiometric coefficient
Burning-to-residence time ratio
Mass flux ratio
Molar mass complex
Relaxation-to-residence time ratio.
Subscripts
Oxide
Combustion chamber
Burning
Binder
Burning surface
Inert gas
Gas phase
Limit
Metal
Metal in binder
Oxidizing gas
Oxidizer particle
Channel outlet
Combustion products
Residence
Initial value.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

Support for this work was provided by the Ministry of Education and Science of Russian Federation, Agreement no. 10.1329.2014/K.