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

VEnzin.for
c:/dna/source/
C
C Parameter variables
C
INTEGER ANTLK, ANTEX,ANTKN, MEDIE(5), ANTPK, ANTM1,
$ ANTME,VARME(3),varel(antst,3),antel(3)
DOUBLE PRECISION X_J(MAXME,ANTST), PAR(2), RES(8), MDOT(3),
: P(3),H(3),E,Q,ZA
CHARACTER*80 KOMTY,ZANAM(1)
C
C Local variables
C
DOUBLE PRECISION Virk, Temp, T2, T3, S, V, X, U
CHARACTER*100 K_PAR(2)
CHARACTER*1000 KOMDSC,K_INP
CHARACTER*500 K_LIG(8), K_BET
CHARACTER*100 KMEDDS(5)
EXTERNAL STATES
INTRINSIC DABS
INCLUDE ’THERPROP.INI’
C========================================================================
GOTO (100,200,1,400,400,200) FKOMP
1 RETURN
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Component name
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
100 CONTINUE
KOMTY = ’ELECTROLYSER’
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Component characteristics
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
200 CONTINUE
KOMTY = ’ELECTROLYSER’
ANTKN = 5
ANTPK = 2
ANTLK = 8
ANTEX = 0
ANTM1 = 3
MEDIE(1) = 97
MEDIE(2) = ANYGAS$
MEDIE(3) = ANYGAS$
MEDIE(4) = power$
MEDIE(5) = heat$
ANTME = 3
VARME(1) = NODE1$
VARME(2) = NODE2$
antel(2) = 1
VARME(3) = NODE3$
antel(3) = 1
varel(1,2) = H2$
varel(1,3) = O2$
IF (FKOMP.EQ.6) GOTO 600
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
Temp = PAR(1)
Virk = PAR(2)
C
CALL STATES(P(2),H(2),T2,V,S,X,U,1,2,MEDIE(2))
CALL STATES(P(3),H(3),T3,V,S,X,U,1,2,MEDIE(3))
C
res(1)=x_j(medie(2),H2$)−1.0
res(2)=x_j(medie(3),O2$)−1.0
res(3)=mdot(2)+mdot(1)/m_mol(H2O_L$)*m_mol(H2$)
res(4)=P(1)−P(2)
res(5)=P(1)−P(3)
res(6)=T2−Temp
res(7)=T2−T3
c res(8)=(E+Q)/E−Virk
res(8)=E*virk+ned_H(H2$)/m_mol(H2$)*MDOT(2)
C
IF (FKOMP.EQ.5) GOTO 500
2/67
19−03−2007