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

VEnzin.for
c:/dna/source/
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
DPA = PAR(1)
DPB = PAR(2)
C
C Pressure losses
C
RES(1) = P(1) − P(2) − DPA
RES(2) = P(3) − P(4) − DPB
C
RES(3) = ZA(1) − MDOT(3)/(MDOT(3)+MDOT(1))
RES(4) = ZA(2) − (MDOT(3)/RM_MOL(MEDIE(3)))/(MDOT(3)
$ /RM_MOL(MEDIE(3))+MDOT(1)/RM_MOL(MEDIE(1)))
RES(5) = ZA(2)−ZC(1)
C
RES(6) = H(1) − H(2)
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
IF (MDOT(3).LT.−1D−10) GOTO 550
IF (MDOT(4).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 the molar fraction of
$ fluid 2 (fluid 1 and 2 are considered as a mixture).’
K_PAR(1) = ’Pressure loss side 1, $\\Delta p_{12}$ [bar]’
K_PAR(2) = ’Pressure loss side 2, $\\Delta p_{34}$ [bar]’
K_BET = ’$\\dot{m}_1 \\gt 0 \\\\ \\dot{m}_2 \\lt 0 \\\\
$\\dot{m}_3 \\gt 0 \\\\ \\dot{m}_4 \\lt 0 $’
KMEDDS(1) = ’Fluid 1 inlet’
KMEDDS(2) = ’Fluid 1 outlet’
KMEDDS(3) = ’Fluid 2 inlet’
KMEDDS(4) = ’Fluid 2 outlet’
KMEDDS(5) = ’molar fraction of fluid 2’
K_INP= ’struc set−x SET_X_REALFLUID 691 693 692 694 908 0 0\\\\
$MEDIA 691 STEAM−HF 692 METHANOL\\\\
$addco p 691 1 p 692 1\\\\
$addco m set−x 691 1 t set−x 691 80 t set−x 692 80\\\\
$addco ZC 908 0.36\\\\
$start t set−x 693 60 m set−x 692 1’
C
GOTO 9999
C
9999 CONTINUE
RETURN
END
C=======================================================================
C***********************************************************************
SUBROUTINE MEASURE_FLOW(KOMTY,ANTLK,ANTKN,ANTPK,ANTM1,MEDIE,ANTME
$ ,VARME,MDOT,P,H,ZC,PAR,RES,KOMDSC,KMEDDS,K_LIG,K_BET,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).
56/67
19−03−2007