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

VEnzin.for
c:/dna/source/
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
IF (MDOT(5).LT.−1D−10) GOTO 550
IF (MDOT(6).GT.1D−10) GOTO 550
IF (MDOT(7).LT.−1D−10) GOTO 550
IF (MDOT(8).GT.1D−10) GOTO 550
GOTO 9999
550 FBETI = .FALSE.
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Write component information
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
600 CONTINUE
KOMDSC = ’Distillation stage for a distillation column. A column
$ can consist of several stages in series. The distillation
$ stage calculates VLE (vapor liquid equilibrium) between water
$ and methanol. The component can be used for distillation of
$ other media by changing 3 constants in the source code. The
$ method is NRTL (Non Random Two Liquid). Can only be used with
$ real fluids. The pressure is assigned in node 10 − not in the
$ seperate nodes’
K_BET =’$\\dot{m}_1 \\gt 0 \\\\ \\dot{m}_2 \\lt 0 \\\\
$\\dot{m}_3 \\gt 0 \\\\ \\dot{m}_4 \\lt 0 \\\\ \\dot{m}_5 \\gt 0 \\
$\\ \\dot{m}_6 \\lt 0 \\\\ \\dot{m}_7 \\gt 0 \\\\ \\dot{m}_8 \\lt 0
$$’
KMEDDS(1) = ’Liquid water inlet’
KMEDDS(2) = ’Liquid water outlet’
KMEDDS(3) = ’Water vapor inlet’
KMEDDS(4) = ’Water vapor outlet’
KMEDDS(5) = ’Liquid methanol inlet’
KMEDDS(6) = ’Liquid methanol outlet’
KMEDDS(7) = ’Methanol vapor inlet’
KMEDDS(8) = ’Methanol vapor outlet’
KMEDDS(9) = ’Heat’
KMEDDS(10) = ’Pressure in system’
KMEDDS(11) = ’Temperature of system’
KMEDDS(12) = ’Molar fraction of methanol in the liquid phase’
KMEDDS(13) = ’Molar fraction of methanol in the gas phase’
K_INP= ’struc Dis_stage_1 DISTILLATION_STAGE 701 707 705 703 702
$708 /\\\\ 706 704 343 906 907 908 909 0 0\\\\
$addco q Dis_stage_1 343 0 ZC 906 1\\\\
$addco tsat 701 80 x Dis_stage_1 701 0 m Dis_stage_1 701 1\\\\
$addco tsat 705 70 x Dis_stage_1 705 1 m Dis_stage_1 705 1\\\\
$addco tsat 702 80 x Dis_stage_1 702 0 m Dis_stage_1 702 1\\\\
$addco tsat 706 70 x Dis_stage_1 706 1 m Dis_stage_1 706 1\\\\
$start m Dis_stage_1 707 −1 m Dis_stage_1 708 −1\\\\
$start x Dis_stage_1 708 0 x Dis_stage_1 703 1 x Dis_stage_1 704 1
$\\\\START ZC 907 83’
C
GOTO 9999
C
9999 CONTINUE
RETURN
END
C=======================================================================
C***********************************************************************
SUBROUTINE EL_MOTOR(KOMTY,ANTLK,ANTKN,ANTPK,
& MEDIE,PAR,Q,E,WK,
& RES,komdsc,k_par,kmedds,k_lig,k_bet,k_inp)
C***********************************************************************
C
C EL−MOTOR is a component that converts electrical power to mechanical
C power and a heat loss. The efficiency is a parameter to the component.
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).
41/67
19−03−2007