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

VEnzin.for
c:/dna/source/
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).GT.1D−10) GOTO 550
IF (Q.GT.1D−10) GOTO 550
IF (E.LT.−1D−10) GOTO 550
GOTO 9999
550 FBETI = .FALSE.
GOTO 9999
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
C Write component information
C−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−
600 CONTINUE
KOMDSC = ’Water elektrolysis with efficiency.’
KMEDDS(1) = ’Water in’
KMEDDS(2) = ’H2 out’
KMEDDS(3) = ’O2 out’
KMEDDS(4) = ’Power in’
KMEDDS(5) = ’Heat loss’
K_PAR(1) = ’Temperature: Temp [C]’
K_PAR(2) = ’Efficiency of elektrolysis: $\\eta$ [−−]’
K_LIG(1) = ’Only H2 in node 2: $y_j(H2,medie(2)) = 1$’
K_LIG(2) = ’Only O2 in node 3: $y_j(O2,medie(3)) = 1$’
K_LIG(3) = ’Mol balance:
$$−\\dot{m}_2 = \\dot{m}_1 M_{mol}(H2)/M_{mol}(H2O)$’
K_LIG(4) = ’Equal pressures: $p_1 = p_2$’
K_LIG(5) = ’Equal pressures: $p_1 = p_3$’
K_LIG(6) = ’Temperature of outlet gas: $Temp = T_2$’
K_LIG(7) = ’Equal temperature of outlet gases: $T_3 = T_2$’
K_LIG(8) = ’Efficiency of elektrolysis:
$$\\dot{E} \\eta=−\\dot{m}_2 LHV_{H2}$’
K_BET = ’$\\dot{m}_1\\gt 0 \\\\ \\dot{m}_2 \\lt 0 \\\\
$\\dot{m}_3 \\lt 0 \\\\ \\dot{Q} \\lt 0 \\\\ \\dot{E} \\gt 0$’
k_inp=’struc Elyse ELECTROLYSER 1 2 3 201 301 90 0.8\\\\
$media 2 H2 3 O2\\\\
$addco p 1 1 m Elyse 1 1 t Elyse 1 15\\\\
$start m Elyse 2 −1 m Elyse 3 −1\\\\
$start t Elyse 2 90 t Elyse 3 90\\\\
$start e Elyse 201 100\\\\
$start y_j H2 H2 1 y_j O2 O2 1\\\\
$start p 2 1 p 3 1 q Elyse 301 −10’
GOTO 9999
C
9999 CONTINUE
RETURN
END
C=======================================================================
C***********************************************************************
SUBROUTINE DRYER_04(KOMTY,ANTLK,ANTEX,ANTKN,ANTPK,ANTM1,MEDIE,
& ANTME,VARME,ANTEL,VAREL,MDOT,P,H,Q,PAR,RES,X_J,CP,HV,HF,ZA,
$ ZANAM,KOMDSC,K_PAR,K_lig,K_bet,KMEDDS,K_inp)
C***********************************************************************
C
C DRYER_02 is a model of a steam fuel dryer.
C The model does not include equations concerning the heat exchange.
C 1−2 is the heat emitting fluid.
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 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.
3/67
19−03−2007