03964.com

文档资料库 文档搜索专家

文档资料库 文档搜索专家

JOURNAL OF PROPULSION AND POWER Vol. 27, No. 1, January–February 2011

Theoretical Investigation of the Dynamics of a Gas–Liquid Coaxial Swirl Injector

Qing-fei Fu? and Li-jun Yang? Beijing University of Aeronautics and Astronautics, 100191 Beijing, People’s Republic of China

DOI: 10.2514/1.B34004 In this paper, a linear dynamic characteristics model of a gas–liquid coaxial swirl injector is developed by dividing the injector into two parts: the recessed chamber and the section before the recessed chamber. The gaseous injector ?ow before the recessed chamber is considered to be steady; only the response of the mass ?ow rate at the exit of coaxial injector to the pressure ?uctuation at the liquid injector inlet has been investigated. The dynamic characteristics of the coaxial swirl injector can now be considered as two processes: the dynamics of the liquid swirl injector and the dynamics of the recessed chamber. The transfer function of the whole injector is obtained by combining the transfer functions of these two processes. The in?uence of the injector con?guration and working condition on the injector dynamic characteristics was also studied. A smaller mixing ratio, larger mass ?ow rate, higher ambient pressure, smaller ratio of recess length to diameter, and larger geometry parameter would decrease the amplitude response of the mass ?ow rate pulsation, but these parameters have little effect on the phase frequency of the coaxial swirl injector.

I. Introduction

IPROPELLANT coaxial injectors are widely used in the liquid propellant rocket engine (LPRE), especially in liquid oxygen (LOX)/hydrogen rocket engines. Coaxial injectors can be classi?ed into shear coaxial injectors and swirl coaxial injectors, according to the spraying mechanism. For a coaxial swirl injector (Fig. 1), LOX is injected with low velocity through a center swirl injector and gaseous hydrogen is injected with high velocity through an annular gap around the center injector. The swirl causes the LOX sheet to have a hollow cone shape, which impinges on the surrounding hydrogen stream; providing better atomization and higher performance. It has been reported that recessing the inner LOX injector (with respect to the exit surface of the outer hydrogen injector), leads to performance enhancement of LPRE [1]. There have been many studies on the spray and combustion characteristics of the coaxial swirl injector with the recess con?guration [2–6], however, the dynamic characteristics of this type of injector have been less studied. The dynamic characteristics of the recessed chamber are of signi?cant importance to the entire coaxial swirl injector. Three different ?ow models were used for the recessed mixing chamber [7]: the outer mixing ?ow condition, the inner mixing ?ow condition and the critical recess ?ow condition. The gas–liquid interaction in the recessed chamber is very complicated, especially under the inner mixing ?ow condition (Fig. 2). Hence, it is necessary to study the dynamic characteristics of the coaxial swirl injector under the inner mixing ?ow condition. Bazarov [8] reported that during ?re tests at 50–75% of nominal thrust level of LOX=LH2 LPRE, high frequency pressure pulsations appeared in the oxidizer collector of the main combustion chamber, with amplitude up to 0.9 MPa, accompanied by destructive vibrations, with acceleration up to 500–800 g. The origin of these pulsations cannot be explained by the acoustic characteristics of the combustion chamber, because these frequencies were much different from those observed; this phenomenon is called self-pulsation. The reason for the occurrence of self-pulsation may be the hydraulic

Received 6 May 2010; revision received 5 September 2010; accepted for publication 13 September 2010. Copyright ? 2010 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved. Copies of this paper may be made for personal or internal use, on condition that the copier pay the $10.00 per-copy fee to the Copyright Clearance Center, Inc., 222 Rosewood Drive, Danvers, MA 01923; include the code and $10.00 in correspondence with the CCC. ? Ph.D. Candidate, School of Astronautics. ? Professor, School of Astronautics; yanglijun@buaa.edu.cn (Corresponding Author). 144

B

instability of the conical liquid sheet and the momentum transfer of kinetic energy between the high speed gas stream and liquid sheet in the recessed chamber [9]. Because the liquid sheet emanates from the center swirl injector with a speci?c spray angle, the annular gas stream is blocked by the conical liquid sheet. The gas stream then pushes the liquid sheet into the centerline, and the liquid and gas phase violently exchange the momentum. When self-pulsation occurs, a loud “scream” and strong ?ow rate oscillation follow. Such self-pulsation motion was discovered in all types of injectors used in LPRE, such as the shear coaxial injector, coaxial swirl injector with central LOX and exterior hydrogen stages, and both swirl/swirl coaxial injector with central liquid ?ow, etc. To understand the mechanism of the unsteady self-pulsation phenomenon and to eliminate the instability, the dynamic characteristics of the gas–liquid coaxial injector were ?rst studied by Andreyev et al. [10,11]. Andreyev studied the dynamic characteristics of the gas–liquid injector and the effect of injector dynamics on the entire LPRE dynamic system. The theoretical study was based on fundamental ?uid mechanics’ equations and the diagram method, which reveals the interrelation of the dynamic process. Experimental methods were also introduced. Nonlinear dynamics, emphasizing the interaction of the gas and liquid phase, the stability boundary of selfpulsation and the effect of injector geometry on the self-pulsation, were also discussed [10]. Bazarov and Yang [12] performed several experimental studies on the in?uences of operating conditions and design parameters on the occurrence of self-pulsation. According to his results, the LOX postrecess length is the most important parameter in determining self-pulsation characteristics. Im and Yoon [9] investigated selfpulsation characteristics of a gas/liquid swirl coaxial injector by measuring spray patterns, spray oscillation characteristics, and the self-pulsation boundary. Using indirect photography, they observed strong and periodic spray oscillations. The spray patterns varied at various injection, geometric, and ambient pressure conditions. The main conclusion of their experiments is that: the increase of the ambient pressure suppresses self-pulsation by decaying the dominant wave of the liquid sheet; recess is found to strengthen self-pulsation because the interaction between the liquid and gas phases becomes more active as the recess length increases. A strong ?ow rate oscillation is accompanied by the instability phenomenon, or self-pulsation, in a swirl coaxial injector; there is a close relationship between the ?uctuation of the injected propellant mass ?ow rate and the ?uctuation of heat release. Therefore, the aim of this study is to develop a simpli?ed method for calculating the dynamic characteristics of a coaxial swirl injector, and to investigate the in?uences of the injector con?guration and working conditions

FU AND YANG

145

Fig. 4 Interaction of dynamic processes in a recessed gas–liquid coaxial swirl injector.

Fig. 1 Schematic of a coaxial swirl injector.

dividing line gas phase area

the injectors (through couplings 6 and 8, respectively). The propellant supply system responds to these disturbances (2 and 4), and generates a ?ow ?uctuation which, when passing through the injector (1-5-9, 3-7-9), results in ?uctuations of droplet mass and size distributions of the combustible mixture spray. Feedback coupling 6 can form its own self-pulsation circuits, causing ?uctuations in the propellant ?ow at the injector exit [12]. In the current study, it is assumed that there is pressure ?uctuation only at the liquid ?ow tangential inlet; before coming into the recessed chamber, the gaseous ?ow is stable. Furthermore, Bazarov [8] reported that chamber pressure pulsation Pc0 in the self-pulsation case were suf?ciently weaker and did not exceed 10%Pc. Therefore, the pressure at the exit of the coaxial swirl injector is assumed to be constant, and only the dynamic route of the liquid phase (1-5-9-6) is studied. To study its dynamics, a recessed coaxial swirl injector can be divided into two components: the liquid swirl injector and the recessed chamber.

A. Linear Dynamic Characteristics of the Recessed Chamber

gas-liquid phase area

Fig. 2 Schematic of the inner mixing ?ow condition.

on the dynamics. This paper will emphasize the ?uctuation of the mass ?ow rate of the gas–liquid injector caused by pressure ?uctuation in the liquid propellant supply line.

II. Linear Dynamics of Recessed Gas–Liquid Coaxial Swirl Injector

The con?guration, as shown in Fig. 3, is a conventional gas–liquid coaxial swirl injector. In Fig. 3, liquid propellant comes to the swirl chamber through the tangential inlet, located at the swirling chamber wall; it forms a swirling ?ow, and is injected to the nozzle post located at the center of the coaxial swirl injector. Gaseous propellant is supplied through the annular gap between the recessed chamber wall and liquid injector. The exit surface of the inner liquid injector is located at a certain length on the inward side from the exit surface of the outer injector. Section 1-1 is at the tangential inlet of the liquid swirl injector, sections 2-2 and 3-3 are de?ned as the inlet and outlet of the recessed chamber, respectively. The diameter and length of the recessed chamber are dR and lR . The diameters of the tangential inlet, swirl chamber and spout ori?ce for the liquid swirl injector are dt , dk , and d0 respectively. Figure 4 shows the dynamic map of a recessed gas–liquid coaxial swirl injector. During internal mixing, the pressure ?uctuations of the combustion chamber affect the recessed chamber through the feedback coupling 10, and consequently, the liquid and gas stages of

1 2 3

lR d0

dk dt 2 1 gas

dR

liquid

3 x

Fig. 3 Schematic of a gas–liquid coaxial swirl injector with recessed chamber.

The following assumptions were made before the study: 1) The breakup length of the liquid sheet is ignored and the gas– liquid two-phase ?ow in the recessed chamber is regarded as homogeneous ?ow, i.e., the droplets and the gas phase are uniformly distributed in the recessed chamber. Under the action of centrifugal force, the liquid spreads out in the form of a conical sheet as soon as it leaves the central swirl injector. In the recessed chamber, the velocity of the gas is much higher than the velocity of the liquid phase. The liquid sheet disintegrates under the action of the gas ?ow. Considering the inner mixing ?ow condition af?rmed by the author [7], the recessed length is large enough so that the swirling liquid sheet will touch the inner surface of the outer gas port, and the gas– liquid phase area will occupy a large section of the recessed mixing chamber. In this case, the outer gas ?ow and the inner liquid sheet mainly interact in the recessed mixing chamber, and the ?ow is referred to as an inner mixing ?ow. Because both the gas ?ow and the inner surface of the outer port will act on the inner swirl liquid ?lm and contribute to atomization, we assume the primary atomization is accomplished in a very short distance and the liquid sheet is disintegrated into numbers of fragments and droplets. The liquid droplets disperse in the gas ?ow uniformly and pass through the recessed chamber with the gas phase. This assumption can be veri?ed to some extent by analyzing the spray pattern observed in experiments [4,13], where it was shown that the spray presents the shape of a fully developed cone when the co?ow gas is introduced, without recess. When the center injector recessed, the liquid sheet blocking the annular passage was blown off in the shape of a mushroom. The spray cone does not form until the gas–liquid mixture moves downstream for a distance. This indicates that the gas phase and liquid phase mix together in the recessed chamber. If the gas and liquid phase separate from each other, the spray cone would form immediately when they exit the spout exactly as occurs without the recess. The literature [14] also states: “For the developed ?ow pattern, mixing of the liquid and gaseous jet were promoted and the mixture ?ow occupied the recessed region. Therefore, only the mixture ?ow was injected from the injector exit.” Based on the preceding analysis, we think the homogeneous ?ow assumption can be used in the dynamic analysis for the recessed chamber, at least in the inner mixing condition.

146

FU AND YANG

2) Because the viscous force is considerably small compared with the inertial force in the present study (Reynolds number is large), the ?ow in the recessed chamber is considered to be inviscid. For example, a back-of-the-envelope calculation for an O=F (LOX=GH2 ) ratio of eight and a total mass ?ow rate of 0:3 kg=s gives an average velocity of 57 m=s and the equivalent density for the mixture is about 126 kg=m3 . Hence, the equivalent Reynolds number is about 4 105 (LOX ? 1:46e 4 kg=m s, H2 ? 8:99e 6 kg=m s [15]). A large Reynolds number means that the viscous force is considerably small and can be ignored. 3) The ?ow is one dimensional. Because the cross-sectional area of the recessed chamber is constant and the radial velocity is much smaller compared with the axial velocity, only the axial movement is considered. 4) Since the viscous can be negligible, there is no friction pressure drop. The pressure drop mainly contributes to the accelerating of the ?uid column. The mass ?ow rate of the two-phase mixture in the recessed chamber is M, and the density of the two-phase mixture is ? g ? l ?1 ?, where l , g are the density of the liquid and gas, respectively, and is the volumetric fraction of the gas phase. Because the ?ow is considered to be inviscid, the viscous term is neglected and the momentum equation in the recessed chamber can be written as @v @v ?v ? @x @t 1 @p @x (1)

the perturbation quantity of the liquid mass ?ow rate, v is the velocity of two-phase uniform ?ow. Then the relative perturbation quantity of mass ?ow rate of the two-phase ?ow can be expressed as M0 ? 0g 1 _ m0 ? Svl l (6)

And the ?ow in the recessed chamber is assumed to be isentropic: 0g ? p0 g p
(7)

where
is the ratio of speci?c heat of two-phase mixture. Because the pressure at the exit of the coaxial swirl injector is assumed to be constant, Eq. (7) can be written in this form [17]: 0g ? p0 g p
(8)

Because there is no heat transfer and friction, the gas phase will comply with pg

g

?C

(9)

where
g is the ratio of speci?c heat of gas phase, C is a constant. If there is a liquid phase of m unit mass in every unit mass of gas phase, this relationship in the two-phase mixture holds: g ? 1 ? m Substitute Eq. (10) to Eq. (9): p

g

where v is the velocity of the two-phase mixture,p denotes the pressure. Equation (1) is actually the Eulerian equation of motion, which applies to both the compressible and incompressible ?uid. For the reason of inviscid ?ow, the axial gradient for the velocity in Eq. (1) can be neglected, i.e., @v can be neglected [15]. Then dv ? @v, @x dt @t Eq. (1) transforms into dv 1 @p ? ?0 @x dt (2)

(10)

? C0

(11)

Integrating Eq. (2) along x direction, the pressure drop in the recessed chamber can be obtained: p2 p3 ? lR dM S dt (3)

It can be seen from Eq. (11) that the ratio of speci?c heat is identical for the mixture and gas phase under isentropic conditions. Combining the Eqs. (4–8), the following is obtained: _ MlR g p p0 m _ (12) s ? l m0l p 0 ? p2
pS Svl Thus, the transfer function of the liquid mass ?ow rate ?uctuation to the pressure drop in the recessed chamber can be expressed as R ? _ p0 sMlR ml
p2 ? ?S
p2 sMlR g ?Sl pv _ m0l (13)

Where the mass ?ow rate M ? Sv, S is the sectional area of recessed chamber. Equation (3) can also be found in section 3.2.2 of Harrje and Reardon’s report [16]. It is easier to study dynamic characteristics with a frequencydomain analysis than with a time-domain analysis method. The frequency response equation can be obtained by transforming the ?uid mechanics equation with a Laplace transform. Linearize Eq. (3) and the ?uctuation parameter is normalized by the mean variables: p 0 ? MlR 0 sM pS (4)

B.

Dynamic Model of a Gas–Liquid Coaxial Swirl Injector

Where s ? j! is the Laplace operator, ! ? 2f is the angular frequency, p ? p2 p3 , overbar denotes the relative quantity, and 0 denotes the perturbation quantity. According to the conservation of mass, the relation M ? Sv holds. Because the sectional area of recessed chamber is constant, the following equation can be obtained: M0 ? Sv0 ? Sv0 (5)

The motivation of this section is to deduce the transfer function for the entire coaxial swirl injector by combining the transfer function of the recessed chamber and the liquid injector. Because the ?ow in the recessed chamber is assumed to be homogeneous (analogous to the two-phase ?ow passing through an ori?ce plate), the relation between the pressure drop and the mixture mass ?ow rate is [18,19] p ? M2 2 2Ct "2 A2 t 0 (14)

Assuming liquid density is constant, 0 ? 0g (the gas phase volumetric fraction is so large that the variation of gas volumetric fraction caused by the variation of the liquid mass ?ow rate can be 1 _ neglected). v0 ? Sl m0l is obtained by the assumption that the gaseous _ ?ow before coming into the recessed chamber is stable, where m0l is

where Ct is the discharge coef?cient when the two-phase mixture ?owing through the ori?ce plate, "t is the expansion factor of the mixture, A0 is the area of the ori?ce. Similar to the mixture ?ow in the recessed chamber, Ct and "t are set to be one, the following equation holds: p???????????????? ~ ~ ~ M ? S 2p (15) Where the tilde indicates the ?uctuation parameter, and it can be expressed as the sum of steady value and perturbation value, i.e.,

FU AND YANG

147

0.04

Amplitude-frequency Phase-frequency

140 120 100 80

0.03

conditions (such as two-phase mass ?ow rate, mixing ratio, and property of the propellant), all in?uence the entire injector transfer function; therefore the effect of these parameters on the dynamics of a gas–liquid coaxial swirl injector will be discussed.

–?φ/°

0.02

Πinj

60 0.01 40 20 0 -0.01 0 1000 2000 3000 4000 -20 5000

III. Results and Discussion

The parameters of the baseline injector are as follows: the length and diameter of the recessed chamber are 24 and 7.3 mm; do , dk , and dt are 5.7, 8, and 1.6 mm, respectively. The mass ?ow rate and mixing ratio of the injector are 0.18 kg/s and 8. The number of the tangential inlet n equals 1. A is de?ned as the geometric parameter of the liquid swirl injector, which can be expressed as [23] ?dk dt ?do 2 ndt

0.00

f/Hz

Fig. 5

Dynamics of a liquid swirl injector.

A?

(23)

9 ~ M ? M ? M0 = ~ p ? p ? p0 ; ~ ? ? 0

(16)

Substituting Eq. (16) to Eq. (15), linearizing it and omitting the highorder quantity, Eq. (15) changes to M0 ? S2 ?p0 ? p0 ?=M (17)

The propellant is LOX, and gaseous hydrogen, whose properties are shown in Table 1. In this calculation, the steady pressure drop of the recessed chamber should be given ?rst. In [25], the pressure drop of the recessed chamber is 0.15 MPa when the mixing ratio O=F is six. In the current research, the order of magnitude of pressure drop in recessed chamber is set to be 0.1 MPa for simplicity.

A. Effect of Mixing Ratio

Because the liquid is assumed to be incompressible, the variation of density is only caused by the variation of gas density, and Eq. (8) is used again, i.e., 0 ? 0g ? p0 g p
(18)

The relationship between the relative perturbation quantity of mixture mass ?ow rate and that of the pressure drop in recessed chamber is M 0 ? p0 ? ? g p=p2
?S2 p=M2 From Eq. (13), the following can be derived: _ p 0 ? R m0l (20) (19)

_ where m0l is the relative perturbation of the ?ow rate of liquid swirl injector. According to the liquid injector dynamics theory [12], _ m0l ? inj p0inj , where inj is the transfer function of liquid swirl injector; p0inj is the relative perturbation of pressure at the liquid swirl injector inlet. For a convergent nozzle swirl injector, inj ? BX c k:cII T R2 a 2T ?k:3II ? k:3III ? ? 1 (21)

Where T is a complex transfer function of the tangential channel, k:cII and c are the transfer functions of the vortex chamber and spout nozzle, respectively, k:3II and k:3III are the complex transfer functions of the vortex chamber, which correspond to the two mechanisms of disturbance propagation from the tangential channel to the vortex chamber. A more thorough description is given by [12]. Many investigators [20–22] have used Eq. (21) to study the dynamics of a liquid swirl injector. Figure 5 shows the typical dynamic characteristics of a liquid swirl injector. It can be seen that the nondimensional amplitude of ?ow pulsation decreases and the phase contrast increases with the increase of pulsation frequency. Moreover, the undulation on the amplitude-frequency curve is caused by the re?ected wave in the swirl chamber. By combining Eqs. (13), (19), and (20), the transfer function of a gas–liquid coaxial swirl injector can be obtained, as follows: ? M0 =p0inj ? R inj ? ? g p=
p2 ?S2 p=M2 (22)

The mixing ratio is de?ned as the ratio of mass ?ow rate of the LOX and hydrogen. The baseline injector parameters are taken as the basic geometry, and four different mixing ratios are calculated. Figure 6 shows the results. Figure 6a is the amplitude-frequency diagram, and Fig. 6b is the phase-frequency diagram. It can be seen from the amplitude-frequency diagram that the amplitude response of the ?uctuation of two-phase mixture mass ?ow rate demonstrates a trend of increasing ?rst, and then decreasing in general, as the pulsation frequency increases. At some frequency point (f1 , f2 , f3 , etc., in Fig. 6), there are spikes of amplitude response. According to the literature [8,9], there exists self-pulsation phenomena in a recessed coaxial swirl injector. The frequency at which a sharp increase in the amplitude response is found is considered to be a selfpulsation frequency. When a weak disturbance occurs in the injector, it is composed of every monochromatic component with different frequencies and wavelengths, according to the dispersed wave theory. Just these components whose frequency equals f1 , f2 , etc., would be ampli?ed and cause the occurrence of self-pulsation. It is believed that this ?uctuation frequency response was generated by the coupling of the recessed chamber and the liquid swirl injector. In Fig. 5 small undulations can be seen at the frequency of about 600 Hz and 1300 Hz. In Bazarov’s theory of injector dynamics [26], these undulations were caused by the re?ected wave in the swirl chamber. The recessed chamber and the liquid swirl injector can form a time lag feedback. Through this feedback, these undulations may be ampli?ed. Hence, the ?uctuating frequency response presents the results, as shown in the calculation results. Moreover, it can be seen from Fig. 6 that the mixing ratio has great effect on the amplitude-frequency characteristics: the mixing ratio is larger, the amplitude response is greater, and the self-pulsation frequencies do not vary with the mixing ratio. Im and Yoon [9], pointed out that the frequency of the spray oscillation by selfpulsation without gas effects matches the frequency of the dominant wave of the liquid sheet, and this frequency is shifted slightly by the

Table 1 Element O2 H2

a

Injector ?ow, ?uid properties Temperature, K 100 150 Density, kg=m3 1119 15.4

Combustion chamber pressure, MPa 10.1 10.1

From Eq. (22) it can be seen that the injector con?guration parameters (such as recessed length, recessed diameter), working

a Critical properties of O2 : pc ? 5:04 MPa and Tc ? 155 K. Critical properties of H2 : pc ? 1:30 MPa and Tc ? 33:2 K [24].

148

0.05

FU AND YANG

600

f2

0.04

f3

f4

f5 f6 f7 f8

f1

0.03

O/F=3 O/F=6 O/F=9 O/F=12

400

O/F=3 O/F=6 O/F=9 O/F=12

??φ/°

0.02 0.01 0.00 0 2000 4000 6000 8000 10000

Π

200

0

-200

0

1000

2000

3000

4000

5000

f/Hz

f/Hz

a) Amplitude-frequency characteristics

b) Phase-frequency characteristics

Fig. 6 Dynamics of coaxial injector with different mixing ratio. C. Effect of Recessed Length

effects of the annular gas. Because the oscillation wave on the liquid sheet depends mainly on the liquid swirl injector, the mixing ratio has little effect on the self-pulsation frequency. It can be seen from the phase diagram that the phase contrast between the ?ow rate pulsation and pressure pulsation increases linearly as the pulsation frequency increases. As the mixing ratio increases, the phase contrast increases. However, the in?uence of the mixing ratio on the phase-frequency characteristics is small. The phase contrast occurs because the propagation speed of a pulsation wave in a two-phase mixture with different mixing ratios is distinct.

B. Effect of Mass Flow Rate

Taking the baseline injector parameter as the basic geometry, the effect of two-phase mass ?ow rate on the coaxial swirl injector dynamics is analyzed. Figure 7 shows the dynamics of the gas–liquid coaxial swirl injector with four different mass ?ow rates. Similarly, the amplitude response presents the trend of increasing ?rst, then decreasing in general, and the phase contrast increases as the pulsation frequency increases. A sharp and narrow increase of the amplitude response occurs at certain frequencies. Furthermore, as the mass ?ow rate increases, the amplitude response and the phase contrast both decrease. The in?uence of the mass ?ow rate on the phase-frequency characteristics is minor. The literature [8], reported that the model injectors tested at atmospheric antipressure showed that injectors with decreased propellant mass ?ow have a self-pulsating ?ow regime in a somewhat wide range of pressure drops of gas and liquid, accompanied by pressure pulsations in a test chamber, gas manifold and liquid feed line, coincident with this analysis. As the injector dynamics theory depicts, the amplitude response increases while the propellants mass ?ow rate decreases, so self-pulsation is easier to observe than pulsation under larger mass ?ow rate.

0.03

M=0.15kg/s M=0.2kg/s M=0.25kg/s M=0.3kg/s

To investigate the effect of the recessed chamber length on the dynamics of the coaxial injector, injectors with four different recessed chamber lengths are calculated. The inner liquid swirl injectors are recessed at 5, 10, 15 and 20 mm, which correspond to 0.7, 1.4, 2.1 and 2:8dR , respectively. Figure 8 shows the dynamics of the coaxial swirl injector with different recessed lengths. It shows that the amplitude response increases and the phase contrast decreases as the recessed length increases. This conclusion also agrees with the experimental results reported in [9]. Their experiments showed that self-pulsation is not observed with a short recess length, but as the recess length increases, self-pulsation is detected. The increased recess length quickens the occurrence of self-pulsation under the same injection conditions.

D. Effect of Ambient Pressure

Figure 9 shows the effect of ambient pressure on the dynamics of a coaxial swirl injector. The ambient pressures are 6, 8, 10.1 MPa. The amplitude-frequency diagram shows that as the ambient pressure increases, the relative amplitude of mass ?ow rate pulsation and the phase contrast both decrease; this complies with the phenomenon observed in the experiments [9]. Self-pulsation can be observed in many cases at the ambient pressure of 0.1 MPa; however, selfpulsation is detected only when gas velocity is high at the ambient pressure of 0.5 MPa. A further increase in ambient pressure leads to the suppression of self-pulsation. Moreover, the ambient pressure does not change the self-pulsation frequency.

E. Effect of Geometric Parameters of the Liquid Swirl Injector

As described in Sec. III.A, there is a relationship between the selfpulsation frequency and the dominant frequency of the unstable waves on the liquid ?lm. Thus, it can be deduced that the variation of

600

M=0.15kg/s M=0.2kg/s M=0.25kg/s M=0.3kg/s

400

0.02

??φ

0.01 0.00 0 2000 4000 6000 8000 10000

Π

200

0

-200

0

1000

2000

3000

4000

5000

f/Hz

f/Hz

a) Amplitude-frequency characteristics

Fig. 7

b) Phase-frequency characteristics

Dynamics of coaxial injector with different mass ?ow.

FU AND YANG

0.020

Lr=5mm Lr=10mm Lr=15mm Lr=20mm

149

600

Lr=5mm Lr=10mm Lr=15mm Lr=20mm

0.015

400

0.010

??φ

0 2000 4000 6000 8000 10000

Π

200

0.005

0

0.000

-200

0

1000

2000

3000

4000

5000

f/Hz

f/Hz

b) Phase-frequency characteristics a) Amplitude-frequency characteristics Fig. 8 Dynamics of coaxial injector with different length recessed chamber.

0.04

pc=6MPa pc=8MPa

500 400 300

pc=6MPa pc=8MPa pc=10.1MPa

0.03

pc=10.1MPa

0.02

??φ/°

200 100 0 -100

Π

0.01

0.00

0

2000

4000

6000

8000

10000

-200

0

1000

2000

3000

4000

5000

f/Hz

f/Hz

a) Amplitude-frequency characteristics

Fig. 9

b) Phase-frequency characteristics

Dynamics of coaxial injector with different ambient pressure.

0.04 f1(A=10) 0.03 f1(A=12.5) 0.02 f1(A=14.5) 0.01

A=10 A=12.5 A=14.5

600 500 400

A=10 A=12.5 A=14.5

??φ

300 200 100 0

Π

0.00

0

1000

2000

3000

4000

5000

-100

0

1000

2000

3000

4000

5000

f/Hz

f/Hz

a) Amplitude-frequency characteristics

b) Phase-frequency characteristics

Fig. 10 Dynamics of coaxial injector with different geometric parameter A.

liquid swirl injector geometry would cause variation of the liquid ?lm in the swirl injector and the waves on the liquid ?lm, further leading to variation of the self-pulsation frequency of the gas–liquid coaxial swirl injector. Here we investigate the effect of liquid swirl injector con?guration on the dynamics of the coaxial injector by changing the geometric parameter of the liquid swirl injector A from 10 to 14.5. Figure 10 shows the calculated results. It can be seen from the amplitude diagram that the amplitude response decreases as the liquid swirl injector’s geometric parameter A increases. This is because A denotes the swirl intensity of the liquid sheet. Increasing A causes the liquid ?ow area to decrease. A smaller liquid ?ow area decreases the collision of the liquid and gas phase, further decreasing the self-pulsation amplitude. Moreover, increasing A would decrease

the amplitude of the liquid mass ?ow rate pulsation according to the liquid swirl injector dynamics theory, further decreasing the amplitude of the two-phase mixture mass ?ow rate. Furthermore, with the increase of A, pulsation frequency decreases. As shown in Fig. 6, the ?rst self-pulsation frequency decreases from 1192 to 719 Hz. The phase-frequency diagram shows that the phase contrast increases somewhat as A increases.

IV. Conclusions

A theoretical model of the linear dynamics of a gas–liquid coaxial swirl injector was developed. The model considers the in?uence of the liquid phase only, and divides the dynamic process of the coaxial

150

FU AND YANG

swirl injector into two parts: the ?ow in the liquid swirl injector and the ?ow in the recessed chamber. The linear transfer function of the recessed chamber was deduced based on fundamental ?uid mechanics equations. The transfer function of the entire gas–liquid coaxial swirl injector was obtained by combining the transfer function of the recessed chamber and the liquid swirl injector. The transfer function depicts the characteristics of the coaxial swirl injector in the frequency domain. Based on the deduced transfer function, the in?uences of the con?guration parameters on the injector dynamics were calculated. The qualitative trend of the predicted dynamics agrees reasonably well with the observations reported in the literature. The effects of the parameters were analyzed and the conclusions were obtained as follows: 1) The amplitude response of the ?uctuation of two-phase mixture mass ?ow rate demonstrates a trend of increasing ?rst, and then decreasing in general, as the pulsation frequency increases. The phase contrast between the relative pulsation of mass ?ow rate and injection pressure increases linearly as the pulsation frequency increases. 2) At some frequency point, spikes of amplitude response occur, which are considered to be the self-pulsation frequency, and the selfpulsation frequency does not vary with other con?guration parameters or working conditions, except the geometric parameter of the liquid swirl injector. With the increased geometric parameter, the self-pulsation frequency decreases. 3) A smaller mixing ratio, larger mass ?ow rate, higher ambient pressure, smaller ratio of recess length-to-diameter, and larger geometric parameter, would all serve to decrease the amplitude response of the mass ?ow rate pulsation. But these parameters have little effect on the phase frequency of the coaxial swirl injector.

[8] [9] [10] [11] [12]

[13] [14]

[15] [16] [17] [18] [19] [20]

Acknowledgments

The ?nancial support of China National Nature Science Funds (Support number: 50406007) and the Innovation Foundation of Beijing University of Aeronautics and Astronautics for Ph.D. Graduates is gratefully acknowledged.

[21]

References

[1] Gill, G. S., “Liquid Rocket Engine Injectors,” NASA SP-8089, 1976. [2] Inamura, T., Miyata, K., Tamura, H., and Sakamoto, H., “Spray Characteristics of Swirl Coaxial Injector and its Modeling,” AIAA Paper 2001-3570, 2001. [3] Cohn, R. K., Strakey, P. A., Bates, R. W., Talley, D. G., Muss, J. A., and Johnson, C. W., “Swirl Coaxial Injector Development,” AIAA Paper 2003-125, 2003. [4] Tharakan, T. J., and Gupta, N. K., “Effect of Recess on the Spray Characteristics of Gas–Liquid Coaxial Swirl Injectors,” Proceedings of the International Conference on Aerospace Science and Technology, Bangalore, India, 2008. [5] Soller, S., Wagner, R., Kau, H., Martin, P., and M?ding, C., “Combustion Stability Characteristics of Coax-Swirl-Injectors for Oxygen/ Kerosene,” AIAA Paper 2007-5563, 2007. [6] Zhang, M., Wang, L., Zhang, Z., and Yang, G., “Flame Characteristic Test Study of Coaxial-Swirl Injector,” Combustion Science and Technology, Vol. 14, No. 1, 2008, pp. 1–5. [7] Yang, L., Ge, M., Zhang, M., Fu, Q., and Cai, G., “Spray Characteristics [22] [23] [24]

[25] [26]

of a Recessed Gas–Liquid Coaxial Swirl Injector,” Journal of Propulsion and Power, Vol. 24, No. 6, 2008, pp. 1332–1339. doi:10.2514/1.23977 Bazarov, V., “Self-Pulsations in Coaxial Injectors with Central Swirl Liquid Stage,” AIAA Paper 1995-2358, 1995. Im, J.-H., and Yoon, Y., “The Effects of the Ambient Pressure on SelfPulsation Characteristics of a Gas/Liquid Swirl Coaxial Injector,” AIAA Paper 2008-4850, 2008. Andreyev, A. V., and Bazarov, V. G., Dynamics of Gas–Liquid Injectors, Mashinostroenie, Moscow, 1991 (in Russian). Andreyev, A. V., and Chepkin, V., “Autovibration of Coax Injector Elements of O2 =H2 Staged Combustion Liquid Rocket Engines,” AIAA Paper 1995-2837, 1995. Bazarov, V., and Yang, V., “Liquid-Propellant Rocket Engine Injector Dynamics,” Journal of Propulsion and Power, Vol. 14, No. 5, 1998, pp. 797–806. doi:10.2514/2.5343 Sasaki, M., Sakamoto, H., Takahashi, M., and Tomita, T., “Comparative Study of Recessed and Non-Recessed Swirl Coaxial Injectors,” AIAA Paper 1997-2907, 1997. Nunome, Y., Tamura, H., Onodera, T., Sakamoto, H., Kumakawa, A., and Inamura, T., “Effect of Liquid Disintegration on Flow Instability in a Recessed Region of a Shear Coaxial Injector,” AIAA Paper 20095389, 2009. Luo, Z., Fluid Network Theory, China Machine Press, Beijing, 1988, p. 207 (in Chinese). Harrje, D., and Reardon, F., “Liquid Propellant Rocket Combustion Instability,” NASA SP-194, 1972, p. 108. Feiler, C. E., and Heidmann, M. F., “Dynamic Response of GaseousHydrogen Flow System and its Application to High-Frequency Combustion Instability,” NASA TN D-4040, 1967. Yan, C., Gas–Liquid Two Phase Flow, Harbin Engineering Univ. Press, Harbin, China, 2010 (in Chinese). Upp, E. L., and LaNasa, P. J., Fluid Flow Measurement, A Practical Guide to Accurate Flow Measurement, Gulf Professional, Woburn, MA, 2002. Fu, Q., Yang, L., and Wang, X., “Theoretical and Experimental Study of the Dynamics of a Liquid Swirl Injector,” Journal of Propulsion and Power, Vol. 26, No. 1, 2010, pp. 94–102. doi:10.2514/1.44271 Seo, S., Kim, S., and Choi, H., “Combustion Dynamics and Stability of a Fuel-Rich Gas Generator,” Journal of Propulsion and Power, Vol. 26, No. 2, 2010, pp. 259–266. doi:10.2514/1.46568 Ismailov, M., and Heister, S., “Nonlinear Modeling of Classical Swirl Injector Dynamics,” AIAA Paper 2009-5402, 2009. Bayvel, L., and Orzechowski, Z., Liquid Atomization, Taylor and Francis, Washington, D.C., 1993. Oefelein, J. C., “Thermophysical Characteristics of Shear-Coaxial LOX-H2 Flames at Supercritical Pressure,” Proceedings of the Combustion Institute, Vol. 30, No. 2, 2005, pp. 2929–2937. doi:10.1016/j.proci.2004.08.212 Sun, J., Zhuang, F., and Wang, J., “Effects of Recess on Coaxial Injector’s Discharge Coef?cient and Performance,” Journal of Propulsion Technology, Vol. 24, No. 5, 2003, pp. 452–455. Bazarov, V., Dynamics of Liquid Injectors, Mashinostroenie, Moscow, 1979 (in Russian).

D. Talley Editor-in-Chief

相关文章:

更多相关标签:

- Linear Stability Analysis of a Conical Liquid Sheet
- A New Nonlinear Liquid Drop Model. Clusters as Solitons on The Nuclear Surface
- Measurement of annular liquid film thickness in an open-end swirl injector
- A MOLECULAR DYNAMICS STUDY OF INTERFACE PHENOMENA OF A LIQUID DROPLET
- satellite formation and merging in liquid jet breakup
- nonlinear breakup of coaxial liquid jet in a swirling gas stream
- An Experimental and Computational Study of Fluid Dynamics in Gas- Liquid Chemical Reactors
- Measurement of annular liquid film thickness in an open-end swirl injector
- CFD Modeling of Gas-Liquid-Solid Fluidized Bed-NIT
- Experimental and theoretical investigation of the dynamics of the SCR - DeNOx reaction
- Bubble formation and dynamics in gas–liquid–solid fluidization—A review
- Two-phase gas–liquid flow characteristics inside a plateheatexchanger
- The New Gas Stripper and Charge State Separator of the Gsi High Current Injector
- The Glass Transition and Liquid-Gas Spinodal Boundaries of Metastable Liquids
- gas PIN injector d5 fns50