27.07.2013 Views

Design og modellering af metanolanlæg til VEnzin-visionen Bilag

Design og modellering af metanolanlæg til VEnzin-visionen Bilag

Design og modellering af metanolanlæg til VEnzin-visionen Bilag

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>VEnzin</strong>.for<br />

c:/dna/source/<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

C Component name<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

100 CONTINUE<br />

KOMTY = ’GASCLE_2’<br />

MMVAR(1) = ANTST<br />

GOTO 9999<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

C Component characteristics<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

200 CONTINUE<br />

KOMTY = ’GASCLE_2’<br />

ANTKN = 4<br />

ANTPK = 1+MMVAR(1)<br />

ANTLK = 4+2*ANTST<br />

ANTM1 = 3<br />

MEDIE(1) = −4<br />

MEDIE(2) = −4<br />

MEDIE(3) = −4<br />

MEDIE(4) = 300<br />

ANTME = 3<br />

VARME(1) = −1<br />

VARME(2) = −2<br />

VARME(3) = −3<br />

ANTEL(1) = 0<br />

ANTEL(2) = ANTST<br />

ANTEL(3) = ANTST<br />

DO I = 1,ANTST<br />

VAREL(I,2) = I<br />

VAREL(I,3) = I<br />

ENDDO<br />

C<br />

IF (FKOMP.EQ.6) GOTO 600<br />

*** FKOMP = 3<br />

GOTO 9999<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

C Component equations. All in residual form.<br />

C Do not include the conservation laws, since these are treated<br />

C automatically by DNA.<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

400 CONTINUE<br />

DELP = PAR(1)<br />

Nr_clean=ANTPK−1<br />

C<br />

DO I = 1,Nr_clean<br />

Compound_clean(I)=PAR(I+1)<br />

ENDDO<br />

C<br />

M_BL1 = 0.D0<br />

DO I=1, ANTST<br />

M_BL1 = M_BL1 + X_J(MEDIE(1),I)*M_MOL(I)<br />

XUD(I) = 0.D0<br />

ENDDO<br />

C<br />

DO I=1, ANTST<br />

N(I) = X_J(MEDIE(1),I)*MDOT(1)/M_BL1<br />

ENDDO<br />

C<br />

NTOT=0D0<br />

M3=0D0<br />

DO I=1, Nr_clean<br />

NTOT = NTOT+N(Compound_clean(I))<br />

M3=M3+N(Compound_clean(I))*M_MOL(Compound_clean(I))<br />

ENDDO<br />

C<br />

DO I=1, Nr_clean<br />

XUD(Compound_clean(I)) = N(Compound_clean(I))/NTOT<br />

N(Compound_clean(I)) = 0D0<br />

ENDDO<br />

C<br />

T_SUM = 0D0<br />

DO I=1, ANTST<br />

T_SUM = T_SUM + N(I)<br />

ENDDO<br />

C<br />

C Find mole ratios in syngas<br />

C<br />

CALL STATES(P(2),H(2),T2,V,S,X,U,1,2,MEDIE(2))<br />

14/67<br />

19−03−2007

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!