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

VEnzin.for
c:/dna/source/
KOMTY = ’SET_M’
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Component characteristics
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
200 CONTINUE
KOMTY = ’SET_M’
ANTKN = 3
ANTPK = 0
ANTLK = 2
ANTM1 = 2
MEDIE(1) = ANYGAS$
MEDIE(2) = ANYGAS$
MEDIE(3) = 999
ANTME = 2
VARME(1) = NODE1$
VARME(2) = NODE1$
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. These are treated automatically
C by DNA.
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
400 CONTINUE
C
RES(1) = P(1) − P(2)
RES(2) = ZC(1)−(X_J(MEDIE(1),H2$)−X_J(MEDIE(1),CO2$))/
$ (X_J(MEDIE(1),CO2$)+X_J(MEDIE(1),CO$))
C
IF (FKOMP.EQ.5) GOTO 500
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Solution check
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
500 CONTINUE
IF (MDOT(1).LT.−1D−10) GOTO 550
IF (MDOT(2).GT.1D−10) GOTO 550
IF (X_J(MEDIE(1),CO$)+X_J(MEDIE(1),CO2$).LT.1D−3) GOTO 550
GOTO 9999
550 FBETI = .FALSE.
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Write component information
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
600 CONTINUE
KOMDSC = ’Utility component for setting the M−factor for a syngas
$ used for production of a liquid fuel. The M−factor is defined
$ in the equations below. If the M−factor is 2 for methanol
$ production it means that all the syngas in theory can be
$ converted to methanol. The M−factor is defined as a control
$ variable (ZC).’
K_LIG(1) = ’Equal pressures: $p_1 = p_2$’
K_LIG(2) =
$ ’M−factor: $ZC(1)=\\frac{y_{H_2}−y_{CO_2}}{y_{CO_2}+y_{CO}}$’
K_BET = ’$\\dot{m_1} \\gt 0 \\\\ \\dot{m_2} \\lt 0 \\\\
$y_{CO}+y_{CO_2} \\gt 0$’
KMEDDS(1) = ’Gas in’
KMEDDS(2) = ’Gas out’
KMEDDS(3) = ’M−factor’
K_INP=’struc set−M SET_M 611 612 900\\\\
$MEDIA 611 gas\\\\
$fluid gas N2 0.1 H2 0.6 CO 0.1 CO2 0.2\\\\
$addco m set−M 611 1 t set−M 611 50 p 611 1’
C
GOTO 9999
C
9999 CONTINUE
RETURN
END
C
C=======================================================================
C***********************************************************************
SUBROUTINE GASCOOL4(KOMTY,ANTLK,ANTKN,ANTPK,ANTM1,ANTM2,MEDIE,
$ ANTME,VARME,ANTEL,VAREL,MDOT,P,H,Q,PAR,RES,X_J,ZA,ZANAM,ZC,
$ ANTEX,KOMDSC,K_PAR,K_lig,K_bet,KMEDDS,K_inp)
C***********************************************************************
29/67
19−03−2007