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

VEnzin.for
c:/dna/source/
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Component name
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
100 CONTINUE
KOMTY = ’MEOH_CONVERTER’
MMVAR(1) = 0
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Component characteristics
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
200 CONTINUE
KOMTY = ’MEOH_CONVERTER’
ANTKN = 3
ANTPK = 1
ANTEX = 6
ANTM1 = 2
ZANAM(1) = ’MULTIPLIER H’
ZANAM(2) = ’MULTIPLIER C’
ZANAM(3) = ’MULTIPLIER O’
ZANAM(4) = ’GIBBS ENERGY’
ZANAM(5) = ’T−COND’
MEDIE(1) = ANYGAS$
MEDIE(2) = ANYGAS$
MEDIE(3) = HEAT$
ANTME = 2
VARME(1) = NODE1$
VARME(2) = NODE2$
C
204 ANTEL(2) = 37
ANTLK = 1+ANTEL(2)
VAREL(1,2) = H2$
VAREL(2,2) = CO$
VAREL(3,2) = CO2$
VAREL(4,2) = H2O_G$
VAREL(5,2) = CH3OH$
VAREL(6,2) = 2
VAREL(7,2) = 3
VAREL(8,2) = 5
DO I=8,36
VAREL(I+1,2) = I
ENDDO
C
IF (FKOMP.EQ.6) GOTO 600
IF (CON) GOTO 404
*** FKOMP = 3
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Component equations. All in residual form.
C Do not include the conservation laws, since these are treated
C automatically by DNA.
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
400 CONTINUE
C
CON=.true.
GOTO 204
404 Continue
C
C Pressure
C
ANTGIBBS=5
RES(1) = P(1) − P(2) − PAR(1)
RES_NR=1
C
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Find the composition of the equilibrium gas
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C
C Calculate mole flow of each species in
C
M_BL(1)=0.D0
DO I=1,ANTST
M_BL(1)=X_J(MEDIE(1),I)*M_MOL(I)+M_BL(1)
ENDDO
DO I=1,ANTST
NIN(I)= MDOT(1)/M_BL(1)*X_J(MEDIE(1),I)
ENDDO
C
C Calculate mole flow of each species out
25/67
19−03−2007