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

VEnzin.for
c:/dna/source/
RES(5) = H3 − H(3)
RES(6) = H4 − H(4)
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
NINST=NIN*X_J(MEDIE(1),7)
NINME=NIN*X_J(MEDIE(1),CH3OH$)
N=NIN−NINST−NINME
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 Outlet composition of gas
C
DO I = 1,6
RES(I+6) = X_J(MEDIE(2),I)*NOUT − X_J(MEDIE(1),I)*NIN
ENDDO
DO I = 8,35
RES(I+5) = X_J(MEDIE(2),I)*NOUT − X_J(MEDIE(1),I)*NIN
ENDDO
C
RES(41)=1.0D0
DO I = 1,ANTST
RES(41) = RES(41)−X_J(MEDIE(2),I)
ENDDO
C
C Check for condensing water or methanol
C
b_2_1=−1062.945621D0
b_1_2=3538.709318D0
c J/mol
alpha=0.2994D0
R_u=8.314D0
c J/(mol*K)
c
T=ZA(4)
X0 = 0.D0
CALL STATES(P_SAT_ST,HD,T,V,S,X0,U,3,6,MEDIE(3))
X0 = 0.D0
CALL STATES(P_SAT_ME,HD,T,V,S,X0,U,3,6,MEDIE(4))
T_K=T+273.15D0
c
tau_2_1=b_2_1/(R_u*T_K)
tau_1_2=b_1_2/(R_u*T_K)
c
x_ME_3=ZA(8)
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)))))
c
P_ST_IN=X_J(MEDIE(1),7)*P(1)
P_ME_IN=X_J(MEDIE(1),CH3OH$)*P(1)
c
P_ME=P_ME_IN*0.999d0
x_ME=P_ME/(gamma_2*P_SAT_ME)
x_ST=1D0−x_ME
P_ST_2=gamma_1*P_SAT_ST*x_ST
P_ST=(P_ME_IN−P_ME+P_ST_IN)−(P_ME_IN−P_ME)/x_ME
if (T.gt.239.1) then
RES(52)=T−238
else
RES(52)=P_ST−P_ST_2
32/67
19−03−2007