Atomization of Gelled Propellant Simulant with Carbon Particles

QIANG Hong-fu, LIU Hu, HAN Qi-long, WANG Guang, HAN Ya-wei

(1. 601 Staffroom, Xi'an Hi-Tech Institution, Xi′an 710025, China; 2. 96263 Unit of the Second Artillery, Luoyang 471000, China)

1 Introduction

As a new kind of rocket propellant, gelled propellants have advantages for both liquid and solid propellants, including high density, combustion energy, security and long reservation period, etc. However, compared with the traditional liquid propellants, the non-Newtonian character makes the gelled propellants difficult to be atomized. Atomization has attracted widespread and lasting attentions since the beginning of the gel propulsion technology[1-2].

The addition of energetic particles like carbon, aluminum or boron is essential to gelled propellants, and it may significantly increase the energy content per unit volume of the gelled propellant and improve the performance of gelled propulsion system[3-4]. As the addition of energetic particles would alter the rheological properties of gelled propellant[5], some researches have been made to investigate the atomization characteristics of gelled propellants with energetic particles. Jayaprakash et al[6]investigated the injection and atomization characteristics of gelled kerosene with 30%(mass fraction) Al particles, they found that Sauter mean diameter(SMD) of the gel spray was more sensitive to the impingement angle and was dependent on the injection pressure in a highly non-linear manner. Kampen et al[7-8]detail studied the influence of Al particles content of gelled Jet A-1 fuels on rheology, atomization and combustion. The gels prepared in their research appeared "solid" at unstressed ambient conditions and a distinct yield stress occurred. With different generalized Reynolds numbers, different atomization modes were observed. Baek et al[9]investigated the atomization behavior of C934 Carbopol gels with and without 15% SUS304 nanoparticles. They found the nanoparticles decreased the gel strength and made breakup length of gel with nanoparticles remarkably shorter than that of the pure Carbopol gel.

Overall, the mechanics of atomization of gelled propellant with energetic particles are still far from being fully understood in these limited researches. In this work, a new simulant containing carbon particles was prepared and used in atomization experiments. The rheological properties were measured and a series of atomization experiments were made. The linear stability theory was adopted to study the breakup characteristics of the liquid sheet.

2 Experiment

2.1 Preparation of Gelled Propellant Simulants

The gelled propellant simulant was prepared by dissolving 1%(mass fraction) high-molecular polymer, 5%(mass fraction) carbon particles with an average diameter of about 5 μm in de-ionized water and mixing with an electric mixer for 20 min at 2000 revolutions per minute. This carbon-loaded simulant is named as S1 in this paper. We also prepared another simulant S2 with 99%(mass fraction) de-ionized water and 1%(mass fraction) high-molecular polymer for comparison. The physical and rheological properties of simulants S1 and S2 are similar to gelled propellant: their densities (ρ) are 1010.1 kg·m-3and 1001.7 kg·m-3, and surface tension coefficients (σ) 0.067 N·m-1and 0.072 N·m-1, respectively. As the shear rate of the gelled propellant is high in atomization. The rheological properties of the simulants were measured by a rotational rheometer when shear rate and by pipe-flowing experiments whenγ≥103s-1.

The relationship of apparent viscosity (η) and shear rate(γ) can be described by power-law equation, Herschel-Bulkley (HB) equation, Herschel-Bulkley Extended(HBE) equation, etc.[8]. The simulants prepared in this paper appear “syrupy” at unstressed ambient conditions. The experimental results show that the yield stresses (τ0) of the two simulants are very low and both below 10 Pa. Therefore, the yield stress was neglected and the most common power-law equation was adopted.

(1)

ηandγof simulants obtained from experiment and the fitted power-law constitutive curves are given in Fig.1, and their physical properties are listed in Table 1. From Fig.1, it can be seen that simulants S1 and S2 are pseudo plastic fluids and theirηdecrease with the increases ofγ. And compared with simulant S2 without carbon particles in Fig.1, simulant S1 shows higher apparent viscosity under the same shear rate than that of S2. The addition of carbon particles increases the consistency coefficient and decreases the flow index of the simulant, which indicates that the addition of carbon particles evidently changes the physical and rheological properties of the gel. As the de-ionized water is main component of the simulants, the densityρand surface tension coefficientσof the two simulants are similar to ones of water (ρ=1000.0 kg·m-3andσ=0.073 N·m-1), as shown in Table 1.

Fig.1 Apparent viscosity and shear rate of gelled propellant simulants

Table 1 Physical property of gelled propellant simulants with and without carbon particles

simulantscarbonparticlesρ/kg·m-3σ/N·m-1power-lawconstitutiveparametersK/Pa·snnS1with1010.10.06716.590.29S2without1001.70.0727.080.37

2.2 Experimental Apparatus

Fig.2 is the schematic diagram of gelled propellant atomization experiment system. At the beginning of the experiment, high pressure gas would be filled into the tank to force the gelled propellant simulant to the jet injector across pipes and values. The impingement angle 2θand the jet velocityvjetwere tuned by adjusting the angle of the doublet injectors and changing the mass flow rate by the control value, respectively. The atomization processes were recorded by a Phantom V12.1 high speed camera with 784×800 image resolution and 5 μs shutter speed. The atomization images were passed to the data acquisition system for further analysis. Meanwhile, important data in the experiment, such as mass flow rate, pressure in the tank, pressure in the injectors, etc. were measured and recorded by the data acquisition system during the whole experiment.

Fig.2 Schematic diagram of gelled propellant atomization experiment system

2.3 Experimental Design

In order to analyze the atomization characteristics of gelled propellant simulant with carbon particles, 10 atomization experiments were designed with different impingement angles 2θ, jet velocitiesvjet, injector orifice diametersdand injector orifice length to diameter ratioL/d0, etc., as shown in Table 2.

The generalized Reynolds number(Regen) was used to describe the flow behavior of the power-law fluid, which is defined as ref.[10]:

(2)

whereρpis the density of the power-law fluid.

Table 2 Conditions for the atomization experiment

case2θ/(°)vjet/m·s-1d/mmL/d0Regen1601018160626015183211360221861824901018160659015183211690221861827601513.532118602213.56182960100.5813131060101.581806

3 Results and Analysis

3.1 Experiment Results and Analysis

In this paper, the atomization quality is evaluated by the atomization angle (β) and atomization patterns. As we know, when two jets impinge with each other, a fan-shaped liquid sheet forms, and the angle between the left and right rims of the liquid sheet is called atomization angle (β, as shown in Fig.3). Generally, largerβmeans better atomization quality. According to former researches[8,11], with different gels, jet velocities, impingement angles, generalized Reynolds numbers, etc., there are different atomization patterns, including close-rim, open-rim, ligament and fully-developed patterns, etc.[8,11]The close-rim and open rim patterns indicate poor atomization qualities, while the other two indicate better atomization qualities.

Fig.3 is atomization images with the impingement angle 2θ=60° and different velocities for cases 1-3. As shown in Fig.3a, with a low jet velocity (10.2 m·s-1) andRegen, a fan-shaped liquid sheet forms, and it has a distinct rim at the upstream and breakups into ligaments and large drops downstream. The atomization pattern in Fig. 3a is the so-called “open-rim pattern” with the atomization angleβof about 70°. Fig.3b also shows an open-rim pattern withβ=85° under a larger jet velocity (15.68 m·s-1) andRegen=3211. As shown in Fig.3c, withvjet=22.14 m·s-1andRegen=3211, the rim of the liquid sheet becomes indistinct, and clearer bow-shaped impact waves occurs and makes the liquid sheet more unstable and breakup into more ligaments and drops. The pattern in Fig.3c is called “ligament pattern”.βin Fig.3c is about 100°. As discussed in Ref.[8], in cases 1-3, the generalized Reynolds numberRegen, which is completely determined by jet velocity, could be adopted as the indication of the atomization quality. LargerRegen(jet velocity) means a larger kinetic energy in jet impingement and would lead to better atomization quality.

Fig.4 shows atomization images with 2θ=90° and different velocities for cases 4-6. Compared with cases 1-3 in Fig.3, theβin Fig.4 are larger, which are about 80°, 110°, 140° in Fig.4a, Fig.4b, Fig.4c, respectively. Fig.4a shows an open-rim pattern, and Fig.4b and Fig.4c show ligament patterns. The atomization qualities at impingement angle 2θ=90° are better than ones with 2θ=60° and the similar jet velocities.

Fig.5 shows that the atomization images with ratio of injector orifice length to diameter (L/d0=3.5)under different jet velocities for cases 7 and 8. Compared with case 2 and case 3 withL/d0=8 shown in Fig.3b and Fig.3c, there are no obvious differences in Fig.5. Atomization patterns in Fig.5a and Fig.5b are almost identical to patterns in Fig.3b and Fig.3c respectively, which indicates that there are no essential differences on the flow characteristics of the jets sprayed from injectors withL/d0=3.5 and 8. As a result, the ratio of injector orifice length to diameter has little influences on the atomization patterns.

Fig.6 shows atomization images with impingement angle 2θ=60° and different injector orifice diameters for cases 9 and 10. As shown in Fig.3a (case 1) and Fig.6, the larger orifice diameters lead to larger mass flow rates and will produce larger liquid sheets. But there are also no essential differences on the atomization patterns of cases 1, 9 and 10, they are all open-rim patterns with almost the same atomization angle about 70°.

In all the cases studied in this paper, the gelled propellant simulant could hardly to be atomized into fine drops, the main atomization products are ligaments and large drops. Within the investigation conditions, the atomization quality increases with the increase of jet velocity and impingement angle, while the changes of injector length to diameter ratio and the injector orifice diameters influence the atomization patterns little.

a.vjet=10.2 m·s-1, b.vjet=15.68 m·s-1, c.vjet=22.14 m·s-1,

β=70°β=85°β=100°

Fig.3 Atomization images of cases 1-3

a.vjet=9.77 m·s-1, b.vjet=15.28 m·s-1, c.vjet=21.71 m·s-1

Fig.4 Atomization images of cases 4- 6

a.vjet=15.35 m·s-1b.vjet=22.22 m·s-1

Fig.5 Atomization images case 7 and case 8

a.vjet=9.80 m·s-1b.vjet=10.49 m·s-1

d0=0.5 mmd0=1.5 mm

Fig.6 Atomization images of case 9 and case 10

3.2 Linear Stability Analysis

Linear stability theory is widely used to evaluate the instabilities of the liquid sheet[12-13]. In this theory, the liquid sheet instability is mainly due to the aerodynamics interactions between the liquid and its surrounding gas. There are two kinds of disturbances that will occur on the liquid sheet: symmetric and anti-symmetric. Squire[14]showed that the anti-symmetric disturbance played a dominant role on breaking the liquid sheet into fragments. Therefore, only the anti-symmetric disturbance is considered here. Schematic of a moving liquid sheet under anti-symmetric disturbance is shown in Fig.7, a two dimensional liquid sheet moves into a quiescent, inviscid, incompressible gas with velocityUs, the thickness of the liquid sheet is 2hs, the surface tension of liquid isσ, the densities of liquid and gas areρ1andρgrespectively, the density ratio of gas and liquid isRg1=ρg/ρ1. Generally, the wave amplitude on the liquid sheet can be expressed as:

(3)

Fig.7 Schematic of a moving liquid sheet under anti-symmetric disturbance

Chojnaki[15]deduced the dispersion relation for a plan liquid sheet based on the power-law constitutive:

(4)

(5)

(6)

If the liquid sheet breakups when the wave amplitude reachesηb, the breakup timeτbcan be obtained as follows:

τb=ln(ηb/η0)/ωi,max

(7)

whereωi,maxis the maximum grow rate, then the breakup length can be calculated as:

Lb=Usln(ηb/η0)/ωi,max

(8)

Here ln(ηb/η0) is set to be 12 according to Ref. [13].

In this paper, the breakup length of the liquid sheet is defined as the axial distance from the impingement point to the point where the liquid sheet along the axis begins to breakup, as shown in Fig.8. We assume the liquid sheet speedUs=0.92vjetaccording to ref. [16]. Solving equation (4) withρ1=1010.1 kg·m-3,ρg=1.225 kg·m-3,K=16.59 Pa·sn,n=0.29,σ=0.067 N·m-1and 2hs=2.0×10-4m under jet velocity of 10 m·s-1(cases 1, 4), 15 m·s-1(cases 2, 5) and 22 m·s-1(cases 3, 6), the effects of sheet velocity on the stability of the liquid sheet can be obtained, as shown in Fig.9. It can be seen that the maximum disturbance wave grow rate increases with the increase of the sheet velocity. It means that the liquid sheet will become more unstable at larger sheet velocities, which agrees well with the experiment results of cases 1-3 and cases 4-6 under the same velocities.

With equations (4), (7) and (8), the breakup length of the liquid sheet can be predicted. Fig.10 is the comparison of breakup lengths of liquid sheets predicted by linear stability theory and measured from experiments. As shown in Fig.10, the variation trend of breakup lengths of liquid sheets calculated from linear stability theory is consistent with the ones measured from experiments. At a low Weber number, the linear stability analysis evidently overestimates the breakup length when compared with the experiment. The relative error between the predicted and measured breakup lengths is about 24% atWes=128. As the Weber number increases, the relative errors decrease to 14.9% atWes=618. The errors of the predicted values are considered as a results of neglect of instabilities caused by jet impingement in the linear stability theory.

Fig.8 Definition of the breakup length of the liquid sheet

Fig.9 Influence of sheet velocityUson the stability of the liquid sheet calculated by linear stability theory

Fig.10 Comparison of breakup lengths of liquid sheets from linear stability analysis and experiments

4 Conclusions

(1) The carbon particles increase the consistency coefficient and decrease the flow index of the simulant S1, which makes simulant S1 show high apparent viscosity under the same shear rate.

(2) The simulant S1 can only be atomized into ligaments and large drops. The atomization quality improves with the increases of jet velocity and impingement angles, while the changes in injector length to diameter ratio and the injector orifice diameters show little influences on the atomization patterns.

(3) There are about 14.9%-24% relative errors between the predicted and measured breakup lengths, but the tendency of the breakup lengths predicted by the linear stability theory agrees well with the experiment ones.

[1] Natan B, Rahimi S. The status of gel propellants in year 2000[J].InternationalJournalofEnergeticMaterialsandChemicalPropulsion, 2002,5(1-6): 172-192.

[2] LIU Hu, QIANG Hong-fu, WANG Guang. Review on Jet Impingement Atomization on Gelled Propellant[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2015, 23(7): 697-708.

[3] Hodge K, Crofoot T, Nelson S. Gelled propellants for tactical missile applications. AIAA 99-2976[R], 1999.

[4] Haddad A, Natan B, Arieli R. The performance of a boron-loaded gel-fuel ramjet[J].ProgressinPropulsionPhysics, 2011,2: 499-518.

[5] Wu Zhijian, Hu Lirong. Performance Research of Metallized Gelled Propellant[J].MissileandSpaceVehicle, 2006, 283(03): 52-55.

[6] Jayaprakash N, Chakravarthy S R. Impingement Atomization of Gelled Fuels. AIAA 2003-316[R], 2003.

[7] von Kampen J, Madlener K, Ciezki H K. Characteristic Flow and Spray Properties of Gelled Fuels with Regard to the Impinging Jet Injector Type. AIAA 2006-4573[R], 2006.

[8] von Kampen J, Alberio F, Ciezki H K. Spray and combustion characteristics of aluminized gelled fuels with an impinging jet injector[J].AerospSciTechnol, 2007,11:77-83.

[9] Baek G, Kim S, Han J, et al. Atomization characteristics of impinging jets of gel material containing nanoparticles[J].JNon-NewtonFluid, 2011,166(21): 1272-1285.

[10] Metzner A B, Reed C J. Flow of non-Newtonian fluids-correlation of the laminar, transition, and turbulent-flow regions[J].AmericanInstituteofChemicalEngineeringJournal, 1955,4(1): 189-204.

[11] Fu Q, Yang L, Zhuang F. Effects of Orifice Geometry on Spray Characteristics of Impinging Jet Injectors for Gelled Propellants[R].AIAA 2013-3704, 2013.

[12] Ryan H M, Anderson W E, Pal S, et al. Atomization characteristics of impinging liquid jets[J].JPropulPower, 1995,11(1): 135-145.

[13] Dombrowski N, Johns W R. The aerodynamic instability and disintegration of viscous liquid sheets[J].ChemEngSci, 1963,18(3): 203-214.

[14] Squire H B. Investigation of the instability of a moving liquid film[J].BritishJournalofAppliedPhysics, 1953,4: 167-169.

[15] Chojnacki K T. Atomization and mixing of impinging non-Newtonian jets[D]. Huntsville: University of Alabama-Huntsville, 1997.

[16] Heislbetz B, Madlener K, Ciezki H K. Breakup Characteristics of a Newtonian Liquid Sheet formed by a Doublet Impinging Jet Injector.AIAA2007-5694[R], 2007.