Reaction Kinetics of Ammonia & Nitric Acid
Reaction Kinetics of Ammonia & Nitric Acid
Reaction Kinetics of Ammonia & Nitric Acid
o , f'
^A't^^.,
.."'- J-
Approyad
August, 1973
f\
fiit-5057
r3
Cop. a
ACKNOWLEDGEMENTS
11
TABLE OF CONTENTS
ACKNOWLEDGEMENTS ii
LIST OF TABLES V
LIST OF FIGURES vi
CHAPTER
I. INTRODUCTION AND PRELIMINARY INVESTIGATIONS 1
II. LITERATURE SURVEY 10
III. THEORY 26
IV. EXPERIMENTAL SYSTEMS 39
The Backmix Reactor System 39
Backmix Reactor Experimental Procedure 47
Fixed-Bed Reactor System 48
Fixed-Bed Reactor Experimental Procedure 51
Chemical Analysis Methods 53
Ammonia Analysis 53
Chromatographic Ammonia Analysis 54
Colorimetric Ammonia Analysis and 55
Calibration
Nitrogen Dioxide Analysis 61
Nitrous Oxide Calibration and Analysis 62
V. EXPERIMENTAL DATA AND NUMERICAL PROCEDURES 65
VI. RESULTS OF EXPERIMENTAL WORK 81
Long-Time Data 81
Short-Time Data 86
Backmix Reactor 98
111
AM^
IV
7. Calibration Apparatus 59
19. Regeneration I 95
20. Regeneration II 96
21. Regeneration V 99
VI
CHAPTER I
INTRODUCTION AND PRELIMINARY INVESTIGATIONS
kcal/mole
Temp. ** R Heat of Reaction K
Reaction I
Reaction II
Reaction III
Reaction IV
Surface Area
Catalyst Material Company Designation
2% manganese dioxide Harshaw Mn-0501 258
supported on silica
0.3% palladium Harshaw Pd-0501 186
supported on high
activity alumina
0.5% platinum Englehard 0.5% Pt 93
supported on alumina supported on alumina
0.5% ruthenium Englehard 0.5% Pt 74
supported on aliomina supported on alumina
4% nickel oxide, 4% Harshaw Ni-1601 78
cobalt oxide, 4% iron
oxide supported on
activated alumina
14% nickel oxide Harshaw Ni-0707 140
supported on high
activity alumina
10% cuprous oxide Harshaw Cu-08 03 137
supported on high
activity alimiina
19% manganese dioxide Harshaw Mn-0201 69
supported on activated
alumina
20% ferric oxide Harshaw Fe-0301 41
supported on activated
alumina
3.5% silver supported Harshaw Ag-0101
on inert alumina
5% cobalt oxide, 5% Harshaw Co-0901 59
copper oxide supported
on activated alumina
Not reported General American Not reported
Transportation
"Gatalyst"
Mixed oxides of Mine Safety Appliance 156
manganese and copper "Hopcalite"
screen matrix. All of the materials tested showed
some activity for the removal of ammonia; however,
platinum, rutheniiim, and Hopcalite seemed most
promising.
4J
4J
c
d
p
(0
o
u
<u
0)
o
p
CO
u
H
+J
M
Id
04
<
1 I I
700 600 500 400 300 350
Temperature (** F)
Figure 1 - Catalytic Activity Comparison
large fraction of highly undesirable nitrogen dioxide
(NO2) product. "Gatalyst", copper oxide, and nickel
oxide all were considerably less active than the
better materials. The iron oxide, silver and
manganese dioxide were not effective showing only
a small amount of reaction at the highest temperatures
used.
While platinum, ruthenium, and Hopcalite all
performed well under dry conditions, only the platinum
was active under humid conditions (4). The activity
of ruthenium was reported by Gully et al^. (4) to be
"poor" in comparison to the activity of the platinum
catalyst. Hopcalite had also performed poorly under
high humidity operation (3). Based on these results,
platinum was chosen for this study.
The purpose of this study is actually threefold.
The first is to provide the experimental data for
the oxidation of ammonia over a platinum catalyst
in the ppm range which is so obviously missing from
the literature. The second is to use these data to
obtain a quantitative mathematical description or
model of the rate of ammonia oxidation with emphasis
on the utilization of this model for design purposes.
The third is to study catalyst stability and be able
to quantify the deactivation.
CHAPTER II
LITERATURE SURVEY
TABLE 3
The Effect of Hydrogen on Activation Energy
NH + O2 ^ HNO2 (2-5)
HNO2 ^ NO + OH (2-6)
NH + NH -> N2 + H2 (2-10)
16
NHO(a) + H20(g)
NH3(g)
02(g)
N2 1 0.5 -1.0
22
Physical Properties:
Temperature = 350F
Pressure = atmospheric
e = 0.3
d = 0.3175 cm
P
2
D = 0.44 cm /sec
G^ = 1.9 7 g mol/hr/cm^
y = 0.024 cp
p = 0.79 X lO"^ g/cm*^
2
a = 10.66 cm /g
^NH = 25 X 10 g mol/hr/g catalyst
Calculated Quantities:
^Re = 21
N^
Sc = 0.437
JD = 0.40
2
K = 1.53 3 g mol/hr/cm /atm
(P,
b - Ps ) = 1.53 X 10"^ atm
P - P
^2- s ^ 0.005
^b
31
Steps 2 and 6 involve intraparticle diffusion.
Physical examination of the catalyst particles reveals
a thin black layer of material on the outer surface
of the catalyst particle. Initial observation leads
one to suspect that this dark ring contains essentially
all of the platinum. Correspondence with the
manufacturer verified that the platinum was deposited
in the thin outer layer (15). Additional reinforcement
is provided from an effectiveness factor of 0.95
calculated for the 0.5 per cent platinum catalyst.
The effectiveness factor is defined as
ef f
where ^^ ^^ is the reaction rate in a porous catalyst
when there is no internal diffusion resistance and
r is the observed rate. The calculation and the
procedure are given in Appendix B.
The high value of the effectiveness factor is
expected since all catalyst was reported to be
essentially on the catalyst surface. In addition,
experimental runs were conducted to see if alumina
could catalyze the reaction. Within the temperature
range of this study, no oxidation of auranonia was
observed when alumina was used as the catalyst. Since
32
alumina cannot catalyze the reaction and all the
platinum is essentially on the catalyst surface,
internal diffusion is judged not to be the rate-
limiting step.
As is frequently done, a kinetic expression of
the Langmuir-type was chosen on a trial-and-error
basis as a possible description of the oxidation
reaction. This Langmuir form was found to fit the
data quite well. Verification of this fit is given
in Chapter VI. However, if the mechanistic steps
given by Reactions 3-5 through 3-10 are assumed the
Langmuir-type equation can be derived.
^6
Un + Pt ^==: (NH-.Pt) (3-5)
^6
^6 (1 + H ""nn/H^
Since the air used during this study contained
approximately 21 per cent oxygen, the difference
between the inlet and outlet oxygen concentration
caused by the ammonia oxidation can be considered
negligible. Consequently, the term C^ can be
"2
considered a constant. Combination of the constants
yields a rate expression of the form
^ ^ ^2 S H 3
- ^g|i = kg C^ Pt - k^ P f X . (3-20)
36
W -1
F in ^ + K4 (C^ - CQ) K3 + K^e"^2^ (3-29)
W (3-30)
F = Kj 1^ C^ -^ ^7 (^I - ^0^
0)
o
H
u
> }-4
<u
01 0)
[o}- (0 i
rd 5
s o
CO
:3
<d
a
/ ^ ^ a
0) (U
M >H o
0) (U r H 4->
u
rd (U
I^ <d o
-P (d
Clj u iS^-l O TS C 0)
H X
c CM c J-l -H
fd
(d en rd
to ^^ <u g
(0 E-i X u
O w Ky u <d w fd
CQ
(U
^ - 1 I
CN
(U
u
H
U4
<d Cu
6
u
o
4-)
rd
>1 rH
41
H
fd
4-)
0)
Q
u
o
4J
U
<d
O
fd
CQ
ro
H
44
B
0)
4->
(0
>1
w
o
-p
u
fd
<u
0^
(U
CQ
I
0)
X
I
in
0)
H
51
The reactor detail is shown in Figure 6. The
reactor was constructed of a 1-inch O.D. type 347
stainless-steel tube approximately 6 inches long.
The inside diameter of the tube was 0.76 5 inches.
All lines, fittings, and valves used in both
reactor systems were constructed of 316 stainless
steel except those specifically mentioned.
V"
% Stainless-Steel Tubing
h" Swagelok
h" Pipe to k" Tubing
1" Stainless Steel
Swagelok Cap
Catalyst
Alumina Pellets
roooi
Feed
F i g u r e 6 - Fixed-Bed R e a c t o r D(>tail
53
Ammonia Analysis
A relatively extensive effort was made in the
early stages of this investigation to develop a
reliable method of analysis for ammonia using a helium-
ionization gas chromatograph (Varian Model 1532-B).
54
TABLE 5
Gas Chromatograph Parameters
4->
0)
0 I
0)
3:
(U
e
w
3
-P
fd
fd
a,
<
CU o
o
H
4-1 0)
0^
3
H
-P
fd
fd x> u
Q) 3 H
(U rH
u u fd
3
u
-P
fd
T B
5-1
Q)
Pu
-p
>-l
3
(D
fd en
m H
-P
fd
4J
j
(0
c
o
u
K h M
0 (U
p ^ TJ
fd H d
rH < -H
3 rH
en
Q)
>i
U
iii__._
o 60
1 1 T" 1 1 i/ 4 1 I 1
/ O
1 O
fN
O
o
fd
c
O -H
00
S O
B -P
X
o U
0)
rH
CO
to
/ (U
/ o 00
^
3
/
H
/ o
CN
'
)
1 1 1 1 1 1 1 1
O o o o o
fN
o 00
e)^ue::^:^TlusueJJ, :^U93 jaj
61
1.500 g sulfanilamide
15.0 g tartaric acid.
100
90
80
70
60
i\
U 50
c N
fd ^
-p
-p 40
H
B
(0
fd N.
u 30
EH
-P
C
(D
U
V
20
\
s
10
0 10 20 30 40
PPM NO2
Figure 9 - Saltzman Calibration Curve
CHAPTER V
F. F
1 I 1 o
V ^
^Ai ^Ao
At steady state a material balance on component A
yields
F. C^. - F C^^ + Wr^ = 0 (5-1)
1 Al o Ao A
where F. = total molar flow rate in
1
F
= total molar flow rate out
o
C,. = inlet concentration of A
Al
C, = outlet concentration of A
Ao
r. = rate of production of A
W = weight of catalyst.
This study involved the reaction of ammonia in
the ppm range in air, consequently, the total molar
flow rate in can be assumed equal to the total molar
65
66
^4 S H - 1 SH^
67
I
O
o
B
fd P vo in CN t^ r^ [^ o CT> cy> in
fd VD ^ 00 r>- 1^ in fsj vo o r* ^
CN X> ^ r^ o CN VD
u
o ro CN fN CO CN ro CO vo in
B
a.
fa
o
O
in
en VD
P o
fd iH
0)
(d -p X
fd
D
s 4J
c fd fH r^ o rH o CN CO in 00 'ii' in cy in
0 in iH a\ CN CO in t^ CN r^ CN r-\ in r^
PQ o o VO
< p +J Cr rH 00 VD 00 'sr CN ro r- o vo i n (T 00 a\
EH o ^' CO rH CVJ CM CN cn CN iH r-{ CN CN
fd u
fd *-%.
CN CN iH
(U 0)
(hr
c
f^
X ffi rH
OUI
H z
;^
o
fd
CQ
CN fvj CN ro *^ ro CN vo <*
in v) vo 00 O 00 CN CM r* iH vo in o rg
rH CM r^ r-
u iH
iH ^ in CN a^ cj> r^ in ro
ro
a in r>j vo VO CN
i n c N ' : : j < i n i n i n o o c N f M i n i n o o o o
CN CN ro
in
'
o
\ in
rH
fa
0
o
in
ro
in -P
(N fd
r-i
B
ro CO
K H
2 3
fd
eCL jr:
O
o
o
ex (U
iH ^ s
\ 3
c
0 0
H H
-P 4J
fd O
>-i fd
\ -p Q)
in
c:
r-
<u IH
u
3 0
0
U -P
to
Ai <u
c t^
fd
H 1
1
o
Figure 10
\ in
\
\
in
\ CN
\
00 VO
o
g_OTXiom/(:;po 6) (J^) (uidd) 'o:^py/ D
70
TABLE 7
Kinetic Data :Eor Fixed -Bed Reactor at 300F
ppm NH3 Air Catalyst Weight
Date l ft3/hr* grams
ppm NH3
Date Time Hours Exit
8-6-72 5:15 P 3.25 31.3
Regeneration II
ppm NH3
Date Time Hours Exit
8-9-72 10:45 P 1 9.8
Regeneration III
ppm NH3
Date Time Hours Exit
8-17-72 5:00 p 1 64.4
ppm NH3
Date Time Hours Exit
8-23-72 1:00 P 1 20.2
ppm NH3
Date Time Hours Feed Exit
8-21-72 1:15 p 1 184.8 8.9
8-21-72 3:15 p 3 206.2 8.9
8-21-72 9:15 p 9 210.0 8.0
8-22-72 10:30 a 22.25 215.0 10.7
Long-Time Data
Application of the rate equation to integral or
fixed-bed data requires using the fixed-bed reactor
design equation that has been integrated from C to
C^. This equation as given in Chapter III for the case
where time is no longer a factor is
in Table 16.
The constant K. for the 400F data was determined
to be zero based on the least-squares analysis. If K.
is zero, the proposed rate expression reduces to the
81
82
form expected of a first-order reaction. This result
was not unexpected because the 400F data indicated
that the catalyst activity was independent of time
at that temperature which caused some suspicion that
a mechanism change might occur between 350F and 400F.
This point is not the objective of this work, because
it lies outside the 300 to 350F temperature range of
this study. The 400F data were taken only to provide
an additional point for an Arrhenius plot.
K = Ae"^/^^ (6-2)
K3 K4 X 10^
Temperature mol/(hr) (g cat) ppm~l
30
20
K4 X 102
10
9 V^
8 VX
7 V T ' -
6
5
3 4 -
1
fd
ro
Ui
(0
p \
3 >
fd
4J
(0
3
O
\^K3
u
<u
+3
fd
1
.9
.8
.7
.6 V
.5 v
.4 \
\
#.
.3
K3 = m o l / ( h r ) ( g ca :)
K4 = ppm" 1
20 21 22 23 24 25
1/T, O
OK
v-l X 10'
Figure 11 - T^rhenius Plot for Constants K3 and K4
86
tration are given in Figures 12, 13, 14, and 15 for
the 300, 320, 350, and 400**F data, respectively. The
points on these plots are the observed data points
while the solid lines are the predicted performance
based on the model represented by Equation 6-4. Visual
inspection of these plots indicates that the model
represents the experimental data quite wel4.. This is
further reinforced when the next three plots. Figures
16, 17, and 18, are viewed. These plots slvOW the
calculated reactor exit concentration plotced against
the experimental outlet concentration. If the model
was an exact representation of the experimental data,
all points would fall on a line with a slope of one
and an intercept of zero. The solid lines drawn on
these figures are such lines.
Short-Time Data
During the experimental work, it was observed
that the platinum catalyst could be regenerated by
passing pure ammonia over the deactivated catalyst
at 350-400F. At 350F and with a humidified air feed,
the catalyst activity was observed to decay quite
rapidly to a fairly constant value over a period of
about 75 hours. After this time the decay was
87
o
in
, Model
ro
u
o 3
\ o fd
erform
ro
o p^
in
ctor
(N
^ \
+J (d Pri
o (U (U
o iH o
CN 3 O
\ * H TJ O
(U ro
1
1 4J
ro TJ fd
o ffi (U
\ . in S X
rH H
B fa
\ a. 1
a
o rv)
o f-K
iH
igure
^ k\
O
in
fa
in o o
o
O o o
r- in in o in
(N CN CN
o
o fa
ro 0
o
CN
Model
ro
+J
o fd
in 0)
rformanc
rsi
o -p Q)
o (D PK
CM r-^
3 u
acto
M
1
ro <D
K Pi^
o 2
in TJ
rH B (U
a
Fixed-
a
O
o 1
ro
3
Cn
H
o fa
in
o o O o o
in o in o in
rvj CN
o
o fa
ro O
o
in
ro
Model
fd
Q)
CM U
3
fd
B
u
o
u
o -p
o
CN
3 U
H O
+J
I U
fd
ro Q)
o
in 13
B Q)
a PQ
a I
13
\ 0)
X
H
fa
o
o I
\
\ 0)
3
Cn
H
fa
in
o o o o o
in o in o in
CN
^ o
o fa
ro 0
from rable 10
Model
o
o
'^
4-)
fd
o
in <u
rformanc
rsj
\ <
<^ o -p 0)
o Q) Cn
rsj rH
\
3 u
acto
< M
1
ro Q)
ffi oi
o 2
I in 13
iH
B 0)
1 Ou CQ
Fixed-
O.
o
o 1
\
in
0)
5-1
3
Cn
H
o fa
in
o o o o o
in o in o in
CM CM
:^ax:vno - ^HN uidd
o 91
m
ro
o
o
ro
fa
0
o
o
ro
o 4J
in fd
(N
rH
fd
\ rH 4J
fd 3
4J 0)
3 e
(D H
o B )-l
^ o -H 0)
rsi U a
0) X
X w
a
\ H CO
\ >
^
4-> 13
o 3 Q)
in O -P
fH fd
ro iH
K 3
S U
fH
B fd
a u
O a 1
o
<-^
Figure 16
9
V
o
in
o o o o o
in o in o in
CNJ CNJ rH iH
o
o
ro
fa
0
o
i i
CM
ro
I o 4-)
^ ^ in fd
r>j
' i iH
I I fd
i rH 4->
fd 3
-p Q)
3 B
(U H
o B U
r \ \ o H (U
(N U a
Q) X
a.
X
w
w >to
^
4-) 13
\ | o 3 (U
in 4J
rH o fd
ro .H
E 3
2 u
B rH
fd
a u
o
a. 1
o
Figure 17
f-{
o
in
o o o o o
in o in o in
CN rI r1
o
o
ro
fa
0
o
in
ro
o 4->
in fd
CM
rH
fd
r-\ 4-)
! fd 3
4J 0)
3 B
I Q) H
o B U
o -H Q)
CN U
I
a
X
i
o fa
a
X
w to
K
>
4- 13
O 3 (U
-^f^; 1 in
iH
O -P
fd
ro iH
K 3
S U
rH
B (d
Cu u
O
a 1
% \ , .
o
Figure 18
iH
\^
O
in
\
\
v
o o o o o
in o in o in
CN CN
3
0
H B fa
4-) CU 0 o
195 ]
rrt O. r^j
35(
Humidified Operc
Temperature
NH3 Inlet
o
o
3
O
H
4J
fd
o CO }^
0)
00 u 3
3 Q)
Mode 1
O
Q)
PC
0)
B I
H
o CT\
0)
V4
\
3
v H
fa
c
N.
\
o
rsj
^ \ ^ ^
*^^^^^^#
o o o o o o
rsj o 00 CM
e
1
a 1
a
1
in
in
3
89.
0 iH
H fa
4J 3 o iH
fd 0 o (U
H in 13
u
Q) 4-> ro 0
Ch fd
O U Q) s o
4J ^ in
13 3 3
(U 0) 4-1
H u fd
4-1 3 S-i
H 0 Q)
13
H u g
a
M
g 4J (1)
in M
3 Q) EH
ffi ,_l CN
ation
3 fH
M
ro 1
l*
ffi
^. to u
U Q)
3 3
O 0 (D
O
a: CP
iH (U
^
<D
B 1
H
EH O
CM
0)
in M
Figu
r--
c
in
in
\ CN
o o o
o o in
CN in o
uidd ' u o T : ^ p j : ; u a o u o D t^ax'^^no HN
97
An attempt was made to determine the length of
time required to regenerate completely the catalyst.
A regeneration time of 15 minutes at 350'*F resulted
in the data given in Table 13. Based on the outlet
concentration at one hour for both Regenerations I
and III, the activity seems to be 50 per cent of the
totally regenerated catalyst. The comparison is
valid since operating conditions are essentially
identical for both runs.
The fourth regeneration was conducted at 400F
for a period of 14 hours. The reactor was then returned
to service but at 400F. After 28 hours, the conversion
was still at 93 per cent which was only slightly lower
than the 9 5 per cent conversion at one hour. Based
on these data, the catalyst activity decay was judged
negligible at 400F.
The last regeneration was to determine if the
elimination of water from the air feed stream would
affect the catalyst decay curve. The same analysis
was applied to the data obtained after Regeneration V
as was applied to the data from Regeneration I and II.
The constants K, and K2 obtained were 3.26 mol/(hr)
(g cat) and 0.016 hrs" , respectively.
98
In order to determine if there was a significant
difference between the rate of decay obtained for the
humidified runs and for the non-humidified runs,
constant K2 was statistically tested as described
by Volk (20). The results of the statistical test
show that there is a significant difference between
the slopes at the 99.9 per cent level. These results
indicate that the humidity in the air does affect the
rate of decay of the catalyst.
If the constants obtained for the 350F non-
humidified run are used in the model, the curve given
in Figure 21 is obtained.
Backmix Reactor
The constants K3 and K. obtained for the long-time
data were used to calculate the reaction rate from the
non-integrated rate equation. The curve obtained is
plotted in Figure 22. This shows how the constants
obtained from integral data agree with differential
data. The agreement seems to be fair, but the rates
predicted at high concentrations seem to be slightly
higher than those obtained in the backmix reactor.
This could be due to physical factors such as not
having a perfectly mixed reactor, leaks, or analytical
problems.
99
o
in
3
O o
Non-Humidified Operati
t
!
207 ppm
350F
o
Temperature
Inlet
o
o
NH
i 3
1
1 O
H
4J
fd
o U u
00
3 <D
O 3
(U
CP
0)
O^
B
i c
H I
CN
QJ
U
3
^ Cn
H
fa
c
\
\ ^
o
CN
o
o o o O
X. o
rsi
00
Model
to o
^ 11
in
13 3
0) H
CQ X 0
1 H Cl^
13 B
\ ^ M fd
\ X O 4J
\ ' - ' fd fd
CQ Q in
\ ^ CN
rH
son
B H
ih 5H
Cu fd
O ^ a.
B
o 3 o
-H 0
-H
u
4-) iH
td (U
u 13
4-)
3
o
dJ
\ U 1
\ % in 3
A V
r^ 0 CN
CN
u
ro (D
\ ffi U
Figu
S
\
o
in
\ <
> in
CN
Nv
o o o
ro
O
CN
in
CONCLUSIONS
j^ = 0.357/(Nj^-^^^e)
k = reaction rate constant
k' = reaction rate constant for reverse reaction
K = empirical constant
N^ = Schmidt number
Sc
N = Reynolds number
Re
P = total pressure (atmospheres)
Greek Symbols
y = viscosity (g/cm/sec)
3
p = density (g/cm )
APPENDIXES
106
k.
APPENDIX A
Platinum Catalyst Data
Calculational Procedure
V = 0.0251 cm^
P
S = 0.47 5 cm
X
Vg = 0 . 3 9 cm^/g
2
D, = 0.025 cm /sec
k
(kn)^Ko = 1-02 (l)/(hr)(g cat)
1 ODS
lb.fiV