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

VEnzin.for
c:/dna/source/
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Write component information
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
600 CONTINUE
KOMDSC = ’Utillity component for converting the pressure to a
$ control signal (ZC). The component can therefore be used for
$ measuring or setting the pressure’
K_LIG(1) = ’Equal pressures: $p_1 = p_2$’
K_LIG(2) = ’Pressure: $p_1 = ZC(1)$’
K_BET = ’$\\dot{m_1} \\gt 0 \\\\ \\dot{m_2} \\lt 0 $’
KMEDDS(1) = ’Fluid in’
KMEDDS(2) = ’Fluid out’
KMEDDS(3) = ’Pressure’
K_INP=’struc set−pres SET_PRES 635 636 990\\\\
$MEDIA 635 STEAM\\\\
$addco t set−pres 635 50 m set−pres 635 1 p 635 1’
C
GOTO 9999
C
9999 CONTINUE
RETURN
END
C***********************************************************************
SUBROUTINE ENTHALPY(J,T,H)
C***********************************************************************
C
C GIBBS finds Gibbs energy at pressure PRES and tempera−
C ture T. The calculation is based on values for enthalpy and entropy.
C
C***********************************************************************
C
CA J − INPUT − For which compound should Gibbs energy be found.
CA TEMP − INPUT − Temperature used finding Gibbs energy.
CA PRES − INPUT − Pressure
CA G − OUTPUT − Gibbs energy.
C
CL T0 Reference temperature (25 C).
CL T Temperature [K].
CL H_MOL Enthalpy [J/mol].
CL S_MOL Entropy [J/(mol K)].
C
CP Programmer : Bent Lorentzen 1994
CP Lab. for Energetics, DTH, Denmark.
C***********************************************************************
C
C Include the common "environment" and thermodynamic properties
C
INCLUDE ’ENVIRO.INI’
CBE INCLUDE ’THERPROP.INI’
INCLUDE ’THERPROP.DEC’
C
C Parameter variables
C
INTEGER J
DOUBLE PRECISION H,T
C
C Local variables
C
DOUBLE PRECISION H_MOL,T0,T00
INTRINSIC DLOG
C
INCLUDE ’THERPROP.INI’
C=======================================================================
T0 = 298.15D0
T00 = 273.15D0
T=T+T00
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Find enthalpy using Cp−polynomia
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
H_MOL = AH(J) +
: T*(A(1,J)+T*(A(2,J)/2D0+T*(A(3,J)/3D0+T*(A(4,J)/4D0+T*
: (A(5,J)/5D0+T*(A(6,J)/6D0+T*A(7,J)/7D0)))))) −
: T0*(A(1,J)+T0*(A(2,J)/2D0+T0*(A(3,J)/3D0+T0*(A(4,J)/4D0+
: T0*(A(5,J)/5D0+T0*(A(6,J)/6D0+T0*A(7,J)/7D0))))))
H=H_MOL/m_mol(J)
C
66/67
19−03−2007