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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

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

c:/dna/source/<br />

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

100 CONTINUE<br />

KOMTY = ’MIXING_TANK’<br />

GOTO 9999<br />

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

C Component characteristics<br />

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

200 CONTINUE<br />

KOMTY = ’MIXING_TANK’<br />

ANTKN = 5<br />

ANTPK = 2<br />

ANTLK = 3<br />

ANTEX = 2<br />

ANTM1 = 2<br />

ANTM2 = 2<br />

MEDIE(1) = REALFL$<br />

MEDIE(2) = REALFL$<br />

MEDIE(3) = REALFL$<br />

MEDIE(4) = REALFL$<br />

MEDIE(5) = HEAT$<br />

ANTME = 4<br />

VARME(1) = NODE1$<br />

VARME(2) = NODE1$<br />

VARME(3) = NODE3$<br />

VARME(4) = NODE3$<br />

PARNAM(1) = ’Pressure loss side 1’<br />

PARNAM(2) = ’Pressure loss side 2’<br />

ZANAM(1) = ’Flu 2 mass−%’<br />

ZANAM(2) = ’Flu 2 vol−%’<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. These are treated automatically<br />

C by DNA.<br />

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

400 CONTINUE<br />

C<br />

DPA = PAR(1)<br />

DPB = PAR(2)<br />

C<br />

C Pressure losses<br />

C<br />

RES(1) = P(1) − P(2) − DPA<br />

RES(2) = P(3) − P(4) − DPB<br />

C<br />

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

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

RES(3) = T2−T4<br />

C<br />

RES(4) = ZA(1) − MDOT(3)/(MDOT(3)+MDOT(1))<br />

RES(5) = ZA(2) − (MDOT(3)/RM_MOL(MEDIE(3)))/(MDOT(3)<br />

$ /RM_MOL(MEDIE(3))+MDOT(1)/RM_MOL(MEDIE(1)))<br />

C<br />

IF (FKOMP.EQ.5) GOTO 500<br />

GOTO 9999<br />

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

C Solution check<br />

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

500 CONTINUE<br />

IF (MDOT(1).LT.−1D−10) GOTO 550<br />

IF (MDOT(2).GT.1D−10) GOTO 550<br />

IF (MDOT(3).LT.−1D−10) GOTO 550<br />

IF (MDOT(4).GT.1D−10) GOTO 550<br />

IF (Q.GT.1D−10) GOTO 550<br />

GOTO 9999<br />

550 FBETI = .FALSE.<br />

GOTO 9999<br />

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

C Write component information<br />

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

600 CONTINUE<br />

KOMDSC = ’Heat exchanger used as a model for a mixing tank. This<br />

$ means that the outlet temperatures are set equal.’<br />

K_PAR(1) = ’Pressure loss side 1, $\\Delta p_{12}$ [bar]’<br />

K_PAR(2) = ’Pressure loss side 2, $\\Delta p_{34}$ [bar]’<br />

K_BET =<br />

37/67<br />

19−03−2007

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

Saved successfully!

Ooh no, something went wrong!