Coker A.K. Fortran Programs for Chemical Process Design, Analysis, and Simulation
Подождите немного. Документ загружается.
15
PROGRAM PROG61
******************************************************************
THIS PROGRAM CALCULATES EITHER POLYTROPIC OR ADIABATIC HEADS,
THE ASSOCIATED GAS HORSEPOWER AND GAS DISCHARGE TEMPERATURE,
AS LONG AS THE POLYTROPIC OR ADIABATIC EFFICIENCY IS USED
WITH THE COMPANION HEAD CALCULATION. THE PROGRAM ALSO CALCULATES
ADIABATIC EFFICIENCY FROM POLYTROPIC EXPONENT "K" FROM THE
MOLAL HEAT CAPACITY (MCp).
Hpoly(PH)
ZAV
R
TS
MW
N
RC
TII
TI2
K
Epoly (PE)
HP
W
Had (AH)
T21
T22
Ead (AE)
MCp (MCP)
P1
P2
CHOICE
= POLYTROPIC HEAD, ft
= AVERAGE COMPRESSIBILITY FACTOR FOR GAS FROM SUCTION TO
DISCHARGE CONDITIONS. A VALUE OF 1.0 WILL GIVE
CONSERVATIVE RESULTS.
= GAS CONSTANT= 1544.
= SUCTION TEMPERATURE, oF
= MOLECULAR WEIGTH OF GAS
= POLYTROPIC EXPONENT
= COMPRESSION RATIO (DISCHARGE PRESSURE/SUCTION PRESSURE)
P2/Pl
= INLET TEMPERATURE, oF
= DISCHARGE TEMPERATURE, oF
= ADIABATIC EXPONENT, Cp/Cv
= POLYTROPIC EFFICIENCY (TRY 0.75 FOR PRELIMINARY WORK)
= GAS HORSEPOWER, bhp.
= GAS FLOW, ib/h
= ADIABATIC HEAD, ft
= SUCTION TEMPERATURE, oF
= DISCHARGE TEMPERATURE, oF
= ADIABATIC EFFICIENCY
= MOLAL SPECIFIC HEAT OF GAS AT CONSTANT PRESSURE,
Btu/(ib mol)(oF).
= SUCTION PRESSURE, psia.
= DISCHARGE PRESSURE, psia.
= CHOOSE WHETHER THE PROCESS IS ADIABATIC OR POLYTROPIC
**********************************************************************
CHARACTER*I4 CHOICE
REAL K,MW,N
COMMON/PROGI/W,PI,P2,TII,TI2,T21,T22,K,MW,N,ZAV,AE,AECAL,RC
COMMON/PROG2/POLYH,ADIAH,HPPOLY,HPADIA,PE,PECAL,BHPP,BHPA
COMMON/PROG3/WORK,PROD
OPEN (UNIT=3,FILE='DATA61.DAT',STATUS='OLD',ERR=I8)
OPEN (UNIT =i, FILE ='PRN')
DETERMINE WHETHER THE COMPRESSOR CALCULATION IS POLYTROPIC
OR ADIABATIC.
READ ( 3, 15, ERR=I9) CHOICE
FORMAT (A)
IF (CHOICE .EQ. 'POLYTROPIC' .OR. CHOICE .EQ. 'polytropic') THEN
457
18
I0
19
20
40
30
60
70
80
90
i00
READ (3, *
READ (3, *
READ (3, *
, ERR=I9) W, PI, P2
, ERR=I9) TII, K, MW
, ERR=I9) ZAV, PE
GO TO 40
ELSEIF (CHOICE .EQ. 'ADIABATIC' .OR. CHOICE .EQ. 'adiabatic') THEN
READ (3, *
READ (3, *
READ (3, *
, ERR=I9) W, PI, P2
, ERR=I9) TII, K, MW
, ERR=I9) ZAV, AE
GO TO 50
ENDIF
WRITE (*, i0)
FORMAT (3X, 'FILE DOES NOT EXIST')
GO TO 999
WRITE (*, 20)
FORMAT (3X, 'ERROR MESSAGE IN THE DATA VALUE')
GO TO 999
WRITE (i, 30)
FORMAT
(///,26X,'POLYTROPIC
COMPRESSOR CALCULATION', /, 78(IH*))
CALL POLYTR
WRITE (i, 60) CHOICE, W, PI, P2, RC
FORMAT (5X,'TYPE OF SERVICE:', T60, A,/,
5X,'GAS FLOWRATE, ib/hr.:',T60, F12.3,/,
5X,'SUCTION PRESSURE, psia:', T60, F12.3,/,
5X,'DISCHARGE PRESSURE, psia:', T60, F12.3,/,
5X,'COMPRESSION RATIO:', T60, F12.3)
WRITE (i, 70) TII, K, MW
FORMAT (5X,'INLET TEMPERATURE, oF:',T60, F12.3,/,
5X,'RATIO OF SPECIFIC HEAT CAPACITIES, k:',T60,FI2.3,/,
5X,'MOLECULAR WEIGHT OF GAS, ib/ib-mole:',T60,Fl2.3)
WRITE (i, 80) ZAV, PE
FORMAT (5X,'AVERAGE COMPRESSIBILITY FACTOR, Z:',T60, F8.4,/,
5X,'POLYTROPIC EFFICIENCY:',T60,F8.2)
WRITE (I, 90) POLYH, TI2, HPPOLY
FORMAT (5X,'POLYTROPIC HEAD, ft:',T60,Fl2.3,/,
5X,'GAS DISCHARGE TEMPERATURE, oF:',T60, F12.3, /,
5X,'GAS HORSEPOWER, hp:', T60, F12.3)
WRITE (i, I00) PROD, BHPP, WORK
FORMAT (5X,'MECHANICAL LOSSES, hp:', T60, F12.3,/,
5X,'TOTAL BRAKE HORSEPOWER, hp:',T60, F12.3,/,
5X,'WORK DONE, Ib-mole/hr.:', T60, F12.3)
WRITE (i, 105)
458
105
C
50
ii0
120
130
140
150
160
170
180
190
C
999
FORMAT (78(IH-))
FORM FEED THE PRINTING PAPER TO THE TOP OF THE NEXT PAGE.
WRITE (i, *) CHAR(12)
GO TO 999
WRITE (i, ii0)
FORMAT (///,26X,'ADIABATIC COMPRESSOR CALCULATION', /, 78(IH*))
CALL ADIAB
WRITE (i, 120) CHOICE, W, PI, P2, RC
FORMAT (5X,'TYPE OF SERVICE:', T60, A,/,
* 5X,'GAS FLOWRATE, ib/hr.:', T60, F12.3,/,
* 5X,'SUCTION PRESSURE, psia:', T60, F12.3,/,
* 5X,'DISCHARGE PRESSURE, psia:', T60, F12.3,/,
* 5X,'COMPRESSION RATIO:', T60, F12.3)
WRITE (i, 130) TII, K, MW
FORMAT (5X,'INLET TEMPERATURE, oF:',T60, F12.3,/,
* 5X,'RATIO OF SPECIFIC HEATS CAPACITIES, K:',T60,FI2.3,/,
* 5X,'MOLECULAR WEIGHT OF GAS, ib/ib-mole:',T60, F12.3)
WRITE (i, 140) ZAV, AE
FORMAT (SX,'AVERAGE COMPRESSIBILITY FACTOR, Z:',T60, F8.4,/,
* 5X,'ADIABATIC EFFICIENCY:',T60, F12.3)
WRITE (I, 150) ADIAH, T22, HPADIA
FORMAT (SX,'ADIABATIC HEAD, ft.:',T60, F12.3,/,
*
5X,'GAS DISCHARGE TEMPERATURE, oF:', T60, F12.3,/,
* 5X, 'GAS HORSEPOWER, hp:', T60, F12.3)
WRITE (i, 160) PROD
FORMAT (5X,'MECHANICAL LOSSES, hp:',T60, F12.3)
WRITE (i, 170) BHPA, WORK
FORMAT (5X,'TOTAL BRAKE HORSEPOWER, hp:',T60, F12.3,/,
*
5X,'WORK DONE, ft/ib-mole:',T60,Fl2.3)
WRITE (I, 180) AECAL
FORMAT (5X,'CALCULATED ADIABATIC EFFICIENCY FROM',/,
* 5X,'POLYTROPIC EFFICIENCY:',T60, F12.3)
WRITE (i, 190)
FORMAT (IH ,78(IH-))
FORM FEED THE PRINTING PAPER TO TOP OF THE NEXT PAGE.
WRITE (i, *) CHAR(12)
CLOSE (3,STATUS='KEEP')
CLOSE (i)
STOP
END
459
THIS PROGRAM CALCULATES THE POLYTROPIC HEAD,THE ASSOCIATED GAS
HORSEPOWER AND GAS DISCHARGE TEMPERATURE IN POLYTROPIC MACHINES
SUBROUTINE POLYTR
REAL K,MW,N,MLOSS
COMMON/PROGI/W,PI,P2,TII,TI2,T21,T22,K,MW,N,ZAV,AE,AECAL,RC
COMMON/PROG2/POLYH,ADIAH,HPPOLY,HPADIA,PE,PECAL,BHPP,BHPA
COMMON/PROG3/WORK, PROD
R = 1545.0
TII = TII + 460.
CALCULATE THE COMPRESSION RATIO (DISCHARGE PRESS./SUCTION PRESS.)
RC = P2/Pl
N = (K*PE)/(K*PE-K+I.0)
CALCULATE THE POLYTROPIC HEAD, POLYH: ft.
VALI = (ZAV*R*TII)/MW
VAL2 = (N/(N-I.0))*(RC**((N-I.0)/N)-I.0)
POLYH = VALI*VAL2
CALCULATE THE GAS HORSEPOWER, HPPOLY: hp.
HPPOLY = (W*POLYH)/(I.98*(10.0**6)*PE)
CALCULATE THE APPROXIMATE MECHANICAL LOSSES AS A PERCENTAGE OF GAS
POWER REQUIREMENT, PROD: hp.
CALL MECHL (HPPOLY, MLOSS)
CALCULATE THE TOTAL BRAKE HORSEPOWER, BHPP: hp.
BHPP = HPPOLY + PROD
CALCULATE THE GAS DISCHARGE TEMPERATURE, TI2:.oF.
VAL3 = (POLYH*MW)/(ZAV*R)
TI2 = VAL3*((N-I.0)/N) + TII
CALCULATE THE WORK DONE ib- mole/hr.
CALCULATE THE BASIS IN moles.
B =
W/MW
CALCULATE THE WORK DONE
WORK = (N/(N-I.O))*(R*TII*((P2/PI)**((N-I.0)/N)-I.0))
WORK = (WORK*B)/(1.98*(10.0**6))
TII = TII - 460.0
TI2 = TI2 - 460.0
RETURN
END
460
THIS PROGRAM CALCULATES THE ADIABATIC HEAD,THE ASSOCIATED GAS
HORSEPOWER AND GAS DISCHARGE TEMPERATURE
*****************************************************************
SUBROUTINE ADIAB
REAL K,MW,N,MLOSS
COMMON/PROGI/W,PI,P2,TII,TI2,T21,T22,K,MW,N,ZAV,AE,AECAL,RC
COMMON/PROG2/POLYH,ADIAH,HPPOLY,HPADIA,PE,PECAL,BHPP,BHPA
COMMON/PROG3/WORK, PROD
PE = 0.75
R = 1545.0
T21 = TII + 460.
CALCULATE THE COMPRESSION RATIO (DISCHARGE PRESS./SUCTION PRESS.)
RC = P2/Pl
CALCULATE THE ADIABATIC HEAD, ADIAH: ft.
VAL4 = (ZAV*R*T21)/MW
VAL5 = (K/(K-I.0))*(RC**((K-I.0)/K)-I.0)
ADIAH = VAL4*VAL5
CALCULATE THE GAS HORSEPOWER, HPADIA: hp.
HPADIA = (W*ADIAH)/(I.98*(10.0**6)*AE)
CALCULATE THE APPROXIMATE MECHANICAL LOSSES AS A PERCENTAGE OF GAS
POWER REQUIREMENT, PROD: hp.
CALL MECHL (HPADIA, MLOSS)
CALCULATE THE TOTAL BRAKE HORSEPOWER, BHPA: hp.
BHPA = HPADIA + PROD
CALCULATE THE GAS DISCHARGE TEMPERATURE, oF.
T22 = T21*(RC**((K-I.0)/K))
T22 = T22 - 460.0
CALCULATE THE ADIABATIC EFFICIENCY FROM THE POLYTROPIC EFFICIENCY
N = (K*PE)/(K*PE-K+I.0)
VAL6 = (RC**((K-I.0)/K))-I.0
VAL7 = (RC**((N-I.0)/N))-I.0
AECAL = VAL6/VAL7
CALCULATE THE WORK DONE, hp
CALCULATE THE BASIS IN ib-moles per hr.
B =
W/MW
CALCULATE THE WORK DONE
461
WORK = (K/(K-I.0))*(R*T21)*((P2/PI)**((K-I.0)/K)-I.O)
WORK = (WORK*B)/(1.98*(10.0**6))
RETURN
END
THIS PROGRAM CALCULATES THE MECHANICAL LOSSES DUE TO
FRICTION IN BEARINGS.
SUBROUTINE MECHL (POLYL, MLOSS)
COMMON/PROG3/WORK,PROD
REAL MLOSS
IF (POLYL .LT. 3000.0 ) THEN
MLOSS = 0.03
PROD = POLYL*MLOSS
ELSEIF (POLYL .GT. 3000.0 .AND. POLYL .LT. 6000.0) THEN
MLOSS = 0.025
PROD = POLYL*MLOSS
ELSEIF (POLYL .GT. 6000.0 .AND. POLYL .LT. i0000.0) THEN
MLOSS = 0.02
PROD = POLYL*MLOSS
ELSEIF (POLYL .GT. i0000.0) THEN
MLOSS = 0.015
PROD = POLYL*MLOSS
ENDIF
RETURN
END
462
PROGRAM PROG62
********************************************************************
THIS PROGRAM CALCULATES ALL PUMP HEAD CALCULATIONS REQUIRED TO WRITE
A PUMP SPECIFICATION SUITABLE FOR A VENDOR INQUIRY OR PURCHASE
ORDER. THE PROGRAM FURTHER CALCULATES THE CENTRIFUGAL PUMP
EFFICIENCY AND PUMP HORSEPOWER USING AN EQUATION DEVELOPED BY
BY C. BRANAN (The Process Engineer's Pocket Handbook, Gulf Publishing
Company, Houston 1976.)
********************************************************************
FLUID = FLUID NAME
TEMP = FLUID OPERATING TEMPERATURE, oF.
VIS = FLUID VISCOSITY, cP.
GPM = FLOWRATE, gpm
SPG = SPECIFIC GRAVITY OF LIQUID
VP = VAPOR PRESSURE OF PUMPED FLUID AT PUMP SUCTION
TEMPERATURE, psia.
VESP = PRESSURE AT EQUIPMENT ON PUMP SUCTION, psia.
SLL = SUCTION LINE LOSSES (PRESSURE DROP IN PIPING), psi.
SSHD = SUCTION STATIC HEAD OF LIQUID ON PUMP SUCTION OR
DISCHARGE, ft.
NPSH = NET POSITIVE SUCTION HEAD ON PUMP, ft.
DELP = PRESSURE AT EQUIPMENT ON PUMP DISCHARGE, psia.
DSHD = DISCHARGE STATIC HEAD, ft.
CVPD = CONTROL VALVE PRESSURE DROP, psi.
EXPD = HEAT EXCHANGER PRESSURE DROP, psi.
ORPD = ORIFICE PRESURE DROP, psi.
FURPD = FURNACE PRESSURE DROP, psi.
DLL = DISCHARGE LINE LOSSES (PRESSURE DROP IN PIPING), psi.
OTLD = OTHER LOSSES (DISCHARGE:-e.g. FURNACE), psi.
OTLS = OTHER LOSSES (SUCTION:- e.g. STRAINER/FILTER), psi
DELTAP = PRESSURE DROP ACROSS PUMP, psi.
DELPHT = PRESSURE DROP ACROSS PUMP, ft.
HBHP = HYDRAULIC BRAKE HORSEPOWER REQUIRED BY PUMP, bhp.
(FOR HBHP CALCULATION, ASSUME PUMP EFFICIENCY=I.0)
EFF = PUMP EFFICIENCY (%).
F = DEVELOPED PUMP HEAD, ft.
G = FLOW RATE THROUGH PUMP, gpm.
BHP = ACTUAL BRAKE HORSEPOWER REQUIRED FOR PUMP OPERATION
psi = ft(0.433 x sp.gr.)
ft = psi/(0.433 x sp.gr.)
CHARACTER*f4 FLUID
REAL NPSH, ORPD, OTLD, OTLS
COMMON/DATA1/ FLUID, TEMP, VIS, SPG, VP
COMMON/DATA2/ GPM, VESP, SLL, SSHD, SSHDP
COMMON/DATA3/ DELP, DSHD, DSHDP, CVPD, EXPD
COMMON/DATA4/ ORPD, FURPD, DLL, OTLS, OTLD
COMMON/DATA5/ SP, SPGAGE, NPSH, DISP, DISPG
COMMON/DATA6/ DELTAP,DELPHT,HBHP,BHP,EFF
OPEN(UNIT=3,FILE='DATA62.DAT',STATUS='OLD',ERR=I8)
OPEN (UNIT =1, FILE='PRN')
463
18
ii
19
12
15
i0
20
30
40
50
60
70
READ (3, 5, ERR=I9) FLUID
FORMAT (A)
READ (3, *, ERR=I9) TEMP, VIS, SPG
READ(3, *, ERR=I9) VP, GPM
READ(3, *, ERR=I9) VESP, SLL, SSHD
READ(3, *, ERR=I9) DELP, DSHD,CVPD
READ(3, *, ERR=I9) EXPD, ORPD, FURPD
READ(3, *, ERR=I9) DLL, OTLS, OTLD
GO TO 15
WRITE(*, Ii)
FORMAT(6X,'FILE DOES NOT EXIST')
GO TO 999
WRITE(*, 12)
FORMAT(6X,'ERROR MESSAGE IN THE DATA VALUE')
GO TO 999
CALL PUMPHA
WRITE(l, i0)
FORMAT(26X,'CENTRIFUGAL PUMP SIZING',/,IH ,78(IH*))
WRITE (i, 20) FLUID, TEMP, VIS, SPG, VP, GPM
FORMAT (5X,'FLUID NAME:', T65, A,/,
*
5X,'FLUID OPERATING TEMPERATURE, oF:', T65, F9.3,/,
* 5X,'FLUID VISCOSITY, cP:', T65, F9.4,/,
* 5X,'FLUID SPECIFIC GRAVITY:', T65, F9.3,/,
* 5X,'FLUID VAPOR PRESSURE, psia:', T65, F9.3,/,
* 5X,'VOLUMETRIC FLOWRATE:, gpm:', T65, F9.3)
WRITE (i, 30) VESP, SLL, OTLS, SSHD
FORMAT (5X,'PRESSURE AT EQUIPMENT, psia:', T65, F9.3,/,
* 5X,'SUCTION LINE LOSSES, psi:', T65, F9.3,/,
* 5X,'OTHER LOSSES (e.g. STRAINER/FILTER), psi:',T65,F9.3,
* /,5X,'SUCTION STATIC HEAD, ft.:', T65, F9.3)
WRITE (i, 40) SSHDP, SPGAGE, SP
FORMAT (5X,'SUCTION STATIC HEAD, psi:', T65, F9.3,/,
5X,'TOTAL PUMP SUCTION PRESSURE, psig:', T65, F9.3,/,
5X,'TOTAL PUMP SUCTION PRESSURE, psia:', T65, F9.3)
WRITE
(i, 50)
NPSH
FORMAT (5X, 'AVAILABLE NET POSITIVE SUCTION HEAD, ft. : ' , T65, F9.3)
WRITE (i, 60) DELP, DSHD, DSHDP, CVPD
FORMAT (5X,'DELIVERY PRESSURE, psia.:', T65, F9.3,/,
*
5X,'DISCHARGE STATIC HEAD, ft.:', T65, F9.3,/,
* 5X,'DISCHARGE STATIC HEAD, psi:', T65, F9.3,/,
* 5X,'CONTROL VALVE PRESSURE DROP, psi.:', T65, F9.3)
WRITE (i, 70) EXPD, ORPD, FURPD, DLL, OTLD
FORM_AT (5X,'EXCHANGER PRESSURE DROP, psi.:', T65, F9.3,/,
5X,'ORIFICE PRESSURE DROP, psi.:', T65, F9.3,/,
5X,'FURNACE PRESSURE DROP, psi.:', T65, F9.3,/,
464
80
90
i00
ii0
120
130
999
5X,'DISCHARGE LINE LOSSES, psi.:', T65, F9.3,/,
5X,'OTHER LOSSES, psi.:', T65, F9.3)
WRITE(l, 80)DISPG,DISP
FORMAT(5X,'TOTAL PUMP DISCHARGE PRESSURE, psig:', T65, F9.3,/,
5X,'TOTAL PUMP DISCHARGE PRESSURE, psia:', T65, F9.3)
WRITE(I, 90)DELTAP,DELPHT
FORMAT(5X,'TOTAL PUMP PRESSURE DROP, psi.:', T65, F9.3,/,
5X,'TOTAL PUMP PRESSURE DROP, ft.:', T65, F9.3)
WRITE(I, 100)HBHP
FORMAT(5X,'HYDRAULIC BRAKE HORSEPOWER REQUIRED BY THE PUMP,'\,
9 bhp:', T65, F9.3)
WRITE(I, II0)BHP
FORMAT(5X,'ACTUAL BRAKE HORSEPOWER REQUIRED FOR PUMP OPERATION',
\,' bhp:', T65, F9.3)
WRITE (i, 120) EFF
FORMAT (5X, 'PUMP EFFICIENCY, (%) : ', T65, F9.3)
WRITE (i, 130)
FORMAT (IH ,78(IH-) )
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 SUCTION PRESSURE AND THE NET POSITIVE
SUCTION HEAD
********************************************************************
SUBROUTINE PUMPHA
CHARACTER*I4 FLUID
REAL NPSH, OTLS, OTLD, ORPD
COMMON~DATA1~
FLUID, TEMP, VIS, SPG, VP
COMMON/DATA2/ GPM, VESP, SLL, SSHD, SSHDP
COMMON/DATA3/ DELP, DSHD, DSHDP, CVPD, EXPD
COMMON/DATA4/ ORPD, FURPD, DLL, OTLS, OTLD
COMMON/DATAS/
SP, SPGAGE, NPSH, DISP, DISPG
COMMON/DATA6/ DELTAP, DELPHT, HBHP, BHP, EFF
VALI = 0.4329*SPG
SSHDP = SSHD*VALI
SP = VESP + SSHDP- SLL- OTLS
465
SPGAGE = SP-14.7
AVAILABLE NET POSITIVE SUCTION HEAD
SLLI = SLL/VALI
VESPI = VESP/VALI
VPI = VP/VALI
DEF = VESPI-VPI
NPSH = DEF + SSHD - SLLI
PUMP DISCHARGE PRESSURE
DSHDP = DSHD*VALI
DISP = DELP+DSHDP+CVPD+EXPD+ORPD+FURPD+DLL+OTLD
DISPG = DISP-14.7
TOTAL PUMP DELTAP = PUMP DISCH. PRES. - PUMP SUC. PRES.
DELTAP = DISP-SP
DELPHT = DELTAP/VALI
HYDRAULIC BRAKE HORSEPOWER REQUIRED BY PUMP
PUMP EFFICIENCY IS ASSUMED = 1.0
HBHP = (DELTAP*GPM)/(1715.0*I.0)
THIS PROGRAM CALCULATES EFFICIENCY AND PUMP HORSEPOWER USING
C. BRANAN'S EQUATION (e.g. C Brannan: The Process Engineer's
Pocket Handbook, Gulf publishing Co.)
F = DELPHT
G = GPM
VAL2=0.2855*F
VAL3=3.78*F*G*(10.0**(-4))
VAL4=2.38*F*(G**2)*(10.0**(-7))
VAL5=5.39*(F**2)*(IO.0**(-4))
VAL6=6.39*(F**2)*G*(10.0**(-7))
VAL7=4.0*(F**2)*(G**2)*(10.O**(-10))
EFF=(80.0-VAL2+VAL3-VAL4+VAL5-VAL6+VAL7)
BHP=(100.0*HBHP)/EFF
RETURN
END
466