Coker A.K. Fortran Programs for Chemical Process Design, Analysis, and Simulation
Подождите немного. Документ загружается.
45
C
50
55
C
999
CALL OUTPUT (L, G, DG, DL, PF, DELP, LVIS)
PRINT THE RESULTS ONTO A PRINTER
CALL PAPER (L, G, DG, DL, PF, DELP, LVIS)
GO TO 999
CALL PROG4 (L, G, DG, DL, PF, LVIS)
PRINT THE RESULTS ON THE SCREEN
CALL OUTPUT (L, G, DG, DL, PF, DELP, LVIS)
PRINT THE RESULTS ONTO A PRINTER
CALL PAPER (L, G, DG, DL, PF, DELP, LVIS)
GO TO 999
CALL PROG5 (L, G, DG, DL, PF, LVIS)
PRINT THE RESULTS ON THE SCREEN
CALL OUTPUT (L, G, DG, DL, PF, DELP, LVIS)
PRINT THE RESULTS ONTO A PRINTER
CALL PAPER (L, G, DG, DL, PF, DELP, LVIS)
GO TO 999
CALL PROG6 (L, G, DG, DL, PF, LVIS)
PRINT THE RESULTS ON THE SCREEN
CALL OUTPUT (L, G, DG, DL, PF, DELP, LVIS)
PRINT THE RESULTS ONTO A PRINTER
CALL PAPER (L, G, DG, DL, PF, DELP, LVIS)
GO TO 999
CLOSE (UNIT = 3, STATUS='KEEP')
CLOSE (UNIT = i)
STOP
END
***************************************************************
PACKED-TOWER SIZING WITH A PRESSURE DROP OF 0.05 inch H20/ft
***************************************************************
SUBROUTINE PROG1 (L, G, DG, DL, PF, LVIS)
REAL L,LVIS
COMMON/TOWER1/ GC, PI
COMMON/TOWER2/ X, Y, Xl, Y1
COMMON/TOWER3/ AREA, DIA
COMMON/TOWER4/ G1, G2, VG
567
COMMON/TOWER5/ AREAF, DIAF
CONSTANTS FOR DELP=0.05
CO = -6.30253
C1 = -0.60809
C2 = -0.11932
C3 = -0.00685
C4 = 0.000320
WHERE 0.001 < Y < 0.15
S~I
= ALOGCX)
SUM2 = (ALOG(X))**2
SUM3 = (ALOG(X))**3
SUM4 = (ALOG(X))**4
SUMT = CO+(CI*SUMI)+(C2*SUM2)+(C3*SUM3)+(C4*SUM4)
Y1 = EXP(SUMT)
CALCULATE VAPOR RATE ib/s.ft^2
VALI = u
VAL2 = PF*(LVIS)**0.1
RATIO = VALI/VAL2
G1 = (RATIO)**0.5
CALCULATE TOWER AREA AND TOWER DIAMETER
AREA = G/(3600.*GI)
DIA = 12.0*((4.0*AREA)/PI)**0.5
CALCULATE Y CO-ORDINATE IN GR-109 R-5 PLOT
VAL3 = (GI**2)*PF*(LVIS)**0.1
VAL4 = DG*(DL-DG)*GC
Y = VAL3/VAL4
CALCULATE THE PACKED TOWER FLOODING CONDITIONS
CALL FLOOD (L, G, DG, DL, PF, LVIS)
RETURN
END
***************************************************************
PACKED-TOWER SIZING WITH A PRESSURE DROP OF 0.i0 inch H20/ft
***************************************************************
SUBROUTINE PROG2 (L, G, DG, DL, PF, LVIS)
REAL L,LVIS
COMMON/TOWER1/ GC, PI
COMMON/TOWER2/ X, Y, XI, u
COMMON/TOWER3/ AREA, DIA
COMMON/TOWER4/ GI, G2, VG
COMMON/TOWER5/ AREAF, DIAF
CONSTANTS FOR DELP=O.I
568
CO = -5.50093
C1 = -0.78508
C2 = -0.13496
C3 = 0.00134
C4 = 0.00174
WHERE 0.001 < Y < 0.15
SUM1 = ALOG(X)
SUM2 = (ALOG(X))**2
SUM3 = (ALOG(X))**3
SUM4 = (ALOG(X))**4
SUMT = C0+(CI*SUMI)+(C2*SUM2)+(C3*SUM3)+(C4*SUM4)
Y1 = EXP(SUMT)
CALCULATE VAPOR RATE, ib/s.ft^2
VALI = YI*DG*(DL-DG)*GC
VAL2 = PF*(LVIS)**O.I
RATIO = VALI/VAL2
G1 = (RATIO)**0.5
CALCULATE TOWER AREA AND TOWER DIAMETER
AREA = G/(3600.*GI)
DIA = 12.0*((4.0*AREA)/PI)**0.5
CALCULATE u CO-ORDINATE IN GR-109 R-5 PLOT
VAL3 = (GI**2)*PF*(LVIS)**0.1
VAL4 = DG*(DL-DG)*GC
Y = VAL3/VAL4
CALCULATE THE PACKED TOWER FLOODING CONDITIONS.
CALL FLOOD (L, G, DG, DL, PF, LvIs)
RETURN
END
***************************************************************
PACKED-TOWER SIZING WITH A PRESSURE DROP OF 0.25 inch H20/ft
***************************************************************
SUBROUTINE PROG3 (L, G, DG, DL, PF, LVIS)
REAL L, LVIS
COMMON~TOWER1~
GC, PI
COMMON/TOWER2/ X, Y, Xl, Y1
COMMON/TOWER3/ AREA, DIA
COMMON/TOWER4/ GI, G2, VG
COMMON/TOWER5/ AREAF, DIAF
CONSTANTS FOR DELP=0.25
CO = -5.00319
C1 = -0.95299
C2 = -0.13930
C3 = 0.01264
569
C4 = 0.00334
WHERE 0.001 < Y < 0.15
SUM1 = ALOG(X)
SUM2 = (ALOG(X))**2
SUM3 = (ALOG(X))**3
SUM4 = (ALOG(X))**4
SUMT = CO+(CI*SUMI)+(C2*SUM2)+(C3*SUM3)+(C4*SUM4)
Y1 = EXP(SUMT)
CALCULATE VAPOR RATE, ib/s.ft^2
VALI = YI*DG*(DL-DG)*GC
VAL2 = PF*(LVIS)**0.1
RATIO = VALI/VAL2
G1 = SQRT(RATIO)
CALCULATE TOWER AREA AND TOWER DIAMETER
AREA = G/(3600.0*GI)
DIA = 12.0*SQRT((4.0*AREA)/PI)
CALCULATE Y CO-ORDINATE IN GR-109 R-5 PLOT
VAL3 = (GI**2)*PF*(LVIS)**0.1
VAL4 = DG*(DL-DG)*GC
Y = VALB/VAL4
CALCULATE THE PACKED TOWER FLOODING CONDITONS.
CALL FLOOD (L, G, DG, DL, PF, LVIS)
RETURN
END
***************************************************************
PACKED TOWER SIZING WITH A PRESSURE DROP OF 0.5 inch H20/ft
***************************************************************
SUBROUTINE PROG4 (L, G, DG, DL, PF, LVIS)
REAL L,LVIS
COMMON/TOWER1/ GC, PI
COMMON/TOWER2/ X, Y, XI, Y1
COMMON/TOWER3/ AREA, DIA
COMMON/TOWER4/ GI, G2, VG
COMMON/TOWER5/ AREAF, DIAF
CONSTANTS FOR DELP=0.5
CO = -4.39918
C1 = -0.99404
C2 = -0.16983
C3 = 0.00873
C4 = 0.00343
WHERE 0.001 < u < 0.15
SUM1 = ALOG(X)
SUM2 = (ALOG(X))**2
570
SUM3 = (ALOG(X))**3
SUM4 = (ALOG(X))**4
SUMT = C0+(CI*SUMI)+(C2*SUM2)+(C3*SUM3)+(C4*SUM4)
Y1 = EXP(SUMT)
CALCULATE VAPOR RATE, ib/s.ft^2
VALI = YI*DG*(DL-DG)*GC
VAL2 = PF*(LVIS)**0.1
RATIO = VALI/VAL2
G1 = (RATIO)**0.5
CALCULATE TOWER AREA AND TOWER DIAMETER
AREA = G/(3600.*GI)
DIA = 12.0*((4.0*AREA)/PI)**0.5
CALCULATE Y CO-ORDINATE IN GR-109 R-5 PLOT
VAL3 = (GI**2)*PF*(LVIS)**0.1
VAL4 = DG*(DL-DG)*GC
Y = VAL3/VAL4
CALCULATE THE PACKED TOWER FLOODING CONDITIONS.
CALL FLOOD (L, G, DG, DL, PF, LVIS)
RETURN
END
***************************************************************
PACKED-TOWER SIZING WITH A PRESSURE DROP OF 1 inch H20/ft.
***************************************************************
SUBROUTfNE PROG5 (L, G, DG, DL, PF, LVIS)
REAL L,LVIS
COMMON/TOWER1/ GC, PI
COMMON/TOWER2/ X, Y, X1, Y1
COMMON/TOWER3/ AREA, DIA
COMMON/TOWER4/ GI, G2, VG
COMMON/TOWERS/ AREAF, DIAF
CONSTANTS FOR DELP=I.0
CO = -4.09505
C1 = -1.00120
C2 = -0.15871
C3 = 0.00797
C4 = 0.00318
WHERE 0.001 < u < 0.15
SUM1 = ALOG(X)
SUM2 = (ALOG(X))**2
SUM3 = (ALOG(X) 7"'3
SUM4 = (ALOG(X))**4
SUMT = C0+(CI*SUMI)+(C2*SUM2)+(C3*SUM3)+(C4*SUM4)
Y1 = EXP(SUMT)
571
C CALCULATE VAPOR RATE, Ib/s.ft^2
VALI = u
VAL2 = PF*(LVIS)**0.1
RATIO = VALI/VAL2
G1 = (RATIO)**0.5
CALCULATE TOWER AREA AND TOWER DIAMETER
AREA = G/( 3600. *G1 )
DIA = 12.0"( (4.0*AREA)/PI)**0.5
CALCULATE u CO-ORDINATE IN GR-I09 R-5 PLOT
VAL3 = (GI**2)*PF*(LVIS)**0.1
VAL4 = DG*(DL-DG)*GC
Y = VAL3/VAL4
CALCULATE THE PACKED TOWER FLOODING CONDITIONS.
CALL FLOOD (L, G, DG, DL, PF, LVIS)
RETURN
END
PACKED-TOWER SIZING WITH A PRESSURE DROP OF 1.5 inch H20/ft.
***************************************************************
SUBROUTINE PROG6 (L, G, DG, DL, PF, LVIS)
REAL L, LVIS
COMMON/TOWER1/ GC, PI
COMMON/TOWER2/ X, Y, XI, Y1
COMMON~TOWER3~
AREA, DIA
COMMON/TOWER4/ GI, G2, VG
COMMON/TOWER5/ AREAF, DIAF
CONSTANTS FOR DELP=I.5
CO = -4.02555
C1 = -0.98945
C2 = -0.08291
C3 = -0.03237
C4 = 0.00532
WHERE 0.001 < Y < 0.15
SUM1 = ALOG(X)
SUM2 = (ALOG(X))**2
SUM3 = (ALOG(X))**3
SUM4
=
(ALOG(X))**4
SUMT = C0+(CI*SUMI)+(C2*SUM2)+(C3*SUM3)+(C4*SUM4)
Y1 = EXP(SUMT)
CALCULATE
VAPOR RATE ib/s.ft^2
VALI = YI*DG*(DL-DG)*GC
VAL2 = PF*(LVIS)**0.1
RATIO = VALI/VAL2
572
G1 = (RATIO)**0.5
CALCULATE TOWER AREA AND TOWER DIAMETER
AREA
= G/(3600.*GI)
DIA = 12.0"( (4.0*AREA)/PI)**0.5
CALCULATE Y CO-ORDINATE IN GR-109 R-5 PLOT
VAL3 = (GI**2)*PF*(LVIS)**0.1
VAL4 = DG* (DL-DG) *GC
Y = VAL3/VAL4
CALCULATE THE PACKED TOWER FLOODING CONDITIONS.
CALL FLOOD (L, G, DG, DL, PF, LVIS)
RETURN
END
***************************************************************
THIS PROGRAM CALCULATES THE VAPOR RATE AT FLOODING CONDITION
***************************************************************
SUBROUTINE FLOOD (L, G, DG, DL, PF, LVIS)
REAL L, LVI S
COMMON/TOWER1/ GC, PI
COMMON/TOWER2/
X, Y, Xl,
Y1
COMMON/TOWER3/ AREA, DIA
COMMON/TOWER4/ G1, G2, VG
COMMON/TOWER5/ AREAF, DIAF
X1 = (L/G)*SQRT(DG/DL)
TI = 62.4/DL
Y1=i0.0.* (-i. 6678-( i. 085*ALOGI0 (Xl)
7-( 0. 29655. (ALOGI0 (Xl) )*.2) )
VAL1 = YI*DL*DG*GC
VAL2 = PF*TI*(LVIS)**2
RATIO = VAL1/VAL2
G2 = SQRT(RATIO)
CALCULATE VAPOR RATE AT 70% OF FLOODING CONDITION
G2 = 0.7"G2
CALCULATE VAPOR VELOCITY AT 70% FLOODING CONDITION
VG=G2/DG
CALCULATE THE TOWER AREA AT 70% FLOODING CONDITION.
SINCE THE FLOODING CORRELATION IS NOT PERFECT, HENCE
TO HAVE A 95% CONFIDENCE, USE A SAFETY FACTOR OF 1.32
FOR THE CALCULATED CROSS-SECTIONAL AREA (Bolles and Fair,
Chem. Eng., 89 (14), 109 (July 12, 19827.
AREAF = 1.32*(G/(3600.0*G2) )
CALCULATE THE TOWER DIAMETER AT 70% FLOODING CONDITION
DIAF = SQRT((4.0*AREAF)/PI)
573
120
130
140
150
160
170
180
190
200
210
220
230
240
DIAMETER IN inches:
DIAF = 12.0*DIAF
RETURN
END
****************************************************************
THIS PROGRAM OUTPUTS THE RESULTS ON THE SCREEN
****************************************************************
SUBROUTINE OUTPUT (L, G, DG, DL, PF, DELP, LVIS)
REAL L,LVIS
COMMON~TOWER2~
X, Y, XI, Y1
COMMON/TOWER3/ AREA, DIA
COMMON/TOWER4/ GI, G2, VG
COMMON/TOWERS/ AREAF, DIAF
WRITE (*, 120)
FORMAT(20X,'PACKED TOWER DESIGN')
WRITE (*, 130)
FORMAT(78(IH*))
WRITE (*, 140) L
FORMAT(5X,'LIQUID FLOW RATE, ib/hr.:', T60, F12.3)
WRITE (*, 150) G
FORMAT(5X,'GAS FLOW RATE, ib/hr.:', T60, F12.3)
WRITE (*, 160) DG
FORMAT(5~,'GAS DENSITY, ib/ft^3: ' , T60, F8.4)
WRITE (*, 170) DL
FORMAT(5X,'LIQUID DENSITY, ib/ft^3: ' , T60, F8.4)
WRITE (*, 180) PF
FORMAT(5X,'TOWER PACKING FACTOR:', T60, FS.I)
WRITE (*, 190) DELP
FORMAT(5X,'PRESSURE DROP, inch H20/ft:', T60, F8.3)
WRITE (*, 200) LVIS
FORMAT(5X,'LIQUID VISCOSITY, cP:', T60, F8.3)
WRITE (*, 210) X
FORMAT(5X,'X CO-ORDINATE IN THE GENERALIZED PLOT:', T60, F8.4)
WRITE (*, 220) Y
FORMAT(5X,'Y CO-ORDINATE IN THE GENERALIZED PLOT:', T60, F8.4)
WRITE (*, 230) G1
FORMAT(5X,'CALCULATED VAPOR RATE, Ib/s.ft^2: ' , T60, F8.4)
WRITE (*, 240) X1
FORMAT(5X,'X CO-ORDINATE AT FLOODING CONDITION:', T60, F8.4)
WRITE (*, 250) Y1
574
250
260
270
280
290
300
310
320
330
340
350
360
370
380
FORMAT(SX,'Y CO-ORDINATE AT FIX)ODING CONDITION:', T60, F8.4)
WRITE (*, 260) G2
FORMAT(5X,'VAPOR RATE AT 70% FLOODING CONDITION, ib/s.ft^2: ' ,\,
* T60, F8.4)
WRITE (*, 270) VG
FORMAT(5X, 'VAPOR VELOCITY AT 70% FLOODING CONDITION, ft./s. : ' ,\,
* T60, F8.2)
WRITE (*, 280) AREA
FORMAT(SX, 'TOWER AREA, ft^2: ' T60, F8.2)
WRITE (*, 290) DIA
FORMAT(5X, 'TOWER DIAMETER, inches: ', T60, F8.2)
WRITE (*, 300) AREAF
FORMAT(5X, 'TOWER AREA AT 70% FLOODING CONDITION, ft^2: ' ,T60,F8.2)
WRITE (*, 310) DIAF
FORMAT(5X, 'TOWER DIAMETER AT 70% FLOODING CONDITION, inches: ' ,\,
* T60, F8.2)
WRITE (*, 320)
FORMAT (78(IH-))
RETURN
END
******************************************************************
THIS PROGRAM PRINTS THE RESULTS ONTO A PRINTER
******************************************************************
SUBROUTINE PAPER (L, G, DG, DL, PF, DELP, LVIS)
REAL L, LVIS
COMMON/TOWER2/ X, Y, XI, Y1
COMMON/TOWER3/ AREA, DIA
COMMON/TOWER4/ GI, G2, VG
COMMON/TOWER5/ AREAF, DIAF
WRITE (i, 330)
FORMAT (///,25X,'PACKED TOWER DESIGN')
WRITE (I, 340)
FORMAT (78(IH*))
WRITE (I, 350) L
FORMAT (5X,'LIQUID FLOW RATE, ib/hr.:', T60, F12.3)
WRITE (i, 360) G
FORMAT (5X,'GAS FLOW RATE, ib/hr.:', T60, F12.3)
WRITE (i, 370) DG
FORMAT (5X,'GAS DENSITY, ib/ft.^3: ' , T60, F8.4)
WRITE (I, 380) DL
FORMAT (5X,'LIQUID DENSITY, Ib/ft.^3: ' , T60, F8.4)
WRITE (i, 390) PF
575
390 FORMAT (5X,'TOWER PACKING FACTOR:', T60, F8.1)
400
410
420
430
440
450
460
470
480
490
500
510
520
530
WRITE (i, 400) DELP
FORMAT (5X,'PRESSURE DROP, INCH-H20:', T60, F8.3)
WRITE (i, 410) LVIS
FORMAT (5X,'LIQUID VISCOSITY, cP:', T60, F8.3)
WRITE (i, 420) X
FORMAT (5X,'X CO-ORDINATE IN THE GENERALIZED PLOT:', T60, F8.4)
WRITE (I, 430) Y
FORMAT (5X,'Y CO-ORDINATE IN THE GENERALIZED PLOT:', T60, F8.4)
WRITE (i, 440) G1
FORMAT (5X,'CALCULATED VAPOR RATE, Ib/s.ft^2: ' , T60, F8.4)
WRITE
(
FORMAT
WRITE (
FORMAT
WRITE (
FORMAT
I, 450) X1
(5X,'X CO-ORDINATE AT FLOODING CONDITION:', T60, F8.4)
i, 460) Y1
(5X,'Y CO-ORDINATE AT FLOODING CONDITION:', T60, F8.4)
i, 470) G2
(5X,'VAPOR RATE AT 70% FLOODING CONDITION, Ib/s.ft^2:',\,
T60, F8.4)
WRITE (I, 480) VG
FORMAT (5X,'VAPOR VELOCITY AT 70% FLOODING CONDITION, ft./sec:',\,
* T60, F8.2)
WRITE (i, 490) AREA
FORMAT (5X,'TOWER AREA, ft^2: ' , T60, F8.2)
WRITE (i, 500) DIA
FORMAT (5X,'TOWER DIAMETER, inches:', T60, F8.2)
WRITE (i, 510) AREAF
FORMAT (5X,'TOWER AREA AT 70% FLOODING CONDITION, ft^2:',\,
* T60,F8.2)
WRITE (i, 520) DIAF
FORMAT (5X,'TOWER DIAMETER AT 70% FLOODING CONDITION, inches:',\,
* T60, F8.2)
WRITE (i, 530)
FORMAT (78(IH-))
FORM FEED THE PRINTING PAPER TO THE TOP OF THE NEXT PAGE.
WRITE (I, *) CHAR(12)
RETURN
END
576