Coker A.K. Fortran Programs for Chemical Process Design, Analysis, and Simulation
Подождите немного. Документ загружается.
33
C
C
44
C
C
55
C
C
66
C
GO TO 44
ELSEIF (TS .EQ. 24.0) THEN
GO TO 55
ELSEIF (TS .EQ. 30.0) THEN
GO TO 66
ELSEIF (TS .EQ. 36.0) THEN
GO TO 77
ENDIF
CALL PROGI
PRINT THE RESULTS ON THE SCREEN
CALL OUTPUT
PRINT THE RESULTS ONTO A PRINTER.
CALL PAPER
GO TO 999
CALL PROG2
PRINT THE RESULTS ON THE SCREEN
CALL OUTPUT
PRINT THE RESULTS ONTO A PRINTER
CALL PAPER
GO TO 999
CALL PROG3
PRINT THE RESULTS ON THE SCREEN
CALL OUTPUT
PRINT THE RESULTS ONTO A PRINTER
CALL PAPER
GO TO 999
CALL PROG4
PRINT THE RESULTS ON THE SCREEN
CALL OUTPUT
557
77
C
999
PRINT THE RESULTS ONTO A PRINTER.
CALL PAPER
GO TO 999
CALL PROG5
PRINT THE RESULTS ON THE SCREEN.
CALL OUTPUT
PRINT THE RESULTS ONTO A PRINTER.
CALL PAPER
CLOSE (3, STATUS='KEEP')
CLOSE (i)
STOP
END
**************************************************************
TOWER SIZING FOR VALVE TRAYS USING A 12 INCH TRAY SPACING.
**************************************************************
SUBROUTINE PROG1
REAL LCFS
COMMON PI
COMMON/MULTI/ X, SF, RD, V, VRD
COMMON/MULT2/ AB, ADC, AT, DIA, RDCA
COMMON/MULT3/ WV, WL, DENL, DENV, FF, REST, TS, SURT
COMMON/MULT4/ LCFS, VCFS
G = 0.77174
H = -0.02964
CALCULATE THE TRAY SPACING FACTOR, TSF
TSF = G+H*ALOG(X)
CALCULATE THE VAPOR VOLUMETRIC RATE, VCFS
VCFS = WV/(DENV*3600.0)
CALCULATE THE LIQUID VOLUMETRIC RATE, LCFS
LCFS = WL/(DENL*3600.0)
IF (VRD .GT. 0.5) THEN
VRD=0.5
ELSE
ENDIF
CALCULATE THE BUBBLING AREA, AB
AB = (VCFS*RD+LCFS)/(TSF*VRD*FF)
558
ADC = (REST*LCFS*I2.0)/(FF*TS)
CALCULATE THE TOWER AREA, AT
AT = (2.0*ADC+AB)/SF
CALCULATE THE TOWER DIAMETER, DIA.
DIA = SQRT((4.0*AT)/PI)
CALCULATE THE DOWNCOMER AREA RATIO, RDCA
RDCA = ADC/(2.*ADC + AB)
RETURN
END
TOWER SIZING FOR VALVE TRAYS USING AN 18 INCH TRAY SPACING.
SUBROUTINE PROG2
REAL LCFS
COMMON PI
COMMON/MULTI/ X, SF, RD, V, VRD
COMMON/MULT2/ AB, ADC, AT, DIA, RDCA
COMMON/MULT3/ WV, WL, DENL, DENV, FF, REST, TS, SURT
COMMON/MULT4/ LCFS, VCFS
G = 0.93655
H = -0.02310
CALCULATE THE TRAY SPACING FACTOR, TSF
TSF = G+H*ALOG(X)
CALCULATE THE VAPOR VOLUMETRIC RATE, VCFS
VCFS = WV/(DENV*3600.0)
CALCITE THE LIQUID VOLUMETRIC RATE, LCFS
LCFS = WL/(DENL*3600.0)
IF ( VRD .GT. 0.5) THEN
VRD = 0.5
ELSE
ENDIF
CALCULATE THE BUBBLING AREA, AB.
AB = (VCFS*RD+LCFS)/(TSF*VRD*FF)
CALCULATE THE DOWNCOMER AREA, ADC
ADC = (REST*LCFS*I2.0)/(FF*TS)
559
ADC = (REST*LCFS*I2.0)/(FF*TS)
CALCULATE THE TOWER AREA, AT.
AT = ( 2. O*ADC+AB)/SF
CALCULATE THE TOWER DIAMETER, DIA.
DIA = SQRT((4.0*AT)/PI)
CALCULATE THE DOWNCOMER AREA RATIO, RDCA
RDCA = ADC/(2.0*ADC+AB)
RETURN
END
***************************************************************
TOWER SIZING FOR VALVE TRAYS USING A 24 INCH TRAY SPACING.
***************************************************************
SUBROUTINE PROG3
REAL LCFS
COMMON PI
COMMON/MULTI/ X, SF, RD, V, VRD
COMMON/MULT2/ AB, ADC, AT, DIA, RDCA
COMMON/MULT3/ WV, WL, DENL, DENV, FF, REST, TS, SURT
COMMON/MULT4/ LCFS, VCFS
CALCULATE THE TRAY SPACING FACTOR, TSF
FOR ALL VALUES OF X, TSF=I. 0
TSF = 1.0
CALCULATE THE VAPOR VOLUMETRIC RATE, VCFS
VCFS = WV/(DENV*3600.0)
CALCULATE THE LIQUID VOLUMETRIC RATE, LCFS
LCFS = WL/(DENL*3600.0)
IF ( VRD .GT. 0.5) THEN
VRD=0.5
ELSE
ENDIF
CALCULATE THE BUBBLING AREA, AB
AB = (VCFS*RD+LCFS)/(TSF*VRD*FF)
CALCULATE THE DOWNCOMER AREA, ADC
ADC = (REST*LCFS*I2.0)/(FF*TS)
CALCULATE THE TOWER AREA, AT.
560
AT = (2.0*ADC+AB)/SF
CALCULATE THE TOWER DIAMETER, DIA.
DIA
=
SQRT((4.0*AT)/PI)
CALCULATE THE DOWNCOMER AREA RATIO, RDCA
RDCA = ADC/((2.0*ADC)+AB)
RETURN
END
***************************************************************
TOWER SIZING FOR VALVE TRAYS USING A 30 INCH TRAY SPACING.
***************************************************************
SUBROUTINE PROG4
REAL LCFS
COMMON PI
COMMON/MULTI/ X, SF, RD, V, VRD
COMMON/MULT2/ AB, ADC, AT, DIA, RDCA
COMMON/MULT3/ WV, WL, DENL, DENV, FF, REST, TS, SURT
COMMON/MULT4/ LCFS, VCFS
G = 0.98067
H = 0.03220
CALCULATE THE TRAY SPACING FACTOR, TSF
TSF = G + H*ALOG(X)
IF ( X .LT. 1.85 ) THEN
TSF = 1.0
ELSE
ENDIF
CALCULATE THE VAPOR VOLUMETRIC RATE, VCFS
VCFS = WV/(DENV*3600.0)
CALCULATE THE LIQUID VOLUMETRIC RATE, LCFS
LCFS = WL/(DENL*3600.O)
IF ( VRD .GT. 0.5) THEN
VRD = 0.5
ELSE
ENDIF
CALCULATE THE BUBBLING AREA, AB
AB = (VCFS*RD+LCFS)/(TSF*VRD*FF)
CALCULATE THE DOWNCOMER AREA, ADC.
ADC =(REST*LCFS*I2.0)/(FF*TS)
561
CALCULATE THE TOWER AREA, AT
AT = (2.0*ADC+ AB)/SF
CALCULATE THE TOWER DIAMETER, DIA.
DIA = SQRT((4.0*AT)/PI)
CALCULATE THE DOWNCOMERAREA RATIO, RDCA
RDCA = ADC/((2.0*ADC)+AB)
RETURN
END
***************************************************************
TOWER SIZING FOR VALVE TRAYS USING A 36 INCH TRAY SPACING.
***************************************************************
SUBROUTINE PROG5
REAL LCFS
COMMON PI
COMMON/MULTI/ X, SF, RD, V, VRD
COMMON/MULT2/ AB, ADC, AT, DIA, RDCA
COMMON/MULT3/ WV, WL, DENL, DENV, FF, REST, TS, SURT
COMMON/MULT4/ LCFS, VCFS
G = 0.96583
H = 0.06162
CALCULATE THE TRAY SPACING FACTOR, TSF.
TSF = G+H*ALOG(X)
IF (X .LT. 1.85) THEN
TSF = 1.0
ELSE
ENDIF
CALCULATE THE VAPOR VOLUMETRIC RATE, VCFS
VCFS = WV/(DENV*3600.0)
CALCULATE THE LIQUID VOLUMETRIC RATE, LCFS.
LCFS = WL/(DENL*3600.O)
IF (VRD .GT. 0.5) THEN
VRD = 0.5
ELSE
ENDIF
CALCULATE THE BUBBLING AREA, AB.
AB = (VCFS*RD+LCFS)/(TSF*VRD*FF)
CALCULATE THE DOWNCOMERAREA, ADC
562
i00
ii0
120
130
140
150
160
ADC = (REST*LCFS*I2.0)/(FF*TS)
CALCULATE THE TOWER AREA, AT.
AT = (2.0*ADC+AB)/SF
CALCULATE THE TOWER DIAMETER, DIA.
DIA = SQRT((4.0*AT)/PI)
CALCULATE THE DOWNCOMER AREA RATIO, RDCA
RDCA = ADC/((2.0*ADC)+AB)
RETURN
END
**************************************************************
THIS PROGRAM PRINTS THE RESULTS ONTO THE SCREEN
SUBROUTINE OUTPUT
REAL LCFS
COMMON/MULTI/ X, SF, RD, V, VRD
COMMON/MULT2/ AB, ADC, AT, DIA, RDCA
COMMON/MULT3/ WV, WL, DENL, DENV, FF, REST, TS, SURT
COMMON/MULT4/
LCFS, VCFS
WRITE (*, i00)
FORMAT (//, 25X, 'TOWER DESIGN FOR VALVE TRAYS',/,78(IH*))
WRITE (*, Ii0) WV, VCFS, WL, LCFS, DENL
FORMAT (5X,'VAPOR FLOW RATE, ib/hr.:', T60, FI0.1,/,
5X,'VAPOR VOLUMETRIC RATE, ft^3/sec.: ' , T60, FI0.3,/,
5X,'LIQUID FLOW RATE, ib/hr.:', T60, FI0.1,/,
5X,'LIQUID VOLUMETRIC RATE, ft^3/sec.: ' , T60, FI0.3,/,
5X,'LIQUID DENSITY, Ib/ft^3.: ~ , T60, FI0.3)
WRITE (*, 120) DENV, FF, REST
FORMAT (SX,'VAPOR DENSITY, ib/ft^3.: ' , T60, FIO.3,/,
5X,'FOAM FACTOR:', T60, FI0.3,/,
5X,'DOWNCOMER RESIDENCE TIME, sec.:', T60, FI0.3)
WRITE (*, 130) TS, SURT
FORMAT (5X,'TRAY SPACING, in.:', T60, F8.0,/,
5X,'SURFACE TENSION, dyne/cm.:', T60, F8.2)
WRITE (*, 140) X, SF, RD
FORMAT (5X,'SURFACE TENSION/VAPOR DENSITY (X):', T60, FI0.3,/,
5X,'SAFETY FACTOR:', T60, FI0.3,/,
5X,'DENSITY RADICAL:', T60, FI0.3)
WRITE (*, 150) AB, ADC, AT
FORMAT (5X,'BUBBLING AREA, ft^2: ' , T60, FI0.3,/,
5X,'DOWNCOMER AREA, ft^2: ' , T60, FI0.3,/,
5X,'TOWER AREA, ft^2: ' , T60, FIO.3)
WRITE (*, 160) DIA, RDCA
FORMAT (5X,'TOWER DIAMETER, ft:', T60, FI0.3,/,
563
170
180
190
200
210
220
230
240
250
5X,'DOWNCOMER AREA RATIO:', T60, FI0.4)
WRITE (*, 170)
FORMAT (78(IH-))
RETURN
END
***************************************************************
THIS PROGRAM PRINTS THE RESULTS ON THE PRINTER.
***************************************************************
SUBROUTINE PAPER
REAL LCFS
COMMON /MULTI/ X, SF, RD, V, VRD
COMMON/ MULT2/ AB, ADC, AT, DIA, RDCA
COMMON/ MULT3/ WV, WL, DENL, DENV, FF, REST, TS, SURT
COMMON/ MULT4/ LCFS, VCFS
WRITE (i, 180)
FORMAT (///, 25X, 'TOWER DESIGN FOR VALVE TRAYS',/,78(IH*))
WRITE (i, 190) WV, VCFS, WL, LCFS, DENL
FORMAT (5X,'VAPDR FLOW RATE, ib/hr.:', T60, FI0.1,/,
* 5X,'VAPOR VOLUMETRIC RATE, ft^3/sec.: ' , T60, FI0.3,/,
* 5X,'LIQUID FLOW RATE, ib/hr.:', T60, FI0.1,/,
* 5X,'LIQUID VOLUMETRIC RATE,
ftA3/sec.:',T60,FlO.3,/,
* 5X,'LIQUID DENSITY, Ib/ft^3.: ' , T60, FI0.3)
WRITE (i, 200) DENV, FF, REST
FORMAT (5X,'VAPOR DENSITY, ib/ft^3.: ' , T60, FI0.3,/,
* 5X,'FOAM FACTOR:', T60, FI0.3,/,
*
5X,'DOWNCOMER RESIDENCE TIME, sec.:', T60, FI0.3)
WRITE (i, 210) TS, SURT
FORMAT (5X,'TRAY SPACING, in.:', T60, F8.0,/,
* 5X,'SURFACE TENSION, dyne/cm.:', T60, F8.2)
WRITE (i, 220) X, SF, RD
FORMAT (5X,'SURFACE TENSION/VAPOR DENSITY (X):', T60, FIO.3,/,
* 5X,'SAFETY FACTOR:', T60, FI0.3,/,
* 5X,'DENSITY RADICAL:', T60, FI0.3)
WRITE (i, 230) AB, ADC, AT
FORMAT (5X,'BUBBLING AREA, ft^2: ' , T60, FI0.3,/,
* 5X,'DOWNCOMER AREA, ft^2: ' , T60, FI0.3,/,
* 5X,'TOWER AREA, ft^2: ' , T60, FI0.3)
WRITE (i, 240) DIA, RDCA
FORMAT (5X,'TOWER DIAMETER, ft:', T60, FI0.3,/,
* 5X,'DOWNCOMER AREA RATIO:', T60, FI0.4)
WRITE (i, 250)
FORMAT (78(IH-))
FORM FEED THE PRINTING PAPER TO TOP OF THE NEXT PAGE.
WRITE ( I, *) CHAR(12)
RETURN
END
564
PROGRAM PROG74
PACKED-TOWER SIZING
THE TOWER SIZING METHOD IS BASED ON THE GENERALIZED PRESSURE
DROP CORRELATION OF SHERWOOD, SHIPLEY AND HOLLOWAY AND
REPLOTTED IN FIGURE GR-109 R-5 BY NORTON CO. THIS CORRELATION
IS A PLOT OF X .VS. Y FOR PRESSURE DROP CURVES RANGING FROM
0.5 TO 1.5 INCH OF WATER. FROM THESE CURVES THE TOWER DIAMETER
CAN BE EVALUATED.
***************************************************************
L
G
DG
DL
X
Y
G1
PF
LVIS
GC
AREA
DIA
DELP
= LIQUID RATE, Ib/h
= VAPOR RATE, Ib/h
= VAPOR DENSITY, ib/ft^3
= LIQUID DENSITY, ib/ft^3
= X CO-ORDINATE IN GR-109 R-5 PLOT
= Y CO-ORDINATE IN GR-109 R-5 PLOT
= VAPOR RATE, ib/s.ft^2
= PACKING FACTOR FOR PACKING USED
= LIQUID VISCOSITY, cP
= GRAVITATIONAL CONSTANT = 32.2 ft/s^2
= TOWER AREA, ft^2
= TOWER DIAMTER, ft
= PRESSURE DROP inch H20/ft.
PI = CONSTANT = 3.1415927
***************************************************************
COMMON/TOWER1/ GC, PI
COMMON/TOWER2/ X, Y, XI, Yl
COMMON/TOWER3/ AREA, DIA
COMMON/TOWER4/ GI, G2, VG
COMMON/T~)WERS/ AREAF, DIAF
REAL L, LVIS
OPEN (UNIT = 3, FILE='DATA74.DAT'
OPEN (UNIT = I, FILE='PRN')
, STATUS='OLD', ERR=IS)
18
i00
19
READ LIQUID FLOW RATE, Ib/hr.
READ GAS FLOW RATE, ib/hr.
READ GAS DENSITY, ib/ft^3.
READ LIQUID DENSITY, ib/ft^3.
READ TOWER PACKING FACTOR
READ PRESSURE DROP, inch-H20/ft.
READ LIQUID VISCOSITY, cP.
READ (3, *, ERR =19) L, G, DG
READ (3, *, ERR =19) DL, PF, DELP, LVIS
GO TO I0
WRITE (*, i00)
FORMAT (5X, 'DATA FILE DOES NOT EXIST')
GO TO 999
WRITE (*, ii0)
= L
= G
= DG
= DL
= PF
= DELP
= LVIS
565
ii0
i0
30
C
35
C
40
FORMAT (5X, 'ERROR MESSAGE IN THE DATA VALUE')
GO TO 999
PI = 3.1415927
GC = 32.2
WHERE 0.01 < X < 2.0
SUMI=DG/( DL-DG )
X = L/G*SQRT(SUMI)
IF (DELP .EQ. 0.05) THEN
GO TO 30
ELSEIF(DELP .EQ. 0.i0) THEN
GO TO 35
ELSEIF(DELP .EQ. 0.25) THEN
GO TO 40
ELSEIF(DELP .EQ. 0.50) THEN
GO TO 45
ELSEIF(DELP .EQ. 1.0) THEN
GO TO 5O
ELSEIF(DELP .EQ. 1.5) THEN
GO TO 55
ENDIF
CALL PROGI (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 PROG2 (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 PROG3 (L, G, DG, DL, PF, LVIS)
PRINT THE RESULTS ON THE SCREEN
566