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

VEnzin.for
c:/dna/source/
MEDIE(2) = ANYFLU$
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
C Pressure balance
C
c RES(1) = P(1) − P(2)
C
IF (FKOMP.EQ.5) GOTO 500
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Solution check
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
500 CONTINUE
IF (MDOT(1).LT.−1D−10) GOTO 450
IF (MDOT(2).GT.1D−10) GOTO 450
GOTO 9999
450 FBETI = .FALSE.
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Write component information
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
600 CONTINUE
KOMDSC = ’Utility component for seperating nodes which both have
$ assigned the pressure. Can also be used as a throttle valve.’
K_BET = ’$\\dot{m}_1 \\gt 0 \\\\ \\dot{m}_2 \\lt 0$’
KMEDDS(1) = ’Fluid in’
KMEDDS(2) = ’Fluid out’
k_INP=’struc dummy PRES_SEP 1 2\\\\
$media 1 STEAM\\\\
$addco m dummy 1 1 t dummy 1 50 p 1 1 p 2 1’
C
GOTO 9999
C
9999 CONTINUE
RETURN
END
C=======================================================================
C***********************************************************************
SUBROUTINE SET_X(KOMTY,ANTLK,ANTKN,ANTPK,ANTM1,MEDIE,ANTME
$ ,VARME,MDOT,P,H,ZC,PAR,RES,X_J,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.
52/67
19−03−2007