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

VEnzin.for
c:/dna/source/
c
c
c
E2=(H(1)−HPH)*MDOT(1)
RES(46) = ZA(9)
else
X0 = 0.D0
CALL STATES(P_SAT_ST,HD,TPL1+ZA(1),V,S,X0,U,3,6,MEDIE(3))
X0 = 0.D0
CALL STATES(P_SAT_ME,HD,TPL1+ZA(1),V,S,X0,U,3,6,MEDIE(4))
T_K=TPL1+ZA(1)+273.15D0
tau_2_1=b_2_1/(R_u*T_K)
tau_1_2=b_1_2/(R_u*T_K)
x_ME_3=ZA(9)
x_ST_3=1D0−x_ME_3
gamma_2=exp((x_ST_3**2*(tau_1_2*((exp(−(alpha
$ *tau_1_2))/(x_ME_3+x_ST_3*exp(−(alpha*tau_1_2))))
$ )**2+(tau_2_1*(exp(−(alpha*tau_2_1))/(x_ST_3+x_ME_3
$ *exp(−(alpha*tau_2_1)))**2)))))
gamma_1=exp((x_ME_3**2*(tau_2_1*((exp(−(alpha
$ *tau_2_1))/(x_ST_3+x_ME_3*exp(−(alpha*tau_2_1))))
$ )**2+(tau_1_2*(exp(−(alpha*tau_1_2))/(x_ME_3+x_ST_3
$ *exp(−(alpha*tau_1_2)))**2)))))
P_ME=gamma_2*P_SAT_ME*x_ME_3
x_ST=1D0−x_ME_3
P_ST=gamma_1*P_SAT_ST*x_ST
NOUTME=P_ME/P(2)*NOUT
NOUTST=P_ST/P(2)*NOUT
NCONDST=NINST−NOUTST
NCONDME=NINME−NOUTME
x_ME_2=NCONDME/(NCONDME+NCONDST)
RES(46)=x_ME_3−x_ME_2
MPH=MDOT(1)−NCONDST*M_MOL(7)−NCONDME*M_MOL(CH3OH$)
CALL STATES(P(3),HST,TPL1+ZA(1),V,S,X0,U,1,3,MEDIE(3))
CALL STATES(P(4),HME,TPL1+ZA(1),V,S,X0,U,1,3,MEDIE(4))
CALL ENTHALPY(7,TPL1+ZA(1),HSTV)
CALL ENTHALPY(CH3OH$,TPL1+ZA(1),HMEV)
HPH=(HPH*MDOT(1)−HSTV*NCONDST*M_MOL(7)−HMEV*NCONDME
$ *M_MOL(CH3OH$))/MPH
E2=H(1)*MDOT(1)−(HPH*MPH+HST*NCONDST*M_MOL(7)+HME
$ *NCONDME*M_MOL(CH3OH$))
endif
RES(48)=E2−E1
RES(2) = (TPL1+ZA(1)) − ZA(3)
if (ZA(1).lt.(TPH−TPL)) then
TPH=TPL1+ZA(1)
TPL=TPL1
endif
ELSE
RES(46) = ZA(9)
RES(48) = ZA(1)
RES(2) = ZA(3)
ENDIF
ELSE
RES(46) = ZA(9)
RES(48) = ZA(1)
RES(2) = ZA(3)
ENDIF
C
C Check if pinch point is at the beginning or the end of the component
C
CALL STATES(P(1),H(1),T1,V,S,X,U,1,2,MEDIE(1))
CALL STATES(P(2),H(2),T2,V,S,X,U,1,2,MEDIE(2))
CALL STATES(P(5),H(5),T5,V,S,X,U,1,2,MEDIE(5))
CALL STATES(P(6),H(6),T6,V,S,X,U,1,2,MEDIE(6))
C
if (((T2−T5).lt.(TPH−TPL)).or.((T1−T6).lt.(TPH−TPL))) then
IF ((T1−T6).LT.(T2−T5)) THEN
TPH = T1
TPL = T6
ELSE
TPH = T2
TPL = T5
ENDIF
ENDIF
C
C At the pinch point DTP is used
34/67
19−03−2007