Chemical Engineering CHFEN6553

Fall 1998

Dr. Terry A. Ring, MEB 3290

 

Computer Homework I - Due 7 October 1998

In the western USA and Canada, there are many natural gas wells that are sour which are treated by the Claus Reaction to remove H2S. The first step involves oxidation to a sufficient fraction of the H2S to SO2

2 H2S(g) + 3 O2(g)Û 2 SO2(g) + 2 H2O(g)

for a subsequent reaction that is performed catalytically. Other components of the natural gas also are oxidized to give CO, CO2 and water. After this oxidation, the stream is cooled and passed to a Claus reactor. A typical gas composition in mole % fed to a catalytic Claus reactor is :

N2 2%

CH4 74%

C2H6 8%

H2S 4%

H2O 4%

SO2 4%

CO2 2%

CO 2%

The Claus reaction is given below and is performed on a g-Al2O3 catalyst.

H2S(g) + 1/2 SO2(g)Û H2O(g) + 3/2 S(g and/or s)

There are two forms of sulphur which can result, gas or solid, depending upon conditions. If solid sulphur results, due to supersaturation, it will precipitate in the pores of the catalyst and decrease the number of catalytically active sites available, thus deactivating the catalyst. The catalyst is porous g-Al2O3 beads with a diameter of 0.2 cm which has a bulk density 1.4 gm/cc, a particle density of 2.1 gm/cc and a specific surface area of 100 m2/gm (the density of g-Al2O3 is 3.6 gm/cc). The forward rate on a fresh catalyst is pseudo-first order with respect to H2S with a rate constant given by :

k1 = 8.9 x 1016 exp[(-40,000 cal/mole)/RT)] .

a) For a 10 m tubular reactor what conversion of the feed is expected? The feed flow rate is 40Kg/hr at a temperature of 150oC and a pressure of 200KPa. Assume that the reactor is well insulated and consists of a packed catalyst bed with an inside diameter of 30 cm. Also assume that the catalyst is not deactivated by the solid sulphur for this calculation.

b) Plot the temperature, total pressure, partial pressure, conversion, equilibrium conversion and reaction rate profile down the length of the reactor at steady state assuming the catalyst is not deactivated by the solid sulphur.

c) Run the same calculations from a and b above accounting for deactiviation of surface sites as the sulfur is condensed in the pores of the catalyst. Pore condensation follows the Modified BET isotherm and is limited to n layers adsorbed. The adsorbed amount, G, is given by:

where Gm is the monolayer coverage, c = 0.1, n = 10 and x is the ratio of partial pressure of sulfur in the gas phase to its equilibrium partial pressure.

d) How long will it take for the reactor conversion to decrease to 10% of its initial value?

Thermodynamic Data at 298.15˚K

Taken from Appendix D, Chemical Thermodynamics by F.T. Wall, 3rd ed. Freeman Press, San Francisco, 1974, Handbook of Chemistry and Physics, 67th edition, 1986, and Perry's Chemical Engineering Handbook 5 th edition 1973.

Substance

?Hf°

(Kcal/mole)

∆Gf°

(kcal/mole)

(cal/deg/mole)

Cp°

(cal/deg/mole)

Al2O3(s)

-400.5

-3787.2

12.17

18.89

C2H6(g)

-20.24

-7.86

54.85

12.58

CH4(g)

-17.88

-12.13

44.492

8.439

H2S(g)

-4.93

-8.02

49.16

8.18

S(s-rhombic)

0

0

7.60

5.41

S2(g)

30.68

18.96

54.51

7.76

SO2(g)

-70.944

-71.748

59.3

9.53

H2O(g)

-57.796

-54.634

45.104

8.025

SO3(g)

-94.58

-88.69

61.34

12.11

H2(g)

0

0

31.208

6.889

N2(g)

0

0

45.77

6.961

CO(g)

-21.416

-32.780

47.219

6.959

CO2(g)

-94.051

-94.254

51.06

8.87

Molar Heat Capacities of Gases at Zero Pressure.

Data for solid sulphur

 

S(s-rhombic) MP =112.8˚C ∆Hfusion= 1.7175 KJ/mole

density = 2.07gm/cc BP = 444.6˚C ∆Hvaporization = 2,200 cal/mole

Cp˚(cal/deg/mole) = a +bT+cT2+dT3 (T˚K)

Taken from Hougen, O. A., Watson, K.M. and Ragatz, R.A., "Chemical Process Principles, Part 2- Thermodynamics" Wiley and Sons, 1947, Table D in Appendix.

Substance

a

b

c

d

N2(g)

6.903

-0.03753x10-2

0.1930x10-5

-0.6861x10-9

C2H6(g)

1.648

4.124x10-2

-1.530x10-5

1.740x10-9

CH4(g)

1.625

6.661x10-2

-3.737x10-5

8.307x10-9

H2(g)

6.952

-0.04576x10-2

0.09563x10-5

-0.2079x10-9

H2S(g)

7.070

0.3128x10-2

.01364x10-5

-0.7867x10-9

S(s-rhombic)

-

-

-

-

S2(g)

6.499

0.5298x10-2

-0.3888x10-5

0.9520x10-9