Design og modellering af metanolanlæg til VEnzin-visionen Bilag
Design og modellering af metanolanlæg til VEnzin-visionen Bilag
Design og modellering af metanolanlæg til VEnzin-visionen Bilag
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
<strong>VEnzin</strong>.for<br />
c:/dna/source/<br />
CALL STATES(P(3),H(3),T3,V,S,X,U,1,2,MEDIE(3))<br />
C<br />
C Pressure<br />
C<br />
RES(1) = P(3) − P(2)<br />
RES(2) = P(1) − P(2) − DELP<br />
C<br />
C Outlet temperatures are identical<br />
C<br />
RES(3) = T2 − T3<br />
C<br />
C Mass flow of removed gas found using M1 and X_J<br />
C<br />
RES(4) = MDOT(3) + M3<br />
C<br />
C Variable mole ratios (first removed gas, then syngas)<br />
C<br />
DO I = 1,ANTST<br />
RES(4+I)=X_J(MEDIE(3),I)−XUD(I)<br />
RES(4+ANTST+I)=X_J(MEDIE(2),I)−N(I)/T_SUM<br />
ENDDO<br />
C<br />
IF (FKOMP.EQ.5) GOTO 500<br />
GOTO 9999<br />
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />
C Solution check<br />
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />
500 CONTINUE<br />
IF (MDOT(1).LT.−1D−10) GOTO 550<br />
IF (MDOT(2).GT.1D−10) GOTO 550<br />
IF (MDOT(3).GT.1D−10) GOTO 550<br />
IF (Q(1).GT.1D−10) GOTO 550<br />
GOTO 9999<br />
550 FBETI = .FALSE.<br />
GOTO 9999<br />
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />
C Write component information<br />
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />
600 CONTINUE<br />
KOMDSC = ’Syngas cleaning. The syntax is struc Componentname<br />
$ GASCLE_2 7 1 2 3 301 0 8 9 10 31 32 38 39. So the first<br />
$ number determines how many compounds is off−seperated. The<br />
$ 7 numbers from 8 to 39 is therefore compound numbers (8<br />
$ is $NH_3$).’<br />
K_PAR(1) = ’Pressure loss: $\\Delta p$ [bar]’<br />
DO I = 2,ANTST+1<br />
K_PAR(I) = ’Compound number’<br />
ENDDO<br />
K_BET = ’$\\dot{m}_1 \\gt 0 \\\\ \\dot{m}_2 \\lt 0<br />
$\\\\ \\dot{m}_3 \\lt 0 \\\\ \\dot{Q} \\lt 0$’<br />
KMEDDS(1) = ’Dirty gas in’<br />
KMEDDS(2) = ’Clean gas out’<br />
KMEDDS(3) = ’Off−separated gas’<br />
KMEDDS(4) = ’Heat loss’<br />
k_inp=’struc Cleaner GASCLE_2 7 13 14 15 308 0 8 9 10 31 32 38 39<br />
$\\\\media 13 GG 14 GG−clean 15 H2S\\\\<br />
$fluid GG H2 0.5 CO 0.3 CO2 0.1 H2O−G 0.05 CH4 0.02 H2S 0.03\\\\<br />
$addco Q Cleaner 308 0\\\\<br />
$addco t Cleaner 13 70 p 13 1 m Cleaner 13 1\\\\<br />
$START M Cleaner 14 −0.93 t Cleaner 14 70\\\\<br />
$START Y_J GG−clean H2 0.5 Y_J GG−clean CO 0.30\\\\<br />
$start Y_J GG−clean CO2 0.1 Y_J H2S H2S 1’<br />
C<br />
9999 CONTINUE<br />
RETURN<br />
END<br />
C=======================================================================<br />
C***********************************************************************<br />
SUBROUTINE STEAM_REFORMER(KOMTY,ANTLK,ANTEX,ANTKN,ANTPK,ANTM1,<br />
$ MMVAR,PARNAM,ZANAM,MEDIE, ANTME,VARME,ANTEL,VAREL,MDOT,P,H,Q<br />
$ ,PAR,ZA, RES,X_J,KOMDSC,KMEDDS,K_PAR,K_LIG ,K_BET,k_inp)<br />
C***********************************************************************<br />
C<br />
C STEREF_0 is a model of a steam reformer. The outlet composition is<br />
C based on an assumption of chemical equilibrium.<br />
C<br />
C***********************************************************************<br />
C<br />
15/67<br />
19−03−2007