Coker A.K. Fortran Programs for Chemical Process Design, Analysis, and Simulation
Подождите немного. Документ загружается.
13
18
20
19
21
C
15
READ (3, 13, ERR=I9)FSS
FORMAT (A)
READ THE
READ THE
READ THE
READ THE
READ THE
READ THE
READ THE
READ THE
READ THE
READ THE
READ THE
READ THE
NUMBER OF SHELL: NS
SHELL SIDE INSIDE DIAMETER, inch: SID
NUMBER OF BAFFLES: NB
BAFFLE SPACING, inch: BSP
SHELL SIDE FOULING FACTOR, SFF
INLET TEMPERATURE, oF: TSI
OULET TEMPERATURE, oF: TS2
SHELL SIDE SPECIFIC HEAT CAPACITY, Btu/ib.oF.: SSHC
FLUID VISCOSITY, cP: SVIS
FLUID DENSITY, ib/ft^3.: SDEN
FLUID THERMAL CONDUCTIVITY, Btu/hr.ft.oF.: SFK
FLUID FLOWRATE, ib/hr.: SW.
READ (3, *
READ (3, *
READ (3, *
, ERRs19) NS, SID, NB, BSP
, ERR=I9) SFF, TSI, TS2, SSHC
, ERR=I9) SVIS, SDEN, SFK, SW
GO TO 15
WRITE (*, 20)
FORMAT (3X, 'ERROR MESSAGE IN THE DATA VALUE')
GO TO 999
WRITE (*, 21)
FORMAT (3X, 'DATA VALUE DOES NOT EXIST')
GO TO 999
TUBE SIDE CALCULATION
CALL TUBE
SHELL SIDE CALCULATION
CALL SHELL
CALCULATE THE HEAT LOAD
Q = TW*TSHC*ABS(TT2-TTI)
CALCULATE THE LOG MEAN TEMPERATURE DIFFERENCE (LMTD)
VALI = (TSI-TT2)-(TS2-TTI)
VAL2 = ALOG((TSI-TT2)/(TS2-TTI))
LMTD = ABS (VALIMAL2)
CALCULATE THE APPROACH FACTOR
R = (TSI-TS2)/(TT2-TTI)
P = (TT2-TTI) / (TSI-TTI)
IF (NS .EQ. 1.0 .AND. R .EQ. 1.0) THEN
PX = P/(NS-NS*P+P)
687
55
C
i00
C
C
SUM1 = (PX*(R*R+I.0)**0.5)/(I.0-PX)
SUM2 = 2.0/PX-I.0-R+(R*R+I.0)**0.5
SUM3 = 2.0/PX-I.0-R-(R*R+I.0)**0.5
SUM4 = ALOG(SUM2/SUH3)
F = SUMI/SUM4
GO TO i00
ELSEIF (NS .EQ. 1.0 .AND. R .GT. 1.0) THEN
SUM5 = ((R*P-I.0)/(P-I.0))**(I.0/NS)
PX
= (I.0-SUMS)/(R-SU~5)
SUM6 = (R*R+I.0)**0.5/(R-I.0)
SUM7 = ALOG((I.0-PX)/(I.0-R*PX))
SUM41 = ((2.0/PX)-I.0-R+(R*R+I.0)**0.5)
SUM42 = ((2.0/PX)-I.0-R-(R*R+I.0)**0.5)
IF (SUM41 .LT. 0.0 .OR. SUM42 .LT. 0.0) THEN
NS = NS+I.
TNP = TNP+2.
GO TO 55
ELSE
SUM43 = ALOG(SUM41/SUH42)
ENDIF
F = (SUM6*SUM7)/SUM43
IF (F .LT. 0.8) THEN
NS = NS+I.0
TNP = TNP+2.
GO TO 55
ELSE
END IF
GO TO I00
ELSEIF (NS .GE. 2.0 .AND. NT .GE. 4.0) THEN
SUM8 = (R*R+I.0)**0.5/(2.0*(R-I.0))
SUM9 = ALOG((I.0-P)/(I.0-P*R))
SUM10 = 2/P-I.0-R+((2/P)*((I.0-P)*(I-P*R))**0.5)+(R*R+I.)**0.5
SUM11 = 2/P-I.0-R+((2/P)*((I.0-P)*(I.-P*R))**0.5)-(R*R+I.)**0.5
F = (SUMS*SUM9)/ALOG(SUMI0/SUMII)
ENDIF
CORRECTED MEAN T~PERATURE DIFFERENCE
CMTD = F*LMTD
CALCULATE THE OVERALL HEAT TRANSFER COEFFICIENT
Btu/hr.ft^2.oF
U = 1.0/(I./HT+I./HS+TFF+SFF)
CALCULATE THE OUTSIDE AREA OF THE HEAT EXCHANGER UNIT ft^2.
AREA = (TOD/12.0*PI*NT*NS*TL)
688
999
Q1 = AREA*U*CMTD
CALL RESULT
FORM FEED THE PRINTING PAPER TO TOP OF THE NEXT PAGE.
WRITE (I, *) CHAR(12)
CLOSE ([]NIT =3, STATUS='KEEP')
CLOSE (UNIT = i)
STOP
END
689
*THIS PROGR2~4 CALCULATES T}{E TUBE SIDE HEAT TRANSFER COEFFICIENT, *
*REYNOLDS NUMBER, PRANDTL NUMBER, FRICTION FACTOR AND PRESSURE DROP.*
SUBROUTINE TUBE
PARA!4ETER (PI : 3.1415927, G:32.174)
CHARACTER*30 FTS, FSS
CHARACTER TYPE
COMMON/TUBE1/ FTS, TYPE
COMMON/TUBE2/ NT, TID, TOD, TP, TL, TNP
COMMON/TUBE3/ TFF, TTI, TT2, TSHC, TVIS, TDEN
CO~4ON/TUBE4/ TFK, TW
COHMON/TUBES/ RET, PRT, TFRF, HT, TV
COM!.'.ON/TUBE6/
TDELP, TRDELP, TOTDP
COMMON/SHELL1/ FSS
COMHON/SHELL2/
NS, SID, NB, BSP, SFF
CO.~[MON/SHELL3/ TSl, TS2, SSHC, SVIS, SDEN
COM/.[ON/SHELL4/ SFK, SW
COMMON/SHELLS/ RES, PRS, SFRF, HS, SV
CO~4MON/SHELL6/
SDELP
COMMON~OUT~
Q, AREA, F, I2,~TD, U, QI, CHTD
REAL NT, NS, NB, LMTD
CALCULATE THE FLOW AREA: TFA, ft ^2
TFA = (PI*TID**2*NT)/(4.0*I44.0*TNP)
CALCULATE THE I.~SS VELOCITY, ib/hr.ft~2
GT : TW/TFA
CALCULATE THE REYNOLDS NUMBER: RET
RET : (GT*TID)/(12.0*2.42*TVlS)
CALCULATE THE PRANDTL NUT[BER: PRT
PRT = (TSHC*TVIS*2.42)/TFK
CHECK TI{A~ THE REYNOLDS NUHBER IS IN THE TRANSITION REGION
( i. e. 2000<RET<I0,000 )
IF (RET .GT. 2000.0 .AND. RET .LT. i0000.0) THEN
X1 : RET/1000.0
X2 = i0.0*TID/TL
x3
= (i.0-(Xi/i0.0)**0.256)
Xa = X2**X3
B! : (-3.08+3.075"XI+(0.32567*X!**2)-(0.02185*XI**3))*X4
HT = (16. I/TOD) * (BI*TFK*PRT**0.33)
GO TO ii0
ELSE IF (RET
.GT.
i0000.0) THEN
690
Ii0
HT = 0.027*TFK*(12.0/TID)*(RET**0.8)*(PRT**0.33)
ENDIF
HT = HT*TID/TOD
CALCULATE THE TUBE FRICTION FACTOR
IF (RET .GT. i000.0) THEN
TFRF = 0.0029/(RET**0.25603)
ELSE
ENDIF
CALCULATE THE FLUID VELOCITY, ft./s.
TV = GT/(3600.0*TDEN)
CALCULATE THE TUBE SIDE PRESSURE DROP, psi.
TDELP = (TFRF*(GT/3600.0)**2*TL*TNP)/(2.0*G*TDEN*TID/12.0)
CALCULATE THE PRESSURE DROP OF RETURN LOSSES, psi.
USING FRANK'S CORRELATION.
TRDELP = (2.5*TNP*TDEN*TV**2)/(2*G*I44.0)
CALCULATE THE TOTAL PRESSURE DROP EXCLUDING THE NOZZLE IN THE
TUBE SIDE, psi.
TOTDP = TDELP+TRDELP
RETURN
END
691
C
120
C
****************************************************************
*THIS PROGRAM CALCULATES THE SHELL SIDE }{EAT TRANSFER *
*COEFFICIENT, REYNOLDS NUMBER, PRANDTL NUMBER, FRICTION FACTOR *
*AND PRESSURE DROP. *
****************************************************************
SUBROUTINE SHELL
PARAMETER (PI = 3.1415927, G=32.174)
CHARACTER*30 FTS, FSS
CHARACTER TYPE
COMMON/TUBE1/ FTS, TYPE
COMMON/TUBE2/ NT, TID, TOD, TP, TL, TNP
COMMON/TUBE3/ TFF, TTI, TT2, TSHC, TVIS, TDEN
COMMON/TUBE4/ TFK, TW
COMMON/TUBE5/ RET, PRT, TFRF, HT, TV
COMMON/TUBE6/ TDELP, TRDELP, TOTDP
COMMON/SHELL1/ FSS
COMMON/SHELL2/ NS, SID, NB, BSP, SFF
COMMON/SHELL3/ TSI, TS2, SSHC, SVIS, SDEN
COMMON/SHELL4/ SFK, SW
COMMON/SHELL5/ RES, PRS, SFRF, HS, SV
COMMON/SHELL6/ SDELP
COMMON/OUT/ Q, AREA, F, LMTD, U, QI, CMTD
REAL NT, NS, NB, LMTD
CHECK WHETHER PITCH TYPE IS TRIANGUIJkR OR SQUARE
IF (TYPE .EQ. 'T' .OR. TYPE .EQ. 't') THEN
DE = ((I.72*TP**2)-(0.5*PI*TOD**2))/(0.5*PI*I2.0*TOD)
GO TO 120
ELSE IF (TYPE .EQ. 'S' .OR. TYPE .EQ. 's') THEN
DE = 4.*((TP**2)-(PI*TOD**2/4.0))/(PI*I2.0*TOD)
ENDIF
CALCULATE THE TUBE CLEARANCE, inch
TC = TP-TOD
CALCULATE THE SHELL FLOW AREA, ft^2.
IF (NS .GE. 2.0 .AND. NT .GE. 4.0) THEN
SFA =
(0.5*(SID*TC*BSP))/(144.0*TP)
GO TO 130
ELSE
SFA = (SID*TC*BSP)/(144.0*TP)
ENDIF
692
C
130
C
140
CALCULATE THE SHELL SUPERFICIAL FLOW RATE, Ib/ft^2.hr.
GS = SW/SFA
CALCULATE THE SHELL VELOCITY, ft./sec.
SV = SW/(3600.0*SDEN*SFA)
CALCULATE THE REYNOLDS NUMBER
RES = (DE*GS)/(2.42*SVIS)
CALCULATE THE PRANDTL NUMBER
PRS = (2.42*SSHC*SVIS)/SFK
CALCULATE THE SHELL SIDE HEAT TRANSFER COEFFICIENT
HS = 0.36/DE*SFK*(RES**0.55)*(PRS**0.333)
CHECK WHETHER THE SHELL'S REYNOLDS NUMBER IS LESS OR GREATER
THAN 500
IF (RES .LT. 500.0) THEN
SFRF = 0.II183/(RES**0.59246)
GO TO 140
ELSEIF (RES .ST. 500.0) THEN
SFRF
=
0.01159/(RES**0.18597)
ENDIF
CALCULATE THE SHELL SIDE PRESSURE DROP, psi.
VALI =
SFRF*(GS/3600.0)**2*SID/12.0*(NB+I.)
VAL2 = 64.4*SDEN*DE
SDELP = (VALI/VAL2)*NS
RETURN
END
693
***************************************************************
THIS PROGRAM PRINTS THE RESULTS OF THE SHELL AND TUBE HEAT
EXCHANGER ONTO A PRINTER
***************************************************************
SUBROUTINE RESULT
CHARACTER*30 FTS,FSS
CHARACTER TYPE
COMMON/TUBE1/ FTS, TYPE
COb~ON/TUBE2/
NT, TID, 'rOD, TP, TL, TNP
COMTi.ON/TUBE3/ TFF, TTI, TT2, TSHC, TVIS, TDEN
COMMON/TUBE4/ TFK, TW
COMMON~TUBES~
RET, PRT, TFRF, HT, TV
COMMON/TUBE6/ TDELP, TRDELP, TOTDP
COMMON/SHELL1/ FSS
COCk, ION/SHELL2/ NS, SID, NB, BSP, SFF
COMMON/SHELL3/ TSI, TS2, SSHC, SVIS, SDEN
COM}!ON/SHELL4/ SFK, SW
COMMON/SHELL5/ RES, PRS, SFRF, HS, SV
COMMON/SHELL6/ SDELP
COMMON/OUT/ Q, AREA, F, LMTD, U, QI, CMTD
REAL NT, NS, NB, LMTD
WRITE (I, 150)
150 FO~4AT (///,18X,'*SHELL AND TUBE HEAT EXCHANGER ~hTING*')
WRITE (i, 160)
160 FOm4AT (78(IH*))
WRITE (i, 170)
170 FORMAT (40X,'TUBE SIDE',IOX,'SHELL SIDE')
WRITE (i, 180)
180 FORMAT (78(IH-))
WRITE (i, 190) FTS, FSS
190 FORM_AT (2X,'FLUID TYPE:',T40, AI5, T60, AI5)
WRITE (i, 200) TW, SW
200 FORMAT (2X,'FLUID FLOW RATE
,
!b/hr.:',T40, FI0.0,T60,FI0.O)
WRITE (I, 210) TSHC, SSHC
210 FORMAT (2X,'SPECIFIC HEAT CAPACITY, Btu/ib.oF:',T40, FS.3,T60,
* F8.3)
WRITE (I, 220)TDEN, SDEN
220 FORMAT
(2X,'FLUID DENSITY,
Ib/ft^3:',T40,FS.3,T60,F8.3)
WRITE (i, 230) TVIS, SVIS
230 FORMAT (2X,'FLUID VISCOSITY, cP:',T40,FS.3,T60,FS.3)
WRITE (i, 240) TFK, SFK
240 FORMAT (2X,'FLUID THERMAL COND. Btu/ft.h.oF:',T40,F8.3,T60,FS.3)
WRITE (i, 250) TTI, TSI
250 FORMAT (2X,'FLUID INLET TEMPERATURE, oF:',T40,FS.0,T60,FS.0)
694
WRITE (i, 260) TT2, TS2
260 FORMAT (2X,'FLUID OUTLET TEMPERATURE, oF:',T40,F8.0,T60,F8.0)
WRITE (i, 270) TID, SID
270 FORMAT (2X,'INSIDE DIAMETER, inch:',T40,F8.3,T60,FS.3)
WRITE (i, 280) TOD
280 FORMAT (2X,'OUTSIDE DIAMETER, inch:',T40,F8.3)
IF (TYPE .EQ. 'T' .OR. TYPE .EQ. 't') THEN
WRITE (I, 290) TP
290 FORMAT (2X,'TUBE PITCH, inch:',T40,FS.3,T60,'TRIANGULAR PITCH')
ELSEIF (TYPE .EQ. 'S' .OR. TYPE .EQ. 's') THEN
WRITE (i, 300) TP
300 FORMAT (2X,'TUBE PITCH, iDch:',T40,F8.3,T60,'SQUARE PITCH')
ENDIF
WRITE (i, 310) TL
310 FORMAT (2X,'TUBE LENGTH, ft.:',T40,F8.3)
WRITE (i, 315) NT
315 FORMAT (2X,'NUMBER OF TUBES:', T40, F6.0)
WRITE (i, 320) TNP, NS
320 FORMAT (2X,'NUMBER OF PASSES:',T40,F4.0,T60,F4.0)
WRITE (I, 340) NB
340 FORMAT (2X,'NUMBER OF BAFFLES:',T60,F4.0)
WRITE (i, 350) BSP
350 FORMAT (2X,'BAFFLE SPACING, inch:',T60,F4.1)
WRITE (i, 360) TV, SV
360 FORMAT (2X,'FLUID VELOCITY, ft/sec.:',T40,FS.3,T60,F8.3)
WRITE (I, 370) RET, RES
370 FORMAT (2X,'REYNOLDS NUMBER:',T40,F9.0,T60,F9.0)
WRITE (I, 380) PRT, PRS
380 FORMAT (2X,'PRANDTL NUMBER:',T40,F8.4,T60,F8.4)
WRITE (i, 390) TFRF, SFRF
390 FORMAT (2X,'FRICTION FACTOR (ft^2/in^2):',T40,F8.5,T60,F8.5)
WRITE (I, 400) HT, HS
400 FORMAT (2X,'HEAT COEFF., Btu/hr.ft^2.oF:',T40,F9.l,T60,F9.1)
WRITE (i, 410) TOTDP, SDELP
410 FORMAT (2X,'PRESSURE DROP, psi:',T40,FS.4,T60,F8.4)
WRITE (i, 420) TFF, SFF
420 FORMAT (2X,'FOULING FACTOR:',T40,F8.4,T60,F8.4)
WRITE (I, 430) F
430 FORMAT (2X,'APPROACH FACTOR, F:', T40, F8.4)
695
440
450
460
470
480
490
500
WRITE (i, 440) LMTD
FORMAT (2X,'LOG MEAN TEMP. DIFF. oF:',T40,F8.2)
WRITE (I, 450) CMTD
FORMAT (2X,'CORRECTED MEAN TEMP. DIFF. oF:',T40,F8.2)
WRITE (I, 460) U
FORMAT (2X,'OVERALL HEAT COEFF., Btu/hr.ft^2.oF:',T40,F9.1)
WRITE (i, 470) AREA
FORMAT (2X,'OUTSIDE AREA OF UNIT, ft^2:',T40,F9.1)
WRITE (I, 480) Q
FORMAT (2X,'HEAT LOAD ON UNIT, Btu/hr:',T40,Fl2.0)
WRITE (i, 490) Q1
FORMAT (2X,'HEAT TRANSFERRED, Btu/hr:',T40,Fl2.0)
WRITE (i, 500)
FORMAT ( IH , 78 (IH-) )
FORM FEED THE PRINTING PAPER TO THE TOP OF THE NEXT PAGE.
WRITE (i, *) CHAR(12)
RETURN
END
696