Elementary Principles 323 - 8.66 (cont'd) c. C* C* 1 C*...

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

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: 8.66 (cont'd) c. C* C* 1 C* PROGRAM FOR PROBLEM 8.66 IMPLICIT REAL (N) READ (5, 1) A1, B1, C1, A2, B2, C2 ANTOINE EQUATION COEFFICIENTS FOR A AND B FORMAT (8F10.4) READ (5, 1) TRA, TRB ARBITRARY REFERENCE TEMPERATURES (DEG.C.) FOR A AND B READ (5, 1) CAL, TBPA, DHVA, CAV1, CAV2 READ (5, 1) CBL, TBPB, DHVB, CBV1, CBV2 CP(LIQ, KS/MBL-DEG.C.), NORMAL BOILING POINT (DEG.C), HEAT OF VAPORIZATION (KJ/MOL), COEFFICIENTS OF CP(VAP., KJ/MOL-DEG.C) = CV1 + CV2*T(DEG.C) READ (5, 1) XF, TF, P MOLE FRACTION OF A IN FEED, FEED TEMP.(DEG.C), EVAPORATOR PRESSURE (MMHG) WRITE (6, 2) TF, XF, P FORMAT (1H0, 'FEEDbATb', F6.1, 'bDEG.CbCONTAINSb', F6.3,' bMOLESbA/MOLEbT *OTAL'//1X'EVAPORATORbPRESSUREb=', E11.4, 'bMMbHG'/) ITER = 0 DT = 0.5 HAF = CAL*(TF TRA) HBF = CBL*(TF TRB) F1 = XF*HAF + (1.0 XF)*HBF F2 = CAL*(TBPA TRA) + DHVA CAV1*TBPA 0.5*CAV2*TBPA**2 F3 = CBL*(TBPB TRB) + DHVB CBV1*TBPB 0.5*CBV2*TBPB**2 T = TF INTER = ITER + 1 IF(ITER 200) 30, 30, 25 WRITE (6, 3) FORMAT (1H0, 'NO CONVERGENCE') STOP PAV = 10.0** (A1 B1/(T + C1)) PAV = 10.0** (A2 B2/(T + C2)) XL = (P PBV)/(PAV PBV) XV = XL*PAV/P NL = (XV XF)/(XV XL) NV = 1.0 NL IF (XL.LE.00.OR.XL.GE.1.0.OR.NL.LE.0.0.OR.NL.GE.1.0) GO TO 45 HAL = CAL*(T TRA) HBL = CBL*(T TRB) HAV = F2 + CAV1*T + 0.5*CAV2*T**2 HBV = F3 + CBV1*T + 0.5*CBV2*T**2 C* C* C* 2 20 25 3 30 8-44 ...
View Full Document

This note was uploaded on 02/27/2012 for the course CHEMICAL E 312 taught by Professor Cheung during the Fall '11 term at The University of Akron.

Ask a homework question - tutors are online