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

VEnzin.for
c:/dna/source/
C
DO I=2,ANTKN
if (I.ne.M_nr) then
RES(RES_nr) = H(1) − H(I)
RES_nr=RES_nr+1
endif
ENDDO
IF (FKOMP.EQ.5) GOTO 500
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Solution check
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
500 CONTINUE
C
IF (MDOT(1).LT.−1D−10) GOTO 550
DO I=1,ANTKN−1
IF (MDOT(I+1).GT.1D−10) GOTO 550
ENDDO
IF (abs(H(M_nr)−H(1)).GT.1D−1) GOTO 550
C
GOTO 9999
550 FBETI = .FALSE.
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Write component information
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
600 CONTINUE
KOMDSC = ’Flow splitter − variable number of outlets
$ for splitting up and later re−uniting the massflow − this
$ means that only the pressure in node 2 is set.This
$ documentation is for 3 outlets.’
K_LIG(1) = ’Equal pressures: $p_1 = p_2$’
K_LIG(2) = ’Equal enthalpy of outlets: $h_1 = h_2$’
K_LIG(3) = ’Equal enthalpy of outlets: $h_1 = h_3$’
K_BET = ’$\\dot{m}_1\\gt 0 \\\\ \\dot{m}_2\\lt 0
$\\\\ \\dot{m}_3\\lt 0 \\\\ \\dot{m}_4\\lt 0$’
KMEDDS(1) = ’Fluid in’
KMEDDS(2) = ’Fluid out’
KMEDDS(3) = ’Fluid out’
KMEDDS(4) = ’Fluid out’
k_inp= ’struc split splitter2 4 1 2 3 4\\\\
$media 1 SIMPLE_AIR\\\\
$addco m split 1 10 t split 1 50 p 1 1\\\\
$addco m split 2 −3 m split 3 −3\\\\
$addco p 3 1 p 4 1\\\\
$start t split 2 50 t split 3 50’
GOTO 9999
C
9999 CONTINUE
RETURN
END
C=======================================================================
C***********************************************************************
SUBROUTINE SPLITTER3(KOMTY,ANTLK,ANTKN,ANTPK,ANTM1,MEDIE,ANTME
$ ,VARME,PAR,MDOT,P,H,RES,komdsc,k_par,kmedds,k_lig,k_bet,k_inp
$ )
C***********************************************************************
C
C SPLITTER is splitting one mass flow in two. The two outlets have the
C same enthalpy.
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 KOMTY − OUTPUT − Component name.
CA ANTLK − OUTPUT − Number of 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.
45/67
19−03−2007