Boron Particle Ignition in Secondary Chamber of Ducted Rocket

Abstract

In the secondary chamber of ducted rocket, there exists a relative speed between boron particles and air stream. Hence, the ignition laws under static conditions cannot be simply applied to represent the actual ignition process of boron particles, and it is required to study the effect of forced convective on the ignition of boron particles. Preheating of boron particles in gas generator makes it possible to utilize the velocity difference between gas and particles in secondary chamber for removal of the liquid oxide layer with the aid of Stoke's forces. An ignition model of boron particles is formulated for the oxide layer removal by considering that it results from a boundary layer stripping mechanism. The shearing action exerted by the high-speed flow causes a boundary layer to be formed in the surface of the liquid oxide layer, and the stripping away of this layer accounts for the accelerated ignition of boron particles. Compared with the King model, as the ignition model of boron particles is formulated for the oxide layer removal by considering that it results from a boundary layer stripping mechanism, the oxide layer thickness thins at all times during the particle ignition and lower the ignition time.

1. Introduction

Systematic advances in missile propulsion systems technology have provided large increases in missile performance capabilities. Replacement of much of solid propellant rocket oxidizer with free stream air, as in the ducted rocket concept, offers over 5 to 1 increase in missile range capability [1]. Boron has been considered for many years as a prime candidate used for increasing the ducted rocket capabilities based on its high potential energy release on both a volumetric and gravimetric basis coupled with a high energy of combustion, high combustion temperature, and low-molecular-weight products [2]. These properties make boron an attractive material for use in ducted rocket propellants [3]. In order for these advantages to be realized, however, the boron particles must ignite and burn completely within a very limited residence time. Since boron particles are generally initially coated with an oxide layer with inhibits combustion and since boron has an extremely high boiling point, which necessitates surface burning subsequent to oxide removal, this can become difficult, particularly under adverse operating conditions [4].

In previous work, two scenarios for boron ignition and combustion have been identified, which differ in the controlling element in the overall process mechanism [5]. Two controversial theories remain to be resolved for describing boron ignition (removal of oxide coating) in terms of the rate limiting steps [4, 6, 7]. Previous experiments of Macek and semple [8] and Li and Williams [7] have demonstrated that when the boron particle temperature is high (e.g., 1800 K), the rate of boron diffusion exceeds the rate of oxygen diffusion across the oxide layer. Recent experimental observations indicating significant agglomeration and liquefaction of boron at 1213 K in both O₂ and Ar environments have also been interpreted as evidence supporting the diffusion of boron instead of oxygen within the oxide coating as the limiting process [3, 9].

Because ignition and combustion of boron particles involve the release of heat, layer evolution assumes the form of a thermo-hydrodynamic process, coupled with the solid boron substrate underneath and with the ambient oxidizing atmosphere, which is assumed to be stagnant [5]. It is important to note that most the current models of boron particles ignition constrain the symmetry to be conserved. That is to say that the oxidizing atmosphere is assumed to be stagnant and uniform if the particle and its initial oxide layer are given to be spherically symmetric. In the secondary combustion chamber of ducted rocket, there always exists a relative speed between the boron particle and air stream [10]. Hence, the ignition laws under static conditions cannot be simply applied to represent the actual ignition process of boron particles, and it is required to study the effect of forced convective on the ignition of the boron particle. Much research has been conducted to study the ignition of the boron particle, both experimentally and theoretically, in pure oxygen and other atmospheres, with relatively little work available on ignition of the boron particle under forced convection conditions [11, 12]. Dirk Meinkohn developed a thin model to allow the oxide flow field to be derived in the creeping-flow approximation and demonstrated that boron particle ignition may then be caused by the Marangoni effect, which is shown to entail the spreading of punctures and ruptures in the oxide layer, leading to layer thinning or even complete layer removal [13, 14]. Preheating of the boron particles in gas generator makes it possible to utilize the velocity difference between gas and particles in secondary chamber for removal of the liquid oxide layer with the aid of Stokes’ forces. Povitsky and Goldman demonstrated there is indication of regime of fast ignition enhanced by Stokes’ forces by experimentally and numerical simulation [15].

The object of this study is to investigate the oxide layer movement and to characterize its influence on ignition of boron particles in the case of high initial relative particle velocity. An ignition model of boron particles is formulated for the oxide layer removal by considering that it results from a boundary layer stripping mechanism.

2. Formulation

2.1. Boundary Layer Stripping Analysis

2.1.1. Assumptions

The movement of the oxide layer on the front half-surface of the boron particle is shown schematically in Figure 1, using a spherical coordinate system with the particle center as origin. 𝑥 is the curvilinear coordinate along the interface separating the two fluids and 𝑦 is the coordinate perpendicular to it. 𝑦=0 represents the interface between the gas layer and liquid oxide one. 𝑦=+∞ represents the value of gas boundary layer. 𝑦=−∞ represents the value of liquid boundary layer. 𝑢_𝑔, 𝑢_𝑙 represent fluid velocity in gas or liquid boundary layer. 𝑅 is the radius of the boron particle. 𝑋 is the thickness of liquid oxide layer. 𝑈_∞ represents gas freestream velocity. 𝑈(𝑥) represents gas fluid velocity around the boron particle.

Figure 1 

Scheme of oxide layer movement relative to surface of flowing boron particles.

In order to solve the problem of liquid oxide layer movement, the following assumptions are adopted.(1)The temperature inside and around the boron particle is uniformity.(2)The flow is steady and incompressible.(3)There is no flow separation at the interface between solid particle and liquid oxide layer.(4)The value of liquid oxide boundary layer is less than the thickness of liquid oxide layer.(5)The boron particle shape can be approximated by a sphere for the simplified analysis here, the oxide layer thickness is much smaller than particle radius, and hence the gas fluid velocity around the boron particle can be represented in the following expression:𝑈(𝑥)=1.5𝑈_∞𝑥sin𝑅. (1)(6)Equating the shear stress in the gas layer to that in the liquid layer at the interface yields the following equation:𝜇_𝑔𝜕𝑢_𝑔𝜕𝑦_𝑦=0=𝜇_𝑙𝜕𝑢_𝑙𝜕𝑦_𝑦=0. (2)(7)The pressure gradient in the gas layer and the liquid layer is equal: 𝜕𝑃𝜕𝑥_𝑙=𝜕𝑃𝜕𝑥_𝑔. (3)

2.1.2. Equations and Solution

An approximate solution to the two-boundary-layer problem can be obtained by assuming arbitrary simple velocity distributions containing a few parameters and the using the momentum differential relations to determine three parameters [16]. If we assume that the flow is steady and incompressible, then it yields the boundary-layer differential equations for the gas and for the liquid.

The basic equations for the gas is given by𝜕𝑢_𝑔𝑟+𝜕𝑣𝜕𝑥^𝑔𝑟𝑢𝜕𝑦=0,_𝑔𝜕𝑢_𝑔𝜕𝑥+𝑣_𝑔𝜕𝑢_𝑔1𝜕𝑦=−𝜌_𝑔𝜕𝑃+𝜇𝜕𝑥𝜌^𝑔_𝑔𝜕²𝑢_𝑔𝜕𝑦². (4) The boundary conditions for equations above is given by𝑦=0∶𝑣_𝑔=0,𝑢_𝑔=𝑢_𝑔(𝑥),𝜇_𝑔𝜕𝑢_𝑔𝜕𝑦=𝜇_𝑙𝜕𝑢_𝑙,𝜕𝑦𝑦=∞∶𝑢_𝑔=𝑈(𝑥), (5) where the pressure gradient in the gas boundary layer is given by−1𝜌_𝑔𝜕𝑃𝜕𝑥=𝑈𝜕𝑈.𝜕𝑥 (6) Then (4) can be given by𝑢_𝑔𝜕𝑢_𝑔𝜕𝑥+𝑣_𝑔𝜕𝑢_𝑔𝜕𝑦=𝑈𝜕𝑈+𝜇𝜕𝑥𝜌^𝑔_𝑔𝜕²𝑢_𝑔𝜕𝑦². (7)

The basic equations for the liquid are given by𝜕𝑢_𝑙𝑟+𝜕𝑣𝜕𝑥^𝑙𝑟𝜕𝑦=0, (8)𝑢_𝑙𝜕𝑢_𝑙𝜕𝑥+𝑣_𝑙𝜕𝑢_𝑙=𝜌𝜕𝑦𝜌^𝑔_𝑙𝑈𝜕𝑈+𝜇𝜕𝑥𝜌^𝑙_𝑙𝜕²𝑢_𝑙𝜕𝑦². (9)

The boundary conditions for equations above are given by𝑦=0∶𝑣_𝑙=0,𝑢_𝑙=𝑢_𝑔(𝑥),𝜇_𝑙𝜕𝑢_𝑙𝜕𝑦=𝜇_𝑔𝜕𝑢_𝑔,𝜕𝑦𝑦=−∞∶𝜇_𝑙𝜕𝑢_𝑙𝜕𝑦=0. (10)

In order for equations solution, the stream function, 𝜓, is defined by𝐿𝑢=𝑟𝜕𝜓𝐿𝜕𝑦,𝑣=−𝑟𝜕𝜓𝜕𝑥, (11) where 𝐿 denotes the latent length.

Some normalized variable parameters are defined by 𝜉=^𝑥₀𝑟²(𝑥)𝑈(𝑥)𝐿²𝑈_∞(𝑥)𝑟(𝑥)𝑑𝑥,𝜂=𝑦𝑈𝐿2𝜇_𝑔𝑈_∞𝜉/𝜌_𝑔,𝜓=𝐺𝜂2𝜇_𝑔𝑈_∞𝜉𝜌_𝑔. (12)

Then the equations for the gas can be simplified as𝐺^′′′+𝐺𝐺^′′+𝛽(𝜉)1−𝐺^′2=0, (13) where 𝛽(𝜉)=2𝜉𝑈_∞((𝐿²/𝑟²)𝑑𝑈/𝑑𝑥(1/𝑈²)).

The boundary conditions for equations above can be simplified as𝜂=0∶𝐺=0,𝐺^′=𝑢(𝑥),𝜂=∞∶𝐺^′=1. (14)

Then the equations for the liquid can be simplified as𝐿^′′′+𝐿𝐿^′′(𝜉)𝜌+𝛽^𝑔𝜌_𝑙−𝐿^′2=0. (15)

The boundary conditions for equations above can be simplified as𝜂=0∶𝐿=0,𝐿^′=𝐺^′,𝐿^′′=𝐺^′′Γ,𝜂=−∞∶𝐿^′′=0, (16) where Γ=𝜌_𝑔𝜇_𝑔/(𝜌_𝑙𝜇_𝑙), gas viscosity 𝜇_𝑔=6.2×10⁻⁵Pa · s, liquid viscosity 𝜇_𝑙=3.2×10⁻⁴Pa · s.

After what may seem like an interminable manipulation of the governing equations, we have finally set up that particular form of the equations that will be most appropriate as well as convenient for the gas and liquid boundary layer flow. In order to solve (13) and (15), 𝐺^′(0) and 𝐺^′′(0) should be first got. According to literature [17], 𝐺^′(0)=𝑢(𝑥)=(𝜌_𝑔𝜇_𝑔/𝜌_𝑙𝜇_𝑙)^1/3. Firstly, 𝐺^′′(0) is given or estimated, the classical four-order Runge-Kutta method is used to solve these ordinary differential equations. If the obtained boundary values satisfy the actual boundary conditions at the outer edge; namely, 𝐺^′(+∞)=1 and 𝐿^′′(−∞)=0, the calculation can be stopped, and the solutions have been found. Otherwise, we have to revise the estimated value of 𝐺^′′(0) and repeat the numerical integrating procedure until all the boundary conditions at the outer edge are satisfied.

2.1.3. The Conditions of Aerodynamic Shattering

During aerodynamic shattering, the tangential friction stress is propitious to sweep the liquid oxide layer from solid particle, and the surface tension is unpropitious to sweep the liquid oxide layer.

The tangential friction stress can be expressed by𝜏=𝜇_𝑔𝜕𝑢_𝑔𝜕𝑦_𝑦=0=𝜇_𝑔𝑈²𝑟𝐺^′′(0)𝐿2𝜈_𝑔𝑈_∞𝜉. (17)

For a spherical particle, the surface tension can be expressed by𝐶𝐹=¹𝜎𝑅. (18) According to literature [15], the constant 𝐶₁=2.0 and the coefficient of surface tension 𝜎=0.01.

We and Re can be expressed by𝜌We=^𝑔𝑢²_𝑟𝑑_𝑝𝜎𝜌,Re=^𝑔𝑢_𝑟𝑑_𝑝𝜇_𝑔, (19) where 𝑢_𝑟=𝑈_∞−𝑊 represents the relative velocity between gas flow and boron particles. 𝑊 represents the velocity of boron particles. 𝑑_𝑝 represents the diameter of boron particles.

While the surface tension is less than the tangential friction stress, aerodynamic shattering may happen. That is to say, part of liquid oxide layer is detached from the particle if the following equation comes into existence:We√≥8Re9𝐶₁𝐺^′′(0). (20)

The critical velocity 𝑢_𝑟min can be reduced by the relation𝑢_𝑟min=64𝐶²₁𝜎²81𝐺^′′(0)𝜌_𝑔𝜇_𝑔𝑑_𝑝^1/3. (21) Beyond 𝑢_𝑟min, part of liquid oxide layer is detached from the particle and it is shifted to the rear. The formula above permits determination of the range of critical velocities for boron particle coated by a thin liquid oxide layer. The critical velocity versus particle diameter for various gas densities is shown in Figure 2. It may be seen that it is difficult for aerodynamic shattering to happen when the particle diameter is less than 15 microns.

Figure 2 

Critical velocity versus particle diameter for various gas densities.

The mass of fluid in the circumferential liquid oxide layer being swept along by the gas steam at a distance 𝑥=𝜋(𝑅+𝑋)/2 from the stagnation point is𝑑𝑚𝑑𝑡=𝜋𝑑_𝑝𝜌_𝑙⁰_−∞𝑢_𝑙𝑑𝑦=𝜋𝑑_𝑝𝜌_𝑙𝑑_𝑝𝑢_𝑟𝜇_𝑙𝜌_𝑙⁰_−∞𝐿^′𝜂𝑑𝜂. (22)

We can assume that the thickness of liquid oxide layer of boron particle is still uniformity along the uniformity after aerodynamic shattering. The following formula can be reduced𝑑𝑚𝑑𝑡=𝜌_𝑙𝜋𝑑²_𝑝𝑑𝑋𝑑𝑡. (23)

The change in the liquid oxide layer of boron particles due to aerodynamic shattering can be reduced:𝑑𝑋=𝑑𝑡𝑢_𝑟𝜇_𝑙𝑑_𝑝𝜌_𝑙^0−∞𝐿^′𝜂𝑑𝜂. (24)

2.2. Ignition Model of Boron Particle in High-Speed Flow

The results obtained above are used for developing an ignition model of boron particles, adopting a spherical form for the particles. For nonspherical particles, the aerodynamic shattering of liquid oxide layer is the same in principle, but the relevant expressions are more complicated.

In the analysis below we resort to King’s model, based on the assumption that the exothermic reaction between boron and oxygen takes place at the boron-boron oxide interface as a result of oxygen diffusion across the melting oxide layer [4, 15]. As the oxygen converts the boron into boric oxide, thickens the oxide layer, and releases heat, the oxide simultaneously vaporizes at the outer surface, this second process is endothermic and reduces the heating rate. Heat is added to or removed from the particle by convection and radiation, while inside the particle it is transferred by conduction from the surface to the center. Ignition occurs when the liquid oxide layer at the leading point of the particle is completely removed by the tangential friction stress effects [4, 15].

Step 1. B(s)+34O₂⟶12B₂O₃(l)+𝑄_𝑅𝑋,B(l)+34O₂⟶12B₂O₃(l)+𝑄_𝑅𝑋2, (25) where 𝑄_𝑅𝑋=6.132×10⁵J/mole, 𝑄_𝑅𝑋2=6.342×10⁵J/mole. 𝑤_B=64.8×10⁻⁸(𝑅+𝑋)²𝑇_𝑝𝑃_O2𝑒^{−22600/𝑇𝑝}𝑋(mole/s), (26) where 𝑤_B is the molar consumption rate of boron.

Step 2. B₂O₃(l)⟶B₂O₃g−Δ𝐻_vap,𝑤_𝐸=1.005×10¹⁰(𝑅+𝑋)²𝛼𝑒^{−44000/𝑇𝑃}𝑇_𝑝^0.51+4.5×10⁷𝛼𝑃(𝑅+𝑋)/𝑇_𝑝(),Numole/s (27)(28) where 𝑤_𝐸 is the molar rate of vaporization of boron, Δ𝐻_vap=3.78×10⁵ J/mole.

Step 3. B₂O₃(l)+H₂O⟶2HBO₂−Δ𝐻_𝐻,E𝑤_𝐻=9.15×10⁻⁷Nu𝑃(𝑅+𝑋)𝑇_𝑝^1/2exp18.11−2100𝑇_𝑝×+−0.150.0225+𝑃_H2O×exp−18.11−2100𝑇_𝑝^0.5moles, (29)(30) where 𝑤_𝐻 is the molar rate of removal of boric oxide, Δ𝐻_𝐻=3.15×10⁵ J/mole. The equation for particle radius and oxide layer thickness dynamic can be written as 𝑑𝑅𝑤𝑑𝑡=−^B𝑀_B4𝜋𝑅²𝜌_B,𝑑𝑋=𝑤𝑑𝑡^B/2−𝑤_𝐸−𝑤_H𝑀_B2O34𝜋𝑅²𝜌_B2O3−𝑢_𝑟𝜇_𝑙𝑑_𝑝𝜌_𝑙^0−∞𝐿^′𝜂𝑑𝜂. (31)(32) The first term on the right-hand side of formula (33) is due to King’s model, and the second refers to our model of liquid oxide layer movement. For 𝑢_𝑟 less than 𝑢_𝑟min, the second term vanishes. The heat balance of liquid oxide layer is described by the equations (note that three different enthalpy balances must be employed depending on whether the particle temperature is less than, equal to, or greater than the melting point of boron, 2450 K): 𝑑𝑇_𝑝=𝑄𝑑𝑡4𝜋𝑅³𝜌_B𝑐pBs/3+4¹𝜋𝑅²𝑋𝜌_B2O3𝑐_𝑝B2O3𝑇_𝑝<2450K,𝑓=0,𝑑𝑓=𝑄𝑑𝑡4𝜋𝑅³𝜌_BΔ¹𝐻_𝑀𝑇/3_𝑝=2450K,0<𝑓<1,𝑑𝑇_𝑝=𝑄𝑑𝑡4𝜋𝑅³𝜌_B𝑐_pBl/3+4²𝜋𝑅²𝑋𝜌_B2O3𝑐_𝑝B2O3𝑇_𝑝>2450K,𝑓=1, (33) where 𝑓 is fraction of boron in the liquid phase, 𝑇₀ is free stream gas temperature, 𝑇_𝑝 is particle temperature, 𝑇_RAD is surrounding radiation temperature, 𝑄₁,𝑄₂ are heat transfer between particle and the surroundings,𝑄₁=𝑤_B𝑄_𝑅𝑋−𝑤_𝐸Δ𝐻_vap−𝑤_𝐻Δ𝐻_𝐻+4𝜋(𝑅+𝑋)²[ℎ(𝑇_∞−𝑇_𝑝)+𝜎𝜀𝛼_𝑅(𝑇⁴_RAD−𝑇⁴_p)], 𝑄₂=𝑤_B𝑄_𝑅𝑋2−𝑤_𝐸Δ𝐻_vap−𝑤_𝐻Δ𝐻_𝐻+4𝜋(𝑅+𝑋)²[ℎ(𝑇_∞−𝑇_𝑝)+𝜎𝜀𝛼_𝑅(𝑇⁴_RAD−𝑇⁴_p)]. The following values were used for other parameters appearing in equations above: 𝑀_B=10.82, 𝜌_B=2.33×10³kg/m³, 𝑀_B2O3=69.64, 𝜌_B2O3=1.85×10³kg/m³, Pr=0.72, 𝜀=0.8, 𝛼_𝑅=1, Δ𝐻_𝑀=2.1×10³ J/mole, ℎ=1.46×10⁻²×Nu×𝑇_∞^0.8/(𝑅+𝑋)J/(m² · s · K),Nu=2+0.6Re^1/2Pr^1/3, 𝑐_𝑝Bs=2130+0.294𝑇_𝑝 J/(kg · K), 𝑐_𝑝Bl=2.84×10³ J/(kg · K), 𝑐_𝑝B2O3=1.84×10³ J/(kg · K).

3. Results and Discussion

The model of boron ignition developed in this study involves nine independent parameters, value of which must be input to the resulting numerical computer program for prediction of particle ignition time and minimum gas temperature required for particle ignition [18]. There variables are relative velocity between gas flow and boron particles, initial oxide thickness, initial particle temperature, initial particle size, ambient temperature (surroundings gas temperature), effective surroundings radiation temperature, ambient pressure, water vapor mole fraction, and oxygen mole fraction. As part of this study, each of the above independent variables has been systematical varied to determine its effect on whether the particle will ignite and, if so, what the ignition time will be.

3.1. Effect of Ambient Temperature

In Figure 3 plots of oxide layer thickness, particle temperature, particle radius, and fraction boron melted versus time presented for a typical case in which the numerical analysis predicts particle ignition. The particle treated is 10 microns in radius with an initial oxide layer thickness of 0.1 micron and initial temperature of 1800 K. The pressure is 0.5 MPa, the oxygen mole fraction is 0.2, the water vapor mole fraction is zero, and the relative velocity between gas flow and boron particles is 100 m/s. The ambient temperature and effective surroundings radiation temperature are both 2100 K.

(a)

(b)

Figure 3 

Predicted time dependence of important variables for particle which ignites at ambient temperature of 2100 K ((a) King’s model (b) our model with film moving).

As may be shown in Figure 3, once the particle temperature exceeds approximately 2000 K, the oxide evaporation and the change due to aerodynamic shattering rise above the oxide generation rate, resulting in thinning of the oxide layer. The boron particle temperature continue to increase as the oxide layer thins while the boron melting point is reached. At this point the oxide thickness and particle temperature remain constant until the boron melts. After the boron particle has melted, the particle temperature resumes its rise and the oxide layer continues to thin until the particle ignite. Compared with King’s model, as the ignition model of boron particles is formulated for the oxide layer removal by considering that it results from a boundary layer stripping mechanism, the oxide layer thickness thins at all times during the particle ignition and lower the ignition time.

In Figure 4, similar results are presented for a case identical except for reduction of the ambient temperature and effective surroundings radiation temperature by 200 K to 1900 K. The initial stages of the ignition process are quite similar. However, heat losses to the surroundings when the particle temperature rises above the 1900 K ambient temperature are sufficient to retard the evaporation rate sufficiently that it drops back to the generation rate with the result that the oxide layer ceases thinning. In this case, a stable quasi-steady-state is reached. The boron particle temperature remains a steady value, and the ignition does not occur.

(a)

(b)

Figure 4 

Predicted time dependence of important variables for particle which ignites at ambient temperature of 1900 K (a) King’s Model (b) our model with film moving).

3.2. Effect of Total Pressure and Oxygen Mole Fraction

The effects of total pressure and oxygen mole fraction in the gas-free stream on minimum gas temperature required for particle ignition and ignition time are shown in Figures 5 and 6. These calculations were all performed for particles of 10 microns radius with an initial temperature of 1800 K and initial oxide thicknesses of 0.1 microns. The calculations were carried out for a 100 m/s gas stream in the surroundings with water vapor mole fraction of zero at the conditions that ambient temperature and effective surroundings radiation temperature are both 2100 K. Similar trends were noted with other particle sizes and initial oxide thicknesses.

Figure 5 

Minimum gas temperature required for ignition versus total pressure for various oxygen mole fractions.

Figure 6 

Ignition time versus total pressure for various oxygen mole fractions.

As shown in Figure 5, decreased total pressure lowers the gas temperature required for ignition at low oxygen mole fractions but raises the required gas temperature at high oxygen mole fractions. These fairly complex dependencies result from interaction of the effect of oxygen mole fraction on the oxide generation rate and the effect of total pressure on the oxide removal rate. Increase of the total pressure intensified oxygen diffusion through the oxide layer but reduced the oxide vaporizing rate. It also leads to increased gas density, with the attendance opposite influence on tangential friction stress [15]. Similar effects on the ignition time at a fixed gas temperature and velocity as functions of total pressure and oxygen mole fraction are shown in Figure 6. At high oxidizer mole fractions, ignition time reduces with reducing total pressure while at low oxidizer mole fractions the pressure dependency is reversed.

3.3. Effect of Particle Radius

The effect of particle radius on ignition time was also studied. In Figure 7, predictions of ignition time versus particle radius for particle which are predicated to ignite are presented for several gas temperatures ranging from 2050 K to 2200 K, and the radiation surroundings temperature is set equal to ambient gas temperature. The particles (radius size range from 5 to 20 microns with an initial oxide layer thickness of 0.1 micron) were assumed to be preheated to 1800 K prior to ejection gas stream. The pressure is 0.5 MPa, the oxygen mole fraction is 0.2, the water vapor mole fraction is zero, and the relative velocity between gas flow and boron particles is 100 m/s. As shown in Figure 7, the ignition time decreases monotonically with decreasing particle radius due to the decrease in particle mass-to-surface area ratio with decreasing radius. For the higher ambient gas temperature, ignition time is shorter.

Figure 7 

Ignition time versus particle radius for various ambient gas temperatures.

3.4. Effect of Water Vapor Mole Fraction

The effect of water vapor mole fraction in the gas-free stream on ignition time is shown in Figure 8. These calculations were all performed for particles of 10 microns radius with an initial temperature of 1800 K and initial oxide thicknesses of 0.1 microns. The pressure is 0.5 MPa, the oxygen mole fraction is 0.2, and the relative velocity between gas flow and boron particles is 100 m/s. The ambient temperature and effective surroundings radiation temperature are both 2100 K. As shown in Figure 8, the ignition time decreases monotonically with increasing water vapor mole fraction due to treating water vapor as reacting by a diffusion-limited reaction with the boric oxide to aid in removal of the oxide layer.

Figure 8 

Oxide layer thickness versus time for various water vapor mole fractions.

4. Conclusions

The shearing action exerted by the high-speed flow causes a boundary layer to be formed in the surface of the liquid oxide layer, and the stripping away of this layer accounts for the accelerated ignition of the boron particle. It is shown that migration of the liquid oxide layer on the front half-surface of the particle can be derived. The rate of disintegration is founded by integrating over the thickness of the liquid boundary layer to determine the mass flux leaving the boron particle surface at its equator. An ignition model of boron particles is formulated for the oxide layer removal by considering that it results from a boundary layer stripping mechanism. Effects of various parameters on boron particle ignition were studied with this model. Compared with King’s model, as the ignition model of boron particles is formulated for the oxide layer removal by considering that it results from a boundary layer stripping mechanism, the oxide layer thickness thins at all times during the particle ignition and lower the ignition time. The present results support a physical mechanism through which preheated boron particles in gas generator of ducted rockets can be ignited in the secondary chamber by ejection into a high-speed stream.

Acknowledgments

The support of the China Sponsorship Council (CSC) (Grant no. 51006118) is gratefully acknowledged. J. X. Hu thanks Benveniste Natan of Technion-Israel Institute of Technology for the help and encouragement.