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

VEnzin.for
c:/dna/source/
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
C Ash characteristics
C
CP(MEDIE(5)) = 1
HF(MEDIE(5)) = −5083.0D0
HV(MEDIE(5)) = X_J(MEDIE(5),28)*NED_H(28)/M_MOL(28)
C
CALCOM(1) = −1
CALCOM(2) = −2
CALCOM(3) = −3
CALCOM(4) = −4
CALCOM(5) = −5
CALCOM(6) = −6
CALCOM(7) = −7
CALCOM(8) = −8
CALCOM(9) = −9
CALCOM(10) = −10
CALCOM(11) = −11
CALCOM(12) = −30
CALCOM(13) = −31
CALCOM(14) = −32
CALCOM(15) = −36
DO I=1,ANTPK−5
DO J=1,15
IF (CALCOM(J).EQ.−NINT(PAR(I))) THEN
CALCOM(J)=−CALCOM(J)
ENDIF
ENDDO
ENDDO
PGAS = PAR(ANTPK−4)
TGAS = PAR(ANTPK−3)+273.15
DELP = PAR(ANTPK−2)
c DMVC = PAR(ANTPK−2)
CC = PAR(ANTPK−1)
cbe 20041711 parameter for "non−equilibrium" methane
METH = PAR(ANTPK)
R = 8.314D0
C
C Pressure
C
RES(1) = P(2) − P(3)
RES(2) = P(3) − P(4) − DELP
RES(3) = P(4) − PGAS
C
C Amount of water relative to coal
C
c RES(4) = MDOT(2) − DMVC*MDOT(1)
RES(4) = MDOT(2) − ZC(1)*MDOT(1)
C
C Ash
C
RES(5) = MDOT(5) +
: MDOT(1)*(X_J(MEDIE(1),38)+X_J(MEDIE(1),28)*(1.0D0−CC))
IF (MDOT(5).EQ.0D0) THEN
RES(6) = X_J(MEDIE(5),38)
ELSE
RES(6) = X_J(MEDIE(5),38) + MDOT(1)*X_J(MEDIE(1),38)/MDOT(5)
ENDIF
RES(7) = X_J(MEDIE(5),28) − (1.0D0−X_J(MEDIE(5),38))
CALL STATES(P(4),H(4),T4,V,S,X,U,1,2,MEDIE(4))
CALL STATES(P(5),H5,T4,V,S,X,U,1,3,MEDIE(5))
RES(8) = H(5) − H5
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Find the composition of the equilibrium gas
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C
C Calculate mole flow of each species in
C
M_BL(3) = 0D0
DO I=1,ANTST
M_BL(3)=M_BL(3)+X_J(MEDIE(3),I)*M_MOL(I)
9/67
19−03−2007