Docsity
Docsity

Prepare-se para as provas
Prepare-se para as provas

Estude fácil! Tem muito documento disponível na Docsity


Ganhe pontos para baixar
Ganhe pontos para baixar

Ganhe pontos ajudando outros esrudantes ou compre um plano Premium


Guias e Dicas
Guias e Dicas

Hydrogen-Oxygen Reaction Kinetics in Rocket Engine: Induction & Postinduction Times, Manuais, Projetos, Pesquisas de Química

An analysis of hydrogen-oxygen chemical reaction kinetics in rocket engine combustion, focusing on induction and postinduction times. The study calculates reaction times and concentration histories for various rocket engine conditions, identifying the importance of temperature and oxidant-fuel weight ratios on induction times. The document also discusses the limitations of the analytical method and the significance of chemical reaction times in rocket combustion.

O que você vai aprender

  • What are the limitations of the analytical method for calculating hydrogen-oxygen reaction times?
  • How do induction and postinduction times impact rocket engine performance?
  • What is the role of oxidant-fuel weight ratios in hydrogen-oxygen reaction times?
  • How does temperature influence hydrogen-oxygen reaction times?
  • What are the key factors affecting hydrogen-oxygen chemical reaction times in rocket engine combustion?

Tipologia: Manuais, Projetos, Pesquisas

2020

Compartilhado em 22/05/2020

rickvenzi
rickvenzi 🇧🇷

2 documentos

1 / 20

Toggle sidebar

Esta página não é visível na pré-visualização

Não perca as partes importantes!

bg1
C
-4
L
NASA TECHNICAL
NOTE
0
Ih
N
P
n
z
c
4
w
4
z
LOF\!l
A
KI
RTI
HYDROGEN-OXYGEN CHEMICAL
REACTION KINETICS
IN
ROCKET ENGINE COMBUSTION
by
Martia Hersch
Lewis
Research Center
Cleveland,
Ohio
NASA
TN
D-42
c,
COPY:
RETURN
TC
,AND
AF5,
N
MEX
.FWL
[WLBL-Z)
50
_-
--I
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION WASHINGTON, D.
C.
DECEMBER
1967
/
https://ntrs.nasa.gov/search.jsp?R=19680002174 2020-05-13T23:31:04+00:00Z
pf3
pf4
pf5
pf8
pf9
pfa
pfd
pfe
pff
pf12
pf13
pf14

Pré-visualização parcial do texto

Baixe Hydrogen-Oxygen Reaction Kinetics in Rocket Engine: Induction & Postinduction Times e outras Manuais, Projetos, Pesquisas em PDF para Química, somente na Docsity!

C^ -4^ L

NASA TECHNICAL N O T E

(^0) Ih N Pn

zc

4 w

4 z

LOF!l

A

KI RTI

HYDROGEN-OXYGEN CHEMICAL

REACTION KINETICS I N

ROCKET ENGINE COMBUSTION

by Martia Hersch

Lewis Research Center

Cleveland, Ohio

N A S A TN D-

c ,

COPY: RETURN TC

,AND AF5, N MEX

.FWL [WLBL-Z)

(^50) _ -

--I

N A T I O N A L A E R O N A U T I C S A N D SPACE A D M I N I S T R A T I O N W A S H I N G T O N , D. C. DECEMBER 1 9 6 7

/

https://ntrs.nasa.gov/search.jsp?R=19680002174 2020-05-13T23:31:04+00:00Z

TECH LIBRARY KAFB, NY

NASA T N D-

HYDROGEN-OXYGEN CHEMICAL REACTION KINETICS

IN ROCKET ENGINE COMBUSTION

By M a r t i n H e r s c h

L e w i s R e s e a r c h C e n t e r C l e v e l a n d , Ohio

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION

For s o l e by the Cleoringhouse for Federal Scientific and Technical Information Springfield, Virginia 22151 - CFSTI price $3.

I I I I I 11.11 I ,^111.^ I.^ ,,...-.,-^ I .I.^. .I^ ..^.^ ..^ ..^.^.^.^ -^.^ .._..^.^ ..^ ..^ ..

the assumption of infinitely fast reaction rates under these conditions may also be

invalid.

Calculations of hydrogen-oxygen (H-0) kinetics for reaction limited ram -jet and

rocket combustor conditions are presented in references 5 to 7. In these studies the

following reaction mechanism was used:

H + O2 = O + OH I

0 + H2 = H + OH 11

H 2 + M = 2 H + M IlI

H + OH + M = H20 + M Iv

0 2 + M = 2 0 + M V

H2 + OH = H + H 2 0 (^) VI

where M is any third body.

In references 5 to 8 the initial temperature was varied from approximately 1000° to

3000 ' K at pressures of 1 to 65 atmospheres (lX105 to 6. 5 ~ 1 0 ~N/sq m). Brokaw, using

a n analytical solution for reactions I, 11, VI, VII, and VIII (ref. 9), shows that the

following reactions are also important, particularly at low temperatures:

HO2 + H2 H202 + H VII

H + 0 2 + M = H 0 2 + M VIII

The importance of these two reactions increases rapidly as the initial temperature de-

creases. It would therefore be desirable to consider a mechanism which includes all

eight reactions.

ical and analytical techniques are used. The effects of the initiation reaction

H2 + O2 - 2 OH is also considered with the use of the analytical technique.

species, and temperature change during reaction. The analytical method is restricted to constant H2 and O2 concentrations, nonreversible reactions, and constant temperature.

Reaction times are presented as functions of initial temperature and oxidant-fuel weight

ratios O/F at a constant pressure. Histories of atom and radial concentrations during

Calculations for this system of eight reactions are presented in this study. Numer-

The numerical method considers reversible reactions, concentration changes of all

the reaction are also shown. Calculations were made for initial temperatures ranging

from 1200' to 2500' K and O/F values of 1 and 10 at a chamber pressure of 20 atmos-

pheres (2. OX106 N/sq m).

C A LC U LAT I ON S

Numerical Solution

The numerical integration program used for this study is that developed under NASA

(1) The mixture is one of thermally perfect gases. (2) Flow is inviscid throughout. (3) Transport properties can be neglected.

(4) Internal energies are in thermal equilibrium.

(5) The law of m a s s action applies throughout. Backward r a t e s a r e calculated, with microscopic reversibility assumed, from the

ratio of forward rate constant to equilibrium constant.

state a r e solved simultaneously with the differential equations for specie concentration changes due to chemical reactions.

tions, forward rate constants, and equilibrium constants. Parameters for the forward

rate constants are given in table I, and parameters for the equilibrium constants in

table II. The equilibrium constants were calculated from the data of reference 10 and

similar unpublished Lewis data.

integration program a r e

contract and described in reference 8. The assumptions made in this program are

Equations of mass, energy, and

Input data for the numerical program are temperature, pressure, species concentra-

The reactions, including those proposed by Brokaw (ref. 9), used for the numerical

H + 0 2 = O + O H In

0 + H2 = H + OH (^) IIn

2 H + M = H 2 + M IIIn

H + O H + M = H 2 0 + M N n

2 0 + M =02 + M (^) Vn

H2 + OH = H + H 2 0 VIn

H 0 2 + H2 = H + H

H + 0 2 + M (^) H 0 2 + M

VIIn

VIIIn

A n a l y t i c a l S o l u t i o n

An analytical solution is obtained with isothermal conditions and constant reactant

(H2 and 02) concentration for the system assumed:

H + 0 2 -0^ K1 + OH

0 + H2& H + OH

H 2 + M - 2 H + M K

0 2 + M - 2 0 + M Kg

la

IIa

IIIa

Va

Ha + OH -HKg + H 2 0 VIa

HOa + H2- K7 H + H202 (^) VIIa

H + 0 2 + M-H02+ K8 M VIILa

H O a + M - H + O Z + M K- 8 -VIIIa

Reactions IIIa and Va are written in the direction of dissociation because during initiation

and induction radical formation is of concern. Reaction lVn is omitted from the analyti-

cal solution because it results in a term containing products of radicals, which results in

a nonlinear equation. This reaction, however, is of little importance during the induction

period. Reaction VLTI is considered in both directions in the analytic solution.

(for the case of reaction VIII in the forward direction)

The differential equations expressing the growth rate of intermediate species are

--dCH - - K I C O p H + K Z C H C O+ 2K3CHpM dt

H02 = - K C 7 H a C (^) H 0 2 + K 8 C 0 2c cM H

dC

d t

  • dCOH = K C C + K C C -KgCH2COH d t 1 0 2 H^ 2 H 2 O

where the C is concentration, K is a reaction rate constant, and t is time.

The general solution for these differential equations is (ref. 9)

hit

C 0 = Aie + a o

i= 1

hit c H = Bie + bH i= 1

hit = C Cie +- 'H H02 i= 1

C

hit COH = Die + dOH i= 1

where Ai, Bi, Ci, and Di are constants with dimensions of concentration, and h is a

constant with dimensions of second-'. The t e r m s a,, bH, cH02, and dOH, which result from initiation reactions, are

given in reference 9. For example,

1111 I 11111111111 l I I I I I

increase, and the system attains equilibrium. This f i n a l period may be referred to as

the postinduction period.

calculations cannot be started at the t r u e beginning of the reaction, with the assumption

that the reaction starts with zero atom and radical concentrations. On the other hand,

the analytical program can be started for the conditions of zero atom and radical con-

centrations at zero time. The analytic method, however, gives no information which

can be used to define the end of induction.

initial temperatures is shown in figure 1. Temperature change is also shown. The

solid and dashed curves represent numerical and analytical results, respectively. The

lower left ends of the curves show the very short period of initial formation of atoms

and radicals. The curves then show the relatively long induction period characterized by constant exponential growth under isothermal conditions. The numerical results show

that the induction period, as defined by constant exponential growth, continues, in some

The numerical and analytical methods are somewhat complementary. The numerical

The OH concentration plotted against time during initiation and induction for various

10

(^0) YIW 0

'I-

Y m? (a) Temperature, 1400° K.

-8 -6 -4 -2 0 2x10- Time, see (b) Temperature, 1600" K. Figure 1. - OH growth and temperature rise.

500-

400

300

200

100- Q

10-

10-

10-

10-

10-l

10-16-10 -8 -6 -4 -2 0 2

3 V

  • VIal 0 - = L z cn mL m E al

I^3 .- c 5-

$ 10-

;10-

al CL I-^ al^ V cO 0 10-

10-

100 10-l

(c) Temperature, Moo" K.

-3 -2 -1 0 1 2 3xW Time, sec (d) Temperature, 2500" K, numerical solution. Figure 1. - Concluded.

cases, well into the region of temperature increase. In this study the induction period is considered ended when the atom and radical concentrations deviate from constant exponential growth. This deviation can be determined, in this study, by the numerical method only. Although the analytical assumptions are violated when the system is

nonisothermal, the limit to which the analytical results can be continued is the point at

which the concentration curves deviate from a straight line, as shown in figure 1. This

point has therefore been selected as the end of induction and is plotted at zero time.

Numerical computation times increased drastically with decreasing temperature.

Therefore, at 1200' K numerical results were obtained only near the end of induction

and thus are not shown in figure 1.

The OH concentration at the start and end of induction are shown in figure 2. The

concentration at the start of induction was determined by the intersection of the ex-

trapolated slopes of the initiation and induction portions of the curves, such as shown

in figure 1 (a). The induction period is often considered ended when the OH concentration reaches the threshold of experimental detectibility, 10" gram-mole per cubic centi-

9

t

meter (ref. 9). This criterion appears to be suitable only for a n initial temperature near 1500' K. The induction time may now be determined by using these limits of OH concentration shown in figure 2.

Since in this study the total reaction time, rather than only induction time, is im-

portant, the postinduction time must also be investigated. Some idea of the post- induction time may be obtained from the curves of figure 3, which shows the temperature

rise for various initial temperatures as a function of time. These results show that for

the temperature range of about 1500' to 2500' K the postinduction time is on the order

of a^ few microseconds.^ Though only^ partial^ results^ are^ shown f o r 1400'^ K,^ apparently

at this temperature, the postinduction time may be on the order of perhaps 30 micro-

seconds.

10-

10- U av), .-^ ai^ E

10-

10- 1200

\
\

1400 1600

. , , ,. I I. , Time.. Sdlution Oxidani-'fuel ratio, 01 F Induction Numerical 1

  • _ _ _ - Induction Analytical 1 Induction Analytical 10 it-- Postinduction Numerical 1

1800 2000 2 m Temperature, "K Figure 4. - Effect of initial temperature on induction and postinduction time.

The induction and approximate postinduction times are shown in figure 4. Numerical

and analytical methods show nearly identical results for the induction time. initial temperature higher than about 1700° K, the induction time is small in comparison with the postinduction time and therefore may be neglected. At lower temperatures the

postinduction time, though, becomes increasingly small in comparison with the induction

time. The induction time increases rapidly with decreasing temperature for temper - a t u r e s below 1600' K. In this region the induction time increases exponentially with the reciprocal of initial temperature. The effect of O/F on induction time was investigated by determining the OH con-

centration at the end of induction f o r a n O/F of 10 by using the numerical program.

Calculations at an O/F of 10 were made for temperatures between 1500' and 2500 ' K.

Detailed calculations at this O/F f o r the entire induction period or f o r lower temper-

atures were not made because of the long computational times. The OH concentration

at the end of induction for a n O/F of 10 is shown in figure 2. At high temperatures the

OH at the end of induction at an O/F of 10 is about 10 times that at a n O/F of 1. At

For a n

temperatures lower than 1700 ' K the OH concentration at the end of induction is approx-

imately the same for both O/F conditions.

The induction time based on these numerically calculated OH concentrations at the

end of induction for a n O/F of 10 is shown in figure 4. The curve was extended to

1200' K by assuming that the OH concentration at the end of induction for an O/F of 10

was the s a m e , a s for a n O/F of 1. These results indicate that at higher temperatures

the induction time decreases when the O/F is increased. At lower temperatures the

effect of O/F on induction time appears to be small.

c 0 c^ mL c al U c 0 U a 5 l CL^ m al

.- L

.- c

Initial temperature, "K (a) Numerical solution.

ized t o OH concentration).

(b) Analytical solution. Figure 5. - Relative product species d u r i n g induction period (normal-

Some idea of the composition during induction is given by figure 5. During initia-

tion relative concentrations of the atom and radical species shift with time. During

induction, however, the relative concentrations remain nearly constant with time.

Relative concentrations of 0, H, H 0 2 , H 2 0 , and H 2 0 2 a r e shown, relative to OH

concentration, during induction as a function of temperature. Almost perfect agreement

is obtained between the numerical and analytical results. The relative concentrations

of 0 and H appear nearly constant during induction over most of the temperature range.

Concentrations of H 0 2 , H 2 0 , and H 2 0 2 a r e much greater than OH at low temperatures.

At higher temperatures the concentrations of 0, H, H 0 2 , and H 2 0 a r e one or two orders

of magnitude greater than OH. At 2500' K the H 2 0 2 concentration during induction is

nearly two orders of magnitude less than that of OH.

with initial concentrations of atoms and radicals determined analytically, because the

numerical program cannot be started with zero atom and radical concentrations. The effect of starting numerical calculations with arbitrary atom and radical concentrations

is shown in figure 6. The dashed lines are concentration histories when the starting

concentrations were calculated analytically. The solid curves a r e concentration

histories for arbitrary initial concentrations. These results show that if the initial

All numerical results presented thus f a r were obtained by starting the computations

ll11llll1l11 I 1 I 1 I I I I 1 II I

0

  • i=E 10-8LLlillllliil 1200 1500 2Mw)

1

Temperature, "K

Figure 7. - Effect of reaction H2 + 02- 20H

o n t i m e to attain a r b i t r a r y OH concentra-

t ion (ana lyt ica I ca I c ulat ion s 1. Oxidant -

f u e l ratio. 1.

induction time would have no effect on the overall reaction time. This reaction might be

of more importance at lower temperatures. It can be shown, however, that for rocket

combustion studies this is not a significant effect.

SIGNIFICANCE OF CHEMICAL REACTION

TIMES IN ROCKET COMBUSTION

Liquid -propellant rocket -engine combustion involves many steps, atomization,

vaporization, gas phase mixing, and chemical reaction. The importance of any steps

in the combustion process may be related to its speed o r completion time.

steps will be more important than the rapid ones. The combustion dead time, a n in-

dication of the time required for the overall process, is reported in reference 2 for a

typical H2-02 rocket engine. It varied from 0.25 to 2.25 milliseconds for a wide range of

operating conditions. High-frequency instability must also be considered. Reference 1 shows that wave times for high-frequency instability range from about 0.1 to 1 milli- second. Thus, the time range of interest in liquid rocket combustion ranges from 0. to several milliseconds. Therefore, if chemical reaction times become greater than about 0.1 millisecond, they must be considered in studying rocket combustion.

ing reactant temperatures if the temperature is below about 1700° K. At 1600' K the

reaction time is about 10 microseconds, which may be considered infinitely fast in com-

parison with other steps in rocket combustion. At about 1400' to 1500' K the reaction

time increases to about

The slower

This study shows that chemical reaction time increases very rapidly with decreas-

second and thus enters the region of interest in rocket

combustion. Above 1700° K the time decreases to several microseconds, and is fairly

insensitive to temperature.

in the combustion chamber. Nor can it be determined here what the gas temperature is

at the start of reaction. Presumably, though, the propellants in actual rocket combus-

tors a r e substantially heated soon after injection by convection, radiation, and recircu- lation. The results of this study indicate that if the unreacted propellants are rapidly

heated to .about 1400' K, then chemical kinetics can be neglected as a rate controlling

mechanism in steady-state combustion for a chamber pressure of 20 atmospheres

(2. OX106 N/sq m).

s u r e gradients. If the instability wave reduces the temperature to below about 1500' K, then reaction times would greatly increase and become approximately equal to the wave period. Under these conditions the reaction kinetics could become an instability driving force.

It is beyond the scope of this study to investigate the heating of unreacted propellants

Combustion instabilities, however, are known to cause large temperature and pres -

S U M M A R Y OF RESULTS

Hydrogen-oxygen reaction times and concentration histories were calculated for rocket combustor conditions. Calculations were made for oxidant-fuel weight ratios

of 1 and 10 and initial reactant temperatures ranging from 1200' to 2500' K at a

chamber pressure of 20 atmospheres (2.0~10 N/sq m).^ A numerical integration pro-

gram and a programmed analytical solution were used. The following principal results

were obtained:

  1. Hydrogen-oxygen reaction times varied inversely with initial reactant temper -

ature, ranging from a few microseconds at 2500' K to approximately 0.01 second at

1200' K. Since processes requiring times less than milliseconds can be neglected in determining performance efficiencies o r combustion characteristics this study indicates

that chemical reactions can be neglected if the injected propellants are heated above

approximately 1400' K.

2. When the initial propellant temperature was higher than about 1700' K, the re-

action time w a s primarily the postinduction or heat-release period. At lower temper-

atures the reaction time was the isothermal induction time.

of the oxidant-fuel mixture ratio on induction time increased with increasing temper- ature.

3. Induction time decreased with increasing oxidant-fuel mixture ratio. This effect

8. Zupnik, T. F. ; Nilson, E. N. ; and Sarli, V. J. : Investigation of Nonequilibrium Flow Effects in High Expansion Ratio Nozzles. Computor Program Manual. Rep. No. UACRL-C910096-11 (NASA CR-54042), United Aircraft Corp., Sept. 15, 1964.

  1. Brokaw, Richard S. : Analytic Solutions to the Ignition Kinetics of the Hydrogen- Oxygen Reaction. NASA TN D-2542, 1964.

Thermodynamic Properties to 6000' K for 210 Substances Involving the First 18 Elements. NASA SP-3001, 1963.

  1. McBride, Bonnie J. ; Heimel, Sheldon; Ehlers, Janet G. ; and Gordon, Sanford:
  2. Westenberg, A. A. : Hydrogen-Air Chemical Kinetic Calculations in Supersonic Flow. Rep. No. CM-1028, Appl. Phys. Lab., Johns Hopkins Univ., Dec. 1962.
  3. Getzinger, R. W. ; and Schott, G. L. : Kinetic Studies of Hydroxyl Radicals in Shock Waves. V. Recombination Via the H+02+M-H02+M Reaction in Lean Hydrogen-Oxygen Mixtures. J. Chem. Phys. , vol. 43, no. 9, Nov. 1, 1965, pp. 3237-3247.
  4. Ripley, Dennis L. ; and Gardiner, W. C., Jr. : Shock-Tube Study of Hydrogen- Oxygen Reaction. II.^ Role of Exchange Initiation.^ J.^ Chem.^ Phys.,^ vol.^ 44, no. 6, Mar. 15, 1966, pp. 2285-2296.

NASA-Langley, 1961 - 21 E-4010 17

"The aeronautical and space activities of the United States shall be conducted so as to contribute... to the expansion of human knowl- edge of phenomena in the atmosphere and space. The Administration shall provide for the widest practicable and appropriate dissemination of information concerning its activities and the results tfiereof." -NATIONALA E R o N A u n c s AND SPACEA C T OF 1958

NASA SCIENTIFIC A N D TECHNICAL PUBLICATIONS

TECHNICAL REPORTS: Scientific and technical information considered important, complete, and a lasting contribution to existing knowledge. TECHNICAL NOTES: Information less broad in scope but nevertheless of importance as a contribution to existing knowledge. TECHNICAL MEMORANDUMS: Information receiving limited distribu- tion because of preliminary data, security classification, or other reasons. CONTRACTOR. REPORTS: Scientific and technical information generated under a NASA contract or,grant and considered an important contribution to existing knowledge. TECHNICAL TRANSLATIONS: Information published in a foreign language considered to merit NASA distribution in English. SPECIAL PUBLICATIONS: Information derived from or of value to NASA activities. Publications include conference proceediags, monographs, data compilations, handbooks, sourcebooks, and special bibliographies. TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- nology used by NASA that may be of particular interest in commercial and other non-aerospace applications. Publications include Tech Briefs, Technology Utilization Reports and Notes, and Technology Surveys.

Details on the availability of these publications may be obtained irom:

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION

Washingon, D.C. PO