27.07.2013 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!