Elementary Principles 459 - TF = SF (N + 1) T = SL (N + 1)...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
10- 11 10.9 (cont’d) 900 * * * FORMAT (A15, F7.4,' mols/s n-pentane', /, 15X, F7.4,' mols/s n-hexane', /, 15X, F7.4,' mols/s n-hephane', /, 15X, F7.2,' K') WRITE (6, 901) Q 901 FORMAT ('Heat Required', F7.2, 'kW') END C SUBROUTINE FLASHIN (SF, SL, SV, N, P, Q) REAL NF, NL, NV, NVP DIMENSION SF(8), SL(8), SV(8) DIMENSION XF(7), XL(7), XV(7), PV(7) DIMENSION A(7), B(7), C(7), CP(7), HV(7) COMMON A, B, C, CP, HV TOL = 1,5 – 6 C Feed Variables NF = 0. DO 100 I = 1, N 100 NF = NF + SF(I) DO 200 I = 1, N 200 XF(I) = SF(I)/NF
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: TF = SF (N + 1) T = SL (N + 1) TC = T 273.15 C Vapor Pressures DO 300 I = 1, N 300 PV(I) = 10.**(A(I) B(I)/(TC + C(I))) C Find NV -- Initial Guess = NF/2 NVP = NF/2 DO 400 ITER = 1, 10 NV = NVP F = 1. FP = 0. DO 500 I = 1, N PPM1 = PV(I)/P 1. F = F + NF*XF(I)/(NF + NV*PPM1) 500 FP = FP PPM1*XF(I)/(NF + NV*PPM1)**2. NVP = NV F/FP IF (ABS((NVP NV)/NVP).LT.TOL) GOTO 600 400 CONTINUE WRITE (6, 900) 900 FORMAT ('FLASHN did not converge on NV') STOP C Other Variables...
View Full Document

Ask a homework question - tutors are online