Coker A.K. Fortran Programs for Chemical Process Design, Analysis, and Simulation
Подождите немного. Документ загружается.
18
12
i0
13
44
I00
ii0
120
130
140
150
55
160
READ(3,*,ERR=I0) P31, P32, F
GO TO 66
WRITE(I,12)
FORMAT(3X,'FILE DOES NOT EXIST')
WRITE(I,13)
GO TO 999
FORMAT(3X,'ERROR MESSAGE IN THE DATA VALUE')
GO TO 999
WRITE (I, i00)
FORMAT (//,20X,'CONTROL VALVE SIZING FOR LIQUID SERVICE',/,
IH
,76(1H*))
CALCULATE THE CONTROL VALVE FOR LIQUID SERVICE
CALL LIQUID
WRITE (i, ii0) CHOICE
FORMAT (5X,'FLUID SERVICE:', T60, A)
IF (CVL .GT. 0.0) THEN
WRITE (I, 1
FORMAT (5X,
5X,
5X,
5X,
5X,
5X,
20) V, SPL, PI, P2, DELPL, CVL
'LIQUID FLOWRATE, gpm:', T60, F9.2,/,
'LIQUID SPECIFIC GRAVITY:', T60, F9.3,/,
'INLET PRESSURE, psia:', T60, F9.3,/,
'OUTLET PRESSURE, psia:', T60, F9.3,/,
'CONTROL VALVE PRESSURE DROP, psi:', T60, F9.3,/,
'CONTROL VALVE CAPACITY COEFFICIENT:', T60, F9.3)
WRITE (I, 130)
FORMAT ( 76 ( IH- ) )
ELSE
WRITE (i, 140) V, SPL, PI, P2, DELPL, CVL
FORMAT (SX,'LIQUID FLOWRATE, gpm:', T60, F9.2,/,
* 5X,'LIQUID SPECIFIC GRAVITY:', T60, F9.3,/,
*
5X,'INLET PRESSURE, psia:', T60, F9.3,/,
* 5X,'OUTLET PRESSURE, psia:', T60, F9.3,/,
*
5X,'CONTROL VALVE PRESSURE DROP, psi:', T60, F9.3,/,
* 5X,'CONTROL VALVE CAPACITY COEFFICIENT:', T60, F9.3)
WRITE ( i, 150)
FORMAT ( 76 ( IH- ) )
ENDIF
GO TO 5
WRITE (I, 160)
FORMAT
(//,20X,'CONTROL
VALVE SIZING FOR GAS SERVICE',/,76(IH*))
CALCULATE THE CONTROL VALVE FOR GAS SERVICE.
397
170
180
190
200
66
210
220
230
240
CALL GAS
WRITE (I, 170) CHOICE
FORMAT (5X,'FLUID SERVICE:', T60, A)
IF (CVG .GT. 0.0) THEN
WRITE (i, 180) Q, SPG, P21, P22, DELPG, T, CVG
FORMAT (5X,'GAS FLOW RATE, ft^3/hr.: ' , T60, F9.3,/,
* 5X,'SPECIFIC GRAVITY:', T60, F9.3,/,
* 5X,'INLET PRESSURE, psia:', T60, F9.3,/,
*
5X,'OULET PRESSURE, psia:', T60, F9.3,/,
* 5X,'CONTROL VALVE PRESSURE DROP, psi:', T60, F9.3,/,
*
5X,'INLET TEMPERATURE, oF:', T60, F9.3,/,
* 5X,'CONTROL VALVE CAPACITY COEFFICIENT:', T60, F9.3)
ELSE
WRITE (i, 190) Q, SPG, P21, P22, DELPG, T, CVG
FORMAT (5X,'GAS FLOW RATE, ft^3/hr.: ' , T60, F9.3,/,
* 5X,'SPECIFIC GRAVITY:', T60, F9.3,/,
* 5X,'INLET PRESSURE, psia:', T60, F9.3,/,
*
5X,'OUTLET PRESSURE, psia:', T60, F9.3,/,
* 5X,'CONTROL VALVE PRESSURE DROP, psi:', T60, F9.3,/,
*
5X,'INLET TEMPERATURE, oF:', T60, F9.3,/,
* 5X,'CONTROL VALVE CAPACITY COEFFICIENT:', T60, F9.3)
ENDIF
WRITE (I, 200)
FORMAT ( 76 ( IH- ) )
GO TO 5
WRITE (i, 210)
FORMAT (//,20X, 'CONTROL VALVE SIZING FOR STEAM SERVICE',
/, 76( IH* ) )
CALCULATE THE CONTROL VALVE FOR STEAM SERVICE.
CALL STEAM
WRITE (i, 220) CHOICE
FORMAT (SX,'FLUID SERVICE:', T60, A)
IF (CVS .GT. 0.0 ) THEN
WRITE (I, 230) W, F, P31, P32, DELPS, CVS
FORMAT (5X,'CALCULATED STEAM FLOW RATE, ib/hr.:', T60, F9.3,/,
*
5X,'DEGREE OF SUPERHEAT, oF:', T60, F9.3,/,
*
5X,'INLET PRESSURE, psia:', T60, F9.3,/,
* 5X,'OUTLET PRESSURE, psia:', T60, F9.3,/,
* 5X,'CONTROL VALVE PRESSURE DROP, psi:', T60, F9.3,/,
* 5X,'CONTROL VALVE CAPACITY COEFFICIENT:', T60, F9.3)
WRITE (I, 240)
FORMAT (76(IH-))
398
250
260
999
FORM FEED THE PRINTING PAPER TO THE TOP OF THE NEXT PAGE.
WRITE (i, *) CHAR(12)
ELSE
WRITE (I, 250) W, F, P31, P32, DELPS, CVS
FORMAT (5X,'STEAM FLOW RATE, ib/hr.:', T60, F9.3,/,
5X,'DEGREE OF SUPERHEAT, oF:', T60, F9.3,/,
5X,'INLET PRESSURE, psia:', T60, F9.3,/,
5X,'OUTLET PRESSURE, psia:', T60, F9.3,/,
5X,'CONTROL VALVE PRESSURE DROP, psi:', T60, F9.3,/,
5X,'CONTROL VALVE CAPACITY COEFFICIENT:', T60, F9.3)
WRITE (i, 260)
FORMAT (76(IH-))
ENDIF
FORM FEED THE PRINTING PAPER TO TOP OF THE NEXT PAGE.
WRITE (i, *)CHAR(12)
CLOSE(3,STATUS='KEEP')
CLOSE (I)
STOP
END
***************************************************************
THIS PROGRAM CALCULATES THE CONTROL VALVE CAPACITY COEFFICIENT
OR THE MAXIMUM FLOWRATE FOR FULLY OPEN CONTROL VALVE AND LIQUID
FLOW CONDITIONS.
***************************************************************
SUBROUTINE LIQUID
COMMON/DATA1/ V, PI, P2, DELPL, SPL, CVL
CHECK WHETHER COMPUTATION REQUIRES CAPACITY COEFFICIENT OR
MAXIMUM LIQUID FLOW RATE.
CALCULATE THE PRESSURE DROP ACROSS THE CONTROL VALVE, psi.
DELPL = P1 - P2
IF (CVL .GT. 0.0) THEN
GO TO i0
ELSE
CALCULATE THE VALVE CAPACITY COEFFICIENT (Cv)
CVL = V*SQRT(SPL)/SQRT(DELPL)
ENDIF
GO TO 30
399
C
i0
30
15
25
i0
CALCULATE THE MAXIMUM LIQUID FLOW RATE, gpm.
V = CVL*SQRT(DELPL)/SQRT(SPL)
RETURN
END
THIS PROGRAM CALCULATES THE CONTROL VALVE CAPACITY COEFFICIENT OR
THE MAXIMUM GAS FLOW RATE FOR GAS FLOW CONDITIONS
SUBROUTINE GAS
COMMON~DATA2~
Q, P21, P22, DELPG, SPG, T, CVG
T = T+460.0
CALCULATE THE CONTROL VALVE PRESSURE DROP, psi.
DELPG = P21 - P22
VALG = P22*DELPG
CHECK WHETHER COMPUTATION REQUIRES CAPACITY COEFFICIENT OR
MAXIMUM GAS FLOW RATE.
IF (CVG .GT. 0.0) THEN
GO TO i0
ELSE
ENDIF
CHECK WHETHER FLOW IS CRITICAL OR SUBCRITICAL
IF ((P21/2.) .GE. P22) THEN
GO TO 15
ELSE
CALCULATE THE VALVE CAPACITY COEFFICIENT FOR SUBCRITICAL FLOW
CVG = Q*SQRT(SPG*T)/(1360.0*SQRT(VALG))
ENDIF
GO TO 2O
CALCULATE THE VALVE CAPACITY COEFFICIENT FOR CRITICAL FLOW.
WRITE (*, 25)
FORMAT (3X,'CRITICAL GAS FLOW')
pause
CVG = Q*SQRT(SPG*T)/(1360.0*O.5*P21)
GO TO 2O
CALCULATE THE MAXIMUM GAS FLOW RATE, ftA3/hr. FOR SUBCRITICAL
FLOW.
IF ((P21/2.0) .GE. P22) THEN
400
C
35
45
20
GO TO 35
ELSE
Q = (CVG*I360.0*SQRT(VALG))/SQRT(SPG*T)
ENDIF
GO TO 2O
CALCULATE THE MAXIMUM GAS FLOW RATE, ftA/hr. FOR CRITICAL FLOW
WRITE (*, 45)
FORMAT(3X,'CRITICAL GAS FLOW')
pause
Q = (CVG*I360.0*0.5*P21)/SQRT(SPG*T)
T = T-
460.0
RETURN
END
*****************************************************************
THIS PROGRAM CALCULATES THE CONTROL VALVE CAPACITY COEFFICIENT
OR THE MAXIMUM STEAM FLOW RATE FOR FULLY OPEN CONTROL VALVE AND
FLOW CONDITIONS.
******************************************************************
SUBROUTINE STEAM
REAL K
COMMON/DATA3/ W, P31, P32, DELPS, F, CVS
K = 1.0+.0007*F
CALCULATE THE CONTROL VALVE PRESSURE DROP, psi.
DELPS = P31 - P32
VALS = P32*DELPS
CHECK WHETHER COMPUTATION REQUIRES CAPACITY COEFFICIENT OR
MAXIMUM FLOW RATE.
IF (CVS .GT. 0. O) THEN
GO TO I0
ELSE
ENDIF
CHECK WHETHER FLOW IS CRITICAL OR SUBCRITICAL
IF ((P31/2.) .GE. P32) THEN
GO TO 15
ELSE
ENDIF
CALCULATE THE CAPACITY COEFFICIENT FOR SUBCRITICAL FLOW
CVS = (W*K)/3.*SQRT(VALS)
, 401
15
25
i0
35
45
20
GO TO 20
WRITE(*,25)
FORMAT(3X,'CRITICAL STEAM FLOW')
pause
CALCULATE THE CAPACITY COEFFICIENT FOR CRITICAL FLOW.
CVS = (W*K)/(I.5*P31)
GO TO 20
CHECK WHETHER FLOW IS CRITICAL OR SUBCRITICAL FLOW AND
CALCULATE THE FLOW.
IF ((P31/2.) .GE. P32) THEN
GO TO 35
ELSE
CALCULATE THE MAXIMUM FLOW RATE FOR CRITICAL FLOW.
W = 3.0*CVS*SQRT(VALS)/K
ENDIF
GO TO 20
WRITE(*, 45)
FORMAT(3X,'CRITICAL STEAM FLOW')
pause
CALCULATE THE MAXIMUM FLOW RATE FOR CRITICAL FLOW.
W=(I.5*P31*CVS)/K
RETURN
END
402
PROGRAM PROG53
RELIEF-VALVE SIZING FOR VAPOR, STEAM, LIQUID, AIR AND FIRE.
THE PROGRAM USES THE EQUATIONS AND CORRELATIONS PRESENTED IN
~4ERICAN PETROLEUM INSTITUTE (API) 520 RECO~DIENDED PRACTICES
FOR THE DESIGN AND INSTALLATION OF PRESSURE-RELIEVING SYSTEMS
IN REFERINERIES FOR THE SIZING OF RELIEF VALVES IN VAPOR,
STEAM, LIQUID, AIR AND FIRE SERVICES.
CHARACTER COUNT(14),REGIME(14),REG(14),REGI(14)
CHARACTER*I4 COUNTI(14),ACOUNT(14),COUNLL(14),ACONT(14)
CHARACTER*I0 CHOICE
REAL OPI, MWI, N1
REAL OP2, KSH, MW2, N2
REAL KV, KVI, KW, KD
REAL OP4
REAL LAT, OP5, MW5, N5
DIMENSION Y(14)
COMMON/NAME/I,Y,IJK
COMMON/DATA/REGI,ACONT
COMMON/DATAI/WGMAX,ADESI,ACALI,CI
COMMON/DATA2/WSMAX,ADES2,ACAL2,C2
COMMON/DATA3/GPIMAX,GP2MAX,ADES3,ACAL3,ACAL31,ADES31,RE
COMMON/DATA4/COUNT,REGIME,COUNT!,COUNLL,REG, ACOUNT
COMMON/DATA5/FORCE
COMMON/DATA6/ADES4,VAMAX,ACAL4,C4
COMMON/DATA7/Q,WF
COMMON/DATA8/KVI
COMMON CHOICE
DATA COUNT/'D' 'E' 'F' 'G 9 'H' 'J' 'K' 'L' 'M' 'N'
s 9 9 n 9 n 9
,
s n
'P','Q','R','T'/
DATA COUNT1/' 1-2 ' 9 ' 1-2 ',
* 9 1.5-2 " , ' 1.5-2 ' ,
* ' 1.5-3 or 2-3 ',' 2-3 or 2.5-4',
* ' 3-4 '9 3-4 or 4-6 ',
* ' 4-6 9 9 ' 4-6 ' ,
* ' 4-6 ' ' 6-8 '
>
~: ' 6-8 or 6-10 ' 9 ' 8-i0 '
OPEN (UNIT=39 FILE='DATA53.DAT'9 STATUS='OLD'9 ERR=I8)
OPEN (UNIT =i9 FILE ='PRN')
READ(3, 59 ERR=I9) CHOICE
FORMAT (A)
IF (CHOICE .EQ. 'VAPOR' .OR. CHOICE .EQ. 'vapor') THEN
GO TO i0
ELSEIF (CHOICE .EQ. 'STEAM' .OR. CHOICE .EQ. 'steam') THEN
GO TO 20
ELSEIF (CHOICE .EQ. 'LIQUID 9 .OR. CHOICE .EQ. 9 9 ) THEN
GO TO 30
~LSEIF (CHOICE .EQ. 9 .OR. CHOICE .EQ. 'air 9 ) THEN
GO TO 4O
ELSEIF (CHOICE .EQ. 'FIRE' .OR. CHOICE .EQ. 'fire') THEN
403
i0
20
GO TO 50
ENDIF
CALCULATE RELIEF VALVE FOR VAPOR CONDITION
INPUT DATA FOR VAPOR CONDITION
WG = VAPOR FLOW RATE, Ib/hr.
SPI = RELIEF VALVE SET PRESSURE, psiq.
OPI = ALLOWABLE OVERPRESSURE, (%).
T1 = RELIEVING TEMPERATURE OF THE INLET VAPOR (oF).
Z1 = COMPRESSIBILITY FACTOR.
MWI : MOLECULAR WEIGHT OF VAPOR.
N1 = RATIO OF THE SPECIFIC HEATS.
READ ( 3, *, ERR:I9 ) WG, SPI
READ (3, *, ERR=Ig) OPI, T1
READ (3, *, ERR:f9) ZI, ~Z<I
READ (3, *, ERR=f9) N1
DETE~4INE THE RELIEF VALVE ORIFICE AREA FOR VAPOR CONDITION
CALL VAPOR (WG, SPI, OPI, TI, ZI, }~i, NI]
PRINT RESULTS
CALL OUT1 (WG, SPI, OPI, TI, ZI, ~Z~I, NI)
GO TO 60
INPUT DATA FOR STEAM RELIEF CONDITION.
WS = STEAM FLOW RATE, !b/hr.
SP2 = RELIEF VALVE SET PRESSURE, psig.
OP2 = ALLOWABLE OVERPRESSURE (%).
KSH = SUPERHEAT STEAM CORRECTION FACTOR.
T2 = RELIEF TEMPERATURE, oF.
HW2 = HOLECUL~R WEIGHT OF STEAM.
N2 = RATIO OF THE SPECIFIC HEATS.
READ (3, *, ERR=f9) WS, SP2
READ (3, *, ERR:Ig) OP2, KSH
READ (3, *, EB_R=Ig) T2, MW2
READ (3, *, ERR=I9) N2
DETERMINE THE RELIEF VALVE ORIFICE AREA FOR STEAM CONDITION
CALL STEAM (WS, SP2, OP2, KSH, T2, M%~2, N2)
PRINT RESULTS
CALL OUT2 (WS, SP2, OP2, KSH, T2, MW2, N2)
GO TO 60
INPUT DATA FOR LIQUID RELIEF
GPM = LIQUID FLOW RATE, U.S. gpm.
SP3 = UPSTREAM RELIEVING SET PRESSURE, psig.
PB = TOTAL BACK PRESSU?.E, psig.
404
30
40
SPL =
VISL =
KV =
KW =
KD =
READ (3, *
READ (3, *
READ (3, *
READ (3, *
SPECIFIC GRAVITY OF THE LIQUID.
VISCOSITY OF LIQUID, cP.
CORRECTION FACTOR DUE TO VISCOSITY.
CORRECTION FACTOR DUE TO BACK PRESSURE.
EFFECTIVE COEFFICIENT OF DISCHARGE.
, ERR =19) GPM, SP3
, ERR = 19) PB, SPL
, ERR=I9) VISL, KV
, ERR=I9) KW, KD
DETERMINE THE RELIEF VALVE ORIFICE AREA FOR LIQUID CONDITION
CALL LIQUID (GPM, SP3, PB, SPL, VISL, KV, KW, KD)
PRINT RESULTS
IF (VISL .LT. i0.0) THEN
CALL OUT3 (GPM, SP3, PB, SPL, VISL, KV, KW, KD)
GO TO 60
ELSE
CALL OUT32 (VISL)
GO TO 60
ENDIF
INPUT DATA FOR AIR RELIEF
VA = REQUIRED AIR CAPACITY, ft^3/min.
T = INLET TEMPERATURE, oF.
P = RELIEVING PRESSURE, psia.
SP4 = RELIEF VALVE SET PRESSURE, psig.
READ (3, *, ERR=I9) VA, T
READ (3, *, ERR=I9) P, SP4
DETERMINE THE RELIEF VALVE ORIFICE AREA FOR AIR CONDITION
CALL AIR (VA, T, P, SP4)
PRINT RESULTS
CALL OUT4 (VA, T, P, SP4)
GO TO 60
DETERMINE THE RELIEF VALVE ORIFICE AREA FOR FIRE CONDITION.
F
A
LAT
SP5
OP5
T
Z
MW5
= ENVIRONMENTAL FACTOR.
= TOTAL WETTED SURFACE AREA, ft^2.
= LATENT HEAT OF VAPORIZATION, Btu/ib.
= RELIEF VALVE SET PRESSURE, psig.
= ALLOWABLE OVERPRESSURE, (%).
= RELIEVING TEMPERATURE, oF.
= COMPRESSIBILITY FACTOR.
= FLUID MOLECULAR WEIGHT.
405
C
C
50
18
750
19
800
60
N5 = RATIO OF THE SPECIFIC HEATS.
INPUT DATA FOR FIRE RELIEF CONDITION.
READ (3, *, ERR=Ig) F, A
READ (3, *, ERR=f9) LAT, SP5
READ (3, *, ERR=f9) OPL, T
READ (3, *, ERR=I9) Z, MW5
READ (3, *, ERR=Ig) N5
DETERMINE THE RELIEF VALVE ORIFICE AREA FOR FIRE CONDITION
CALL FIRE (F, A, LAT, SPL, OPL, T, Z, MWL, NL)
PRINT RESULTS
CALL OUT5 (F, A, LAT, SPL, OPL, T, Z, MWL, NL)
GO TO 60
WRITE (*, 750)
FORMAT (3X, 'DATA FILE DOES NOT EXIST')
GO TO 60
WRITE (*, 800)
FORMAT (3X, 'ERROR MESSAGE IN THE DATA VALUE')
GO TO 60
CLOSE (3, STATUS = 'KEEP')
CLOSE (i)
STOP
END
THIS PROGRAM SIZES THE RELIEF VALVE FOR VAPOR CONDITIONS
SUBROUTINE VAPOR (WG, SPI, OPI, TI, El, MWI, NI)
CHARACTER COUNT(14),REGIME(14), REG(14)
CHARACTER*I4 COUNT1(14), COUNLL(14), ACOUNT(14)
REAL K,KB,MWI,NI,OPI
COMMON/NAME/I,Y,IJK
COMMON/DATAI/WGMAX,ADESI,ACALI,CI
COMMON/DATA4/COUNT,REGIME,COUNTI,COUNLL,REG, ACOUNT
COMMON/DATAL/FORCE
DIMENSION Y(14)
CALCULATE THE GAS COEFFICIENT (CI) AND THE ORIFICE AREA.
K=0.975
KB=I.0
TI=TI+460.0
PI=SPI+(SPI*OPI/100.0)+I4.7
VALI=NI*(2.0/(NI+I.O))**((NI+I.0)/(NI-I.0))
CI=520.0*SQRT(VALI)
VAL2=CI*K*PI*KB*SQRT(MWI)
406