Elementary Principles 467 - HR = HR * EXT C Product Heat...

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

View Full Document Right Arrow Icon
10- 19 10.12 (cont’d) 900 * * * * FORMAT ('Product Stream', F7.3, ' mols/s carbon monoxide', /, 15X, F7.3, 'mols/s oxygen', /. 15X, F7.3, 'mols/s nitrogen', /. 15X, F7.3, 'mols/s carbon dioxide', /, 15X, F7.2, 'C') END C SUBROUTINE REACTAD (SF, SP, NU, N, X, IX) DIMENSION SF(8), SP(8), NU(7), ACP(7), BCP(7), CCP(7), DCP(7), HF(7) COMMON ACP, BCP, CCP, DCP, NF TOL = 1.E-6 C Extent of Reaction EXT = –SF(IX)*X/NU(IX) C Solve Material Balances DO 100 I = 1, N 100 SP(I) = SF(I) + EXT*NU(I) C Heat of Reaction HR = 0 DO 200 I = 1, N 200 HR = HR + HF(I) * NU(I)
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: HR = HR * EXT C Product Heat Capacity AP = 0. BP = 0. CP = 0. DP = 0. DO 300 I = 1, N AP = AP + SP(I)*ACP(I) BP = BP + BP(I)*BCP(I) CP = CP + SP(I)*CCP(I) 300 DP = DP + SP(I)*DCP(I) C Find T TIN = SF (N + 1) TP = TIN D0 400 ITER = 1, 10 T = TP F = HR FP = 0. F = F +T*(AP + T*(BP/2. + T*(CP/3. + T*DP/4.))) * TIN*(AP + TIN*(BP/2. + TIN*(CP/3. + TIN*DP/4.))) FP = FP + AP + T *(BP + T*(CP + T*DP)) TP = T F/FP IF(ABS((TP T)/T).LT.TOL) GOTO 500 400 CONTINUE WRITE (6, 900) 900 FORMAT ('REACTED did not converge') STOP...
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