Data Book on Hydrocarbons

chemical engineering physical property handbook.Descripción completa

Views 177 Downloads 4 File size 60MB

Report DMCA / Copyright

DOWNLOAD FILE

Recommend stories

Citation preview

DATA BOOK ON

HYDROCARBONS APPLICA TION TO PROCESS ENGINEERING

by

J. B. MAXWELL NINTH PRINTING

ROBERT E. KRIEGER PUBLISHING COMPANY MALABAR, FLORIDA



ORIGINAL EDITION 1960 REPRINTED 1977 FROM NINTH PRINTING 1968

Printed and Published by ROBERT E. KRIEGER PUBLISHING COMPANY, INC. KRIEGER DRIVE MALABAR, FLORIDA 32950

© Copyright 1950 by

STANDARD OIL DEVELOPMENT COMPANY Reprinted by Arrangement with VAN NOSTRAND REINHOLD CoMPANY

All rights reserved. No reproduction in any form of this book, in while or in part (except for brief quotation in critical articles or reviews), may be made without written authorization frOm the publisher.

PRINTED IN THE UNITED STATES OF AMERICA

Library of Congress Cataloging in Publication Data

Maxwell, J B 1902Data book on hydrocarbons. Reprint of the 9th printing published in 1968 by Van Nostrand, Princeton, N. J., in The Esso series. Includes bibliographies. 1. Hydrocarbons. I. Title. TP690.M35 1975 661'.81 74-30163 ISBN 0-88275-257-X

PREFACE The primary purpose of this book is to provide (1) basic data on hydrocarbons and petroleum fractions, (2) methods of applying these data to process engineering, including illustrative examples and some fundamental theory, and (3) applications of a few of the unit operations of chemical engineering uscd extensively in the petroleum industry. Earlier editions of the present volume have been used in the Standard Oil Development Company and other affiliates of the Standard Oil Company (New .Jersey). Because this book has proved to be quite valuable to technical personnel, the Standard Oil Development Company has decided to make it available for practicing engineers and students of petroleum technology. The author is very much indebted to many associates in the preparation of thi s book and, in particular, to W. H. Hatch for invaluable assistance in editing the text and preparing the charts for publication, to C. O. Rbys, Sr., for the derivation of the .Mollier diagrams and other charts, to C. J. Robrecht (or constructive criticism and advice during the preparation of the manuscript. Furthermore, any list of acknowledgments would be incomplete without mentioning R. S. Piroomov who was responsible for the early development of a company data book. J. B. MAXWELL Standard Oil Development Company Linden, New Jersey



CONTENTS PHYSICAL DATA SElCTIOl\

PAGE

1. PHYSICAL CONSTANTS.......................................

1

Hydrocarbons, 2-Miscellaneous Organic Compounds, 6--MisceIlaneous Gases, 9 2. CHARACTERISTICS OF PETROLEUM FRACTIONS... . . . . . . . ..

10

Average Boiling Point, 14-Characterization Factor, Hi-Gravity, 18 3. MOLECULAR WEIGHT. . . . . . . . . . . . . . . . . .. . . . . .. ... ... . . . . . . . . .

19

Paraffins, 20-Petroleum Fractions, 21 4. VAPOH PRESSURE

.

24

Paraffins and Olefins, 27-Diolefins and Acetylenes, 35-Aromatics, 37 -Cycloparaffins, 39-Hydrocarbons, 40-Gasolines, 44 5. FUGACITY

.

45

Fugacity Function of Individual Hydrocarbons, 49-Fugacity Function of Hydrogen, Ol-Fugacity of Hydrocarbon Vapors, 62-Relative Volatility of LigM Hydrocarbons, 6~-Fugacity Correction Factor for Light Hydrocarbons in Absorber Oib, 67 6. CRITICAL PROPERTIES.......................................

68

Critical Temperature of Pure Hydrocarbons, 69-Critical Temperature of Light Hydrocarbon Mixtures, 'i'O-Critical Pressure of K ormal Paraffins, 71-Critical Temperature and Pressure of Petroleum Fractions, 72 7. THEHMAL PHOPERTlES

.

Specific Heats of Gases and Vapors, 88-Enlhalpy-Presoure Relationship for Hydrocarbon Vapors, 92-Bpecifjr Heats of Liquid Hydrocarbons and Petroleum Fractions, 93-Latenl Ileat of Vaporization of Light Hydrocarbons and Normal Paraffins, 94-Enthalpy of Individual Hydrocarbons, 98-Enthalpy of Petroleum Fractions, 114-Mollier Diagrams for Light Hydrocarbons, 128 VII

75

viii

CONTENTS PAGE

SECTION

8. DENSITy........................... . . . . . . . . . . . . . . . . . . . . . . . . . . .

136

Conversion Charts for 0 API Gravity, 138-Specific Gravity of Saturated Hydrocarbon Liquids, 14o-Thermal Expamlion of Liquid Petroleum Fractions, 143-P-V-T Relations of Hydrocarbon Vapors, 148 9. VISCOSITY. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

155

Conversion Charts, 158-Viscosity of Hydrocarbons and Crude Fractions, 161-Viscosity-Temperature Charts, 166--Viscosity Index of Lubricating Oils, 168-Viscosity Blending Index, 173-Viscosity of Hydrocarbon Vapors and Miscellaneous Gases, 174 10. COMBUSTION. . .

178

Heat of Combustion of Petroleum Fractions and Hydrocarbon Gases, 18o-Enthalpy of Flue Gas Components, 182-Heat Available from the Combustion of Refinery Gases and Fuel Oils, 184-Properties of Flue Gases, 189 UNIT OPERATIONS 11. FLOW OF FLUIDS.............................................

193

Friction Factor for Fluid Flow, 19B-Pressure Drop in Commercial Pipes, 199-Equivalent Length of Fittings, 202-Friction Loss Due to Contraction and Enlargement, 204-Discharge Characteristics of Weirs, 205-Pressure Drop Across Tube Banks, 206 12. FLOW OF HEAT. . . . . . . . . . . . . . . . . . . . . . . . . . .. . .

.. .

207

Heat Loss by Radiation and Natural Convection, 209-Heat Transfer to Fluids Inside Tubes, 211-Heat Transfcr to Fluids Outside Tubes, 212-Thermal Conductivity of Petroleum Fractions, Water, and Gases, 213-Logarithmic Mean Temperatme Difference, 217 13. EQUILIBRIUM FLASH VAPORIZATION.. . . . . . . . . . . . . . . . . . . . . ..

222

14. FRACTIONATING TOWERS.. ..

230

Minimum Reflux Ratio and Theoretical Steps, 23O-Correlation of Theoretical Steps with Reflux Ratio, 244-0verall Plate Efficiency, 245-Packed Towcrs, 246 CONVERSIOX FACTORS. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..

249

INDEX............................................................

253

Section I

PHYSICAL CONSTANTS In the following tables the more common physical constants are given for hydrocarbons, certain other organic series, and miscellaneous gases. While these constants, in general, are based upon reliable data, estimated "alues were included in a few instances where available data were considered questionable. Where no reasonably good basis was available for either estimating or calculating the constants, they are omitted. The density, boiling point, melting point, and heat of eo'mbustion for most of the hydrocarbons are those given in the Burea1t of Standards Circular C461. 1 GENERAL REFERENCES Annual Tables of Physical Constants, Nat. Research Council (19-11, 1942). Beattie, Poffenberger and Hadlock, J. Chem. Phys. 3, 96 (1935). Beattie, Simard and Su, J. Alii. Chem. Soc. 61, 24 (1939); 61,924 (1939). Cole and Cole, J. Chem. Phys. 9, 341 (1941). Doss, "Physical ~nstants of the Principal Hydrocarbons," 4th Edition, The Texas Co., New York, N.Y. (1943). Ginnings, J. Am. Cltell/. Soc. 62, 1923 (1940). Ginnings and Baum, J. Am. Chem. Soc. 59, 1111 (1937). Ingersoll, Thesis, ~Iass. Inst. Tech. (1930). International Critical Tables, Vols. I and III. Kay, Ind. Eng. Chem. 30, 459 (1938). Kharasch, J. Research Nat. Bur. Standards 2,359 (1929). Krase and Goodman, Ind. Eng. Chelll. 22, 13 (1930). Meyers, Scott, Brickwede and RAnds, Unpublished Data, Nat. Bur. Standards. Pickering, Bur. Standards Sci. Paper 511 (1926). Rintelen, Gross and Saylor, J. Am. Chelll. Soc. 62, 1923 (19-10). Tables anntlelles de wnstantes et dunnee nUllteriqlte, Vols. VII to XIII (1925-1939). I'

"Sclcdcd Values of Propertips of Hydrocarbons," Nal. Bur.

(947).

1

~lalldards

Circular Cl,61

~

PHYSICAL CONSTANTS OF HYDROCARBONS

. FOR~tULA

MOLEC. WT.

BOILING POINT of

:.fELTING POINT of

NORMAL PARAFFINS Methane ................... Ethane .................... Propane ................... Bu~no ....................

CH, C,H. C,H, C,H,o

Pentane ................... Hexane .................... Heptane ................... Octane ...... '" ...........

C,H 12 C,H" C,H 16 C,H,s

72.1 86.2 100.2 114.2

96.9 -201. 5 155.7 -139.5 209.2 -131. 1 258.2 - 70.3

Nonane .................... Decane .................... Undecane .................. Dodecane .... , .............

C,H,o C1oH" C"H,. C12 H "

128.2 142.3 156.3 170.3

303.4 345.2 384 .4 421.3

ISO-PARAFFINS Isobutane ..................

C,H 1o

58.1

10.9 -255.0

C,H 12

72.1

82.2 -255.5

C,H"

72.1

49.0

2-Methylpentanc (Isohexane) , 3-Methylpcnt,ane, ........ 2,2-Dimethylbutane ();eohexane) . ' .. , ......... , ... 2,3-Dimethylbutane (Diisopropyl) ...............

C,H" C,H"

86.2 86.2

140.5 -245 145.9 -180

C,H"

86.2

C.H"

86.2

2-Met,hylhexane (Isoheptane) . 3-Methylhexanc ............ 3-Ethylpentane .... ...... 2,2-Dimethylpentane ....... '

C,H I6 C,H" (',H" (',H"

2,3·Dimethylpentane ........ 2,4-Dimethylpentane ........ 3,3-Dimethylpentane ....... ,

C,H" C,H" C,H"

2-Methylbutane (Isopentane), 2,2-Dimeth:dpropane (~eopentane). . ..............

16.0 -258.9 -296.5 30.1 -128.0 -297.8 44.1 - 43.8 -305.7 58.1 31.1 -216.9

DENSITY

°API

Sp Gr 60°/60° Lb/gal

t

P Atm

G/ml

206.3 306

45.8 48.2 42.0 37.4 32.6 29.4 26.8 24.6

of

D

HEAT OF COMBUSTION

@

60°F-BTU /Ib

Gross

Net

0.162 .203 .226 .225

23,860" 22,300" 21,650" 21,290'

21,500" 20,42019,930" 19,670"

.232 .234 .234 .233

21,070" 20,780 20,670 20,590

19,500" 19,240 19,160 19,100

20,530 20,480 20,450 20,420

19,050 19,020 19,000 18,980

0.30 .374 .508 .584

2.50 3.11 4.23 4.86

92.7 81.6 74.2 6ti,{j

.631 .664 .688 .707

5.25 5.53 5.73 5.89

386.5 455.0 512.5 565

64.5 61.3 58.7 56.4

. -.).) ,- ... .734 .744 , 753

6.01 6.11 6.19 6.27

612' 654' 695' 731"

.563

4.69

275

36

.234

21,240'

19,610"

.625

5.20

369.5

32.4

.234

21,030"

19,450"

.597

4.97

329"

35'

-

20,960'

19,330"

o;.:-

83.5 80.0

.658 , 669

5.48 5.57

437' 443'

31' 30"

-

-

20,750 20,760

19,210 19,220

tp

121.5 -147 6

84.9

.654

5.44

415'

31'

-

20,700

19,160

136.4 -198.8

81.0

.666

5.54

441

31

20,740

19,200

100.2 100.2 100.2 100.2

194.1 197.5 200.2 174.6

-180.8 -182.9 -181.5 -190.8

75.i 73.0 69.ti 77.'2

. 68:l .692 .703 .678

5.68 5.76 5.85 5.64

496 504 508' 475'

28" 28.5" 28.5 28.5

20,650 20,660 20.670 20,600

19,140 19,150 19,160 19,090

100.2 100.2 100.2

193.6 176.9 -183.1 186.9 -211.0

-') , I ._

70.6

.700 .678 .698

5.83 5.54 5.81

498' 477 487'

29 28.5' 28'

20,540 20,620 20,620

19,130 19,110 19,110

+

+

64.5 21.5 14.1 14.7

340 247 147 111

CRITICAL CONSTANTS

-116.3

+ 90.1

23" 22"

20" lb'

-

-

d

;.:>-3

;.:tp

o o

~

o Z

:I:

+

2.1

120 94.9 105

71.2

..

. 241

-

-

-

-

-

>


Z

~ lfJ

-

-

20,150 Q

-

-

20,320Q -

19,210Q

-

-

20,060"

18,950Q

...

;> t"'i

19,180Q

19,040" 19,04QQ

* Mixture of cis- and ** Sublimes.

-


~ c; o o ~

o

Z ::t: ~

-

-114

9.8 -153

+

47.7 -188.5 80.4 - 26.0

104.4 -159 132.8 -148

209

60~

Pentyne-1 (Propylacetylene). Pentyne-2 .................. 3-Methylbutyne-1 (Isopropylacetylene) ................

CsH s CsH s

68.1 68.1

CsHs

68.1

79.7

.670

5.58

410"

Hexyne-1 (Butylacetylene) ... Hexyne-2 .................. Hexyne-3 ...............•..

CsH lo CSH lO CsH lO

82.1 82.1 82.1

160.9 -205.6 184.1 -126.4 179.2 -149.8

65.0 60.8 63.1

.720 .736 .727

5.99 6.13 6.05

-

-

4-Methylpentyne-1 .......... 4-Methylpentyne-2 .......... 3,3-Diroethylbutyne-1 .......

CSH lO CSH lO CaRlO

82.1 82.1 82.1

142.1 -157.1 162 100.0 -114.2

67.5 65.3 78.7

.711 .719 .673

5.92 5.98 5.60

-

-

C.H.

52.1

73.9

.689

5.73

365"

OLEFINS-ACETYLENES Buten-3-yne-1 (VinylAcetylene) .................•..

D G/rol

HEAT OF COMBUSTION @ 60°F-BTU lib

82

42

-

-

20,550

20,74~

19,46()Q

-

-

-

-

-

75.;

-

-

-

-

-

-

tj

~

o o > ~ c;

C

Z

en

-

58.7

0.744

6.19

49.4

.782

6.51

-

56.4 32.8

.753 .861

6,27 7.17

-

-

176.2

41.9

28.6

.884

7.36

231.1

-139.0

30.8

.872

7.26

106.2 106.2 106.2 106.2

292.0 - 13.3 282,4 - 54.2 281.0 55.9 277 .1 -138.9

28.4 31.3 31.9 30.8

.885 .869 .866 .872

C,Hu

120.2

349.0 -

13.8

25.7

C,Hu

120.2

336.5 -

47.3

C,H"

120.2

C9H l2 C9Hn C,H" C,H" C,R"

120.2 120.2 120.2 120.2 120.2

318.6 306.3 329.2 322.7 324.5

CYCLOPARAFFINS Cyclopropane ... ............

C,H.

42.1

-

Cyclobutane . ...............

C.H.

56.1

+ 54.7

Cyclopentane . ...•.......... Methylcyc1opentane . ........ 1, I-Dimethylcyclopentane ....

CeHlo CtHn C,H u

l,2··)imethylcyc1opentane-cis. 1,2-Dimethylcyclopentanetrans . ................... 1,3-Dimethylcyclopentanetrans . ................... Ethylcyclopentane ..........

Cyclohexane ................ Methylcyclohexane . .........

Penten-l-yne-3 . ............ (AllylacetyPenten-l-yne-4 lene) . ................... 2-Methylbuten·l-yne·3 . .....

C.He

66.1

138.6

C,H e CoHo

66.1 66.1

107 90

Hexen-l-yne-3 . ............. Hexen-l-yne-5 .. ............ 2-Methylpenten-l-yne-3 ..... 3--Methylpenten-3-yne-l* .....

C.H. CoHo CoHo CoHo

80.1 80.1 80.1 80.1

185 158 169 156

AROMATICS Benzene ... .................

CoHo

78,1

Toluene ..... ...............

C,H.

92.1

o-Xylene .... ............... m-Xylene .........•....... . ;p-Xylene . .............. , ... Ethylbenzene .... ...........

C.H IO CaHlo CaHlo CaH lo

1,2,3-Trimethylbenzene . ..... 1,2,4-Trimethylbenzene (Pseu· documene) . .............. 1,3,5·Trimethyl benzene (Mesitylene) . ................ Propylbenzene . ............. lsopropylbenzene (Cumene) .. I-Methyl-2-Ethylbenzene .... I-Methyl-3-Ethylbenzene .... I-Methyl-4-Ethylbenzene ....

-

-

-

-

-

-

-

-

-

-

-

-

---

-

--

-

551.3

47.9

0.304

17,990

17.270

609.1

41.6

.292

18,270

17,450

7,37 7.24 7.21 7.26

675 655· 652 655

37 36' 35' 38

.288· .288· .270'

-

18.500 18,500 18,430 18,490

17,610 17,610 17,540 17.600

.900

7.49

72C/'

32'

-

-

-

29.1

.881

7.34

708'

33

-

18,570

17.620

328.3 - 48.6

31.1

.870

7.24

700'

-

18,620

17,670

-147,1 -140.8 -126,6

.866 .866 .883 .870 .868

7,21 7.21 7.35 7.24 7.23

690 68C/' 702< 695' 696'

-

18,660 18.670

17,710 17,720

o

34' 34'

82.7

31.9 31.9 28.7 31.1 31.5

33 /' 34' 35'

27.0 -196.6

98.6

.615

5.12

256

54

-

-

58.0

74.8

.686

5.71

385'

50'

-

cB

70.1 84.2 98.2

120,7 -136.7 161.3 -224.4 189.5 -105

56,9 56.2 54.7

.751 .754 .760

6.25 6.28 6.33

470' 520' 550'

46'

C7H u

98.2

210.7 -

62

50.4

.778

6.48

570'

40'

C,H u

98.2

197.4 -182

65.4

.757

6.30

560'

41'

-

C 7H 14 C 7H u

98.2 98.2

195.4 -213 218.2 -217

57.2 52.0

.750 .771

6.24 6.42

555' 580'

41' 40'

-

20,110

18.760

CeRn C 7H 14

84.2 98.2

177.3 44 213.6 -195.6

49.0 61.3

.784 .774

6,53 6.44

538 575

40.4 40'

.273

20.030 20.000

18.680 18.650

+

, Heat of combustion &8 a gas-otherwise &8 a liquid. • Estimated.

-

-

-

+

-

-

-

-

Critical temperature-boiling point correlatioil. ., Vapor pressure curve or correlation.

C

34'

42' 42'

-

-

-

-

-

-

--

-

-

-

-

-

20,350' 20,110

18.760

20,020

18,670

20.020

18,670

-

-

"'d

~ ~

......

o

~

19.‫סס‬OO

-

-

• Mixture of cis- and trans' isomers. •• Sublimes.

o:n

PHYSICAL CONSTANTS OF ORGANIC COMPOUNDS 0:.

FORMULA

MOLEC. WT.

BOILING POINT of

MELTING POINT of

CRITICAL CONSTANTS

DENSITY

Sp Gr

600 j600 Lb jgal

t

of

HEAT OF COMBUSTION HEAT OF VAPORIZ.

@B.P. P D Atm Gjml BTU lib

@ 60°F-BTU /lb Gross

Net

-- -ALCOHOLS Methanol (Methyl Alcohol) .. CHaOH Ethanol (Ethyl Alcohol) ..... CH aCH20H Propanol-1 (Normal Propyl Alcohol) ................. CH aCH2CH20H Propanol-2 (Isopropyl Alcohoi) ........ , ............ (CHa),CHOH Butanol-1 (Normal Butyl Alcohol) ................. Butanol-2 (Sec. Butyl Alcohol) 2-Methylpropanol-l (Isobutyl Alcohol) ................. 2-Methylpropanol-2 (Tert. Butyl Alcohol) ........... Pentanol-1 (Normal Amyl Alcohol) ................. Pentanol-2 (Sec. Amyl Alcohoi) ..................... Pentanol-3 (Diethyl Carbinol) 2-Methylbutanol-l (Sec. Butyl Carbinol) ................ 2-Methylbutanol-2 (Tert. Amyl Alcohol) ........... 3-Methylbutanol-1 (Isoamyl Alcohol) ........ , ........ 3-Methylbutanol-2 (Methyl Isopropyl Carbinol) ....... 2.2-Dimethylpropanol-l (Tert. Butyl Carbinol) ..........

474

9760

8580

.794

6.61 469.6 63.1

.275

361

12,780

11,550

207.0 -195

.808

6.73 506.7 50.0 .273

296

14,450

13,190

180.2 -129

.789

6.57

289

14.350

13.090

> ~

148.1 -143.7 0.796

46.1

173.0 -174

60.1 60.1

tj

254 242

15.500

14,220

-

-

-

-

-

249

15.450

14.170

o

(6.60)

-

-

-

235

15,290

14,010

.819

6.82

-

-

223*

16,220

14,930

.814 .826

6.78 6.88

-

-

-

213* 211*

-

...

.820 6.83 - .825 -6.87

-

-

-

218*

-

-

77.9

(.793)

280.4 -109.8

(CH ahCHCH 2OH

74.1

226.4 -162

(CHa)aCOH

74.1

180.7

CH a(CH 2) aCH 20H

88.1

-

-

CH a(CH2hCH(OH)CH a (CH.CH 2hCHOH

88.1 88.1

247.1 240

-

CH aCH2CH(CH a)CHIOH

88.1

264

-

to

-

6.71

-129.6

-

-

.806

243.9 211.1

-

o o

6.78 549 6.75

74.1 74.1 ,

-

.814 .811

CH a(CH 2hCH2OH CH aCH 2CH(OH)CH a

48

-

~

Z

:I:

~ ~

o

§o Z

00

-

-

203*

16,030

14,740

-

216

16.150

14,860

-

-

-

-

209*

-

-

210*

-

-

9.31

-

-

-

344

8250

7340

CH aCH 2C(OH) (CHah

88.1

215.8

15

.815

6.79

-

(CHahCHCH2CH t OH

88.1

269.2 -179

.814

6.78

-

(CHahCHCH(OmCH.

88.1

233

.825

6.87

(CH a)aCCH20H

88.1

236

-

62.1

387.5

GLYCOLS AND GLYCEROL Ethanediol-l,2 (Ethylene Glycol) ..................... CH2(OH)CH 2OH

-

6.63 464.0 78.7 0.272

32.0

120-125

9

1.118

..

Propanediol-1,2 (Propylene Glycol) ................. CH3CH(OH)CH 2OH Propanediol-1,3 (Trimethy- CH 2 (OH)CH2CH 2(OH) lene Glycol)

76.1 76.1

371 850 (appr .)

Propanetriol-1,2,3 (Glycerol). CH2(OH)CH(OH)CH 2OH

92.1

554

ETHERS Methyl Ether .............. CH,OCH 3 Ethyl Ether ............... CH,CH 2OCH 2CH,

46.1

-

94.1

-

-

273* 266*

-

-

-

260

52

5.99 381

35

8.68

65.0 1.265

10.53

-11.5 -217

74.1

-

1.042 -

-

-

-177 .3 0.719

-

-

-

10,350 10,450

9350 9450

-

7760

6940

0.271

187

13,570u

12,340u

.262

151

15,840

14,560

-

129 120

16,930 16,870

15,630 15,570

17,560

16,250

-

-

Propyl Ether .............. CH 3(CH 2)20(CH 2)2CH 3 Isopropyl Ether ........ '.' .. (CH')2CHOCH(CH,)2

102.2 102.2

194.2 -188 155.3 -122

.752 .729

6.26 6.07

-

Butyl Ether ............... CH 3(CH 2)30(CH2)3CH 3 Sec. Butyl Ether ........... [CH3CH 2CH (CH,) 120

130.2 130.2

288.0 -144 250 -

.773 .760

6.44 6.33

-

-

-

115* 109*

-

-

-

-

-

320*

.786

6.54

-

-

-

257*

11,400

10,540

ALDEHYDES Methanal (Formaldehyde) ... HCHO

30.0 -

3

-180

-

-

8050U

-

7420U

~

P::

~ 00

......

(1

> ~

Ethanal (Acetaldehyde) ..... CHaCHO

44.0

Propanal lPropionaldehyde)

58.1

120

-114

.812

6.76

-

-

-

215*'

13,400

12,420

72.1

167.2 -144

.809

6.74

-

-

-

189*

14,640

13,590

00

72.1

142

.799

6.65

-

-

-

180*

14,600

13,550

58.1

133.0 -138.8

.795

6.62

-

-

-

220

13,260

12,280

> Z ;1

72.1

175.5 -123.5

.810

6.74

-

-

-

190

14,540

13,490

86.1 86.1

216.1 -108.0 215.2 - 40

.812 .820

6.76 6.83

-

15,430 15,380

14,330 14,280

200.7 -134

.820

6.83

-

168* 168*

86.1

-

165*

15,350

14,250

100.2

240.6 -119

.806

6.71

-

-

152*

15,980

14,840

CH,CH2CHO

Butanal (Butyraldehyde) .... CH 3CH 2CH2CHO 2-Methylpropanal (Isobutyraldehyde) ............... (CH 3)2CHCHO KETONES Propanone (Acetone) ....... CH 3COCH a Butanone (Methyl Ethyl Ketone) ................. CH 3COCH 2CH, Pentanone-2 (Methyl Propyl Ketone) ................. CH 3COCH 2CH 2CH, Pentanone-3 (Diethyl Ketone) (CH 3CH 2hCO 3-Methylbutanone-2 (Methyl Isopropyl Ketone) ........ CH 3COCH (CH a) 2 4-Methyl Pentanone-2 (Methyl Isobutyl Ketone) .

CH,COCH~H(CH,h

68.5 -190.3

-

-

-

87

-

-

* Calculated or estimated with a probable accuracy of ±2%.

(1

o

~

~

"'-J

u Heat of combustion as a gas-otherwise as a liquid.

..

PHYSICAL CONSTANTS OF GASES

FOn~ULA

MOLEC. WT.

BOILING POINT

MELTING POINT

of

of

CRITICAL CONSTANTS

HEAT OF COMBUSTION

@

60°F-BTU j1b

,~~

P Atm

D G /ml

Gross

270.3

111.5

0.235

9670

8000

88.0

73.0

.460

-

-

-

-220.4

34.5

.301

4345

4345

291

76

.57

-

-

7.51

369

51.6

.33

-

-

-400

12.8

.031

61,100

51,600

124.5

81.6

.42

-

-

°API

Sp Gr Lb/gal 60°/60°

97.5

0.617

5.15

'42.0

.815

6.78

t of

NH 3

17.0

-

Carbon Dioxirle ......

COz

44.0

-109.3*

-

Carbon Monoxirle ....

CO

28.0

-312.7

-::137.0

-

-

Chlorine .............

Ch

70.9

-

30

-151

-

-

Ethyl Chloride .......

CzH~CI

64.5

54.1

-214

Hydrogen ...........

Hz

2.0

-423.0

-434.5

-

-

-

Hydrogen Chloride ...

HCl

36.5

-121.0

-173.6

-

-

-

Hydrogen Sulfide .....

H 2S

34.1

-

76.5

-122.0

46.0

.797

6.64

212.7

88.9

-

7100

6550

Methyl Chloride .....

CH 3Cl

50.5

-

11.6

-143.8

20.3

.931

7.76

289.6

65.8

.37

-

-

Nitrogen ............

N2

28.0

-320.5

-346.0

-

-

-

-232.8

33.5

.31

-

-

Oxygen .............

Oz

32.0

-297.4

-362.0

-

-

-

-181. 9

49.7

.43

-

-

Sulfur Dioxide .......

S02

64.1

14.0

-

-

315.0

77.7

.52

-

-

Ammonia ............

28.1

-107.9

DENSITY

69.9

98.9

25.5

.901

1.394 11.62

"'tl ~

~ ~

o > ~ o o

z

U1

"-3

>

Z

"-3

U1

*.Sublimes.

~

..

Section 2

CHARACTERISTICS OF PETROLEUM FRACTIONS Average Boiling Point of Petroleum Fractions

,.

Many physical properties of pure hydrocarbons can be correlated with specific gravity and normal boiling point as independent variables. However, for use in the petroleum industry, these correlations must also be applicable to petroleum fractions which are mixtures of a large number of components, usually having a wide variation in boiling points. While the average specific gravity is a property of the petroleum fraction which can be measured directly, just as in the case of pure compounds, there is not an analogous average normal boiling point for a mixture. By integrating or averaging its distiUation curve (temperature vs. liquid volume percent distilled), a volume average boiling point can be determined for the mixture. However, as Watson and Nelson! and Smith and Watson 2 have pointed out, this has no special significance as a true average boiling point and many physical properties can be better correlated by the use of some other average boiling point, i.e., weight average, molal average, etc. Consequently, in all correlations involving boiling points of petroleum fractions, the proper average should be used. For the following physical properties, these are: Average Boiling Point

1 2

Physical Property

Volume average

Viscosity Liquid specific heat

Weight average

True critical temperature

Molal average

Pseudo-critical temperature Thermal expansion of liquids

Mean average

Molecular weight Characterization factor Specific gravity Pseudo-critical pressure Heat of combustion

Watson and Nelson, Ind. Eng. Chem. 26, 880 (1933). Smith and Watson, Ind. Eng. Chem. 29,1408 (1937). 10

CHARACTERISTICS OF PE'l'ROLEUM FRACTIONS

11

Since a distillation curve is usually available and a volume average boiling point is readily obtained therefrom, the other average boiling points are given as a function of these data. The chart on page 14 is based on an assay (True Boiling Point) distillation 3 of the whole crude, while the chart on page 15 refers to the 1070 (or ASTM) distillation of the fraction itself. The chart on page 14 was derived empirically from crude assay fractions of a number of crudes. For narrow boiling fractions, all of the average boiling points approach each other and the volume average boiling point may be used for any of the others. Then, by appropriately combining the volume average boiling points of the narrow cuts, the various average boiling points of wider cuts were determined. The weight and molal average boiling points of the wider cuts were calculated directly by combining the narrow cuts on the basis of their weight and mole fractions, respectively. The mean average boiling point could not be calculated in the same manner since it is not a direct average or integral of its fractional parts. As used herein, mean average boiling point is defined as the boiling point which best correlates the molecular weight of petroleum fractions. Consequently, the mean average boiling point for wider cuts was determined indirectly from the generalized molecular weight chart on page 21. Although Smith and Watson proposed a cubic average boiling point for the correlation of characterization factor, specific gravity-boiling point relations forthe different crudes indicate that the present mean average boiling point can be used for correlating gravity, and consequently characterization factor. Smith and Watson also used cubic average boiling point for correlating viscosity, but the present data indicate that the volume average is the proper boiling point. Since these boiling point correlations were developed directly from crude assay distillations, this chart should always be used 4 if an assay is available. Otherwise, the 10% (or ASTM) distillation of the fraction may be used in conjunction with the other chart. The latter was derived from the crude assay chart and an empirical correlation between the two types of distillation curves. The difference between the two sets of curves at zero slope represents the thermometer stem corrections for the 10% distillations. In the case of light hydrocarbon mixtures, where the analysis is known, the volume, weight, and molal average boiling points can be calculated directly from the boiling points of the components and their volume, weight, and mole fractions, respectively. On the oth~r hand, the mean average boiling point must be determined indirectly from the average molecular weight of the mixture. Up to an Approximately 15 theoretical plates and 5 to 1 reflux ratio. Below slopes of 2°F/% for low boiling fractions (V.A.B.P. < 500°F) and 3°F/% for high boiling fractions (V.A.B.P. > 500°F), the volume average may be used for the other average boiling points with very little error. 3 4

DATA BOOK ON HYDROCARBONS

12

average molecular weight of 80, the molecular weight-boiling point relation for normal paraffins (page 20) may be used for this purpose, but for higher molecular weights the generalized chart on page 21 should be employed. Characterization Factor Watson and Nelson 1 introduced characterization factor as an index of the chemical character of pure hydrocarbons and petroleum fractions. The characterization factor of a hydrocarbon is defined as the cube root of its absolute boiling point in oR divided by its specific gravity (60°F/60°F), or Characterization Factor

=

yt T B/Sp Gr

Characterization factor is given on page 16 as a function of gravity in °API and boiling point in of for hydrocarbons and petroleum fractions. That characterization factor is only an approximate index of the chemical nature of hydrocarbons is indicated by its variation with boiling point both for members of a homologous series and for fractions from the same crude (page 17). However, it has considerable value in that it can be applied to the entire boiling range of a crude and it has been generally accepted by the petroleum industry. Typical Crude Fractions For approximate use when there are insufficient data, several correlations have been developed for typical crude fractions grouped according to characterization factor and viscosity index. 5 These groups are numbered in order of decreasing paraffinicity and each may be considered representative of the crude fractions within its characterization factor or viscosity index range. The five groups were arbitrarily selected as follows:

Group I II III

IV

v

Characterization Factor . 12.1-12.6 . 11.9-12.2 . 11.7-12.0 . 11.5-11.8 . 11.3-11. 6

Viscosity Index of Lube Fraclions6 80-100 60-80 40-60 20-40 0-20

Fractions from some of the more common crudes are cla5sifil'd in the following table: lS

6

See page 156. Dewaxed to +20°F pour.

CHARACTERISTICS OF PETROLEUM FRACTIONS

13

TYP.lCAL GROUP

White Products Pennsylvania I Rodessa. . . . . . . . . . . . . . . . . . . . . . . . . . . .. I Panhandle . . . . . . . . . . . . . . . . . . . . . . . . . . . II Mid-Continent . . . . . . . . . . . . . . . . . . . . . .. II Kuwait I-II CRUDE

Gas Oils and Heavier I I I II II-III

Iraq Iranian East Texas South Louisiana. . . . . . . . . . . . . . . . . . . . .. Jusepin

II II III III III

II-III II-III II II III

West Texas . . . . . . . . . . . . . . . . . . . . . . . . .. Tia Juana (Med. and 102) Colombian Lagunillas . . . . . . . . . . . . . . . . . . . . . . . . . ..

III III IV V

III IV IV V

Since, in the case of some crudes, the lower boiling fractions belonged in a different group than the higher boiling fractions, they were classified separatelythat is, into white prorlucts having an average boiling point less than 500°F, and gas oils and heavier having an average boiling point greater than 500°F.

+40

+ 30

2

3

4

5

6

WEIGHT ·AVERAGE·

+20 +'0

0 - 10

7

-

10



Ii

iJ

9

~

AVERAGE BOILING POINT OF PETROLEUM FRACTIONS ~ CRUDE ASSAy DISTILLATION

- 20

i

MEAN AVERAGE

-30 -40 -50 -60 -70

2

4

0

5

6

7

8

THE SLOPE AND THE 50% POINT FOR THE VOL. AV. B.P. UNLESS THE DISTILLATION FOR THE FRACTION DEVIATES APPRECIABLY FROM A STRAIGHT LINE. !N THE LATTER EVENT THE FOLLOWING FORMUI AS SHOULD BE USED:

-40 -

_ t7O-t10 S 60 ty = to+4t.50+tIOO 6

-60 In

..

cr....

('~

10

* THE CUT RANGE MAY BE USED FOR

MOLAL AVERAGE

-20

V>

9

-80

FOR WHOLE CRUDES:

t y = ho

-'00 -120 -'40

3

4

5

7

6

14

8

9

10

t

t,50+ teo

t40 WEIGHT AVERAGE +20 ,.L

0

-20

,~ 4

2

5

t . t·

6

7

8 i:J...~t-L_; of

MEAN AVERAGE

i....... .ly.......H-t~.

AVERAGE BOILING POINT

+20

,.

OF PETROLEUM FRACTIONS 0

10 % (A.5.tM.) DISTILLATION

-20 IF AVAILABLE, THE CRUDE ASSAY DISTILLATION SHOULD BE USED FOR DETERMINING AVERAGE BOILING POINTS.

-40

r r-

l/}


2). Vaughan and Graves, Ind. Eng. Chem. 32, 12.;>2 (1940). Wiebe and Brcevoort, J. Am. Chem. Soc. 52, 622 (1930). Wiebe, Hubbard and Breevoort, J. Am. Chem. Soc. 52, 611 (1930) .



_~:

.-- -

4-

::l

HH-

'!' --

80 70 60

50

40 -_.~

30

20

-;-

-300

o

-200

27

100

200

-2'50

-200

-150

.8

-

.1 .6 .5

l--:r.

~

--

-100

~

.. - , ~ ~.-

VAPOR PRESSURE OF ETHANE AND ETHYLENE

.4 .3

.2

200

.08

80

.07 .06

60

.()5

50

70

.04

40

.-

.O~

30

.02

20

all', , II

.1 ~

_.

...-

gr

4 3

2

-200

-tOO

o

100

28

200

300

400

-100 - ,-;;.,,",

i-

-50

::T:-E.=:io€-~:ff-'

~ ::ri..:r--,~f.':'=j-:',

le"",

-,

..6

,,,="

=

-'

.1

'~7,k

,c.:.~3':!..J

,

= VAPOR PRESSURE OF PROPANE AND PROPYLENE

.5 .4

+- ,

'+:

.3 --H

.2

.1 .09 .08 .07

200

100

~

r=

90

_

'''''

" !CA_~{~!

.06

-t-t=f

60

-,.!~

.05

70

''f'

50

.04

40

.03

30

.02

20

II -'I=.:f.: "

.009 ,008

-'

.001

10 9 8 7

IXXJ

6

.005

5

=r

4

T

3

2

o

100

200

300

29

400

500

50

·5

-

.4

VAPOR PRESSURE OF BUTANES AND 8UTENES

:II+-_~,

..3

.2

200

30

I .

II II II

.01 ~-====47

-

1

I I

/I

I

o

V V 1/

I

tOO

200

300

30

400

500

-roo

t.O .9

o

50

p.

: ::r=tfel. ·-:-r-·: ~-.g:

--

VAPOR PRESSURE OF .5 ~

PENTANE

~

-I'

...

::T'

AND ISOPENTANE

.3

.2

'200

11

I

.€ ==':.f"_ -=I::

:t=:!-~:.

90

_ !:::f

.07 .06

1=l;.T:,:

80 70 60

.05

50

.()4

40 30 20

.02

.OOQ

Q'

I

I

.01

.- -

, .

-',--.

-:.

,008 .007 .006

.:

-

-

:~

±

.§'

- ":'::. - .

--

8

7 6

,=:

l:±-

H

.005

5

.004

4

.003

_.

- I::±l

2

.002

It

.

.00IU-J....l....L.LLL.JLJ...,U.I-......AJ...J...Iu1t.U./-l..'~~I~J.:..z.J,;.W-l...l-J...J.-l...J..J..1..J..J...J...L..J...J..J..J..J...J...IW-L...J.-l....J.-l...J..J..u..J-U...J...J..J..J..1..J..J...J...L-U 100

200

300

400

500

31

l--.

.-.---------

600

-50

-100

0

100

- 1 I f t i j - -.:.~

.9

.8

..,

';-..::

T.~:

.

':

="=±

.5

-

50

: .....

'.':~:

ISO

VAPOR PRESSURE OF HEXANE

.4

.a

200 I~

.1 .09

-t...

.

-fl

:.~,"'f> :. 'c: :::'.-:x':~ ::~-,.

-17

.07

.

n-

"

_.

"c't._',_ - _"':4'_ .•••.-.t

'--J-;,'

:: t~

:'':'

--. :r.' ,-- '-,-, , 'f--X:==, T.:i::-'~ ~,

.-

..

-'-+-

.

~.

=- __.---7l~:. ==1':::.

,.. -,-V-FFl r-- ~-;:=::;=:I -I:' ' =l:

'

, -

100 90

80 70

.06

60

.05

50

.04

40

.03

30

.02

20 .'

r

I

/

;=;=. ~ 11.,,':.. '.,~':'1.

£

.008 .OOf .006

.'f=:£:i=E j~ --!'

.''=8, !E::-l,=.:l3:: ' --:=;::1:'. :t+= 'r:

- -,

1-..

-- -

,~.

,~f:~T ~.~:~I­ _:. ,_ ,. ~ -

,-J.'

~t:.:~

'

-, -. . ,.. =.-1-=!iiC:l 10 9

= -

9 7 6

5

.00

4

3

.003

2

I

I I I

I

II .• ;.,

32

-50 .8

.5

0

. =-=

::::.C

50

100

150

·=--==~~-:'~';'_~c~"=§=:l~F-=~. .

_--


...

.4.0 0 F/'Y.

0 COORDINATING RESEARCH

COUNCIL (CRC) HANDBOOK. PP. 244-254 (1946)

44

Q.

~ L&J

Q.

0

0..

0:

It

It

0:

10

REFERENCE:

400

0

TRUE VAPOR PRESSURES CORRESPONDING TO A R.V.P. OF 9.0 LBS./ so. IN. AND A SLOPE OF 4.0°F/'>'. ARE READ FROM THE CHART AS FOLLOWS: TEMPERATURE TRUE VAPOR PRESS. OF LBs/sa. IN. 32 2.8 100 9.9 21.4 150

10

It (/)

l&J

SLOPEc 160-120

20

300 w ;:) (/)

EXAMPLE: DETERMINE THE TRUE VAPOR PRESSURES AT 32°F. 100°F AND 15Q°F OF THE FOLLOWING GASOLINE: REID VAPOR PRESS. - 9.0 LBS./SO. IN. DISTILLATION + LOSS, 5'>'. ~ 120 OF 15'>'. (! 160 OF

30

~

:::>

12 14 16 18 20

40

~

(J) (J)

II

50

0

~ en 5.0 CD

4

;;: 4. •0 5 2 6 UJ 7 Q:" 8 :> en 9 en 10

L&J 0:

80

~ d

...J

l&.

y

*

SLOPE

120

90

50 60

160

100

40

II

12 13 14 15 16 17 18 19 20

30

500 ~ us ;:) 600 It I-

700 800 900 1000

Section 5

FUGACITY Raoult's Law If two or more compounds form an ideal solution in the liquid phase, and if the saturated vapors of the individual components are perfect gases, the system has been termed an ideal system.! For such a system the partial vapor pressure of any component may be calculated from the composition of the liquid phase by Raoult's Law and from the composition of the vapor phase by Dalton's Law. An equation of these two expressions gives the liquid-vapor equilibrium relation for any component, i where i = 1, 2, ... , n:

or

Pi = PiXi = 7rYi

(1 )

== P i /7r = Ki

(2)

yi/Xi

where Pi

=

partial pressure of i

Pi = saturated vapor pressure of i Xi = mole fraction of i in the liquid phase

Yi

= mole fraction

7r = Ki

of i in the vapor phase

total (vapor) pressure of the system

= vapor-liquid equilibrium constant for i at the temperature and pressure of the system

The above equation, usually referred to as the Raoult's Law relation, is true only for ideal systems, ,as defined above. However, it is usually a good approximation for mixtures of homologues and, in general, for mixtures of chemically similar compounds, if none of the saturated 'vapors at the equilibrium temperature deviate too greatly from a perfect gas. Up to moderate pressures (several atmospheres) hydrocarbon mixtures frequently fall within the scope of the Raoult's Law relation. However, its application to these mixtures is rather limited because of the wide differences usually encountered between the boiling points of the most volatile and least volatile components. This results in equilibrium temperatures at which the saturated vapors of the lowest boiling components deviate considerably from a perfect gas, even though the equilibrium pressure of the system may be relatively low. I

Gamson and Watson, Nat. Petroleum News} Technical Section 36, R-258 (1944).

45

DATA BOOK ON HYDROCARBONS

46

Fugacity Functions In order to improve the accuracy in predicting vapor-liquid cquilibrium constants for hydrocarbons at higher pressures, Lewis and Luke 2 and other investigators replaced the pressures in equations (1) and (2) by -analogous fugacities for any component, i, whereby:

or

Ii = fpiXi = fnYi

(3)

Ki

(4)

Yi/Xi

where fi

=

=

fpdf.Tri

=

fugaci·ty of i in either phase of the system

fpi = fugacity of i as a pure saturated liquid (or vapor) at its vapor pressure

corresponding to the equilibrium temperature of the system f-rri

= fugacity of i as a pure vapor at the equilibrium temperature and pressure of the system

Generalized correlations have been developed for the r-atio of fugacity to pressure for pure hydrocarbons as a function of reduced temperature and reduced pressure. A correlation of this type (pages 62 and 63) was used in conjunction with the vapor pressure charts to develop the fugacity function charts for individual hydrocarbons. 3 The fugacity function given by these charts, 7rfp/!-rr, may be considered a corrected vapor pressure and used in place of the latter in ~ny equation pertaining to liquid-vapor equiiibrium such as equations (1) and (2). These simple fugacity relations greatly extend the pressure range for which liquid-vapor equilibria for hydrocarbon systems may be predicted with confidence, and can be used up to equilibrium pressures of 20 to 25 atm with a fair degree of accuracy. Beyond these pressures and especially as the critical point of the mixture is approached, serious deviations from true equilibrium conditions are encountered. Under these circumstances, the assumptions of ideal mixtures no longer hold and the fugacities of the individual compounds are dependent upon the compositions of the liquid and vapor phases as well as temperature and pressure. In the region where the simple fugacity relations no longer apply and consequently beyond the scope of the present charts, there are data in the literature on a number of specific binary and multicomponent hydrocarbon systems. Also, The M. W. Kellogg 00. 4 has published an excellent correlation for lig~1t paraffin and olefin hydrocarbons in which the fugacities of the individual compounds are given as a function of the molal average boiling points of the liquid and vapor 2 Lewis

and Luke, Trans. Am. Soc. M echo Engrs. 54, 55 (1932).

a This method was actually used only up to the critical temperature of each compound. Beyond this point values were calculated from more general fugacity correlations developed by The M. W. Kellogg Co. to avoid using extrapolated vapor pressure curves. 4"Liquid-Vapor Equilibria in Mixtures of Light Hydrocarbons," The M. W. Kellogg Co., New York, N. Y. (1950).

FUGA CITY

47

phases in additi on to the equilib rium tempe rature and pressu re. The Kellogg correl ation was derive d from the applic ation of exact therm odyna mic relatio ns to a comprehensive equati on of state for pure hydro carbon vapors and liquids and their mixtur es. 5 If, in additi on to hydro carbon v-apors, other gases (air, H , CO , etc.) are 2 presen t in the vapor phase, it is recom mende d that an effective pressu 2 re, equal to the produ ct of the total pressu re multip lied by the square root of the mole fractio n of the entire hydro carbon portio n of the vapor, or 7r-VV;;;, be used in determ ining the fugacities of the indivi dual hydroc arbons . Fragm entary data have indioa ted that this effective pressu re gives better results than either the total pressure, 7r, or partia l hydro carbon pressure, 7r'YHC, for determ ining indivi dual fugacities. Then, after the fugacities or fugaci ty functio ns have been read from the charts , the total pressure is again used as a basis for all equilibl:ium calcul ations . The following examp le illustr ates the applic ation of the fugaci ty functi on charts when other gases are presen t in the vapor phase: Examp le 1. Determ ine the pressu re and composition of the liquid phase in equilib rium with a vapor of the following composition at gO°F: 1st Trial Component

Vapor. Mole Fract.

1l'

0.040 .220 .280 .175 .160 .125 1.000

CH 4 C2H6 C3Hs C 4H 1O

= 25 atm

1I"e = 21.5 atm F, atm

Air H2

2nd Trial

* * 180 38.0 13.5 4.9

x

-

0.039 .115 .296 .637 1.087

Interpo lation

11" = 20 I1tm 1I"e = 17.2 atm

F, atm

11"

x

* *

-

180 36.0 12.7 4.4

0.031 .097 .252 .[;68 0.948

= 21.8 atm x

0.034 .104 .269 .593 1.000

• In this example, the fugacity functions of air and H2 are conside red to be infinite.

where 7r = total equilib rium pressu re 7re = 7rVO.740 = effective pressu re used to determ ine fugaci ty functio ns F = 7rfp/!1I" = fugaci ty functi on for pure hydroc arbons

x ~

= 7rY/ F

Relati ve Volati lity Since relativ e volati lity is quite useful in fractio nation problems, curves for the relativ e volatil ities of light unsatu rates and isoparaffins to the corresponding norma l paraffins are given on pages 64 to 66. The curves for the C 4 unsatu rates IS

Benedict, Webb and Rubin, J. Chem. Phys. 8, 334 (1940); 10, 7474 (1942).

48

DATA BOOK ON HYDROCARBONS

inay also be used in conjunction with the normal butane fugacity chart to predict fugacity functions for these compounds. Except for butadiene and the normal butenes, these relative volatility curves were derived from the Kellogg fugacity correlation. Composition was indirectly taken into account to some extent since the fugacities for each pair of compounds were read at the same liquid and vapor molal average boiling points as well as at the same temperatures and pressures. In general, the relative volatility charts may be considered to have a somewhat greater range of applicability than the simple fugacity charts. They may be used up to 25 atm, irrespective of the composition of the liquid and vapor phases of the mixture; beyond this pressure their application is limited to systems in which there is a difference of at least 75°F between the molal average boiling points of the two phases, but under no circumstancel? should the curves be extrapolated. While all of the curves may be considered to be accurate within 25% for the relative volatility minus one (0; - 1), deviations from the solid curves rarely exceed 15% for this difference. Chemical Structure and Liquid Activity Coefficients When components in a hydrocarbon mixture are quite dissimilar chemically, the liquid phase may deviate appreciably from an ideal solution. This effect of chemical structure is not taken into account in any of the fugacity correlations heretofore considered. It has been mentioned that in correlations of the Kellogg type, fugacity is a function of the liquid and vapor compositions, but only with respect to components of similar chemical structure. To correct for chemical dissimilarity in solutions of light hydrocarbons in absorber oils, liquid activity coefficients are given for these light hydrocarbons on page 67. Within the range of the data these activity coefficients were practically independent of temperature (100°F and 220°F) and pressure (500 psia and 1000 psia) .

GENERAL REFERENCES Brown, Souders and Smith, Ind. Eng. Chem. 24, 513 (1932). Dean and Tooke, Ind. Eng. Chem. 38, 389 (1946). Hadden, Chem. Eng. Progress 44, 37 (1948). Kay, Chem. Revs. 29, 501 (1941). Lewis, Ind. Eng. Chem. 28, 257 (1936). Lewis and Kay, Oil and Gas J. 32, 40 (1934). . Lewis and Randall, "Thermodynamics," pp. 190-198, McGraw-Hill Book Co. (1923). Nelson and Bonnell, Ind. Eng. Chem. 36, 204 (1943). Sage and Lacey, Ind. Eng. Chem. 30, 1296 (1938).

~

..

·~-t-toj-or

~T:~

FUGACITY FUNCTION

1m

OF METHANE

_. t

-l- .

200

r v.

'I i . ':.;=

~:-:.:F-~:-~-

-;':~i:£

-.

J:-

10

-::-±

_

I

.• -~

100

90

I::f.;=+

80 70

60

50 40 30

20

-

I

: :.-;;!

~:.,~

=~~~ -'- :.=;-

j. p-L~-t.: :~- .

_

~/J::.1~·_-,2? ~E§= -



.j...

-

=:=.

10 9

8 7

6

5 ,

4

..::~

3

-300

-200

o

-100

49

100

200

300

&I .

..

F~~~CITY FUN~TI~N OF ETHYLENE

I

r +-

-

..

I-- ~

1 00

90

80 70

60 50 40 30

20

-

••

_I~ 8 7 6

5 4

-200

o

-100 50

100

200

300

~

-~=

~

-

--- --

-

-

---

=

= = = = = -

FUGACITY FUNCTION OF ETHANE --

-

~

00



--

--

-

-

-

--

-

- - --

-

--

-

-

-

- ~t a -

--

I 00

-- -

-

90

---

-

80

-

70 60

--

50

40 30

20

~

-

-

~

/ -

--

- -rl=Fl

10

-'

- -

9

8 7 6

5 4

3

2

-zoo

o

-100

51

100

zoo

300

1.0 0.9 0.8 0.7

0.6

--200

o

-100

.=a=::-

-

/M

0.5

.

-

-c

FUGACITY FUNCTION OF PROPYLENE

0.4

-

-

0.3

0.2 I~

/1

AI If I

II

0.1

I 00

90 80 70 60

...

50

40 30

20 I-

-, I I

IJ

iii

II

10 9 8

7 6 5 4 3

2

-100

o

100

52

200

300

400

-200

1.0 0.9

0.8 0.7 0.6

-

o

-100

, 0:;

FUGACITY FUNCTION. OF PROPANE ::

i

0.5

- . -.

0.4

-

-+-+-l

-

0.3

0.2

..

0.1

~

I 00

i

90 80 70 60 50

40

-

30

20

11 -(

If

V I

iW

II

'.



, •

I

8

7 6

4

~

l.

-

9

5

-+-

-l-

10

3

.

,

,

2

I

"'T

I

J

I

i,""'" I

-ir-

,

I

I

!..ll -j(jO

__.x. ~

L

If

.-L~

c.:.'_ _,-

o

100

53

I I

200

300

400

__

:T"-~

,- FUGACITY FUNCTIO:lI "OF ISOBUTANE :: ,

~

-

r

t

0.2

"1J;:,f-;t--l+++H Ifl+I-I-/HI-H+-I--+ ,

I

'I

I

ElIOO

~l--

~

=3...J .

~l='-

90

80 70 60

,-

50 40 ___ t:::t:::t:=t=J:

30

20

I.-

t;o

~

r+- -H-lH+I--t+H+t++-H+-H+-l-H-++++H+H-I-+-1-++I-I4-IA+:I4-jll--W-+-I-l--!-l-l--+-W-!-W~-I--l 1

I ~:-:.!-

3:,-'=i-- -i-:-(=F" -...t.~ ~-- .~~.l

'7

--

""--

It

=-

-

-~Ig

..;

--t.-

8 7 6 5 4

3 -I- -

1/111 1/

I

-tOO

o

tOO

54

200

300

400

o

100

-!!III

·i-tll- f:.Tfl ~aii--~f~ "':::-l_

~---.-.

--r=..t=

:.

-~

8+

-,

I-

-f-

FUGACITY FUNCTION [~I--

r:

§i=-I

OF BUTANE

- ffP-i' -:-.:~. '- Ii § ~~ H-".IEI:- ~=8::=

_,.-,_-r-F

--=+:~-

-

--'-+-

'--4-=1= 0.2

1111

I

--.

100 90 80 70 60

j:-t-;'

±:'--; -I--l.-r·~·

7

.. -f~- 3f~

,.J~ 51~ =t--+=!'I !_-l ~~:

- ..: ->-.-

t-

50

,~

f---

...::;C :-I.,

~

40

i7

30

.-

20

1/ j

=. -.

-,..

---:::-:~-_.

~r_ ~_r=-- ~ ~ -~_

1=I j:::f

...

, .-1'\)

10 9 8 7 6

-f-

5 ..

.

. ..•

~E:::

r-

4

-

:;=. :7t;- ~ 3

I I

I

H--i--H-+-t-l---H-++-+-+-+-IH-+++-t--HH-++-+-+--HH- 'I - - LIl-lH--H..j-+-HH--H-++-HH--H-++-H-+-H-++-H-+-~--H-j I

I

-100

o

100

55

200

300

400

-100

a

100

0.9 0.6 0.7

~

0.6

0.5

~.

FUGACITY FUNCTION OF ISOPENTANE ~ :

0.4

0.3

0.2

t I I

40

30

20

I I

4 3

I

IjjjT±-;1-±±ttttjjjttt::t:t:tttttlltJzt~~:tt:t1tttt:!i!~~~~:t;;'tr~t'ttt:r~;r~11!!!t:!=1

o

100

200

56

300

400

500

1 H-++++ -H-++-- H-t-f-+- t--H 't+t-i-:+t+'

1

., -"1f-f-+I H--H-+ t-H-t-f- +-H-+t +t-H-t- f-+-H-+ t+t-Hf- f-t-I-t-t --H-H- Hf-f-t --:

50

=

40

... .

17'



1

I~

:¥T 1

17'

20

If

11

It

30

, .+-H-H-H-H-i+H-H-H-H-H-H-H-H-H-+++++

. I

-

+++++-hI"'H-T+-H'-H' '.s the latent heat of vaporization of benzene at 1 atm. Example 2. Determine the latent heat of vaporization of the following ga9 oil at 500°F. 10% Distillatioll 10% @ 430°F 50% @ 540°F 70% @ 605°F 90% @ 680°F

Gravity 35°API

Vol. Av. B.P. = 547°F; Slope = 2.9°F/% Mean Av. B.P. = 547 - 9 = 538°F Molec. wt. = 211 Vapor pressure (538°F normal B.P.) = 0.63 atm. at 500°F Pseudo-critical pressure = 266 psia co 18.1 atm I'l'lolee. wt. of normal paraffin (538°F normal B.P.) = 222 Critical pressure of normal paraffin = 15.0 atm Vapor pressure of normal paraffin = (l5.0/18.1)0.63 = 0.52 !!.tm Latent heat of normal paraffin = 104 BTU/lb Latent heat of vaporization of the gll3 oil!!.t 500°F =

22Z X 104 = 108 B1'Ulib 214

78

DATA BOOK ON HYDROCARBONS

Enthalpy of Light Hydrocarbons The enthal py 8 or heat content of low-boiling paraffins, olefins, and aromatics is given by the chart.s on pages 98 to 113. These charts can be applied to mixtures of light hydrocarbons on the basis of the following assumptions:

1. The entha!pies of individual components of a mixture are additive in the !iquid phase, that is, the mola! heat content of the mixtttre equals the sum of the products of the mo!al heat contents of the components by their mo!e fractions. 2. The entha!pies of individua! components are additive in the vapor phase at !ow pressures (0-1 atm). 3. The change in enthalpy of the vapor with pressure at constant temperature is the same for a mixture as for a sing!e compound having the same mo!ecular weight as the mixture. The first assumption is substantially true for hydrocarbon mixtures (especially for homologous series) at temperatures below the critical regions of all components. At temperatures near to or above the critical temperatures of any of the components, the liquid mixture is no longer an ideal solution of its components and there is some deviation from the rule of additive heat contents. However, since these deviations arc not too serious, and since no other simple method has been developed for determining the heat content of a liquid mixture, the rule of additive enthalpies should be used for all hydrocarbon mixtures irrespective of the critical temperatures and chemical composition of the components. The second assumption is strictly true only for vapor mixtures at infinite dilution (0 atm) but is a very close approximation for pressures up to 1 atm. The third assumption is empirical but has been shown indirectly to give quite accurate rcsults for mixtures of homologous series and petroleum fractions. Also, the usc of the average molecular weight to determine the change of enthalpy with pressure is the simplest average which can be used. Above thc critical temperature a dashed line is shown for the heat content of the gas in solution. This line was based on the assumption that thc gas in solution at any temperature would have the same partial density and enthalpy as the pure compound at a pressure corresponding to an extrapolation of its vapor pressure curve above the critical point. Obviously, this is only a rough approximation since both a vapor prcssure curve and an ideal liquid solution arc meaningless in this regIOn.

E:rample 3. Determine the difference in enthalpy bctween the liquid at 100°F and the vapor at iiOO°F and 20 atm for a mixture having the following composition: 8

Based on on entholpy of zero for the saturated liquid

01,

-200"F.

THERMAL PROPERTIES Component

Mole Fraclioll

C,H. CsH. C,H lO C,H, CsH.

0.100 .500 .100 .050 .250 1.000

79

The enthalpy of the mixture as a liquid at 100°F and as a vapor at 500°F and 0-1 atm is computed from the individual components as tabulated below:

Com-

Mole

poncnt

Fract.

CzH s

0.100 .500 .\00 .050 .250

CaH s C 4 H 1O

C2H 4

C3 H,

Enthalpy of Liquid lOO°F

l\lolcr, Wt.

Enthalpy of Vapor 5OQ°F and 0-1 elm

Lb/Mole of Mixture

BTU 1~lole

BTU lIb

3.0 22.0 5.8 1.4 10.5 42.7

of Mixture

239 171 159 223 169

li y (500°F, 0-1 atm) - li L

BTU lIb

720 3760 920 310 1770 7480

BTU ll\lolc

of J\'lixture

553 530 525 506 508

1660 11660 3040 710 5330 22400

= 22,400 - 7480 = 14,920 BTU/mole

The change of enthalpy of the vapor at 500°F betwecn 0-1 atm. and 20 atm. is computed by interpolating between C,H. and CsH.: C,H.: 1f,,(500°F, 20 atm) - li y (500°F, 0-1 atm) = 30(546 - 553) = -210 BTU/mole CsH.: H y (5OO°F, 20 atm) - li y (500°F, 0-1 atm)

=

44(522 - 530)

=

-350 BTU/mole

Mixture: li y (500°F, 20 atm) - li y (500°F, 0-1 atm) = -210 =

- 30 + 42.7 H _ 30 {-350 -

(-21O)J

-340 BTU/mole

Therefore,

li y (500°F, 20 atm) - liL(100°F) = -340

or

14,580 42.7

= 849

+ 14,920 =

14,580 BTU/mole

BTU/Ib

The foregoing procedure can hc simplificd, with a loss of accuracy which does not usually exceed 5%, by interpolating on a basis of molecular wcight and total

so

DATA BOOK ON HYDROCARBONS

olefin content between the initial and final states:

CaH s: H v (500°F, 20 atm) - Ih(IOO°F) = 44(522 - 171) = 15,440 BTU/mole

C ZH 4 : H v (500°F, 20 atm) - Ih(lOO°F) = 28(500 - 223) = 7750 BTU/mole

CaHe:

H v (500°F, 20 atm) - HL(lOO°F) = 42(500 - 169) = 13,900 BTU/mole

Since the average molecular weight of the paraffin portion of the mixture is 44, the propane values can be used directly, making interpolation unnecessary. The average molecular weight of the olefin portion is 39.7; hence the enthalpy difference between the initial and final states will be: 7750

+ 3:;7 _- 2~8 (13,900 -

7750) = 12,880 BTU/mole

Interpolating between the paraffin and olefin portions, Hv(500°F, 20 atm) - HL(lOO°F)

or

X 15,440 + 0.30 X 12,880 = 14,670 BTU/mole

= 0.70

14,670 42.7 = 344 BTU/lb vs, 342 BTU /lb by the longer method.

Enthalpy of Petroleum Fractions The enthalpyO of petroleum fractions is given by the charts on pages 114 to 127 for both paraffinic stocks, having a characterization factor of 12.0, and nonparaffinic stocks, having a characterization factor of 11.0 over a mean average boiling point range from 200°F to 800°F. Theoretically, these charts represent pure hydrocarbons of the designated characterization factor and boiling point, but they may be applied to petroleum fractions if the following assumption is made in addition to the three previous ones pertaining to light hydrocarbon mixtures: 4. 'The avemge difference between the enthalpy of the vapor at low preSSU1'es (0-1 atm) and the enthalpy of the liquid, at constant temperature, is the same f01' a rni.'l:ture of chemically similar hyd1'Ocarbons as for a single compound of the sct'1ne molecular weight (or mean avemge boiling point). \\'hile this assumption is empirical, it is accurate \\'ithin a few percent excrpt in the region of the pseudo-critical temprrature \\'here the enthalpy of the liquid is subject to variation depending upon the true criticaltemperatUl'e of the mixtUl'e. Since the dashed line starting at the pseudo-critical point applies only to a pure compound in solution above its critical point, another dashed line was arbitrarily drawn for mixtUl'es, joining the satUl'aled liquid line below the pseudo-critical 9 Based on all enthalpy of zero for the saturated liquid at 00 F .



THERMAL PROPERTIES I

81

point with the pure compound line about 50°F above the pseudo-critical tempera_ ture. This is more representative of a mixture and should be used in preference to the pure compound line. These charts may be interpolated and extrapolated linearly with both characterization factor and mean average boiling point. Occasionally, in inter·

polating between two adjacent boiling point chart~ the pressure and temperature of the vapor will be such that they fall inside of the "dome" of the higher boiling point chart.. Since it is impossible to use the charts in this region, it is recommended that the two adjacent lower boiling point char'" be extrapolated upward. Following are two examples illustrating the use of these charts: Exan,ple 4. Determine the ditTcrence in enthalpy between the liquid at 500°F and the vapor at 775°F and 25 psig for the following refined oil fraetion: Grauity 40 0 API

Crude Assay DistillatiOlt I.RP. 300°F 50% 440°F F.B.I'. 580°F Vol. Av. RP. = 440°F

· '11 . 81 ope 0 f t}Ie dIstl atlOn eurve

=

580 - 300 100

=

2.8 0 F/%

Mean Av. B.P. = 440 - 6 = 434°F Characterization Factor = 11.65 hI' = Enthalpy of the vapor at 775°F and 2.7 atm (25 psig) h L = Enthalpy of the liquid at 500°F Mean Ch. Cb. Cb.

Au. B.P. -400°F Factor = 12: h" - hL = 567 - 286 = 281 BTU/lb Factor ~ II: h v :- h L = 538 - 263 ~ 275 BTU/lb Faetor = 11.65: hI' - h L = 275 + 0.65(281 - 275) = 279 BTU/lb

Mean Ch. Ch. Ch.

Au. B.P. -500°F Faetor = 12: hI' - h L = 556 - 273 = 283 BTU/lb Factor = II: hI' - h L ~ 534 - 255 ~ 279 BTU/lb Faetor = 11.65: hI' - hi = 279 + 0.65(283 - 279) ~ 282 BTU/lb

Mean Au. B.P. - 434°F Cb. Factor = 11.65: hI' - hL = 279

+ 1.'.. (282

- 279) =.280 BTU/lb

If the char'" for 300°F and 400°F Mean Av. B.P.'s had been extrapolated, the result would have been essentially the same, 281 BTU/lb.

Example 5. Determine the difference in enthalpy between the liquid at 425°F and the vapor at 925'F and 350 psig for the following gas oil:

82

DATA BOOK ON HYDROCARBONS 10% 10% 50% 70% 90%

VoI. Av. B.P. = 455

Distillation @ 455°F @ 560°F @ 620°F @ 695°F

Gravity

15.5°API

+ 2 X 560 + 695 = 5670F 4

. SI ope 0 f di Bt I'II ation curve

- 455 = 620 69 = 2.8°F/% 5 = 562°F

Mean Av. B.P. = 567 Characterization Factor = 10.48

h v = Enthalpy of the vapor at 925°F and 24.8 atm (350 psig) h L = Enthalpy of the liquid at 425°F Av. B.P.-400°F Factor = 12: hv - h L = 662 Factor = 11: h v - h L = 622 Factor = 10.48: h v - h L = 406 Mean Av. B.P.-500°F Ch. Factor = 12: hv - hL = 642 Ch. Factor = 11: h v - lt L = 606 Ch. Factor = 10.48: hv - h L = 398 Mean Av. B.P.-562°F Ch. Factor = 10.48: ltv - h L = 388 Mean Ch. Ch. Ch.

233 = 429 BTU/lb 216 = 406 BTU/lb 0.52(429 - 406) = 394 BTU/lb 224 = 418 BTU/lb 208 = 398 BTU/Ib 0.52(418 - 398) = 388 BTU/lb M(394 - 388) = 385 BTU/Ih

MollieI' Diagrams The MollieI' diagrams for the individual light hydrocarbons are of essentially the same type as the familiar one for steam. To minimize confusion and to make the charts as easily usable as possible, lines of constant volume are omitted and lines of constant temperature replace lines of constant superheat in the superheated vapor region. These charts will be used principally for adiabatic compressions and expansions. In applying the MollieI' diagrams to hydrocarbon mixtures, the mixture should be treated as a single compound of the average molecular weight. An empirical study of the diagrams indicates that successive charts of the same series (paraffin or olefin) may be interpolated (or extrapolated) by assuming a linear relation exists between melecular weight and (1) isentropic change of molal enthalpy with pressure and (2) the product of the square root of the molecular weight and the isentropic ehange of temperature with pressure. If both paraffins and olefins are present in the mixture, the charts of each



,

THERMAL PROPERTIES

83

I

series are interpolated (or extrapolated) to the average molecular weight of t.he total mixture. These values corresponding, respectively, to a 100% paraffin mixture and a 100% olefin mixture are used for linear interpolation to the actual olefin content of !lie mixture. The following example illustrates the application of the MollieI' diagrams to a hydrocarbon mixture: Exmnple 6. Determine the work of compression 1 0 and final temperature when

the following mixture is compressed adiabatically from atmospheric pressure and 60°F to 50 psig: Average Molee. WI.

M ok Fraction

Component

CH, . ....... . C,H, ........ C2 H• ....... . C3 H. C3 H s ..... - .. C,H s C.H, 0 C.H, .......

0.050 .100 .150 .100 .200 .100 .200 .100 -1000

,

0.8 2.8 4.5 4.2 8.8 5.6 11.6 7.2 -45.5

Values corresponding to adiabatic compression from 1 aIm and 60°F to 4.4 atm were read from the individual charls and arc tabulated below: BTU/lb Compound

C,B, C~HIO

C2H.. C3H e

0.763 .680 .935 .780

By interpolation,

lif{ =

6fl M(h, - h,)

I,

S h,

112

301.5 295 300.5 303

338 321 363.5 3!2

'F 1M 135

1610 1510 1760 16iO

221

164

1600 BTU/mole and litVM

61VM

BTU/mole

=

625

570 850 675

620 for a saturated hydro-

carbon mixture of 45.5 malec. wt.

By extrapolation, li.H ~ 1610 BTU/mole and li/V,lf hydrocarbon mixture of 45.5 molec. wI.

=

632 for an unsaturated

By interpolation, li.H ~ 1603 BTU/mole and litVM

=

6z.t for a hydrocarbon

mixture of 45.5 malec. wt. containing 30% unsaturates. 10 Change in enthalpy which includes the difference between the work of expulsion and work of admission.

DATA BOOK ON HYDROCARBONS

84

:. The theoretical work of compression is 35.2 BTU/lb and the final tempera.ture is 152°F. If other gases (H 2 , O2 , H 2 0, etc.) are present in a mixture, it is recommended that effective IJressures equal to 7rVYlfC be used to determine the total work of compression and final temperature of the hydrocarbon portion of the mixture. The inert g&ses usually may be assumed to be ideal and the w'ork of comprestlion and final temperature for this portion of the mixture calculated by the adiabatic compression formulas for perfect gases. The work of compression for the mixture is then evaluated by combining the change of heat content for the hydrocarbon portion with that for the inert gases on the basis of their mole fractions. In determining the final temperature of the mixture, it is assumed that the ch.ange in enthalpy of each portion from its final temperature to that of the mixture is equal and opposite in sign to the other. This method is illustrated by the following example: Example 7. Determine the work of compression and final temperature when

the following mixture is compressed adiabatically from 25 psig and 0°1" to 150 psig: Hydrocarbon Portion Average Molec. Wt.

Mole Fraction

Component

Average Molec. Wt.

-

-

1.0 1.6 4.5 11.0

0.500 .100

H, CH, C,H. C,H.

Mole Fraction

.ISO

.250 1.000

0.200 .300 .500 1.000

18.1

3.2 9.0 22.0 34.2

The effective pressuras to be used for the hydrocarbon portion of the mixture are: 7rei

=

7r.2

=

25

+ 14.7 14.7

150

vO.500 = 1.91 atm

+ 14.7 vO.500 -=

14.7

7.91 atm

Values read from the ethane and propane charts are tabulated below: BTU/lb Compound

S

C,H. CaR.

0.837 0.686

t,

h,

h,

294 278

340 307

6.H

OF

M(h, - hi)

MvM

121 92

1385 1280

663 610

THERMAL PROPERTIES

85

By interpolation, UTi = l35~ BTU/mole and t:.tVM = 6-1.7 for a saturated hydrocarbon mixture of 3-1.2 molee. wt. The corresponding final temperature for the hydrocarbon portion of the mixture is 111°F. For the H 2 portion of the mixture, the work of compression and final

tcm.pera~

ture arc calculated as follows:

MCp

=

UTi

=

2.016 X 3.46

~

K

6.96;

-.

=

6.97 99 6.97 - I.

1.40

~I R1'[("2)K~' - IJ

J( -

"1

[(1647)1.40-' I J

I. 40 X 1.99 X 460 __ . 1.40 1.40 - I 39.7 = 3200(1.502 - I) = 1610 BTU/mole =

I.fO -1

7'2 = (164.7)1:'lE: !4+ a.. ~ "« Hill' I 'H. l. ~ ';iit Z

It

l'

I'

t-!! I '

~

~l..

f

t ~!.

••

I"'"

t

;1'

j

'11"_

~



-I

I

, 1

••

1."

••••

Ih·

-.

~

'j

:.;..

, ......,

.; "" '"

; l i .•~'~ • •

I

....

.'!~

"..!~~

''t'

;',1

. ..• - I: .".

~;:; Ii .. '''' ...

t-yl:

I

I'"

'.'

t

1

J:.

I'

~::: ~:'



-..

"

'I'

'HI .. I.

".

.~.

.

!'t:

I.j. ~

" r • j.,....

.

. ,. . . . ,.;10

I. I

j..



,

1

I.t'r

·'t'

I • Ilp •



; ... , I', I· ,

'I'

t··

• j "" , , . ...I., • , • •••••.

,J: oil:

';

Htill ill "I+H: lj:Hl++hr; tii:nr !~ri ;LJ \ :: '-:~ ..!i;! jl~l L1':': :!;j ~;:l £::1 I';q :~!; ,I;1.n.. t; f. II r~ '; ~\' .~\., d . -.-.. • ;t; . ,.It.tl-.l n:j I,; .. , .pl ;'\,.; . .'~~~ .Li~'. I,U .. 1. ;,.~ I. ... ,t. r" 'jl ,~'. W ' 't. !lfl;' I'll r .II!' r 111't '1. li1 :t;" t,'~~n;,-,,! 1'1 '"" ,.W .'

350

56

340

s.,

330

300

5,3

"'

~

290

~

280

'"0z

270 260 250

100

55

320 310

110

59

5f! >
>iii

'">-

230

'"

70

8

..J

!!! >

0

'">-

80

w

0

8w

240

90

60

-

in 0

tJ:;

UI

....

(j'

80

f--

t

z

850

~

....

f-!:::

T

2300+ 2200t 2100

f T

2000

0

0'" ":;'"

'" :;

70

60

0

al

>ow .-(/) -ct

.-2:e> O(/)

8.,

::l::l

°0 Zw

Oz Oct

...J...J ct...J

::E~ wJ:::E

a:(/)

.o o

0 0

t')'s

122 177 262 350 443 538 636 752 (905) -

129 177 273 369 465 561 656 752 848.

-7

0.40

-

.36 .34 .34 .34 .33

-11

-19 -22 -23 -20

-

-

.33

+57

Flash Vaporization (F) t>t'

FRL

Curve

-3

233 265 329 393 457 521 585 649 713

230 265 325 387 450 513 578 649 .732

-

-4 -6 -7 -8 -7 -

+19

The flash reference line and the atmospheric flash curve of the original crude are-shown on Figure 1. Proceeding from (1), the flash curve of the original crude, the atmospheric flash curves of the stripped and unstripped reduced crudes are derived by the method outlined in the text: (2) Vol. Av. B.P. of whole crude

=

262

+ 538 + 905 = 568°F 3

226

DATA BOOK ON HYDROCARBONS Mean Av. B.P. of whole crude = 568 - 70 = 498°F (Section 2) Molec. wt. of whole crude = 197 (Section 3) Vol. Av. B.P. of 65% overhead = Slope of DRL (65% overhead) =

203

+ 373 + 3

558

= 378°F

495 - 139 60 = 5.9°F/%

Mean Av. B.P. = 378 - 38 = 340°F Molec. wt. of 65% overhead = 139

(3)

(4)

(5)

(6)

Per 100 Gal of Crude Moles of crude = (6.98 X 100)/197 = 3.55 Moles of overhead = (6.57 X 65)/139 = 3.07 Moles of reduced crude 0.48 Partial pressure of reduced crude at the dew point of the original crude 0.48 = X 1 = 0.135 atm. 3.55 The 40% point on the reduced crucle flash curve corresponds to 65 + 0.40 X 35 = 79% or 722°F on the flash curve of the original crude. By extrapolation from 0.135 atm. to 1 atm., the 40% point on the atmospheric flash curve of the reduced crude is 900?F. The atmospheric flash curve of the stripped reduced crude is drawn through the extrapolated point parallel to the 65-100% portion of the flash curve of the original crude. This reduced crude flash curve may be converted to percent on reduced crude by proportioning the 65-100% yield on original crude to 0-100% on reduced crude. Both curves are shown in Figure l. The front end of the atmospheric flash curve on the unstripped reduced crude is constructed by drawing a smooth curve from the 65% point on the flash curve of the odginal crude to the 20% point on the flash curve of the stripped reduced crude as shown in Figure 1. This curve is also given on the basis of 0-100% reduced crude.

GENERAL REFERENCES Edmister and Pollock, Chem. E7l{J. Progress 44, 905 (1948). Katz and Brown, Ind. E7l{J. Chern. 26, 1373 (1933). Packie, Trans. Am. Inst. Chern. Engrs. 37, 51 (1941).

EQUILIBRIUM FLASH VAPORIZATION

227

1100

1000

900

800

700

600

!lOO

400

300

200

10

20

30

40

50 FIGURE

60 1

70

80

90

100

7

PREDICTION OF FLASH REFERENCE LINE FROM DISTILLATION REFERENCE LINES

--JfJ1 .•r

~

5

.. ~.

7

g

w

4

2

. FLASH AND DISTILLATION REFERENCE LINES (FRL AND DRLl ARE STRAIGHT LI NES THROUGH THE 10% AND 70% , POINTS. THE TEMPERATURES AT THE

I -

o

i li

If

3

50% POINTS REFER TO THESE REFERENCE LINES.

3

6

7

8

9

10

II

,,'

2

m I .

12

.~ I~

40

~ ~

E : ... RE~G~l~it!fI:~:':I~-'~ 1

.... ;~~ 2

3

o

60

t



T

i ~If

U4

6

S

7

!'iii R,f;f,Et EN' ~...

!; - .

8

9

10

-

II

12

11- -

40

20

..

II

BI

O -20

-40

-60

-40

2

3

5

6

228

tl±I

.

~

.

, ,

~

,

.

;

, tmI,

I

: t 1.

..

,

Mm c

1

,

, •

too

1.00 .80

lm!l

100

I..

PREDICTION OF FLASH CURVE FROM ITS REFERENCE LINE

I



1'-'-

fmII

mIl tmn

,

II

11UlI



_ 11$ CRUDE , .

.60

ASSAY (T.B.P.) DIS TlLLATION

'gmj

,20

11m!

-,,",p .. I:J±

t"

10

,

1.00

20

30

II

40

50

11m

60

70

90

100

.

·lmJU

.JI#J:

"* lIV IS THE DEPARTURE

n5

.60 f:fk,'ll

I

OF THE ACTUAL FLASH AND D1STI LLATION CURVES FROM THEIR RESPECTIVE REFERENCE LINES. WHILE THE INDIVIDUAL (lIl')'S MAY BE EITHER PLUS OR MINUS, THE RATIO IS ALWAYS POSITIVE.

,

~ , ' '. ..

~l'mIJ

Iflill

.20

o

80

lEE ' 10% (A.s:tM.) DISTILLATION

.80

!l0

•• • II

10

W

30

·w_ •

40

~

60

ro

229

~

90

100

.1

Section 14

FRACTIONATING TOWERS In order to simplify the work involved in making stepwise calculations for the rectification of binary and multieomponent systems, Gilliland' has presented an empirical correlation between theoretical steps and reflux ratio. To use the Gilliland correlation to predict the number of theoretical plates for a given reflux ratio, the minimum number of steps at total reflex and the minimum reflux ratio are required. Minimum Number of Theoretical Steps When a separation is specified with respect to only two components of a multicomponent mixture, the lower boiling of these two components is designated the light key component and the higher boiling the heavy key component, and the minimum number of steps can be calculated by the well-known Fenske equation 2 as follows: a log

8 Af=

(XXLKlV LKD) (X HKlV) XHKD log aLK

[aLK]SM =

or

(XXLKlV LK D) (X HKlV) XHKV

(1)

(1 a)

After equation (1) is solved for 8M , the latter may be substituted in this equation along with the distribution of either key component to prcdict 4 the distribution of the other components, or

(X

,.D) X HKlV) = 8 M log aL XX log ( -

(2)

LK log (X lilY) (X LKD) = 8 M log (a ) XlID XLKlV all

(3)

!-IV

11 K D

Likewise,

In any of the above equations, moles per 100 moles of feed may be replaced by total mo)es, or volume or weight units since in any of these conversions the multiplying factors cancel out. Gilliland, Ind. Eng. Chern. 32, 1220 (1940). Fenske, Ind. Eng. Chern. 24, 482 (1932). A table o( nomenclnture is given on page 243. • This equation may be used (or any pair o( component8.

1

2 S

230

FRACTIONATING TOWERS

2H1

When the dcgree of separation is specified for more than two components, equation (1) must be applied to all critical combinations of these components and the maximum SJ/ determincd for the most difficult case. If the separation is specificd with respect to the total quantity of two or more components, as in the case of Examplc 1, trial and eITor is required for thc solution of SjJ. It should be pointed out that the concentrations calculated by equations (2) and (3) actually apply only to the separation at total rcflux and, with the

exception of the two key components, there will be some variation of thc degrec of separation with the reflux ratio. As the rcflux ratio decreases, there is some improvement in separation betwecn light and heavy componcnts boiling outside the range of the kcy components and some deterioration in the separation of components boiling intermediate bctween the kcy components. However, in so far as the present procedure is concerned, the distillate and bottoms compositions for other reflux ratios are assumed to be thc Same as those calculated for total reflux. Minimum Reflux Ratio

Gilliland 5 has proposed several diffcrent formulas for predicting minimum reflux ratio and all have the disadvantage of being composcd of a number of complex terms in addition to requiring trial and error for solution. Although all these equations appeal' to give satisfactory rcsults, the tcrms are so complcx that it is difficult to bc ccrtain that therc arc no numerical crrors in thcir application. In order to apply the Gilliland method with greater facility, the following equation was developcd for predicting the minimum rcflux ratio of a multicomponent system: (O/Dhf

+ 1=

(aLKTf.FC + ('(LK - 1

1) (XLKD ILK

XIIKD)

(4) (O/D)M can be calculated for two arbitrary states of feed vaporization: 1. "Liquid" feed, cOITesponding to vaporization of the feeu equivalcnt to the

fraction of the feed lighter than the light key component. For the components lighter than the light key, h = ZL/aL and for the light key and heavier components, ILK = ZLK, and III = Zfl.6 Gilliland, Ind. Eng. ehe",. 32, HOI (1940). ]£ components, intermediate between the two key components, are present, they are considered ei her light or heavy componen'" depending upon which key their volatility more nearly approaches. In the case of "liquid" feed, I L = Zl. and I II = ZH for these intermediate components; in the case of "vapor" feed, I L = ZL/aL and I II = ZuaH/aLK. 5

8

232

DATA BOOK ON HYDROCARBONS

2. "Vapor" feed, corresponding to vaporization of the feed equivalent to the fraction of the fecd consisting of the hcavy key component and lighter. For the components lighter than the heavy key, If_ = ZL/CtL and ILK = ZLK/aLK and for the components heavier than the heavy key, 1/1 = ZI/.6 After the minimum reflux ratios have been calculated for "liquid" and "vapor" feeds, the minimum rcflux ratio for the actual vaporization of the feed can be calculated by direct interpolation or extrapolation. However, extrapolation beyond 50% of the difference between "liquid" and "vapor" feed may lead to serious deviations. The first term of the right-hand side of equation (4) is the same as for binary mixtures, and the equation reduces to the cquivalcnt of a binary mixture whon all light components other than the light key have infinite volatility and all heavy component other than the heavy key have zero volatility. Under these circumstances the equation is exact when hK is taken as the ratio of the two components in the liquid phase of the feed. That is, if the feed is introduced as a liquid at its bubble point, hK = ZLK, which is the ratio of the two components in the feed; if the feed is introduced as a vapor at its dewpoint, hK = Z',K/OtLK, which is the ratio of the two components in the cquilibrium liquid. For intermediate stages of vaporization hK can be calculated from the flash vaporization formula, although direct intcrpolation of the minimum reflux ratio on the basis of percentage vaporization between thc sat urated liquid and saturated vapor feeds gives values only slightly in error on the conservative side. In the case of multicomponent mixturcs, equation (4) is semi-empirical since it was necessary to make simplifying approximations in its derivation. Furthermore, the exact values of the various 1's cannot be calculated directly from the composition and state of vaporization of the feed, since the liquid on the feed platc is not identical to the liquid phase of the fecd as in the case of a binary mixture. As a result, it was necessary to define the 1's empirically for two states of fced vaporization, arbitrarily choscn to simulate a binary mixture of the two key components, and then intcrpolatc or extrapolatc to the minimum reflux ratio corresponding to the actual vaporization of the fced. Equation (4j has been checked for a number of multicomponent systems on which the minimum reflux ratio was determined by stepwise trial and error calculations. Generally, unusual systems were chosen with respect to composition and relative volatility in order to reveal the maximum deviations ever likely to be encountercd in practice. The agrecment was quite satisfactory as the average deviation was less than -+-5% and the maximum about lOra. The latter occurred at the limit of extrapolation relative to the arbitrary feed states.



FRACTIONATING TOWERS

233

Also, the minimum reflux ratio was calculated for these same systems by the Colburn method 7 with about the same degree of accuracy. It should be pointed out that the latter gave better results than equation (4) when the relative volatilities and compositions were not so abnormal as the systems selected. However, under these circumstances both methods were quite accurate as the deviations seldom exceeded a few percent, and the present equation has a distinct advantage in that it is explicit and does not require trial and error. Both methods are quite sensitive to the selection of key components, and the selection of the wrong key components can lead to a much greater error than is inherent in either method. If the desired separation is between adjacent components, there is usually no doubt about selecting these as the key components. However, if there are additional specifications relative to other components, it may be necessary to try two or more combinations of key components to make sure that the minimum reflux ratio is sufficient to fulfill all specified conditions. Correlation of Theoretical Steps with Reflux Ratio As mentioned at the beginning of this section, Gilliland correlated the results of a large humber of stepwise calculations on various binary and multicomponent mixtures by plotting [S-SM) I[ S + 1] ~ (S) against [(OlD) - (OID))f)/[ OlD + 1)-F (OlD)

and found that all points could be represented by a single curve irrespective of the type or degree of separation. These points, along with about half again as many additional points, were replotted, and the best curve through them was essentially the same as Gilliland's original correlation. In arriving at the coordinates for the additional points the minimum reflux ratio was calculated by equation (4); therefore these points are a criterion of the present method as well as the curve itself. In no case did the deviations exceed either 3 theoretical steps or 15%, and the average deviation was less than 1 theoretical step and also less than +50/0. To take care of the maximum deviation it is recommended that in any design the number of theoretical steps predicted ~ from the correlation on page 244 be increased by either 3 theoretical steps or 100/0, whichever is greater. Plate Efficiency Because of the large number of factors which undoubtedly influence the plate efficiency of a fractionating tower, any fundamental formula accounting for even the most important variables must necessarily be quite involved. For this reason, a simple empirical correlation of the limited data on hydrocarbon mixtures seemed to be the most promising method of predicting plate efficiency. 7

Colburn, Trans. Am. Inst. ehem. Engrs. 37, 805 (1941).

DATA BOOK ON HYDROCARBONS

234

Gunness 8 correlated the results of several tests on petroleum mixtures on the basTs of vapor pressure of the liquid. As he points out, this is a method of indirectly correlating plate efficiency with liquid viscosity since viscosity of pure hydrocarbons and narrow boiling fractions is an approximate function of vapor pressure over a fai,rly wide range of vapor pressures. In view of the consistent results obtained by Gunness, pla.te efficiency was plotted directly against fluidity (reciprocal viscosity) for a number of tests on commercial towers including those upon which Gunness based his curve. The curve on page 245 represents this correlation. While the overall plate efficiency exceeds 10070 at fluidities greater than 9 Cp-1, this is not inconsistent as the flow of the liquid across the plates results in concentration gradients which may achieve a greater degree of fractionation than predicted by stepwise calculations in which the liquid is assumed to leave the plate in equilibrium with the composite vapor. Lewis 9 has shown theoretically that different combinations of liquid and vapor concentration gradients across the plate may give overall plate efficiencies as high as 200--300'10 when based on stepwise calculations. There is no reason to believe that this correlation applies to mixtures other than hydrocarbons, and with the exception of alcohol-water mixtures there are too little data available to afford a comparison. Although there is considerable variation in the alcohol-water data, there is some indication that plate efficiencies are somewhat greater than for hydrocarbons of the same viscosity. Location of the Feed Plate As a simple approximation for locating the feed plate, it may be assumed that the proportion of actual plates above the feed will be the same as that required to effect the same separation between the key components at total reflux. That is, the number of theoretical steps at total reflux is calculated for the concentration change in the key components between the feed and distillate compositions. It is then assumed that the ratio of this to the total number of theoretical steps at an infinite reflux ratio is the same as the ratio of actual plates above the feed is to the total number of plates. Application of this method is illustrated by Example 1. In some cases where there are oritical components other than the two key components, it may. be necessary to check the total reflux steps above and below the feed on the basis of these components, since the optimum location of the feed plate will be different with each pair of components. Usually the separation of components other than the key components is so complete that only the latter need be considered. 8 9

Gunness, Sc.D. Thesis, Mass. Inst. Tech. (1936). Lewis, Ind. Eng. Chern. 28, 399 (1936).

FRACTIONATING TOWERS

235

Packed Towers The charts on pages 246 to 248 giving the H.E.T.P., capacity and pressure drop in packed towers are self-explanatory. Since practically all of the H.E.T.P. data were on towers less than 12 in. in diameter, caution should be used in the design of larger towers. One of the greatest sources of inefficiency in a packed tower is poor liquid distribution. If good distribution can be achieved by efficient distributors, the extrapolations may be used for larger towers with reasonable assurance. Example 1. At an operating pressure of 100 psig determine the number of plates and reflux ratio required to separate the mixture given below so that the bottoms contain at least 90ro of the butenes-2 present in the feed and at the same time have an isobutene content not greater than 5%:

Component

i-C,H lo i-C,H s C,H s-1

C.H IO t-C,H s-2 c-C,H s-2

Feed (Mole %) 40.0 20.0 15.0 5.0 10.0 10.0 100.0

(1) Dewpoint of Distillate and Bubble Point of Bottoms

In order to calculate the average volatilities, the dewpoint of the distillate and bubble point of the bottoms must be found by trial and error using assumed compositions. These are tabulated below.

Moles Per 100 Moles of Feed Component i-C.H LO i-C.H, C.H...1 C.H,o t-C.H,·2 .,.C.H...2

Mole Fraction

Feeti

Distillate

Bottoms

Distillate

Bottoms

40.0 20.0 15.0 5.0 10.0 10.0 100.0

39.3 18.7 13.0 1.0 1.5 0.5 74.0

0.7 1.3 2.0 4.0 8.5 9.5 26.0

0.530 .253 .176 .014

0.027 .050 .077 .154 .327 .365 1.000

.OW .007 1.000

As a first trial, assume the dewpoint of the distillate is 14tl°F at 7.8 atm (114.7 psia) .

DATA BOOK ON HYDROCARBONS

236

First Trial Component

YD

':"C,H IO i-C,H, C,H,-I C,H,o t-C,H...2 c-C,H ...2

a'D·

Pt

:i;

HO°F

140°F

"yiP

1.29 1.155 1.13 1.00 0.97 0.91

8.4 7.5 7.35 6.5 6.3 5.9

0.493 .263 .187 .017 .025 .009 0.994

0.530 .253 .176 .014 .020 .007 1.000

• Relative volatilities to C 4H 10 or (0")'8 aTC used as a matter of convenience; then, the (a'.,.)'s are converted to (a••)'s, the relative volatilities to t-C,H,-2, which will be seleetcd as the heavy key component. t Computed from the fugacity function of butane multiplied by the relative volatilities.

Since the sum of the x's is 0.994 instead of 1.000, the assumed temperature should be lowered slightly, but the difference would be so small (less than l°F) that the change in the (a'Jl) 's would be imperceptible. Consequently, 140°F will be used as the dewpoint of the distillate. The bubble point of the bottoms is assumed to be 165°F at 8.0 atm 10 for the first trial.

Second Trial

First Trial

Component

i-C,H,o i-C,H, C,H ... I C,H,o t-C,H ...2 c-C,H,-2

XIV

0.027 .050 .077 .154 .327 .365 1.000

a' w· 165°F

pt

Y

165°F

Pxl..

1.26 1.14 1.115 1.00 0.97 0.915

10.7 9.7 9.5 8.5 8.25 7.8

0.036 .061 .091 .164 .337 .356 1.045

a'w· 160°F

160°F

Y

1.265 1.145 1.12 1.00 0.97 0.915

10.25 9.3 9.1 8.1 7.85 7.4

0.035 .058 .088 .156 .321 .338 0.996

Pt

• Relative volatilities to C~HIO or {«')'8 are used ns a mnttcr of convenience; then, the (a'av)'S are converted to (aU\')'s, the relative volatilities to t-C IH:;-2, which will be selected as

the heavy key component. t Computed from the fugacity funet:on of hutane multiplied by the relative volatilities.

The bubble point of the bottoms wiil be taken as 160°F. The relativ~ volatilities are averaged and converted to t-C~ H s -2 as the heavy key in the following table: \0

After allowing 3 Ib/sq in. as the approximate pressure drop through the tower.

FRACTIONATING TOWERS

,

,

,

aD

aw

aA

Component

140°F 7.8 at.m

160°F 8.0 atm

150°F 7.9 atm

i-C,H ID i-C,H. C,H.-l C,H ID t-C,H8-2 c-C,Hg-2

1.29 1.155 1.13 1.00 0.97 0.91

1.265 1.145 1.12 1.00 0.97 0.915

1.275 1.15 1.125 1.00 0.97 0.91

237

,

a a.

(a' Da'wa'.A.)~!l

1.275 1.15 1.125 1.00 0.97 0.91

aa.

1.315 1.185 1.16 1.03 1.00 0.94

(2) Minimum Theoretical Steps (Total Reflu:t) The minimum number of theoretical steps by which the desired separation can be accomplished is calculated as follows:

Let

t = moles of t-C 4 H s-2 in the distillate per 100 moles of feed 10 - t = moles of t-C 4 H s-2 in the bottoms per 100 moles of feed

Since 90% of the butenes-2 must be.retained in the bottoms, the cis-butcne-2 content of the distillate and bottoms will be: (2 - t) moles in the distillate per 100 moles of feed (8 + t) moles in the bottoms per 100 moles of feed

and

Using the previously assumed values of 18.7 moles of isobutene in the distillate and 1.3 moles in the bottoms, the following equations must be satisfied:

C1~;) CO t- t) 18.7) (~) ( 1.3 2- t

= (1.185)8 M =

(1.185)8 0.94

M

A trial and error solution of these equations shows that they are satisfied by SM = 25.5 and t = 1.62. The distribution of the other components can be calculated from SM and the distribution of t-C 4 H s -2. i-C 4 H IO : Let u = moles of i-C 4 H IO in bottoms (

40 u

1') (8.38) = (1.315)25.5 = 1075 1.62

= 0.19 moles of i-C4 H IO in the bottoms C 4H s-1: Let v = moles of C4 H s-1 in the bottoms (

15 v

v) (838) 1.62

=

(1.16)25.5

=

44

v = 1.58 moles of C 4 H s-2 in the bottoms

DATA BOOK ON HYDROCARBONS

238

C.H IO : Let w = moles of C.H IO in the bottoms (

5W

W)(8.38) = (1.03)25.5 = 2.12 1.62

W

= 3.55 moles of C.H 10 in the bottoms

The percentage of i-C.Hs in the bottoms will be: (0.19

+ 1.3 + 1.58 ~33.55 + 8.38 + 9.62) 100 = 5.3%

In order to meet a maximum of 5.0ro i-C 4 H s specified for the bottoms, it is necessary to reduce the 1.3 moles to 1.22 moles in the bottoms. This would require an increase in SM to 25.8 which would modify the distribution of the other components. However, the latter change is so slight that it can be neglected. The composition of the overhead and bottoms will then be: Mole Fraction

Moles Per 100 Moles of Feed Component

i·C,H,o i·C,H, C,H..1 C,H IO /·C,H,-2 .,.C,H..2

Fecd

Distillate

Bottoms

Distillate

Bottoms

40.0 20.0 15.0 5.0 10.0 10.0

39.81 18.78 13.42 1.45 1.62 0.38 75.46

0.19 1.22 ·1.58 3.55 8.38 9.62 24.54

0.528 .249 .178 .019 .021 .005 1.000

0.008 .050 .064 .145 .342 .391 1.000

(3) Minimum Reflux Ratio Since the critical separation is between isobutene and the butenes-2, the former is naturally selected as the light key component and trans-butene-2, since it is more volatile than the cis-butene-2, as the heavy key component. Butene-l is considered a light intermediate component because of the proximity of its relative volatility to that of isobutcne; normal butane is considered a heavy intermediate component since its relative volatility is nearer to the heavy key than the light key. The following tabulation gives the necessary information for calculating the minimum reflux ratios for the two arbitrary states of feed vaporization: Mole Fraction Component

':-C,H,o i-C,H, C,H ..1 C,H IO t-C,H,.2 .,.C,H..2

Type

L LK L H HK H

"'BV

Feed

Distillate

Bottoms

0.400 .200 .150 .050 .100 .100 1.000

0.528 .249 .178 .019 .021 .005 1.000

0.008 .050 .064 .145 .342 .391 1.000

1.315 1.185 1.16 1.03 1.00 0.94

HLiquid tl

lIVapor"

Feed

Feed

3.04 2.00 1.50 4.00

3.04 1.69 1.29 3.48

2.00

2.00

FRACTIONATING TOWERS

239

"Liquid" jeed-40% vaporized (O/D)

+1 = M

1.185 X 2.00 + 1.0 (0.249 _ 1.185 _ 1.0 2.00

+ 01. 31 5 (0.528 .315

1.03

(O/Dhf

° 1) . 02

3.04 X 0.021)

(0.249

+ 1.16 (0.178 0.16

)

0.94

+ 1.185 - 1.03 4.00 - 0.019 + 1.185 - 0.94 = 1.88 + 1.94 + 1.07 + 0.29 + 0.46 - 1 = 4.64

1.50 X 0.021)

(0.249 ) 2.00 - 0.005

"Vapor" jeed-90% vaporized (O/D) M

+1 =

1.185 X 1.69 + 1.0 (0.249 - - - 0.021 ) 1.69 1.185 - 1.0 1.16 ( + 1.94 + 0.16 0.178 1.03

+ 1.185 (O/D)M

1.03

1.29 X 0.021)

(0.249 ) 3.48 - 0.019

= 2.06 + 1.94 + 1.10 + 0.35 + 0.46

- 1

+ 0.46

= 4.91

Assume that the feed is sufficiently preheated to vaporize a percentage equivalent to the distillate or 75.4670. By interpolation, the minimum reflux ratio corresponding to this feed vaporization is: (O/D)M

75.46 - 40)

= 4.64 + ( 90 _ 40

(4.91 - 4.64) = 4.83

(4) Theoretical Steps vs. Reflux Ratio

Using the values determined in preceding sections for minimum theoretical steps,. 25.8, and for minimum reflux ratio, 4.83, the number of theoretical steps for various reflux ratios can be predicted from the correlation on page 244: OlD

F(OID)

4.83 5.25 5.75 6.50 7.50

0.067 .136 .223 .314

.

-

-

(8)

-

0.570 .502 .430 .366

-

8

Theoretical Platea-

61.3 52.7 46.0 41.3 25.8

60.3 51.7 45.0 40.3 24.8

.

.

- The reboiler i. considered the equivalent of one theoretical step. With a partial instead of a total condenser, a second theoretical step also could have been deducted.

DATA BOOK ON HYDROCARBONS

240

(5) Number of Actual Fractionating Plates

To predict the number of actual plates it is necessary to determine the average viscosity of the liquid on the plates. Since the temperature difference between the top and bottom of the tower is so small, the average viscosity may be taken as the viscosity at the average temperature. For this purpose the viscosity of butane at 150°F will be used. Viscosity of C.H IO @ 150°F = 0.216 cs "" 0.216 X 0.523 = 0.113 cp Fluidity = 1/0.113 = 8.9 Cp-l j Plate efficiency = 99% Using a plate efficiency of 99% the number of actual plates is computed for various reflux ratios: OlD 4.83 5.25 5.75 6.50 7.50

.,

S

Theoretical Steps

Actual Plates

.,

.,

.,

61.3 52.7 46.0 41.3 25.8

60.3 51.7 45.0 40.3 24.8

60.9 52.2 45.5 40.7 25.0

The number of actual plates is plotted against reflux ratio in Figure 1. A reflux ratio of 6.50 to 1, or 1.35 times the minimum, is selected. The number of actual plates corresponding to this reflux ratio is 45.5 so that a 50-plate tower would be required. (6) Location of the Feed Plate

The number of plates above the feed is based on the proportion of theoretical steps at total reflux which would be required to effect the change in concentration of the key components between the feed and distillate. This proportion is applied to the actual number of plates (including the reboiler) to determine the number above the feed plate. In order to take into account any appreciable difference in relative volatility above and below the feed, the relative volatility used for calculating the steps at total reflux between feed and distillate is the geometric mean of aD and a" or, an

=(

1.15)~i 1.155 0.97 X 0.97

=

1.19

The number of total reflux steps which would be required between the feed and distillate is calculated by the following equation: l8.78) (~) ( 20 1.62

=

1.19n

= 5.79' '

n = 10.1

FRACTIO ATING TOWERS

241

50

40

30

5

4

7

6 FIGUllE

8

1

Number of actual plates above the feed would then be: 10.1 (50 25.8

+ 1)

=

20

The vaporization of the feed can be taken into account by adding the fraction vaporized to n since 10010 vaporization would be equivalent to a theoretical step at total reflux. This would change the proportion of plates above the feed as follows: (

10.1 + I • 0.75) ( 00" +) 1 = 21.4 pates above the feed 20.8

Feed lines would probably be installed above the 2'!th, the 28th and 32nd plates from tile bottom of the tower.

242

DATA BOOK ON HYDROCARBONS GENERAL REFERENCES

Atkins and Franklin, Refiner Natural Gasoline Mfgr. (Jan. 1936). Brown, Sanders, Nyland and Hesler, Ind. Eng. Chem. 27, 383 (1935). Brown and Souders, Oil and Gas J. 31, 34 (1932). Chilton llnd Colburn, Trans. Am. Inst. Chern. Engrs. 26, 178 (1931). Elgin and Weiss, Ind. Eng. Chem. 31, 435 (1939). Fenske, Lawroski llnd Tongberg, Ind. Eng. Chern. 30, 227 (1938). Fenske, Unpublished data, Pennsylvania State College. Gilliland, Ind. Eng. Chent. 32, 918, 1101, 1220 (1940). Gunness, Ind. Eng. Chern. 29, 1092 (1937). Lewis and Wilde, Trans. Am. Inst. Chern. Engrs. 21, 99 (1928). Perry, "Chemical Engineers' Handbook," pp. 829-832, McGraw-Hill Book Co., New York, N.Y. (1941). Sherwood, Shipley and Holloway, Ind. Eng. Chem. 30, 765 (1938). White, Tram. Am. Imt. Chem. Engrs. 31, 390 (1935).

FRACTIONATING TOWERS

243

.vomenclature X x

moles of any component in distillate or bottoms per 100 moles of feed mole fraction of any component in liquid y mole fraction of any component in vapor D moles of distillate per 100 moles of feed o moles of reflux per 100 moles of feed OlD reflux ratio (OIDhl minimum reflux ratio corresponding to S = 00 S number of steps from still to distillate 8,1/ minimum number of steps corresponding to OlD = 00 P number of theoretical plates; with a partial reboiler and partial condenser, P = S - 2, and with a partial reboiler and total condenser,

P=S-l

ZH aD alV

LK HK L

H D W n m

ratio of mole fraction of any light component to heavy key component in the feed ratio of mole fraction of light key component to any heavy component in feed relative volatility of any component to heavy key at the dew point of the distillate relative volatility of any component at the bubble point of the bottoms relative volatility of any component at the arithmetic average temperature of the dew point of the distillate and the bubble point of the bottoms mean relative volatility of any component, (aD' alV . a.4)fi used as a subscript to refer to the light key component used as a subscript to refer to the heavy key component used as a subscript to refer to any light component used as a subscript to refer to any heavy component used as a subscript to refer to the distillate used as a subscript to refer to the bottoms used as a subscript to refer to the plates above the feed used as a subscript to refer to the plates below the feed

1IIIriCORRELATION OF THEORETICAL STEPS WITH REFLUX RATIO MULTICOMPONENT AND BINARY MIXTURES

.9

.8

.6

.5 .4 .3

.2 .1

244

OVERALL PLATE EFFICIENCY vs. FLUIDITY OF LIQUID ON PLATES

120

1.11111

120

110 100

100

90

90

80·

80

70

ONLY DATA 00 HYDROCARBON MIXTURES WERE USED IN THIS CffiRELATION, AND THERE WERE INSUFFICIENT DATA ON OTHER TYPES TO JUSTIFY A MORE GENERAL USE. HOWEVER, THERE WERE SOME EVIDENCE THAT THE CURVE IS A LITTLE CONSERVATIVE FOR ALCOHOL - WATER

60

40_ 50

70 60 50

MIXTURES.

40

30 20

: :• • ll!IIflE:IJffi 10

_ _ _ _ _ '0 2

3

4

5

6

7

245

8

9

10

II

12

13

14

2

3

45678910

20 .:_.

30 :t';

HEIGHT EQUIVALENT TO A THEORETICAL PLATE

4a!f11i1 2

(I) WHILE THIS CORRELATION WAS DE'

VELOPED fROM DATA ON RASHIG RINGS AND 8ERL SADDLES, IT PR08ABLY APPLIES TO OTHER SIMILAR TYPES OF HOLLOW PACKING. (2) VALUES OF H.E.lP. FROM THIS CHART CORRESPOND TO THE MAXIMUM TOWER CAPACITIES GIVEN BY THE CHART ON THE OPPOSITE PAGE. FOR THE VALUES OF HE.T.P. AT CAPACITIES BETWEEN 80% AND 100% OF THE MAXIMUM, DIVIDE H.E.T.P. FROM CURVES BY THE FRACTION OF ULTIMATE CAPACITY (.80-1.00) AT WHICH THE TOWER WILL OPERATE.

mIR~12ioll~3IoI14~oll~ 60

• 246

7'08090 I

·1

* USE VALUES OF S/F 3 FROM CUR\IE

~-d

FOR RASCHIG RINGS. BERL SADDLES W PACKING UP 10 2 INOiES IN SIZE. FOR SIZES GREATER'THAN 2 INCHES, USE INDIVIDUAL

.2

.,.,, ~li

VALUES OF SAND F.

.1

lilt! 1)

r r 1111

'r-

.08

~

jj

.06 .05 .04

.

e:l3

.f+

"t~

I 1 I

ltllllltltrTl1111 ;~':+fTo"",

1000 800

.

,I';r.j

.....···'11:;::;:1-0-

~~i:-:'l:::.Jx·

""'-:1= ••I"4T'

·r.... '· ....

tI

600 500 400

~

~

=

300

'J="~J';+: -j:,l

r,

rt~'·-I--t"""'

'it

L-

..



n

"llOUIO-

m.

"

80 - reNSlTY OF VAPOR - LeS/CU· FT. 8L fl L1QUIO" .. Uo-SUPERFICIAL VAPOR VELOCITY AT INITIAL FLOODING-FT/SEC.

+

S*-SURFACE AREA OF PACKING- SO. FT./CU.FT. TOWER VOLUME F*-F'RACTION OF FREE VOLUME IN PACKING

.002,

OF LIQUID - CENTIPOISES -GRAVITY CONSTANT-32.2 FT./SEC~

...

SHERWOOu

_......

~MlrLEY

..mI·;,·;:; I ;'i-;t:;:-..J:, AND HOLLOWAY. IND . ENG. DiEM. 30. 765 (1930) ~ rF!¥:tJ~ lii r:'~ •., !

·/+l!;I"·'

."

001

.01

.OZ

.03.04

.06 .08 .fO

.2

.3

I'

I

.,

.4.5.6.7.8.910

"'.

'n

100 80 60

.M -VISCOSITY 9

" .200

f"'";

SUPERFICIAL MASS vaoCITY OF VAPOR-L8SISEC/SO.FT.

2

3

·

~

AP/H