Jawaharlal Nehru Technological University Kakinada University College of Engineering (A) Dept. of PetroleumKakinada Engi
Views 44 Downloads 0 File size 521KB
Jawaharlal Nehru Technological University Kakinada University College of Engineering (A) Dept. of PetroleumKakinada Engineering & Petrochemical Engineering
B. Tech. Petrochemical Engineering Multicomponent Distillation PROGRAM ON MULTICOMPONENT-Lewis Matheson Method using Bonner's Feed match technique
Prof. K. V. Rao Programme Director Petroleum Courses JNTUK
Lewis – Matheson Method: Feed Plate Matching The correction method of Bonner can be used to establish a new value of D, di, and xiD,
Where (xF)T and (xF)B = Liquid compositions calculated from the top and bottom of the column, respectively di, bi = the former moles of individual component in the distillate and bottoms diN, biN = the new values Repetitive calculations are carried out by readjusting the mole of components in the distillate and bottoms until the agreement at the feed plate is met.
REM PROGRAM ON MULTICOMPONENT-Lewis Matheson Method using Bonner's Feed match technique DIM
Y(10,
YHK(100),
100),
XF(10),
X(10, A(10),
100), D1(10),
DELD(10), ALPHA(10) READ NC, NT, NAF, NBF, HK, LK FOR I = 1 TO NC READ Y(I, 1), X(I, NT), ALPHA(I) NEXT I READ D, B, RV, RL, SV, SL, R, Q NF = NAF + 1 NIT = 1
XHK(100), BI(10),
9 I = 1 10 YHK(I) = Y(HK, I) SUM = 0 FOR J = 1 TO NC SUM = SUM + Y(J, I) / (ALPHA(J) * YHK(I)) NEXT J XHK(I) = 1 / SUM FOR J = 1 TO NC X(J, I) = Y(J, I) * XHK(I) / (ALPHA(J) * YHK(I)) NEXT J K = I + 1 FOR J = 1 TO NC Y(J, K) = X(J, I) * (R / (R + 1)) + Y(J, 1) / (R + 1): NEXT J IF K > NF THEN 100 I = K: GOTO 10
100 FOR J = 1 TO NC XF(J) = X(J, NF) NEXT J M = NT 110 XHK(M) = X(HK, M) SUM = 0 FOR J = 1 TO NC SUM = SUM + ALPHA(J) * X(J, M) / XHK(M) NEXT J YHK(M) = 1 / SUM FOR J = 1 TO NC Y(J, M) = ALPHA(J) * YHK(M) * X(J, M) / XHK(M) NEXT J L = M - 1 FOR J = 1 TO NC X(J, L) = SV * Y(J, M) / SL + B * X(J, NT) / SL NEXT J
IF L = NF THEN 200 M = L: GOTO 110 200 FOR J = 1 TO NC A(J) = ABS(XF(J) - X(J, NF)) IF A(J)>10 E-4 THEN 300 NEXT J: GOTO 400 300 NIT = NIT + 1: PRINT "NIT="; NIT IF NIT > 100 THEN 1000 SD = 0: SB = 0 FOR J = 1 TO NC DI(J) = D * Y(J, 1): BI(J) = B * X(J, NT)
DELD(J) = (X(J, NF) - XF(J)) / (XF(J) / DI(J) + X(J, NF) / BI(J)) DI(J) = DI(J) + DELD(J): BI(J) = BI(J) - DELD(J): SD = SD + DI(J): SB = SB + BI(J) NEXT J: D = SD: B = SB FOR J = 1 TO NC Y(J, 1) = DI(J) / D: X(J, NT) = BI(J) / B: NEXT J: SL = R * D + Q * F: SV = D * (R + 1) - F * (1 - Q): GOTO 9
400 PRINT NT, NAF, NBF, NIT FOR J = 1 TO NT FOR I = 1 TO NC PRINT J, X(I, J), Y(I, J) NEXT I, J FOR I = 1 TO NC PRINT I, DI(I), BI(I) NEXT I PRINT D, B: END 1000 PRINT "LIMIT ON NUMBER OF ITERATIONS EXCEEDED" END
DATA 4,16,6,9,3,2 DATA .19,.0005,2.75 DATA .729,.0135,2.11 DATA .069,.449,1 DATA .012,.537,.8 DATA 25,75,100,75,100,175,2.9125,1
Thank You