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

27.07.2013 Views

VEnzin.for c:/dna/source/ CA ANTPK − OUTPUT − Number of parameters for the component. CA ANTLK − OUTPUT − Number of equations in the component. CA ANTEX − OUTPUT − Number of independent equations in the component. CA ANTED − OUTPUT − Number of differential independent equations. CA ANTKN − OUTPUT − Number of nodes connected to the component. CA ANTM1 − OUTPUT − Number of massflows in the first conservation of CA mass equation. CA MEDIE − IN/OUT − Media (fluid) of the connected nodes. CA The values mean : CA 2 : Coal. CA 8 : Water (liquid). CA 27 : Oxygen rich gas. CA 25 : Raw gas. CA 300 : Heat. CA ANTME − OUTPUT − Number of fluids with variable composition. CA VARME − OUTPUT − Fluid numbers (with variable composition). CA ANTEL − OUTPUT − Number of compounds in these variable fluids. CA VAREL − OUTPUT − Compound numbers in variable fluids. CA RES − OUTPUT − Residuals for the component. C CL M4 Mass flow of raw gas. CL DELP Pressure drop through the plant. CL PGAS Gasifier pressure. CL TGAS Gasifier temperature. CL DMVC Amount of water relative to amount of coal. CL XRAW Composition of raw gas. CL K_PAR Parameter description. CL K_LIG Equation description. CL K_BET Condition description. CL K_MED Media description. C C Subroutines : COMINF C REAC C CP Programmer : Brian Elmegaard 2000 (Bent Lorentzen 1994, Niels Emsholm 1991) CP Dept. Energy Engr., DTU, Denmark. C*********************************************************************** C C Include the common "environment" C INCLUDE ’ENVIRO.INI’ INCLUDE ’THERPROP.DEC’ INCLUDE ’GASI.DEC’ C C Parameter variables C INTEGER ANTLK, ANTEX, ANTKN, MEDIE(7), ANTPK, : ANTM1, ANTME, VARME(4), ANTEL(4), : VAREL(ANTST,4),MMVAR(MAXMM) DOUBLE PRECISION X_J(MAXME,ANTST), PAR(20), RES(30), : MDOT(5),P(5),Q(1),CP(MAXME),HV(MAXNM),HF(MAXME), : H(5),ZA(7),ZC(1) CHARACTER*80 KOMTY CHARACTER*80 ZANAM(7),PARNAM(20) C C Local variables C INTEGER K_MED(6),I,J,CALCOM(ANTST) DOUBLE PRECISION DELP, PGAS, TGAS, DMVC, : T4,H5,V,S,X,U,CC,NIN(ANTST+1),NOUT(ANTST+1), : M_BL(5),G(15),R,METH,XEQ(ANTST) CHARACTER*100 K_PAR(5),K_STAT(1) CHARACTER*500 K_LIG(30), K_BET,K_GRAF $ ,warnstring CHARACTER*1000 KOMDSC,K_INP CHARACTER*100 KMEDDS(7) EXTERNAL COMINF,GIBBS,STATES INTRINSIC DABS,DLOG,NINT,EXP INCLUDE ’THERPROP.INI’ INCLUDE ’GASI.INI’ C======================================================================= GOTO (100,200,1,400,400,200,350,1,250) FKOMP 1 RETURN C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− C Component name C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 100 CONTINUE KOMTY = ’GASIFI_3_VENZIN’ 7/67 19−03−2007

VEnzin.for c:/dna/source/ MMVAR(1) = 15 GOTO 9999 C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− C Component characteristics C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 200 CONTINUE KOMTY = ’GASIFI_3_VENZIN’ ANTKN = 7 ANTPK = 5+MMVAR(1) ANTLK = 23 ANTEX = 7 ANTM1 = 5 ZANAM(1) = ’MULTIPLIER H’ ZANAM(2) = ’MULTIPLIER C’ ZANAM(3) = ’MULTIPLIER N’ ZANAM(4) = ’MULTIPLIER O’ ZANAM(5) = ’MULTIPLIER S’ ZANAM(6) = ’MULTIPL Ar’ ZANAM(7) = ’GIBBS ENERGY’ DO I=1,ANTM1 PARNAM(I) = ’CAL COMPOUND’ ENDDO PARNAM(ANTM1+1) = ’EQ PRESSURE’ PARNAM(ANTM1+2) = ’EQ TEMPERAT’ PARNAM(ANTM1+3) = ’PRESSURELOSS’ c PARNAM(ANTM1+4) = ’STEAM FLOW’ PARNAM(ANTM1+4) = ’UNCONV CARBO’ PARNAM(ANTM1+5) = ’METHANE PERC’ MEDIE(1) = −5 MEDIE(2) = 97 MEDIE(3) = −4 MEDIE(4) = −4 MEDIE(5) = −5 MEDIE(6) = 300 MEDIE(7) = 999 ANTME = 4 VARME(1) = −1 VARME(2) = −3 VARME(3) = −4 VARME(4) = −5 ANTEL(1) = 0 ANTEL(2) = 0 ANTEL(3) = 15 ANTEL(4) = 2 VAREL(1,3) = 1 VAREL(2,3) = 2 VAREL(3,3) = 3 VAREL(4,3) = 4 VAREL(5,3) = 5 VAREL(6,3) = 6 VAREL(7,3) = 7 VAREL(8,3) = 8 VAREL(9,3) = 9 VAREL(10,3) = 10 VAREL(11,3) = 11 VAREL(12,3) = 30 VAREL(13,3) = 31 VAREL(14,3) = 32 VAREL(15,3) = 36 VAREL(1,4) = 28 VAREL(2,4) = 38 IF (FKOMP.EQ.6) GOTO 600 *** FKOMP = 3 GOTO 9999 C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− C Initialization of algebraic variables. C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 250 CONTINUE ZA(7) = −.1D+1 C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− C Specific heat, heat of formation, heating value of ashes C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 350 CONTINUE CP(MEDIE(5)) = 1 HF(MEDIE(5)) = −5083.0D0 HV(MEDIE(5)) = X_J(MEDIE(5),28)*NED_H(28)/M_MOL(28) GOTO 9999 C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 8/67 19−03−2007

<strong>VEnzin</strong>.for<br />

c:/dna/source/<br />

MMVAR(1) = 15<br />

GOTO 9999<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

C Component characteristics<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

200 CONTINUE<br />

KOMTY = ’GASIFI_3_VENZIN’<br />

ANTKN = 7<br />

ANTPK = 5+MMVAR(1)<br />

ANTLK = 23<br />

ANTEX = 7<br />

ANTM1 = 5<br />

ZANAM(1) = ’MULTIPLIER H’<br />

ZANAM(2) = ’MULTIPLIER C’<br />

ZANAM(3) = ’MULTIPLIER N’<br />

ZANAM(4) = ’MULTIPLIER O’<br />

ZANAM(5) = ’MULTIPLIER S’<br />

ZANAM(6) = ’MULTIPL Ar’<br />

ZANAM(7) = ’GIBBS ENERGY’<br />

DO I=1,ANTM1<br />

PARNAM(I) = ’CAL COMPOUND’<br />

ENDDO<br />

PARNAM(ANTM1+1) = ’EQ PRESSURE’<br />

PARNAM(ANTM1+2) = ’EQ TEMPERAT’<br />

PARNAM(ANTM1+3) = ’PRESSURELOSS’<br />

c PARNAM(ANTM1+4) = ’STEAM FLOW’<br />

PARNAM(ANTM1+4) = ’UNCONV CARBO’<br />

PARNAM(ANTM1+5) = ’METHANE PERC’<br />

MEDIE(1) = −5<br />

MEDIE(2) = 97<br />

MEDIE(3) = −4<br />

MEDIE(4) = −4<br />

MEDIE(5) = −5<br />

MEDIE(6) = 300<br />

MEDIE(7) = 999<br />

ANTME = 4<br />

VARME(1) = −1<br />

VARME(2) = −3<br />

VARME(3) = −4<br />

VARME(4) = −5<br />

ANTEL(1) = 0<br />

ANTEL(2) = 0<br />

ANTEL(3) = 15<br />

ANTEL(4) = 2<br />

VAREL(1,3) = 1<br />

VAREL(2,3) = 2<br />

VAREL(3,3) = 3<br />

VAREL(4,3) = 4<br />

VAREL(5,3) = 5<br />

VAREL(6,3) = 6<br />

VAREL(7,3) = 7<br />

VAREL(8,3) = 8<br />

VAREL(9,3) = 9<br />

VAREL(10,3) = 10<br />

VAREL(11,3) = 11<br />

VAREL(12,3) = 30<br />

VAREL(13,3) = 31<br />

VAREL(14,3) = 32<br />

VAREL(15,3) = 36<br />

VAREL(1,4) = 28<br />

VAREL(2,4) = 38<br />

IF (FKOMP.EQ.6) GOTO 600<br />

*** FKOMP = 3<br />

GOTO 9999<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

C Initialization of algebraic variables.<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

250 CONTINUE<br />

ZA(7) = −.1D+1<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

C Specific heat, heat of formation, heating value of ashes<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

350 CONTINUE<br />

CP(MEDIE(5)) = 1<br />

HF(MEDIE(5)) = −5083.0D0<br />

HV(MEDIE(5)) = X_J(MEDIE(5),28)*NED_H(28)/M_MOL(28)<br />

GOTO 9999<br />

C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−<br />

8/67<br />

19−03−2007

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

Saved successfully!

Ooh no, something went wrong!