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

VEnzin.for
c:/dna/source/
** 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
C Calculate mole flow in
C
M_BL1 = 0.D0
DO I=1,ANTST
M_BL1=M_BL1+X_J(MEDIE(1),I)*M_MOL(I)
ENDDO
NIN = MDOT(1)/M_BL1
NSTIN = X_J(MEDIE(1),7)*NIN
N = NIN−NSTIN
C
C Calculate mole flow out
C
M_BL2 = 0.D0
DO I=1,ANTST
M_BL2=M_BL2+X_J(MEDIE(2),I)*M_MOL(I)
ENDDO
NOUT = (−MDOT(2))/M_BL2
C
C Same temperature out
C
X0 = 0.D0
CALL STATES(P(2),H(2),T2,V,S,X,U,1,2,MEDIE(2))
X0 = 0.D0
CALL STATES(PSTOUT,HD,T2,V,S,X0,U,3,6,MEDIE(3))
CALL STATES(P(3),H2,T2,V,S,X,U,1,3,MEDIE(3))
RES(1) = H2 − H(3)
C
C Molar flow and ratio of steam out
C
XSTOUT = PSTOUT/P(2)
IF (XSTOUT.GT.X_J(MEDIE(1),7)) XSTOUT = X_J(MEDIE(1),7)
NSTOUT = XSTOUT*(NIN−NSTIN)/(1−XSTOUT)
C
C Molar flow of condensate
C
NCOND = NSTIN − NSTOUT
c if (NCOND.lt.2D−7) NCOND=1D−7
RES(2) = MDOT(3) + NCOND*M_MOL(7)
C
C Pressure out
C
RES(3) = P(1) − P(2) − PAR(1)
RES(4) = P(2) − P(3)
RES(5) = P(4) − P(5) − PAR(2)
C
C Outlet composition of gas
C
DO I = 1,6
RES(I+5) = X_J(MEDIE(2),I)*NOUT − X_J(MEDIE(1),I)*NIN
ENDDO
RES(7+5) = X_J(MEDIE(2),7) − XSTOUT
DO I = 8,ANTST
RES(I+5) = X_J(MEDIE(2),I)*NOUT − X_J(MEDIE(1),I)*NIN
ENDDO
C −−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Pinch point analysis
C −−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C
C Check if pinch point is where the water condenses
C
DTP = PAR(3)
EVAP=.false.
if (MDOT(3).lt.−1d−5) then
EVAP=.true.
X0 = 0.0D0
CALL STATES(P(1)*X_J(MEDIE(1),7),HD,TPH,V,S,X0,U,1,6,
$ MEDIE(3))
CALL STATES(P(1),HSAT0,TPH,V,S,X,U,1,3,MEDIE(1))
HPH = H(5) − MDOT(1)*(H(1)−HSAT0)/MDOT(4)
21/67
19−03−2007