source: trunk/BioModel/separators/flash_w.mso @ 1008

Last change on this file since 1008 was 1008, checked in by Argimiro Resende Secchi, 21 months ago

Adding BioModel? to the MSO library.

File size: 5.1 KB
Line 
1#*-------------------------------------------------------------------
2* Biorrefinaria Petrobras
3*--------------------------------------------------------------------
4* Nome do arquivo: flash_w.mso
5* Projeto: Modelo integrado de producao de etanol 1G/2G
6* Conteudo: tanque flash somente para a agua
7*--------------------------------------------------------------------*#
8
9#*-------------------------------------------------------------------
10*
11* Versao 2.2
12* Data:    03/2016
13* Autores:   Anderson R. A. Lino e Felipe F. Furlan
14*
15*--------------------------------------------------------------------
16*Descricao: modelo do tanque flash que sera empregado na biorrefinaria
17*--------------------------------------------------------------------
18
19*--------------------------------------------------------------------
20*Hipoteses assumidas: 1 - mistura perfeita em ambas as fases
21*                     2 - equilibrio termodinamico
22*                     3 - modelo estatico
23*--------------------------------------------------------------------
24
25*--------------------------------------------------------------------
26*Notas: Foi feito o flowsheet teste para averiguar o modelo. Dependendo
27* das especificacoes, o flash pode funcionar como um evaporador tambem
28*--------------------------------------------------------------------*#
29
30using "water_stream";
31using "energy_stream";
32
33Model flash_wR
34
35        ATTRIBUTES
36        Pallete         = false;
37        Icon            = "icon/flash";
38        Brief           = "Flash with 2 Outlet Streams";
39        Info =
40"== GENERAL ==
41        Model of a steady-state flash tank.
42       
43== ASSUMPTIONS ==
44* Both phases are perfectly mixed;
45* No solid is dragged with in the vapour outlet stream.
46       
47== SPECIFY ==
48* The inlet stream:
49  flow rate
50  temperature
51  pressure
52  stream composition;
53* The outlet pressure (OutletV.P);
54* The outlet temperature OR the heat supplied.
55
56== SET ==
57* Basic composition (mass or molar).
58";
59
60#*-------------------------------------------------------------------
61#Parametros
62*--------------------------------------------------------------------*#
63
64        PARAMETERS
65        propterm as Plugin (Brief = "IAPWS 97 properties of water", Type = "water", Symbol = "{}");
66       
67#*-------------------------------------------------------------------
68* Declaracao de variaveis
69*--------------------------------------------------------------------*#
70
71VARIABLES
72
73in      Inlet   as water_stream;
74out OutletV     as water_stream;
75out OutletL     as water_stream;
76        Q               as power;
77        v               as fraction;
78        P               as pressure;
79        T               as temperature;
80       
81        Hl as enth_mass (Brief = "Saturated Water Enthalpy", Hidden = true);
82        Hv as enth_mass (Brief = "Saturated Steam Enthalpy", Hidden = true);
83        Sl as entr_mass (Brief = "Saturated Water Entropy", Hidden = true);
84        Sv as entr_mass (Brief = "Saturated Steam Entropy", Hidden = true);
85       
86#*-------------------------------------------------------------------
87* Equacoes do modelo
88*--------------------------------------------------------------------*#
89
90        EQUATIONS
91        "Global Mass Balance"
92        Inlet.Fw = OutletV.Fw + OutletL.Fw;
93       
94        "Energy Balance"
95        Inlet.Fw * Inlet.H + Q = OutletV.Fw * OutletV.H + OutletL.Fw * OutletL.H;
96       
97        "Mechanical Equilibrium 1"
98        OutletV.P = OutletL.P;
99       
100        "Thermal Equilibrium 1"
101        OutletV.T = OutletL.T;
102       
103        "Vapour Fraction"
104        OutletV.Fw = v * Inlet.Fw;
105
106        "Steam Mass Entropy and Enthalpy, array = [Sv, Hv]"
107        [Sv, Hv] = propterm.propPTv(P, T);
108       
109        "Water Mass Entropy and Enthalpy, array = [Sl, Hl]"
110        [Sl, Hl] = propterm.propPTl(P, T);
111       
112        "Liquid Stream Vapour Fraction"
113        OutletL.v = 0;
114       
115        "Steam Stream Vapour Fraction"
116        OutletV.v = 1;
117       
118        "Steam Mass Entropy and Enthalpy, array = [Sv, Hv]"
119        [OutletV.S, OutletV.H] = propterm.propPTv(OutletV.P, OutletV.T);
120       
121        "Water Mass Entropy and Enthalpy, array = [Sl, Hl]"
122        [OutletL.S, OutletL.H] = propterm.propPTl(OutletL.P, OutletL.T);
123       
124        "Saturation Condition"
125        T = propterm.Tsat(P);
126       
127        "Outlet Pressure"
128        OutletV.P = P;
129       
130        "Outlet Temperature"
131        OutletL.T = T;
132
133end
134
135FlowSheet teste_flash_w
136       
137        DEVICES
138        W101 as water_sourceR;
139        F101 as flash_wR;
140       
141#*-------------------------------------------------------------------
142* Especifica as conexoes entre os modelos
143*--------------------------------------------------------------------*#
144       
145        CONNECTIONS
146        W101.Outlet to F101.Inlet;
147       
148#*-------------------------------------------------------------------
149#Parametros
150*--------------------------------------------------------------------*#
151       
152        PARAMETERS
153        propterm as Plugin (Brief = "IAPWS 97 properties of water", Type = "water", Symbol = "{}");
154       
155#*-------------------------------------------------------------------
156* Define o valor dos parametros declarados no modelo
157*--------------------------------------------------------------------*#
158       
159        SET
160        W101.ValidPhases = "Liquid-Only";
161       
162#*-------------------------------------------------------------------
163* Especifica variaveis definidas no modelo
164*--------------------------------------------------------------------*#
165       
166        SPECIFY
167        W101.Fw = 100 * 'kg/h';
168        W101.P = 5 * 'atm';
169        W101.T = 420 * 'K';
170       
171        F101.P = 4.5 * 'atm';
172        F101.Q = 1000 * 'W';
173
174#*-------------------------------------------------------------------
175* Condicoes iniciais e opcoes de Solver
176*--------------------------------------------------------------------*#
177       
178        OPTIONS
179        Dynamic = false;
180
181end
Note: See TracBrowser for help on using the repository browser.