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 The values mean: CA 99 : Water. CA 999 : Control. CA ANTME − OUTPUT − Number of fluids with variable composition. CA RES − OUTPUT − Residuals for the component. C CL CV Valve friction coefficient. CL V1 Specific volume for water entering the valve. CL T Temperature. CL S Entropy. CL X Quality. CL U Internal energy. CL K_PAR Parameter description. CL K_LIG Equation description. CL K_BET Condition description. CL K_MED Media description. C C Subroutines : STATES C COMINF C CP Programmer : Bent Lorentzen 1994 CP Lab. for Energetics, DTU, Denmark. C*********************************************************************** C C Include the common "environment" C INCLUDE ’ENVIRO.INI’ C C Parameter variables C INTEGER ANTLK, ANTEX, ANTKN, MEDIE(4), ANTPK, & ANTM1, ANTME, VARME(2) DOUBLE PRECISION MDOT(2), P(2), H(2), RES(3), ZC(2), PAR(2) CHARACTER*80 KOMTY C C Local variables C INTEGER K_MED(3) DOUBLE PRECISION MAXF,maxp,T,V,S,X,U CHARACTER*100 K_PAR(1),K_STAT(2) CHARACTER*500 K_LIG(2), K_BET CHARACTER*1000 KOMDSC,K_INP CHARACTER*100 KMEDDS(4) EXTERNAL COMINF,STATES C======================================================================= GOTO (100,200,1,400,400,200) FKOMP 1 RETURN C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− C Component name C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 100 CONTINUE KOMTY = ’SET_TEMP’ GOTO 9999 C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− C Component characteristics C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 200 CONTINUE KOMTY = ’SET_TEMP’ ANTKN = 3 ANTPK = 0 ANTLK = 2 ANTM1 = 2 MEDIE(1) = ANYFLU$ MEDIE(2) = ANYFLU$ MEDIE(3) = 999 ANTME = 2 VARME(1) = NODE1$ VARME(2) = NODE1$ IF (FKOMP.EQ.6) GOTO 600 ** FKOMP = 3 GOTO 9999 C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− C Component equations. All in residual form. C Do not include the conservation laws. These are treated automatically C by DNA. C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 400 CONTINUE C 61/67 19−03−2007

VEnzin.for c:/dna/source/ RES(1) = P(1) − P(2) C CALL STATES(P(1),H(1),T,V,S,X,U,1,2,MEDIE(1)) C RES(2) = ZC(1) − T C IF (FKOMP.EQ.5) GOTO 500 GOTO 9999 C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− C Solution check C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 500 CONTINUE IF (MDOT(1).LT.−1D−10) GOTO 550 IF (MDOT(2).GT.1D−10) GOTO 550 GOTO 9999 550 FBETI = .FALSE. GOTO 9999 C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− C Write component information C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−− 600 CONTINUE KOMDSC = ’Utillity component for setting or measuring the $ temperature by use of a control signal (ZC).’ K_LIG(1) = ’Equal pressures: $p_1 = p_2$’ K_LIG(2) = ’Temperature: T = ZC(1)’ K_BET = ’$\\dot{m_1} \\gt 0 \\\\ \\dot{m_2} \\lt 0 $’ KMEDDS(1) = ’Fluid in’ KMEDDS(2) = ’Fluid out’ KMEDDS(3) = ’Temperature’ K_INP=’struc set_temp SET_TEMP 1 2 900\\\\ $media 1 STEAM\\\\ $addco p 1 1 m set_temp 1 60 t set_temp 1 60’ C GOTO 9999 C 9999 CONTINUE RETURN END C*********************************************************************** SUBROUTINE SET_TEMP2(KOMTY,ANTLK,ANTKN,ANTPK,ANTM1,MEDIE,ANTME $ ,VARME,MDOT,P,H,ZC,PAR,RES,KOMDSC,K_PAR,K_lig,K_bet,KMEDDS $ ,K_inp) C*********************************************************************** C C SETFLOW1 is a model of a control valve. The valve controls massflow C using an error input from a controller. C C*********************************************************************** C CA FKOMP − INPUT − Flag with the value: CA 1: Initialize the component. CA 2: Initialize with actual system. CA 3: Fluid composition calculation (constant). CA 4: Find residuals. CA 5: Find residuals and check variables. CA 6: Output information about component. CA MDOT − INPUT − Massflows from nodes. CA P − INPUT − Pressure in nodes. CA H − INPUT − Enthalpy of massflows. CA ZC − INPUT − Control variables. CA PAR − INPUT − Parameters of the component. CA KOMTY − OUTPUT − Component name. CA ANTPK − OUTPUT − Number of parameters. CA ANTLK − OUTPUT − Number of equations. CA ANTEX − OUTPUT − Number of algebraic 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 of the connected nodes. CA The values mean: CA 99 : Water. CA 999 : Control. CA ANTME − OUTPUT − Number of fluids with variable composition. CA RES − OUTPUT − Residuals for the component. C CL CV Valve friction coefficient. CL V1 Specific volume for water entering the valve. CL T Temperature. 62/67 19−03−2007

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

c:/dna/source/<br />

CA The values mean:<br />

CA 99 : Water.<br />

CA 999 : Control.<br />

CA ANTME − OUTPUT − Number of fluids with variable composition.<br />

CA RES − OUTPUT − Residuals for the component.<br />

C<br />

CL CV Valve friction coefficient.<br />

CL V1 Specific volume for water entering the valve.<br />

CL T Temperature.<br />

CL S Entropy.<br />

CL X Quality.<br />

CL U Internal energy.<br />

CL K_PAR Parameter description.<br />

CL K_LIG Equation description.<br />

CL K_BET Condition description.<br />

CL K_MED Media description.<br />

C<br />

C Subroutines : STATES<br />

C COMINF<br />

C<br />

CP Pr<strong>og</strong>rammer : Bent Lorentzen 1994<br />

CP Lab. for Energetics, DTU, Denmark.<br />

C***********************************************************************<br />

C<br />

C Include the common "environment"<br />

C<br />

INCLUDE ’ENVIRO.INI’<br />

C<br />

C Parameter variables<br />

C<br />

INTEGER ANTLK, ANTEX, ANTKN, MEDIE(4), ANTPK,<br />

& ANTM1, ANTME, VARME(2)<br />

DOUBLE PRECISION MDOT(2), P(2), H(2), RES(3), ZC(2), PAR(2)<br />

CHARACTER*80 KOMTY<br />

C<br />

C Local variables<br />

C<br />

INTEGER K_MED(3)<br />

DOUBLE PRECISION MAXF,maxp,T,V,S,X,U<br />

CHARACTER*100 K_PAR(1),K_STAT(2)<br />

CHARACTER*500 K_LIG(2), K_BET<br />

CHARACTER*1000 KOMDSC,K_INP<br />

CHARACTER*100 KMEDDS(4)<br />

EXTERNAL COMINF,STATES<br />

C=======================================================================<br />

GOTO (100,200,1,400,400,200) FKOMP<br />

1 RETURN<br />

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

C Component name<br />

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

100 CONTINUE<br />

KOMTY = ’SET_TEMP’<br />

GOTO 9999<br />

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

C Component characteristics<br />

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

200 CONTINUE<br />

KOMTY = ’SET_TEMP’<br />

ANTKN = 3<br />

ANTPK = 0<br />

ANTLK = 2<br />

ANTM1 = 2<br />

MEDIE(1) = ANYFLU$<br />

MEDIE(2) = ANYFLU$<br />

MEDIE(3) = 999<br />

ANTME = 2<br />

VARME(1) = NODE1$<br />

VARME(2) = NODE1$<br />

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

** FKOMP = 3<br />

GOTO 9999<br />

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

C Component equations. All in residual form.<br />

C Do not include the conservation laws. These are treated automatically<br />

C by DNA.<br />

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

400 CONTINUE<br />

C<br />

61/67<br />

19−03−2007

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

Saved successfully!

Ooh no, something went wrong!