Context Navigation

column.mso @ 775

Last change on this file since 775 was 775, checked in by gerson bicca, 13 years ago
changes in column model
File size: 109.8 KB

Line
1	#*-------------------------------------------------------------------
2	* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
3	*
4	* This LIBRARY is free software; you can distribute it and/or modify
5	* it under the therms of the ALSOC FREE LICENSE as available at
6	* http://www.enq.ufrgs.br/alsoc.
7	*
8	* EMSO Copyright (C) 2004 - 2007 ALSOC, original code
9	* from http://www.rps.eng.br Copyright (C) 2002-2004.
10	* All rights reserved.
11	*
12	* EMSO is distributed under the therms of the ALSOC LICENSE as
13	* available at http://www.enq.ufrgs.br/alsoc.
14	*
15	*----------------------------------------------------------------------
16	* File containg models of columns: distillation, stripping, absorbers
17	* rectifier, ....
18	*
19	* The default nomenclature is:
20	* Type_Column_reboilertype_condensertyper
21	*
22	* where:
23	* Type = refluxed or reboiled or section
24	* Column = Stripping, Absorption, Rectifier, Distillation
25	* Reboiler type (if exists) = kettle or thermosyphon
26	* Condenser type (if exists) = with subccoling or without subcooling
27	*
28	*-----------------------------------------------------------------------
29	* Author: Paula B. Staudt
30	* $Id: column.mso 511 2008-05-12 17:25:33Z paula $
31	---------------------------------------------------------------------#
32
33	using "tray";
34	using "reboiler";
35	using "condenser";
36	using "mixers_splitters/splitter";
37	using "tank";
38	using "pressure_changers/pump";
39
40	Model ControlSection
41
42	ATTRIBUTES
43	Pallete = false;
44	Brief = "Control variables for column section.";
45
46	PARAMETERS
47	outer NComp as Integer (Brief="Number of components");
48
49	Tindicator_TrayNumber as Integer (Brief="Temperature Indicator in the tray");
50	Pindicator_TrayNumber as Integer (Brief="Pressure Indicator in the tray");
51
52	end
53
54	Model InitializeSection
55
56	ATTRIBUTES
57	Pallete = false;
58	Brief = "Initial conditions for column section.";
59
60	PARAMETERS
61	outer NComp as Integer (Brief="Number of components");
62
63	TopTemperature as temperature (Brief = "Tray Temperature at Column Top", Default = 300);
64	BottomTemperature as temperature (Brief = "Tray Temperature at Column Bottom", Default = 330);
65	LevelFraction as fraction (Brief = "Tray Level Fraction", Default = 0.5);
66
67	TopComposition(NComp) as fraction (Brief = "Component Molar Fraction at Column Top", Default = 0.30);
68	BottomComposition(NComp) as fraction (Brief = "Component Molar Fraction at Column Bottom", Default = 0.30);
69
70	end
71
72	Model InitializeStage
73
74	ATTRIBUTES
75	Pallete = false;
76	Brief = "Initial conditions for Packed column section.";
77
78	PARAMETERS
79	outer NComp as Integer (Brief="Number of components");
80
81	TopStageTemperature as temperature (Brief = "Tray Temperature at Column Top", Default = 300);
82	BottomStageTemperature as temperature (Brief = "Tray Temperature at Column Bottom", Default = 300);
83
84	TopStageComposition(NComp) as fraction (Brief = "Component Molar Fraction at Top");
85	BottomStageComposition(NComp) as fraction (Brief = "Component Molar Fraction at Bottom");
86	LiquidMolarHoldup as mol (Brief="Molar liquid holdup", Default=0.01);
87
88	end
89
90	Model Section_ColumnBasic
91
92	ATTRIBUTES
93	Pallete = false;
94	Brief = "Model of a basic column section.";
95	Info =
96	"Model of a basic column section containing a vetor of TRAYS numbered from the top-down.";
97
98	PARAMETERS
99	outer PP as Plugin (Brief="External Physical Properties", Type="PP");
100	outer NComp as Integer (Brief="Number of components");
101	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
102	# Section Column Feed Tray - Side Streams Location and Numbering
103	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
104	NumberOfTrays as Integer (Brief="Number of trays", Default=8);
105	FeedTrayIndex(NumberOfTrays) as Integer (Brief="Feed Tray Index", Default=0,Hidden=true);
106	LiqSideTrayIndex(NumberOfTrays) as Integer (Brief="Liquid Side Tray Index", Default=0,Hidden=true);
107	VapSideTrayIndex(NumberOfTrays) as Integer (Brief="Vapour Side Tray Index", Default=0,Hidden=true);
108	FeedTrayLocation as Integer (Brief="Feed tray Location", Default=2);
109	LiquidSideStreamLocation as Integer (Brief="Liquid Side Stream Location", Default=2);
110	VapourSideStreamLocation as Integer (Brief="Vapour Side Stream Location", Default=2);
111	g as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
112	Mw(NComp) as molweight (Brief="Component Mol Weight",Hidden=true);
113	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
114	# Section Column Flow Model for Liquid and Vapour
115	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
116	VapourFlowModel as Switcher (Valid = ["Reepmeyer", "Feehery_Fv", "Roffel_Fv", "Klingberg", "Wang_Fv", "Elgue"], Default = "Reepmeyer");
117	LiquidFlowModel as Switcher (Valid = ["default", "Wang_Fl", "Olsen", "Feehery_Fl", "Roffel_Fl"], Default = "default");
118	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
119	# Section Column Tray Geometry and Auxiliar Parameters
120	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
121	VolumeOfTray as volume (Brief="Total Volume of the tray");
122	HeatSupply as heat_rate (Brief="Rate of heat supply");
123	PlateArea as area (Brief="Plate area = Atray - Adowncomer");
124	HolesArea as area (Brief="Total holes area");
125	WeirLength as length (Brief="Weir length");
126	WeirHeight as length (Brief="Weir height");
127	FeeheryCoeff as Real (Brief="Feeherys correlation coefficient", Unit='1/m^4', Default=1);
128	ElgueCoeff as Real (Brief="Elgues correlation coefficient", Unit='kg/m/mol^2', Default=1);
129	OlsenCoeff as Real (Brief="Olsens correlation coefficient", Default=1);
130	TrayLiquidPasses as Real (Brief="Number of liquid passes in the tray", Default=1);
131
132	V as volume (Brief="Total Volume of the tray",Hidden=true);
133	Q as heat_rate (Brief="Rate of heat supply",Hidden=true);
134	Ap as area (Brief="Plate area = Atray - Adowncomer",Hidden=true);
135	Ah as area (Brief="Total holes area",Hidden=true);
136	lw as length (Brief="Weir length",Hidden=true);
137	hw as length (Brief="Weir height",Hidden=true);
138	beta as fraction (Brief="Aeration fraction");
139	alfa as fraction (Brief="Dry pressure drop coefficient");
140	w as Real (Brief="Feeherys correlation coefficient", Unit='1/m^4', Default=1,Hidden=true);
141	btray as Real (Brief="Elgues correlation coefficient", Unit='kg/m/mol^2', Default=1,Hidden=true);
142	fw as Real (Brief="Olsens correlation coefficient", Default=1,Hidden=true);
143	Np as Real (Brief="Number of liquid passes in the tray", Default=1,Hidden=true);
144
145	VapourFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true);
146	LiquidFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true);
147
148	SET
149	FeedTrayIndex(FeedTrayLocation) =1;
150	VapSideTrayIndex(FeedTrayLocation) =1;
151	LiqSideTrayIndex(FeedTrayLocation) =1;
152	Mw = PP.MolecularWeight();
153
154	V=VolumeOfTray;
155	Q=HeatSupply;
156	Ap=PlateArea;
157	Ah=HolesArea;
158	lw=WeirLength;
159	hw=WeirHeight ;
160	w=FeeheryCoeff;
161	btray=ElgueCoeff;
162	fw=OlsenCoeff;
163	Np=TrayLiquidPasses;
164
165	VARIABLES
166	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
167	# Tray Initialization For the Whole Column Section
168	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
169	INITIALIZATION as InitializeSection (Brief = "Column Model Initialization");
170	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
171	# Control Ports
172	#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
173	CONTROL as ControlSection (Brief = "Control");
174	out TCI as control_signal (Brief="Temperature Indicator", Protected = true, PosX=1, PosY=0.90);
175	out PCI as control_signal (Brief="Pressure Indicator", Protected = true, PosX=1, PosY=0.80);
176	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
177	# Column Feed - Side Stream Flow and Murphree Efficiency
178	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
179	in FeedTray as stream (Brief="Feed stream", PosX=0, PosY=0.55);
180	TRAYS(NumberOfTrays) as tray (Brief="Number of trays in the Column Section");
181	VapourDrawOffFlow as flow_mol (Brief="Stream Molar Flow Rate");
182	LiquidDrawOffFlow as flow_mol (Brief="Stream Molar Flow Rate");
183	MurphreeEff as Real (Brief="Murphree efficiency for All Trays");
184
185	CONNECTIONS
186	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
187	# Connecting Intermediate Trays
188	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
189	TRAYS([2:NumberOfTrays]).OutletVapour to TRAYS([1:NumberOfTrays-1]).InletVapour;
190	TRAYS([1:NumberOfTrays-1]).OutletLiquid to TRAYS([2:NumberOfTrays]).InletLiquid;
191
192	INITIAL
193	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
194	# Tray Initialization
195	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
196	for i in 1:NumberOfTrays do
197
198	"The initial temperature of the TRAYS"
199	TRAYS(i).OutletLiquid.T = INITIALIZATION.TopTemperature+(INITIALIZATION.BottomTemperature-INITIALIZATION.TopTemperature)*((i-1)/(NumberOfTrays-1));
200
201	"The initial Level of the TRAYS"
202	TRAYS(i).Level = INITIALIZATION.LevelFraction*hw;
203	end
204
205	for i in 1:NComp-1 do
206
207	for j in 1:NumberOfTrays do
208
209	"The initial composition of the TRAYS - Normalized"
210	TRAYS(j).OutletLiquid.z(i) = INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) +(INITIALIZATION.BottomComposition(i)/sum(INITIALIZATION.BottomComposition)-INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) )*((j-1)/(NumberOfTrays-1));
211	end
212
213	end
214
215	EQUATIONS
216	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
217	# Control Ports
218	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
219
220	"Tray Temperature Indicator"
221	#TCI*'K' = TRAYS(min([NumberOfTrays, CONTROL.Tindicator_TrayNumber])).OutletVapour.T;
222	TCI*'K' = TRAYS(CONTROL.Tindicator_TrayNumber).OutletVapour.T;
223
224	"Tray Pressure Indicator"
225	#PCI*'atm' = TRAYS(min([NumberOfTrays, CONTROL.Pindicator_TrayNumber])).OutletVapour.P;
226	PCI*'atm' = TRAYS(CONTROL.Pindicator_TrayNumber).OutletVapour.P;
227
228	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
229	# Equating Feed Tray Variables to Trays Variables
230	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
231	"Feed Tray Inlet Flow"
232	FeedTray.F*FeedTrayIndex= TRAYS.Inlet.F;
233
234	"Feed Tray Inlet Temperature"
235	FeedTray.T = TRAYS.Inlet.T;
236
237	"Feed Tray Inlet Pressure"
238	FeedTray.P = TRAYS.Inlet.P;
239
240	"Feed Tray Inlet Composition"
241	FeedTray.z = TRAYS.Inlet.z;
242
243	"Feed Tray Inlet Vapour Fraction"
244	FeedTray.v = TRAYS.Inlet.v;
245
246	"Feed Tray Inlet Enthalpy"
247	FeedTray.h = TRAYS.Inlet.h;
248
249	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
250	# Trays Equations
251	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
252
253	for i in [1:NumberOfTrays] do
254
255	"Murphree Efficiency"
256	TRAYS(i).OutletVapour.z = MurphreeEff * (TRAYS(i).yideal - TRAYS(i).InletVapour.z) + TRAYS(i).InletVapour.z;
257
258	"Level of clear liquid over the weir"
259	TRAYS(i).Level = TRAYS(i).ML*TRAYS(i).vL/Ap;
260
261	"Geometry Constraint"
262	V = TRAYS(i).ML* TRAYS(i).vL + TRAYS(i).MV*TRAYS(i).vV;
263
264	"Energy Holdup"
265	TRAYS(i).E = TRAYS(i).MLTRAYS(i).OutletLiquid.h + TRAYS(i).MVTRAYS(i).OutletVapour.h - TRAYS(i).OutletLiquid.P*V;
266
267	"Energy Balance"
268	diff(TRAYS(i).E) = ( TRAYS(i).Inlet.FTRAYS(i).Inlet.h + TRAYS(i).InletLiquid.FTRAYS(i).InletLiquid.h + TRAYS(i).InletVapour.FTRAYS(i).InletVapour.h- TRAYS(i).OutletLiquid.FTRAYS(i).OutletLiquid.h - TRAYS(i).OutletVapour.F*TRAYS(i).OutletVapour.h
269	-TRAYS(i).VapourSideStream.FTRAYS(i).VapourSideStream.h - TRAYS(i).LiquidSideStream.FTRAYS(i).LiquidSideStream.h + Q );
270
271	switch LiquidFlow
272	case "on":
273	switch LiquidFlowModel
274	case "default":
275	"Francis Equation"
276	TRAYS(i).OutletLiquid.FTRAYS(i).vL = 1.84'1/s'lw((TRAYS(i).Level-(beta*hw))/(beta))^2;
277
278	case "Wang_Fl":
279	TRAYS(i).OutletLiquid.FTRAYS(i).vL = 1.84'm^0.5/s'lw((TRAYS(i).Level-(beta*hw))/(beta))^1.5;
280
281	case "Olsen":
282	TRAYS(i).OutletLiquid.F / 'mol/s'= lwNpTRAYS(i).rhoL/sum(MwTRAYS(i).OutletVapour.z)/(0.665fw)^1.5 * ((TRAYS(i).MLsum(MwTRAYS(i).OutletLiquid.z)/TRAYS(i).rhoL/Ap)-hw)^1.5 * 'm^0.5/mol';
283
284	case "Feehery_Fl":
285	TRAYS(i).OutletLiquid.F = lwTRAYS(i).rhoL/sum(MwTRAYS(i).OutletLiquid.z) * ((TRAYS(i).Level-hw)/750/'mm')^1.5 * 'm^2/s';
286
287	case "Roffel_Fl":
288	TRAYS(i).OutletLiquid.F = 2/3TRAYS(i).rhoL/sum(MwTRAYS(i).OutletLiquid.z)lw(TRAYS(i).MLsum(MwTRAYS(i).OutletLiquid.z)/(Ap1.3)/TRAYS(i).rhoL)^1.5sqrt(2*g/
289	(2(1 - 0.3593/'Pa^0.0888545'abs(TRAYS(i).OutletVapour.Fsum(MwTRAYS(i).OutletVapour.z)/(Ap*1.3)/sqrt(TRAYS(i).rhoV))^0.177709)-1)); #/'(kg/m)^0.0888545/s^0.177709';
290	end
291	when TRAYS(i).Level < (beta *hw) switchto "off";
292
293	case "off":
294	"Low level"
295	TRAYS(i).OutletLiquid.F = 0 * 'mol/h';
296	when TRAYS(i).Level > (beta * hw) + 1e-6*'m' switchto "on";
297	end
298
299	switch VapourFlow
300	case "on":
301	switch VapourFlowModel
302	case "Reepmeyer":
303	TRAYS(i).InletVapour.FTRAYS(i).vV = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)/(TRAYS(i).rhoValfa))*Ah;
304
305	case "Feehery_Fv":
306	TRAYS(i).InletVapour.F = TRAYS(i).rhoV/Ap/w/sum(MwTRAYS(i).OutletVapour.z) sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-(TRAYS(i).rhoVgTRAYS(i).ML*TRAYS(i).vL/Ap))/TRAYS(i).rhoV);
307
308	case "Roffel_Fv":
309	TRAYS(i).InletVapour.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.921e5' = (TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)1e5 - (betasum(TRAYS(i).MMw)/(Ap1.3)g1e5) (TRAYS(i).rhoVAh/sum(MwTRAYS(i).OutletVapour.z))^1.08 * 'm^1.08/mol^1.08';
310
311	case "Klingberg":
312	TRAYS(i).InletVapour.F * TRAYS(i).vV = Ap * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoLgTRAYS(i).Level)/TRAYS(i).rhoV);
313
314	case "Wang_Fv":
315	TRAYS(i).InletVapour.F * TRAYS(i).vV = Ap * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoLgTRAYS(i).Level)/TRAYS(i).rhoV*alfa);
316
317	case "Elgue":
318	TRAYS(i).InletVapour.F = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)/btray);
319	end
320	when TRAYS(i).InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
321
322	case "off":
323	TRAYS(i).InletVapour.F = 0 * 'mol/s';
324	when TRAYS(i).InletVapour.P > TRAYS(i).OutletVapour.P + TRAYS(i).LevelgTRAYS(i).rhoL + 1e-1 * 'atm' switchto "on";
325	end
326
327	end
328
329	end
330
331
332
333
334
335	Model Section_Column
336
337	ATTRIBUTES
338	Pallete = true;
339	Icon = "icon/SectionColumn";
340	Brief = "Model of a column section.";
341	Info =
342	"== Model of a column section containing ==
343	* NumberOfTrays TRAYS.
344
345	== Specify ==
346	* the feed stream of each tray (Inlet);
347	* the Murphree eficiency for each tray Emv;
348	* the InletLiquid stream of the top tray;
349	* the InletVapour stream of the bottom tray.
350
351	== Initial Conditions ==
352	* the TRAYS temperature (OutletLiquid.T);
353	* the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F);
354	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray.
355	";
356
357	PARAMETERS
358	outer PP as Plugin (Brief="External Physical Properties", Type="PP");
359	outer NComp as Integer (Brief="Number of components");
360	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
361	# Section Column Feed Tray - Side Streams Location and Numbering
362	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
363	NumberOfTrays as Integer (Brief="Number of trays", Default=8);
364	FeedTrayIndex(NumberOfTrays) as Integer (Brief="Feed Tray Index", Default=0,Hidden=true);
365	LiqSideTrayIndex(NumberOfTrays) as Integer (Brief="Liquid Side Tray Index", Default=0,Hidden=true);
366	VapSideTrayIndex(NumberOfTrays) as Integer (Brief="Vapour Side Tray Index", Default=0,Hidden=true);
367	FeedTrayLocation as Integer (Brief="Feed tray Location", Default=2);
368	LiquidSideStreamLocation as Integer (Brief="Liquid Side Stream Location", Default=2);
369	VapourSideStreamLocation as Integer (Brief="Vapour Side Stream Location", Default=2);
370	g as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
371	Mw(NComp) as molweight (Brief="Component Mol Weight",Hidden=true);
372	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
373	# Section Column Flow Model for Liquid and Vapour
374	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
375	VapourFlowModel as Switcher (Valid = ["Reepmeyer", "Feehery_Fv", "Roffel_Fv", "Klingberg", "Wang_Fv", "Elgue"], Default = "Reepmeyer");
376	LiquidFlowModel as Switcher (Valid = ["default", "Wang_Fl", "Olsen", "Feehery_Fl", "Roffel_Fl"], Default = "default");
377	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
378	# Section Column Tray Geometry and Auxiliar Parameters
379	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
380	VolumeOfTray as volume (Brief="Total Volume of the tray");
381	HeatSupply as heat_rate (Brief="Rate of heat supply");
382	PlateArea as area (Brief="Plate area = Atray - Adowncomer");
383	HolesArea as area (Brief="Total holes area");
384	WeirLength as length (Brief="Weir length");
385	WeirHeight as length (Brief="Weir height");
386	FeeheryCoeff as Real (Brief="Feeherys correlation coefficient", Unit='1/m^4', Default=1);
387	ElgueCoeff as Real (Brief="Elgues correlation coefficient", Unit='kg/m/mol^2', Default=1);
388	OlsenCoeff as Real (Brief="Olsens correlation coefficient", Default=1);
389	TrayLiquidPasses as Real (Brief="Number of liquid passes in the tray", Default=1);
390
391	V as volume (Brief="Total Volume of the tray",Hidden=true);
392	Q as heat_rate (Brief="Rate of heat supply",Hidden=true);
393	Ap as area (Brief="Plate area = Atray - Adowncomer",Hidden=true);
394	Ah as area (Brief="Total holes area",Hidden=true);
395	lw as length (Brief="Weir length",Hidden=true);
396	hw as length (Brief="Weir height",Hidden=true);
397	beta as fraction (Brief="Aeration fraction");
398	alfa as fraction (Brief="Dry pressure drop coefficient");
399	w as Real (Brief="Feeherys correlation coefficient", Unit='1/m^4', Default=1,Hidden=true);
400	btray as Real (Brief="Elgues correlation coefficient", Unit='kg/m/mol^2', Default=1,Hidden=true);
401	fw as Real (Brief="Olsens correlation coefficient", Default=1,Hidden=true);
402	Np as Real (Brief="Number of liquid passes in the tray", Default=1,Hidden=true);
403
404	VapourFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true);
405	LiquidFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true);
406
407	SET
408	FeedTrayIndex(FeedTrayLocation) =1;
409	VapSideTrayIndex(FeedTrayLocation) =1;
410	LiqSideTrayIndex(FeedTrayLocation) =1;
411	Mw = PP.MolecularWeight();
412
413	V=VolumeOfTray;
414	Q=HeatSupply;
415	Ap=PlateArea;
416	Ah=HolesArea;
417	lw=WeirLength;
418	hw=WeirHeight ;
419	w=FeeheryCoeff;
420	btray=ElgueCoeff;
421	fw=OlsenCoeff;
422	Np=TrayLiquidPasses;
423
424	VARIABLES
425	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
426	# Tray Initialization For the Whole Column Section
427	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
428	INITIALIZATION as InitializeSection (Brief = "Column Model Initialization");
429	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
430	# Control Ports
431	#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
432	CONTROL as ControlSection (Brief = "Control");
433	out TCI as control_signal (Brief="Temperature Indicator", Protected = true, PosX=1, PosY=0.90);
434	out PCI as control_signal (Brief="Pressure Indicator", Protected = true, PosX=1, PosY=0.80);
435	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
436	# Section Inlet Liquid and Outlet Vapour Material Ports
437	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
438	in LiquidInlet as stream (Brief="Liquid Inlet in the section", PosX=0.70, PosY=0);
439	out VapourOutlet as vapour_stream (Brief="Vapour Outlet in the section", PosX=0.30, PosY=0);
440
441	out VapourDrawOff as vapour_stream (Brief="Vapour Draw Off Port in the section", PosX=1, PosY=0.35,Protected = true);
442	out LiquidDrawOff as liquid_stream (Brief="Liquid Draw Off Port in the section", PosX=1, PosY=0.65,Protected = true);
443	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
444	# Column Feed - Side Stream Flow and Murphree Efficiency
445	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
446	in FeedTray as stream (Brief="Feed stream", PosX=0, PosY=0.55);
447	TRAYS(NumberOfTrays) as tray (Brief="Number of trays in the Column Section");
448	VapourDrawOffFlow as flow_mol (Brief="Stream Molar Flow Rate");
449	LiquidDrawOffFlow as flow_mol (Brief="Stream Molar Flow Rate");
450	MurphreeEff as positive (Brief="Murphree efficiency for All Trays");
451	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
452	# Section Inlet Vapour and Outlet Liquid Material Ports
453	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
454	in VapourInlet as stream (Brief="Vapour Inlet in the section", PosX=0.30, PosY=1);
455	out LiquidOutlet as liquid_stream (Brief="Liquid Outlet in the section", PosX=0.70, PosY=1);
456
457	LiquidConnector as stream (Brief="Liquid connection in the middle TRAYS", PosX=0.75, PosY=1,Hidden=true);
458	VapourConnector as stream (Brief="Vapour connection in the middle TRAYS", PosX=0.55, PosY=0,Hidden=true);
459
460	CONNECTIONS
461	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
462	# Connecting Intermediate Trays
463	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
464	TRAYS([2:NumberOfTrays]).OutletVapour to TRAYS([1:NumberOfTrays-1]).InletVapour;
465	TRAYS([1:NumberOfTrays-1]).OutletLiquid to TRAYS([2:NumberOfTrays]).InletLiquid;
466
467	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
468	# Connecting Connectors
469	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
470	LiquidConnector to TRAYS(1).InletLiquid;
471	VapourConnector to TRAYS(NumberOfTrays).InletVapour;
472
473	INITIAL
474	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
475	# Tray Initialization
476	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
477	for i in 1:NumberOfTrays do
478
479	"The initial temperature of the TRAYS"
480	TRAYS(i).OutletLiquid.T = INITIALIZATION.TopTemperature+(INITIALIZATION.BottomTemperature-INITIALIZATION.TopTemperature)*((i-1)/(NumberOfTrays-1));
481
482	"The initial Level of the TRAYS"
483	TRAYS(i).Level = INITIALIZATION.LevelFraction*hw;
484	end
485
486	for i in 1:NComp-1 do
487
488	for j in 1:NumberOfTrays do
489
490	"The initial composition of the TRAYS - Normalized"
491	TRAYS(j).OutletLiquid.z(i) = INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) +(INITIALIZATION.BottomComposition(i)/sum(INITIALIZATION.BottomComposition)-INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) )*((j-1)/(NumberOfTrays-1));
492	end
493
494	end
495
496	EQUATIONS
497	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
498	# Control Ports
499	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
500
501	"Tray Temperature Indicator"
502	#TCI*'K' = TRAYS(min([NumberOfTrays, CONTROL.Tindicator_TrayNumber])).OutletVapour.T;
503	TCI*'K' = TRAYS(CONTROL.Tindicator_TrayNumber).OutletVapour.T;
504
505	"Tray Pressure Indicator"
506	#PCI*'atm' = TRAYS(min([NumberOfTrays, CONTROL.Pindicator_TrayNumber])).OutletVapour.P;
507	PCI*'atm' = TRAYS(CONTROL.Pindicator_TrayNumber).OutletVapour.P;
508
509	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
510	# Equating Feed Tray Variables to Trays Variables
511	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
512	"Feed Tray Inlet Flow"
513	FeedTray.F*FeedTrayIndex= TRAYS.Inlet.F;
514
515	"Feed Tray Inlet Temperature"
516	FeedTray.T = TRAYS.Inlet.T;
517
518	"Feed Tray Inlet Pressure"
519	FeedTray.P = TRAYS.Inlet.P;
520
521	"Feed Tray Inlet Composition"
522	FeedTray.z = TRAYS.Inlet.z;
523
524	"Feed Tray Inlet Vapour Fraction"
525	FeedTray.v = TRAYS.Inlet.v;
526
527	"Feed Tray Inlet Enthalpy"
528	FeedTray.h = TRAYS.Inlet.h;
529
530	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
531	# Trays Equations
532	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
533	for i in [1:NumberOfTrays] do
534
535	"Murphree Efficiency"
536	TRAYS(i).OutletVapour.z = MurphreeEff * (TRAYS(i).yideal - TRAYS(i).InletVapour.z) + TRAYS(i).InletVapour.z;
537
538	"Level of clear liquid over the weir"
539	TRAYS(i).Level = TRAYS(i).ML*TRAYS(i).vL/Ap;
540
541	"Geometry Constraint"
542	V = TRAYS(i).ML* TRAYS(i).vL + TRAYS(i).MV*TRAYS(i).vV;
543
544	"Energy Holdup"
545	TRAYS(i).E = TRAYS(i).MLTRAYS(i).OutletLiquid.h + TRAYS(i).MVTRAYS(i).OutletVapour.h - TRAYS(i).OutletLiquid.P*V;
546
547	"Energy Balance"
548	diff(TRAYS(i).E) = ( TRAYS(i).Inlet.FTRAYS(i).Inlet.h + TRAYS(i).InletLiquid.FTRAYS(i).InletLiquid.h + TRAYS(i).InletVapour.FTRAYS(i).InletVapour.h- TRAYS(i).OutletLiquid.FTRAYS(i).OutletLiquid.h - TRAYS(i).OutletVapour.F*TRAYS(i).OutletVapour.h
549	-TRAYS(i).VapourSideStream.FTRAYS(i).VapourSideStream.h - TRAYS(i).LiquidSideStream.FTRAYS(i).LiquidSideStream.h + Q );
550
551	switch LiquidFlow
552	case "on":
553	switch LiquidFlowModel
554	case "default":
555	"Francis Equation"
556	TRAYS(i).OutletLiquid.FTRAYS(i).vL = 1.84'1/s'lw((TRAYS(i).Level-(beta*hw))/(beta))^2;
557
558	case "Wang_Fl":
559	TRAYS(i).OutletLiquid.FTRAYS(i).vL = 1.84'm^0.5/s'lw((TRAYS(i).Level-(beta*hw))/(beta))^1.5;
560
561	case "Olsen":
562	TRAYS(i).OutletLiquid.F / 'mol/s'= lwNpTRAYS(i).rhoL/sum(MwTRAYS(i).OutletVapour.z)/(0.665fw)^1.5 * ((TRAYS(i).MLsum(MwTRAYS(i).OutletLiquid.z)/TRAYS(i).rhoL/Ap)-hw)^1.5 * 'm^0.5/mol';
563
564	case "Feehery_Fl":
565	TRAYS(i).OutletLiquid.F = lwTRAYS(i).rhoL/sum(MwTRAYS(i).OutletLiquid.z) * ((TRAYS(i).Level-hw)/750/'mm')^1.5 * 'm^2/s';
566
567	case "Roffel_Fl":
568	TRAYS(i).OutletLiquid.F = 2/3TRAYS(i).rhoL/sum(MwTRAYS(i).OutletLiquid.z)lw(TRAYS(i).MLsum(MwTRAYS(i).OutletLiquid.z)/(Ap1.3)/TRAYS(i).rhoL)^1.5sqrt(2*g/
569	(2(1 - 0.3593/'Pa^0.0888545'abs(TRAYS(i).OutletVapour.Fsum(MwTRAYS(i).OutletVapour.z)/(Ap*1.3)/sqrt(TRAYS(i).rhoV))^0.177709)-1)); #/'(kg/m)^0.0888545/s^0.177709';
570	end
571	when TRAYS(i).Level < (beta *hw) switchto "off";
572
573	case "off":
574	"Low level"
575	TRAYS(i).OutletLiquid.F = 0 * 'mol/h';
576	when TRAYS(i).Level > (beta * hw) + 1e-6*'m' switchto "on";
577	end
578
579	switch VapourFlow
580	case "on":
581	switch VapourFlowModel
582	case "Reepmeyer":
583	TRAYS(i).InletVapour.FTRAYS(i).vV = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)/(TRAYS(i).rhoValfa))*Ah;
584
585	case "Feehery_Fv":
586	TRAYS(i).InletVapour.F = TRAYS(i).rhoV/Ap/w/sum(MwTRAYS(i).OutletVapour.z) sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-(TRAYS(i).rhoVgTRAYS(i).ML*TRAYS(i).vL/Ap))/TRAYS(i).rhoV);
587
588	case "Roffel_Fv":
589	TRAYS(i).InletVapour.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.921e5' = (TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)1e5 - (betasum(TRAYS(i).MMw)/(Ap1.3)g1e5) (TRAYS(i).rhoVAh/sum(MwTRAYS(i).OutletVapour.z))^1.08 * 'm^1.08/mol^1.08';
590
591	case "Klingberg":
592	TRAYS(i).InletVapour.F * TRAYS(i).vV = Ap * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoLgTRAYS(i).Level)/TRAYS(i).rhoV);
593
594	case "Wang_Fv":
595	TRAYS(i).InletVapour.F * TRAYS(i).vV = Ap * sqrt(((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)-TRAYS(i).rhoLgTRAYS(i).Level)/TRAYS(i).rhoV*alfa);
596
597	case "Elgue":
598	TRAYS(i).InletVapour.F = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P)/btray);
599	end
600	when TRAYS(i).InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
601
602	case "off":
603	TRAYS(i).InletVapour.F = 0 * 'mol/s';
604	when TRAYS(i).InletVapour.P > TRAYS(i).OutletVapour.P + TRAYS(i).LevelgTRAYS(i).rhoL + 1e-1 * 'atm' switchto "on";
605	end
606
607	end
608
609	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
610	# Equating Inlet Liquid Connector Variables
611	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
612	"Liquid Inlet Flow"
613	LiquidConnector.F = LiquidInlet.F;
614
615	"Liquid Inlet Temperature"
616	LiquidConnector.T = LiquidInlet.T;
617
618	"Liquid Inlet Pressure"
619	LiquidConnector.P = LiquidInlet.P;
620
621	"Liquid Inlet Composition"
622	LiquidConnector.z = LiquidInlet.z;
623
624	"Liquid Inlet Vapour Fraction"
625	LiquidConnector.v = LiquidInlet.v;
626
627	"Liquid Inlet Enthalpy"
628	LiquidConnector.h = LiquidInlet.h;
629	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
630	# Equating Inlet Vapour Connector Variables
631	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
632	"Vapour Inlet Flow"
633	VapourConnector.F = VapourInlet.F;
634
635	"Vapour Inlet Temperature"
636	VapourConnector.T = VapourInlet.T;
637
638	"Vapour Inlet Pressure"
639	VapourConnector.P = VapourInlet.P;
640
641	"Vapour Inlet Composition"
642	VapourConnector.z = VapourInlet.z;
643
644	"Vapour Inlet Vapour Fraction"
645	VapourConnector.v = VapourInlet.v;
646
647	"Vapour Inlet Enthalpy"
648	VapourConnector.h = VapourInlet.h;
649	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
650	# Equating Outlet Liquid Variables
651	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
652	"Liquid Outlet Flow"
653	LiquidOutlet.F = TRAYS(NumberOfTrays).OutletLiquid.F;
654
655	"Liquid Outlet Temperature"
656	LiquidOutlet.T = TRAYS(NumberOfTrays).OutletLiquid.T;
657
658	"Liquid Outlet Pressure"
659	LiquidOutlet.P = TRAYS(NumberOfTrays).OutletLiquid.P;
660
661	"Liquid Outlet Composition"
662	LiquidOutlet.z = TRAYS(NumberOfTrays).OutletLiquid.z;
663	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
664	# Equating Outlet Vapour Variables
665	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
666	"Vapour Outlet Flow"
667	VapourOutlet.F = TRAYS(1).OutletVapour.F;
668
669	"Vapour Outlet Temperature"
670	VapourOutlet.T = TRAYS(1).OutletVapour.T;
671
672	"Vapour Outlet Pressure"
673	VapourOutlet.P = TRAYS(1).OutletVapour.P;
674
675	"Vapour Outlet Composition"
676	VapourOutlet.z = TRAYS(1).OutletVapour.z;
677	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
678	# Equating Outlet Vapour Side Stream Variables
679	#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
680	"Vapour Draw Off Stream - Flow"
681	VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F;
682
683	"Vapour Draw Off Stream - Temperature"
684	VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T;
685
686	"Vapour Draw Off Stream - Pressure"
687	VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P;
688
689	"Vapour Draw Off Stream - Composition"
690	VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z;
691
692	"Vapour Draw Off Stream"
693	VapourDrawOffFlow = VapourDrawOff.F;
694	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
695	# Equating Outlet Liquid Side Stream Variables
696	#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
697	"Liquid Draw Off Stream - Flow"
698	LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F;
699
700	"Liquid Draw Off Stream - Temperature"
701	LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T;
702
703	"Liquid Draw Off Stream - Pressure"
704	LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P;
705
706	"Liquid Draw Off Stream - Composition"
707	LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z;
708
709	"Liquid Draw Off Stream"
710	LiquidDrawOffFlow = LiquidDrawOff.F;
711
712	end
713
714	Model Distillation_kettle_cond as Section_ColumnBasic
715
716	ATTRIBUTES
717	Pallete = true;
718	Icon = "icon/DistillationKettleCond";
719	Brief = "Model of a distillation column with dynamic condenser and dynamic reboiler.";
720	Info =
721	"== Specify ==
722	* the feed stream of each tray (Inlet);
723	* the Murphree eficiency for each tray Emv;
724	* the pump pressure difference;
725	* the heat supllied in reboiler and condenser;
726	* the condenser vapor outlet flow (OutletVapour.F);
727	* the reboiler liquid outlet flow (OutletLiquid.F);
728	* both splitter outlet flows OR one of the splitter outlet flows and the splitter frac.
729
730	== Initial Conditions ==
731	* the TRAYS temperature (OutletLiquid.T);
732	* the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F);
733	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray;
734
735	* the condenser temperature (OutletLiquid.T);
736	* the condenser liquid level (Level);
737	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions;
738
739	* the reboiler temperature (OutletLiquid.T);
740	* the reboiler liquid level (Level);
741	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
742	";
743
744	PARAMETERS
745	CondenserVapourFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true);
746
747	VARIABLES
748	out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.375,Protected = true);
749	out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.42,Protected = true);
750	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
751	# Column Devices
752	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
753	CONDENSER as condenser;
754	REBOILER as reboiler;
755	SPLITTER as splitter2;
756	PUMP as pump;
757	alfaTopo as Real;
758	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
759	# Heat Ports
760	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
761	in HeatToReboiler as power (Brief="Heat supplied to Reboiler",Protected = true, PosX=1, PosY=0.90);
762	in HeatToCondenser as power (Brief="Heat supplied to Condenser", Protected = true, PosX=1, PosY=0.034);
763	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
764	# Column Connectors
765	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
766	ConnectorHeatReboiler as power (Brief="Heat supplied to Reboiler", Hidden=true);
767	ConnectorHeatCondenser as power (Brief="Heat supplied to Condenser", Hidden=true);
768	in ConnectorCondenserVout as stream (Brief="Connector for Vapour outlet stream From Condenser", Hidden=true);
769	in ConnectorReboilerLout as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true);
770	in ConnectorSplitterOut as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true);
771	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
772	# Column products
773	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
774	out VapourDistillate as vapour_stream (Brief="Vapour outlet stream From Condenser", PosX=0.66, PosY=0);
775	out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.20);
776	out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.67, PosY=1);
777
778	EQUATIONS
779
780	switch CondenserVapourFlow
781
782	case "on":
783	CONDENSER.InletVapour.FTRAYS(1).vV = alfaTopo Ah * sqrt(2*(TRAYS(1).OutletVapour.P -
784	CONDENSER.OutletLiquid.P + 1e-8 * 'atm') / (alfa*TRAYS(1).rhoV));
785	when CONDENSER.InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
786
787	case "off":
788	CONDENSER.InletVapour.F = 0 * 'mol/s';
789	when TRAYS(1).OutletVapour.P > CONDENSER.OutletLiquid.P + 1e-1 * 'atm' switchto "on";
790
791	end
792
793	# Condenser Connector Equations
794	ConnectorCondenserVout.T = VapourDistillate.T;
795	ConnectorCondenserVout.P = VapourDistillate.P;
796	ConnectorCondenserVout.F = VapourDistillate.F;
797	ConnectorCondenserVout.z = VapourDistillate.z;
798
799	# Splitter Connector Equations
800	ConnectorSplitterOut.T = LiquidDistillate.T;
801	ConnectorSplitterOut.P = LiquidDistillate.P;
802	ConnectorSplitterOut.F = LiquidDistillate.F;
803	ConnectorSplitterOut.z = LiquidDistillate.z;
804
805	# Reboiler Connector Equations
806	ConnectorReboilerLout.T = BottomProduct.T;
807	ConnectorReboilerLout.P = BottomProduct.P;
808	ConnectorReboilerLout.F = BottomProduct.F;
809	ConnectorReboilerLout.z = BottomProduct.z;
810
811	VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F;
812	VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T;
813	VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P;
814	VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z;
815
816	LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F;
817	LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T;
818	LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P;
819	LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z;
820
821	VapourDrawOffFlow = VapourDrawOff.F;
822	LiquidDrawOffFlow = LiquidDrawOff.F;
823
824	HeatToReboiler = ConnectorHeatReboiler;
825	HeatToCondenser = ConnectorHeatCondenser;
826
827	CONNECTIONS
828	#vapor
829	REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour;
830	TRAYS(1).OutletVapour to CONDENSER.InletVapour;
831
832	#liquid
833	CONDENSER.OutletLiquid to SPLITTER.Inlet;
834	SPLITTER.Outlet2 to PUMP.Inlet;
835	PUMP.Outlet to TRAYS(1).InletLiquid;
836	TRAYS(NumberOfTrays).OutletLiquid to REBOILER.InletLiquid;
837
838	#Connectors
839	CONDENSER.OutletVapour to ConnectorCondenserVout;
840	SPLITTER.Outlet1 to ConnectorSplitterOut;
841	REBOILER.OutletLiquid to ConnectorReboilerLout;
842
843	ConnectorHeatReboiler to REBOILER.InletQ;
844	ConnectorHeatCondenser to CONDENSER.InletQ;
845
846	end
847
848	Model Distillation_thermosyphon_subcooling as Section_ColumnBasic
849	ATTRIBUTES
850	Pallete = true;
851	Icon = "icon/DistillationThermosyphonSubcooling";
852	Brief = "Model of a distillation column with steady condenser and steady reboiler.";
853	Info =
854	"== Specify ==
855	* the feed stream of each tray (Inlet);
856	* the Murphree eficiency for each tray Emv;
857	* the pump head;
858	* the condenser pressure drop;
859	* the heat supllied in top and bottom tanks;
860	* the heat supllied in condenser and reboiler;
861	* the Outlet1 flow in the bottom splitter (spbottom.Outlet1.F) that corresponds to the bottom product;
862	* both top splitter outlet flows OR one of the splitter outlet flows and the splitter frac.
863
864	== Initial Conditions ==
865	* the TRAYS temperature (OutletLiquid.T);
866	* the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F);
867	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray;
868
869	* the top tank temperature (OutletLiquid.T);
870	* the top tank liquid level (Level);
871	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions;
872
873	* the bottom tank temperature (OutletLiquid.T);
874	* the bottom tank liquid level (Level);
875	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
876	";
877
878	PARAMETERS
879
880	CondenserVapourFlow as Switcher(Valid = ["on", "off"], Hidden=true, Default = "on");
881
882	VARIABLES
883	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
884	# Column Vapour and Liquid Draw Sides
885	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
886	out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.41,Protected = true);
887	out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.45,Protected = true);
888	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
889	# Column Devices
890	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
891	CONDENSER as condenserSteady (Brief="steady state condenser with subcooling");
892	ACCUMULATOR_TOP as tank_cylindrical (Brief="vessel drum (layed cilinder)");
893	SPLITTER_TOP as splitter2 (Brief="splitter which separate reflux and distillate");
894	PUMP as pump (Brief="pump in reflux stream");
895	REBOILER as reboilerSteady (Brief="steady state reboiler (thermosyphon)");
896	ACCUMULATOR_BOTTOM as tank (Brief="vessel in the bottom of column");
897	SPLITTER_BOTTOM as splitter2 (Brief="splitter who separate the bottom product and the stream to reboiler");
898	alfaTopo as Real;
899	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
900	# Heat Ports
901	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
902	in HeatToCondenser as power (Brief="Heat supplied to Condenser",PosX=1, PosY=0.065, Protected=true);
903	in HeatToReboiler as power (Brief="Heat supplied to Reboiler", PosX=1, PosY=0.87,Protected=true);
904	in HeatToBottomAccumulator as power (Brief="Heat supplied to Bottom Accumulator",PosX=0, PosY=0.925,Protected=true);
905	in HeatToTopAccumulator as power (Brief="Heat supplied to Top Accumulator", PosX=1, PosY=0.19,Protected=true);
906	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
907	# Column Top and Bottom product
908	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
909	out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.24);
910	out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Bottom Splitter", PosX=0.18, PosY=1);
911	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
912	# Column Model Connectors
913	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
914	in ConnectorSplitterBottom as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true);
915	in ConnectorSplitterTop as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true);
916
917	ConnectorHeatReboiler as power (Brief="Connector for Heat supplied to Reboiler", Hidden=true);
918	ConnectorHeatCondenser as power (Brief="Connector for Heat supplied to Condenser", Hidden=true);
919	ConnectorHeatTopAccumulator as power (Brief="Connector for Heat supplied to Top Accumulator", Hidden=true);
920	ConnectorHeatBottomAccumulator as power (Brief="Connector for Heat supplied to Bottom Accumulator", Hidden=true);
921
922	EQUATIONS
923	# Heat Connector Equations
924	HeatToCondenser = ConnectorHeatCondenser;
925	HeatToReboiler = ConnectorHeatReboiler;
926	HeatToBottomAccumulator = ConnectorHeatBottomAccumulator;
927	HeatToTopAccumulator = ConnectorHeatTopAccumulator;
928
929	# Top Splitter Connector Equations
930	ConnectorSplitterTop.T = LiquidDistillate.T;
931	ConnectorSplitterTop.P = LiquidDistillate.P;
932	ConnectorSplitterTop.F = LiquidDistillate.F;
933	ConnectorSplitterTop.z = LiquidDistillate.z;
934
935	# Bottom Splitter Connector Equations
936	ConnectorSplitterBottom.T = BottomProduct.T;
937	ConnectorSplitterBottom.P = BottomProduct.P;
938	ConnectorSplitterBottom.F = BottomProduct.F;
939	ConnectorSplitterBottom.z = BottomProduct.z;
940
941	VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F;
942	VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T;
943	VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P;
944	VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z;
945
946	LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F;
947	LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T;
948	LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P;
949	LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z;
950
951	VapourDrawOffFlow = VapourDrawOff.F;
952	LiquidDrawOffFlow = LiquidDrawOff.F;
953
954	switch CondenserVapourFlow
955
956	case "on":
957	CONDENSER.InletVapour.FTRAYS(1).vV = alfaTopo Ah * sqrt(2*(TRAYS(1).OutletVapour.P -
958	CONDENSER.OutletLiquid.P + 1e-8 * 'atm') / (alfa*TRAYS(1).rhoV));
959	when CONDENSER.InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
960
961	case "off":
962	CONDENSER.InletVapour.F = 0 * 'mol/s';
963	when TRAYS(1).OutletVapour.P > CONDENSER.OutletLiquid.P + 1e-1 * 'atm' switchto "on";
964
965	end
966
967	CONNECTIONS
968	#vapor
969	REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour;
970	TRAYS(1).OutletVapour to CONDENSER.InletVapour;
971
972	#liquid
973	CONDENSER.OutletLiquid to ACCUMULATOR_TOP.Inlet;
974	ACCUMULATOR_TOP.Outlet to SPLITTER_TOP.Inlet;
975	SPLITTER_TOP.Outlet2 to PUMP.Inlet;
976	PUMP.Outlet to TRAYS(1).InletLiquid;
977	TRAYS(NumberOfTrays).OutletLiquid to ACCUMULATOR_BOTTOM.Inlet;
978	ACCUMULATOR_BOTTOM.Outlet to SPLITTER_BOTTOM.Inlet;
979	SPLITTER_BOTTOM.Outlet2 to REBOILER.InletLiquid;
980
981	#Connectors
982	ConnectorHeatCondenser to CONDENSER.InletQ;
983	ConnectorHeatReboiler to REBOILER.InletQ;
984	ConnectorHeatBottomAccumulator to ACCUMULATOR_BOTTOM.InletQ;
985	ConnectorHeatTopAccumulator to ACCUMULATOR_TOP.InletQ;
986	SPLITTER_TOP.Outlet1 to ConnectorSplitterTop;
987	SPLITTER_BOTTOM.Outlet1 to ConnectorSplitterBottom;
988
989	end
990
991	Model Distillation_thermosyphon_cond as Section_ColumnBasic
992	ATTRIBUTES
993	Pallete = true;
994	Icon = "icon/DistillationThermosyphonCond";
995	Brief = "Model of a distillation column with dynamic condenser and steady reboiler.";
996	Info =
997	"== Specify ==
998	* the feed stream of each tray (Inlet);
999	* the Murphree eficiency for each tray Emv;
1000	* the pump head;
1001	* the condenser vapor outlet flow (OutletVapour.F);
1002	* the heat supllied in bottom tank;
1003	* the heat supllied in condenser and reboiler;
1004	* the Outlet1 flow in the bottom splitter (spbottom.Outlet1.F) that corresponds to the bottom product;
1005
1006	== Initial Conditions ==
1007	* the TRAYS temperature (OutletLiquid.T);
1008	* the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F);
1009	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray;
1010
1011	* the condenser temperature (OutletLiquid.T);
1012	* the condenser liquid level (Level);
1013	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions;
1014
1015	* the bottom tank temperature (OutletLiquid.T);
1016	* the bottom tank liquid level (Level);
1017	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
1018	";
1019
1020	PARAMETERS
1021	CondenserVapourFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true);
1022
1023	VARIABLES
1024	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1025	# Column Vapour and Liquid Draw Sides
1026	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1027	out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.379,Protected = true);
1028	out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.425,Protected = true);
1029	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1030	# Column Devices
1031	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1032	CONDENSER as condenser (Brief="dynamic condenser without subcooling");
1033	SPLITTER_TOP as splitter2 (Brief="splitter which separate reflux and distillate");
1034	PUMP as pump (Brief="pump in reflux stream");
1035	ACCUMULATOR_BOTTOM as tank (Brief="vessel in the bottom of column");
1036	SPLITTER_BOTTOM as splitter2 (Brief="splitter who separate the bottom product and the stream to reboiler");
1037	REBOILER as reboilerSteady (Brief="steady state reboiler (thermosyphon)");
1038	alfaTopo as Real;
1039	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1040	# Heat Ports
1041	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1042	in HeatToCondenser as power (Brief="Heat supplied to Condenser",PosX=1, PosY=0.035, Protected=true);
1043	in HeatToReboiler as power (Brief="Heat supplied to Reboiler", PosX=1, PosY=0.865,Protected=true);
1044	in HeatToBottomAccumulator as power (Brief="Heat supplied to Bottom Accumulator",PosX=0, PosY=0.92,Protected=true);
1045	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1046	# Column Top and Bottom product
1047	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1048	out VapourDistillate as vapour_stream (Brief="Vapour outlet stream From Top Condenser", PosX=0.665, PosY=0);
1049	out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.20);
1050	out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Bottom Splitter", PosX=0.18, PosY=1);
1051	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1052	# Column Model Connectors
1053	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1054	in ConnectorSplitterBottom as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true);
1055	in ConnectorSplitterTop as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true);
1056	in ConnectorCondenserVout as stream (Brief="Connector for Vapour outlet stream From Top Condenser", Hidden=true);
1057
1058	ConnectorHeatReboiler as power (Brief="Connector for Heat supplied to Reboiler", Hidden=true);
1059	ConnectorHeatCondenser as power (Brief="Connector for Heat supplied to Condenser", Hidden=true);
1060	ConnectorHeatBottomAccumulator as power (Brief="Connector for Heat supplied to Bottom Accumulator", Hidden=true);
1061
1062	EQUATIONS
1063	# Heat Connector Equations
1064	HeatToCondenser = ConnectorHeatCondenser;
1065	HeatToReboiler = ConnectorHeatReboiler;
1066	HeatToBottomAccumulator = ConnectorHeatBottomAccumulator;
1067
1068	# Condenser Connector Equations
1069	ConnectorCondenserVout.T = VapourDistillate.T;
1070	ConnectorCondenserVout.P = VapourDistillate.P;
1071	ConnectorCondenserVout.F = VapourDistillate.F;
1072	ConnectorCondenserVout.z = VapourDistillate.z;
1073
1074	# Top Splitter Connector Equations
1075	ConnectorSplitterTop.T = LiquidDistillate.T;
1076	ConnectorSplitterTop.P = LiquidDistillate.P;
1077	ConnectorSplitterTop.F = LiquidDistillate.F;
1078	ConnectorSplitterTop.z = LiquidDistillate.z;
1079
1080	# Bottom Splitter Connector Equations
1081	ConnectorSplitterBottom.T = BottomProduct.T;
1082	ConnectorSplitterBottom.P = BottomProduct.P;
1083	ConnectorSplitterBottom.F = BottomProduct.F;
1084	ConnectorSplitterBottom.z = BottomProduct.z;
1085
1086	VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F;
1087	VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T;
1088	VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P;
1089	VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z;
1090
1091	LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F;
1092	LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T;
1093	LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P;
1094	LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z;
1095
1096	VapourDrawOffFlow = VapourDrawOff.F;
1097	LiquidDrawOffFlow = LiquidDrawOff.F;
1098
1099	switch CondenserVapourFlow
1100	case "on":
1101	CONDENSER.InletVapour.FTRAYS(1).vV = alfaTopo Ah * sqrt(2*(TRAYS(1).OutletVapour.P -
1102	CONDENSER.OutletLiquid.P + 1e-8 * 'atm') / (alfa*TRAYS(1).rhoV));
1103	when CONDENSER.InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
1104
1105	case "off":
1106	CONDENSER.InletVapour.F = 0 * 'mol/s';
1107	when TRAYS(1).OutletVapour.P > CONDENSER.OutletLiquid.P + 1e-1 * 'atm' switchto "on";
1108	end
1109
1110	CONNECTIONS
1111	#vapor
1112	REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour;
1113	TRAYS(1).OutletVapour to CONDENSER.InletVapour;
1114
1115	#liquid
1116	CONDENSER.OutletLiquid to SPLITTER_TOP.Inlet;
1117	SPLITTER_TOP.Outlet2 to PUMP.Inlet;
1118	PUMP.Outlet to TRAYS(1).InletLiquid;
1119	TRAYS(NumberOfTrays).OutletLiquid to ACCUMULATOR_BOTTOM.Inlet;
1120	ACCUMULATOR_BOTTOM.Outlet to SPLITTER_BOTTOM.Inlet;
1121	SPLITTER_BOTTOM.Outlet2 to REBOILER.InletLiquid;
1122
1123	#Connectors
1124	ConnectorHeatCondenser to CONDENSER.InletQ;
1125	ConnectorHeatReboiler to REBOILER.InletQ;
1126	ConnectorHeatBottomAccumulator to ACCUMULATOR_BOTTOM.InletQ;
1127	CONDENSER.OutletVapour to ConnectorCondenserVout;
1128	SPLITTER_TOP.Outlet1 to ConnectorSplitterTop;
1129	SPLITTER_BOTTOM.Outlet1 to ConnectorSplitterBottom;
1130
1131	end
1132
1133	Model Distillation_kettle_subcooling as Section_ColumnBasic
1134
1135	ATTRIBUTES
1136	Pallete = true;
1137	Icon = "icon/DistillationKettleSubcooling";
1138	Brief = "Model of a distillation column with steady condenser and dynamic reboiler.";
1139	Info =
1140	"== Specify ==
1141	* the feed stream of each tray (Inlet);
1142	* the Murphree eficiency for each tray (Emv);
1143	* the pump pressure difference;
1144	* the heat supllied in reboiler and condenser;
1145	* the heat supllied in the top tank;
1146	* the condenser pressure drop;
1147	* the reboiler liquid outlet flow (OutletLiquid.F);
1148	* both splitter outlet flows OR one of the splitter outlet flows and the splitter frac.
1149
1150	== Initial Conditions ==
1151	* the TRAYS temperature (OutletLiquid.T);
1152	* the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F);
1153	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray;
1154
1155	* the top tank temperature (OutletLiquid.T);
1156	* the top tank liquid level (Level);
1157	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions;
1158
1159	* the reboiler temperature (OutletLiquid.T);
1160	* the reboiler liquid level (Level);
1161	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
1162	";
1163
1164	PARAMETERS
1165	CondenserVapourFlow as Switcher(Valid = ["on", "off"], Default = "on",Hidden=true);
1166
1167	VARIABLES
1168	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1169	# Column Vapour and Liquid Draw Sides
1170	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1171	out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.41,Protected = true);
1172	out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.45,Protected = true);
1173	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1174	# Column Devices
1175	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1176	CONDENSER as condenserSteady (Brief="steady state condenser with subcooling");
1177	ACCUMULATOR as tank_cylindrical (Brief="vessel drum (layed cilinder)");
1178	SPLITTER as splitter2 (Brief="splitter which separate reflux and distillate");
1179	PUMP as pump (Brief="pump in reflux stream");
1180	REBOILER as reboiler (Brief="kettle reboiler");
1181	alfaTopo as Real;
1182	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1183	# Heat Ports
1184	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1185	in HeatToCondenser as power (Brief="Heat supplied to Condenser", PosX=1, PosY=0.065, Protected=true);
1186	in HeatToReboiler as power (Brief="Heat supplied to Reboiler", PosX=1, PosY=0.905, Protected=true);
1187	in HeatToAccumulator as power (Brief="Heat supplied to Top Vessel", PosX=1, PosY=0.195, Protected=true);
1188	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1189	# Column Model Connectors
1190	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1191	in ConnectorSplitterOut as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true);
1192	in ConnectorReboilerLout as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true);
1193	ConnectorHeatReboiler as power (Brief="Connector for Heat supplied to Reboiler", Hidden=true);
1194	ConnectorHeatCondenser as power (Brief="Connector for Heat supplied to Condenser", Hidden=true);
1195	ConnectorHeatAccumulator as power (Brief="Connector for Heat supplied to TopVessel", Hidden=true);
1196	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1197	# Column Top and Bottom product
1198	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1199	out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.24);
1200	out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.67, PosY=1);
1201
1202	EQUATIONS
1203
1204	#Heat Connectors
1205	HeatToCondenser = ConnectorHeatCondenser;
1206	HeatToReboiler = ConnectorHeatReboiler;
1207	HeatToAccumulator = ConnectorHeatAccumulator;
1208
1209	# Splitter Connector Equations
1210	ConnectorSplitterOut.T = LiquidDistillate.T;
1211	ConnectorSplitterOut.P = LiquidDistillate.P;
1212	ConnectorSplitterOut.F = LiquidDistillate.F;
1213	ConnectorSplitterOut.z = LiquidDistillate.z;
1214
1215	# Reboiler Connector Equations
1216	ConnectorReboilerLout.T = BottomProduct.T;
1217	ConnectorReboilerLout.P = BottomProduct.P;
1218	ConnectorReboilerLout.F = BottomProduct.F;
1219	ConnectorReboilerLout.z = BottomProduct.z;
1220
1221	VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F;
1222	VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T;
1223	VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P;
1224	VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z;
1225
1226	LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F;
1227	LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T;
1228	LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P;
1229	LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z;
1230
1231	VapourDrawOffFlow = VapourDrawOff.F;
1232	LiquidDrawOffFlow = LiquidDrawOff.F;
1233
1234	switch CondenserVapourFlow
1235	case "on":
1236	CONDENSER.InletVapour.FTRAYS(1).vV = alfaTopo Ah * sqrt(2*(TRAYS(1).OutletVapour.P -
1237	CONDENSER.OutletLiquid.P + 1e-8 * 'atm') / (alfa*TRAYS(1).rhoV));
1238	when CONDENSER.InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
1239
1240	case "off":
1241	CONDENSER.InletVapour.F = 0 * 'mol/s';
1242	when TRAYS(1).OutletVapour.P > CONDENSER.OutletLiquid.P + 1e-1 * 'atm' switchto "on";
1243	end
1244
1245	CONNECTIONS
1246	#vapor
1247	REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour;
1248	TRAYS(1).OutletVapour to CONDENSER.InletVapour;
1249
1250	#liquid
1251	CONDENSER.OutletLiquid to ACCUMULATOR.Inlet;
1252	ACCUMULATOR.Outlet to SPLITTER.Inlet;
1253	SPLITTER.Outlet2 to PUMP.Inlet;
1254	PUMP.Outlet to TRAYS(1).InletLiquid;
1255	TRAYS(NumberOfTrays).OutletLiquid to REBOILER.InletLiquid;
1256
1257	#Connectors
1258	ConnectorHeatCondenser to CONDENSER.InletQ;
1259	ConnectorHeatReboiler to REBOILER.InletQ;
1260	ConnectorHeatAccumulator to ACCUMULATOR.InletQ;
1261
1262	SPLITTER.Outlet1 to ConnectorSplitterOut;
1263	REBOILER.OutletLiquid to ConnectorReboilerLout;
1264
1265	end
1266
1267	Model Rectifier as Section_ColumnBasic
1268
1269	ATTRIBUTES
1270	Pallete = true;
1271	Icon = "icon/RefluxedCond";
1272	Brief = "Model of a rectifier column with dynamic condenser.";
1273	Info =
1274	"== Specify ==
1275	* the feed stream of each tray (Inlet);
1276	* the Murphree eficiency for each tray Emv;
1277	* the InletVapour stream of the bottom tray unless its flow;
1278	* the pump pressure difference;
1279	* the heat supllied in the condenser;
1280	* the condenser vapor outlet flow (OutletVapour.F);
1281	* both splitter outlet flows OR one of the splitter outlet flows and the splitter frac.
1282
1283	== Initial Conditions ==
1284	* the TRAYS temperature (OutletLiquid.T);
1285	* the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F);
1286	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray;
1287
1288	* the condenser temperature (OutletLiquid.T);
1289	* the condenser liquid level (Level);
1290	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions;
1291	";
1292
1293	PARAMETERS
1294	CondenserVapourFlow as Switcher(Valid = ["on", "off"], Default = "on",Hidden=true);
1295
1296	VARIABLES
1297	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1298	# Column Vapour and Liquid Draw Sides
1299	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1300	out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.45,Protected = true);
1301	out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.50,Protected = true);
1302	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1303	# Column Devices
1304	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1305	CONDENSER as condenser (Brief="dymamic condenser without subcooling");
1306	SPLITTER as splitter2 (Brief="splitter which separate reflux and distillate");
1307	PUMP as pump (Brief="pump in reflux stream");
1308	alfaTopo as Real;
1309	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1310	# Column Top product
1311	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1312	out VapourDistillate as vapour_stream (Brief="Vapour outlet stream From Top Condenser", PosX=0.66, PosY=0);
1313	out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.24);
1314	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1315	# Heat Port
1316	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1317	in HeatToCondenser as power (Brief="Heat supplied to Condenser",PosX=1, PosY=0.04, Protected = true);
1318	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1319	# Column Bottom Outlets
1320	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1321	in VapourInlet as stream (Brief="Vapour Inlet in the section", PosX=0.07, PosY=1);
1322	out LiquidOutlet as liquid_stream (Brief="Liquid Outlet in the section", PosX=0.32, PosY=1);
1323	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1324	# Column Model Connectors
1325	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1326	in ConnectorSplitterOut as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true);
1327	in ConnectorCondenserVout as stream (Brief="Connector for Vapour outlet stream From Top Condenser", Hidden=true);
1328	VapourConnector as stream (Brief="Vapour connection at the middle TRAYS", PosX=0.55, PosY=0,Hidden=true);
1329	ConnectorHeatCondenser as power (Brief="Connector for Heat supplied to Condenser", Hidden=true);
1330
1331	EQUATIONS
1332
1333	switch CondenserVapourFlow
1334
1335	case "on":
1336	CONDENSER.InletVapour.FTRAYS(1).vV = alfaTopo Ah * sqrt(2*(TRAYS(1).OutletVapour.P -
1337	CONDENSER.OutletLiquid.P + 1e-8 * 'atm') / (alfa*TRAYS(1).rhoV));
1338	when CONDENSER.InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
1339
1340	case "off":
1341	CONDENSER.InletVapour.F = 0 * 'mol/s';
1342	when TRAYS(1).OutletVapour.P > CONDENSER.OutletLiquid.P + 1e-1 * 'atm' switchto "on";
1343
1344	end
1345
1346	HeatToCondenser = ConnectorHeatCondenser;
1347
1348	# Condenser Connector Equations
1349	ConnectorCondenserVout.T = VapourDistillate.T;
1350	ConnectorCondenserVout.P = VapourDistillate.P;
1351	ConnectorCondenserVout.F = VapourDistillate.F;
1352	ConnectorCondenserVout.z = VapourDistillate.z;
1353
1354	# Splitter Connector Equations
1355	ConnectorSplitterOut.T = LiquidDistillate.T;
1356	ConnectorSplitterOut.P = LiquidDistillate.P;
1357	ConnectorSplitterOut.F = LiquidDistillate.F;
1358	ConnectorSplitterOut.z = LiquidDistillate.z;
1359
1360	LiquidOutlet.F= TRAYS(NumberOfTrays).OutletLiquid.F;
1361	LiquidOutlet.T = TRAYS(NumberOfTrays).OutletLiquid.T;
1362	LiquidOutlet.P = TRAYS(NumberOfTrays).OutletLiquid.P;
1363	LiquidOutlet.z = TRAYS(NumberOfTrays).OutletLiquid.z;
1364
1365	VapourConnector.F= VapourInlet.F;
1366	VapourConnector.T = VapourInlet.T;
1367	VapourConnector.P = VapourInlet.P;
1368	VapourConnector.z = VapourInlet.z;
1369	VapourConnector.v = VapourInlet.v;
1370	VapourConnector.h = VapourInlet.h;
1371
1372	VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F;
1373	VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T;
1374	VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P;
1375	VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z;
1376
1377	LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F;
1378	LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T;
1379	LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P;
1380	LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z;
1381
1382	VapourDrawOffFlow = VapourDrawOff.F;
1383	LiquidDrawOffFlow = LiquidDrawOff.F;
1384
1385	CONNECTIONS
1386	#vapor
1387	TRAYS(1).OutletVapour to CONDENSER.InletVapour;
1388
1389	#liquid
1390	CONDENSER.OutletLiquid to SPLITTER.Inlet;
1391	SPLITTER.Outlet2 to PUMP.Inlet;
1392	PUMP.Outlet to TRAYS(1).InletLiquid;
1393
1394	#Connectors
1395	ConnectorHeatCondenser to CONDENSER.InletQ;
1396	VapourConnector to TRAYS(NumberOfTrays).InletVapour;
1397	SPLITTER.Outlet1 to ConnectorSplitterOut;
1398	CONDENSER.OutletVapour to ConnectorCondenserVout;
1399
1400	end
1401
1402	Model Rectifier_subcooling as Section_ColumnBasic
1403
1404	ATTRIBUTES
1405	Pallete = true;
1406	Icon = "icon/RefluxedSubcooling";
1407	Brief = "Model of a rectifier column with steady condenser.";
1408	Info =
1409	"== Specify ==
1410	* the feed stream of each tray (Inlet);
1411	* the Murphree eficiency for each tray Emv;
1412	* the InletVapour stream of the bottom tray unless its flow;
1413	* the pump head;
1414	* the condenser pressure drop;
1415	* the heat supllied in the top tank;
1416	* the heat supllied in condenser;
1417	* both top splitter outlet flows OR one of the splitter outlet flows and the splitter frac.
1418
1419	== Initial Conditions ==
1420	* the TRAYS temperature (OutletLiquid.T);
1421	* the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F);
1422	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray;
1423
1424	* the top tank temperature (OutletLiquid.T);
1425	* the top tank liquid level (Level);
1426	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions;
1427	";
1428
1429	PARAMETERS
1430	CondenserVapourFlow as Switcher (Valid = ["on", "off"], Default = "on",Hidden=true);
1431
1432	VARIABLES
1433	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1434	# Column Vapour and Liquid Draw Sides
1435	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1436	out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.485,Protected = true);
1437	out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.535,Protected = true);
1438	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1439	# Column Devices
1440	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1441	CONDENSER as condenserSteady (Brief="steady state condenser with subcooling");
1442	ACCUMULATOR as tank_cylindrical (Brief="vessel drum (layed cilinder)");
1443	SPLITTER as splitter2 (Brief="plitter which separate reflux and distillate");
1444	PUMP as pump (Brief="pump in reflux stream");
1445	alfaTopo as Real;
1446	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1447	# Column Bottom outlets
1448	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1449	in VapourInlet as stream (Brief="Vapour Inlet in the section", PosX=0.07, PosY=1);
1450	out LiquidOutlet as liquid_stream (Brief="Liquid Outlet in the section", PosX=0.32, PosY=1);
1451	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1452	# Heat Ports
1453	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1454	in HeatToCondenser as power (Brief="Heat supplied to Condenser", PosX=1, PosY=0.070, Protected=true);
1455	in HeatToAccumulator as power (Brief="Heat supplied to Top Accumulator", PosX=1, PosY=0.23, Protected=true);
1456	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1457	# Column Top Product
1458	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1459	out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.285);
1460	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1461	# Column Connectors
1462	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1463	in ConnectorSplitterTop as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true);
1464	VapourConnector as stream (Brief="Vapour connection at the middle TRAYS", PosX=0.55, PosY=0,Hidden=true);
1465	ConnectorHeatCondenser as power (Brief="Connector for Heat supplied to Condenser", Hidden=true);
1466	ConnectorHeatAccumulator as power (Brief="Connector for Heat supplied to Top Accumulator", Hidden=true);
1467
1468	EQUATIONS
1469	#Heat Connectors
1470	HeatToCondenser = ConnectorHeatCondenser;
1471	HeatToAccumulator = ConnectorHeatAccumulator;
1472
1473	LiquidOutlet.F= TRAYS(NumberOfTrays).OutletLiquid.F;
1474	LiquidOutlet.T = TRAYS(NumberOfTrays).OutletLiquid.T;
1475	LiquidOutlet.P = TRAYS(NumberOfTrays).OutletLiquid.P;
1476	LiquidOutlet.z = TRAYS(NumberOfTrays).OutletLiquid.z;
1477
1478	VapourConnector.F= VapourInlet.F;
1479	VapourConnector.T = VapourInlet.T;
1480	VapourConnector.P = VapourInlet.P;
1481	VapourConnector.z = VapourInlet.z;
1482	VapourConnector.v = VapourInlet.v;
1483	VapourConnector.h = VapourInlet.h;
1484
1485	# Splitter Connector Equations
1486	ConnectorSplitterTop.T = LiquidDistillate.T;
1487	ConnectorSplitterTop.P = LiquidDistillate.P;
1488	ConnectorSplitterTop.F = LiquidDistillate.F;
1489	ConnectorSplitterTop.z = LiquidDistillate.z;
1490
1491	VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F;
1492	VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T;
1493	VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P;
1494	VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z;
1495
1496	LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F;
1497	LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T;
1498	LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P;
1499	LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z;
1500
1501	VapourDrawOffFlow = VapourDrawOff.F;
1502	LiquidDrawOffFlow = LiquidDrawOff.F;
1503
1504	switch CondenserVapourFlow
1505
1506	case "on":
1507	CONDENSER.InletVapour.FTRAYS(1).vV = alfaTopo Ah * sqrt(2*(TRAYS(1).OutletVapour.P -
1508	CONDENSER.OutletLiquid.P + 1e-8 * 'atm') / (alfa*TRAYS(1).rhoV));
1509	when CONDENSER.InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
1510
1511	case "off":
1512	CONDENSER.InletVapour.F = 0 * 'mol/s';
1513	when TRAYS(1).OutletVapour.P > CONDENSER.OutletLiquid.P + 1e-1 * 'atm' switchto "on";
1514
1515	end
1516
1517	CONNECTIONS
1518	#vapor
1519	TRAYS(1).OutletVapour to CONDENSER.InletVapour;
1520
1521	#liquid
1522	CONDENSER.OutletLiquid to ACCUMULATOR.Inlet;
1523	ACCUMULATOR.Outlet to SPLITTER.Inlet;
1524	SPLITTER.Outlet2 to PUMP.Inlet;
1525	PUMP.Outlet to TRAYS(1).InletLiquid;
1526
1527	#Connectors
1528	VapourConnector to TRAYS(NumberOfTrays).InletVapour;
1529	ConnectorHeatCondenser to CONDENSER.InletQ;
1530	SPLITTER.Outlet1 to ConnectorSplitterTop;
1531	ConnectorHeatAccumulator to ACCUMULATOR.InletQ;
1532
1533	end
1534
1535	Model Reboiled_Stripping_kettle as Section_ColumnBasic
1536
1537	ATTRIBUTES
1538	Pallete = true;
1539	Icon = "icon/ReboiledKettle";
1540	Brief = "Model of a reboiled stripping column with dynamic reboiler.";
1541	Info =
1542	"== Specify ==
1543	* the feed stream of each tray (Inlet);
1544	* the Murphree eficiency for each tray Emv;
1545	* the vapour flow leaving the top of the column;
1546	* the InletLiquidiquid stream of the top tray;
1547	* the heat supllied in the reboiler;
1548	* the reboiler liquid outlet flow (OutletLiquid.F);
1549
1550	== Initial Conditions ==
1551	* the TRAYS temperature (OutletLiquid.T);
1552	* the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F);
1553	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray;
1554
1555	* the reboiler temperature (OutletLiquid.T);
1556	* the reboiler liquid level (Level);
1557	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
1558	";
1559
1560	VARIABLES
1561	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1562	# Column Vapour and Liquid Draw Sides
1563	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1564	out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.28,Protected = true);
1565	out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.33,Protected = true);
1566	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1567	# Column Devices
1568	#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1569	REBOILER as reboiler (Brief="Kettle Reboiler");
1570	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1571	# Heat Port
1572	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1573	in HeatToReboiler as power (Brief="Heat supplied to Reboiler",PosX=1, PosY=0.865,Protected=true);
1574	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1575	# Column Top Outlets
1576	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1577	in LiquidInlet as stream (Brief="Liquid Inlet in the section", PosX=0.30, PosY=0);
1578	out VapourOutlet as vapour_stream (Brief="Vapour Outlet in the section", PosX=0.07, PosY=0);
1579	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1580	# Column Bottom Product
1581	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1582	out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.68, PosY=1);
1583	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1584	# Column Connectors
1585	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1586	in ConnectorReboilerLout as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true);
1587	LiquidConnector as stream (Brief="Liquid connection at the middle TRAYS", PosX=0.75, PosY=1,Hidden=true);
1588	ConnectorHeatReboiler as power (Brief="Connector for Heat supplied to Reboiler", Hidden=true);
1589
1590	CONNECTIONS
1591	#vapor
1592	REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour;
1593
1594	#liquid
1595	TRAYS(NumberOfTrays).OutletLiquid to REBOILER.InletLiquid;
1596
1597	#Connectors
1598	REBOILER.OutletLiquid to ConnectorReboilerLout;
1599	LiquidConnector to TRAYS(1).InletLiquid;
1600	ConnectorHeatReboiler to REBOILER.InletQ;
1601
1602	EQUATIONS
1603	#Heat Connectors
1604	HeatToReboiler = ConnectorHeatReboiler;
1605
1606	LiquidConnector.F= LiquidInlet.F;
1607	LiquidConnector.T = LiquidInlet.T;
1608	LiquidConnector.P = LiquidInlet.P;
1609	LiquidConnector.z = LiquidInlet.z;
1610	LiquidConnector.v = LiquidInlet.v;
1611	LiquidConnector.h = LiquidInlet.h;
1612
1613	VapourOutlet.F= TRAYS(1).OutletVapour.F;
1614	VapourOutlet.T = TRAYS(1).OutletVapour.T;
1615	VapourOutlet.P = TRAYS(1).OutletVapour.P;
1616	VapourOutlet.z = TRAYS(1).OutletVapour.z;
1617
1618	# Reboiler Connector Equations
1619	ConnectorReboilerLout.T = BottomProduct.T;
1620	ConnectorReboilerLout.P = BottomProduct.P;
1621	ConnectorReboilerLout.F = BottomProduct.F;
1622	ConnectorReboilerLout.z = BottomProduct.z;
1623
1624	VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F;
1625	VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T;
1626	VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P;
1627	VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z;
1628
1629	LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F;
1630	LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T;
1631	LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P;
1632	LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z;
1633
1634	VapourDrawOffFlow = VapourDrawOff.F;
1635	LiquidDrawOffFlow = LiquidDrawOff.F;
1636
1637	end
1638
1639	Model Reboiled_Stripping_thermosyphon as Section_ColumnBasic
1640
1641	ATTRIBUTES
1642	Pallete = true;
1643	Icon = "icon/ReboiledThermosyphon";
1644	Brief = "Model of a reboiled stripping column with steady reboiler.";
1645	Info =
1646	"== Specify ==
1647	* the feed stream of each tray (Inlet);
1648	* the Murphree eficiency for each tray (Emv);
1649	* the vapour flow leaving the top of the column;
1650	* the InletLiquidiquid stream of the top tray;
1651	* the heat supllied in bottom tank;
1652	* the heat supllied in the reboiler;
1653	* the Outlet1 flow in the bottom splitter (spbottom.Outlet1.F) that corresponds to the bottom product;
1654
1655	== Initial Conditions ==
1656	* the TRAYS temperature (OutletLiquid.T);
1657	* the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F);
1658	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray;
1659
1660	* the bottom tank temperature (OutletLiquid.T);
1661	* the bottom tank liquid level (Level);
1662	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
1663	";
1664
1665	VARIABLES
1666	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1667	# Column Vapour and Liquid Draw Sides
1668	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1669	out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.28,Protected = true);
1670	out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.33,Protected = true);
1671	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1672	# Column Devices
1673	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1674	REBOILER as reboilerSteady (Brief="steady state reboiler (thermosyphon)");
1675	SPLITTER as splitter2 (Brief="splitter which separate the bottom product and the stream to reboiler");
1676	ACCUMULATOR as tank (Brief="vessel in the bottom of column");
1677	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1678	# Heat Ports
1679	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1680	in HeatToReboiler as power (Brief="Heat supplied to Reboiler",PosX=1, PosY=0.84,Protected=true);
1681	in HeatToAccumulator as power (Brief="Heat supplied to Bottom Vessel",PosX=0, PosY=0.91,Protected=true);
1682	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1683	# Column Top Outlets
1684	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1685	in LiquidInlet as stream (Brief="Liquid Inlet in the section", PosX=0.32, PosY=0);
1686	out VapourOutlet as vapour_stream (Brief="Vapour Outlet in the section", PosX=0.07, PosY=0);
1687	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1688	# Column Bottom Product
1689	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1690	out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.19, PosY=1);
1691	LiquidConnector as stream (Brief="Liquid connection at the middle TRAYS", PosX=0.75, PosY=1,Hidden=true);
1692	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1693	# Column Connectors
1694	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1695	in ConnectorSplitterBottom as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true);
1696	ConnectorHeatReboiler as power (Brief="Connector for Heat supplied to Reboiler", Hidden=true);
1697	ConnectorHeatAccumulator as power (Brief="Connector for Heat supplied to Bottom Accumulator", Hidden=true);
1698
1699	CONNECTIONS
1700	#vapor
1701	REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour;
1702
1703	#liquid
1704	TRAYS(NumberOfTrays).OutletLiquid to ACCUMULATOR.Inlet;
1705	ACCUMULATOR.Outlet to SPLITTER.Inlet;
1706	SPLITTER.Outlet2 to REBOILER.InletLiquid;
1707
1708	#Connectors
1709	LiquidConnector to TRAYS(1).InletLiquid;
1710	ConnectorHeatReboiler to REBOILER.InletQ;
1711	ConnectorHeatAccumulator to ACCUMULATOR.InletQ;
1712	SPLITTER.Outlet1 to ConnectorSplitterBottom;
1713
1714	EQUATIONS
1715	#Heat Connectors
1716	HeatToReboiler = ConnectorHeatReboiler;
1717	HeatToAccumulator = ConnectorHeatAccumulator;
1718
1719	ConnectorSplitterBottom.T = BottomProduct.T;
1720	ConnectorSplitterBottom.P = BottomProduct.P;
1721	ConnectorSplitterBottom.F = BottomProduct.F;
1722	ConnectorSplitterBottom.z = BottomProduct.z;
1723
1724	LiquidConnector.F= LiquidInlet.F;
1725	LiquidConnector.T = LiquidInlet.T;
1726	LiquidConnector.P = LiquidInlet.P;
1727	LiquidConnector.z = LiquidInlet.z;
1728	LiquidConnector.v = LiquidInlet.v;
1729	LiquidConnector.h = LiquidInlet.h;
1730
1731	VapourOutlet.F= TRAYS(1).OutletVapour.F;
1732	VapourOutlet.T = TRAYS(1).OutletVapour.T;
1733	VapourOutlet.P = TRAYS(1).OutletVapour.P;
1734	VapourOutlet.z = TRAYS(1).OutletVapour.z;
1735
1736	VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F;
1737	VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T;
1738	VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P;
1739	VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z;
1740
1741	LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F;
1742	LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T;
1743	LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P;
1744	LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z;
1745
1746	VapourDrawOffFlow = VapourDrawOff.F;
1747	LiquidDrawOffFlow = LiquidDrawOff.F;
1748
1749	end
1750
1751	Model ReactiveDistillation
1752
1753	ATTRIBUTES
1754	Pallete = true;
1755	Icon = "icon/DistillationReac";
1756	Brief = "Model of a reactive distillation column with dynamic condenser and reboiler.";
1757	Info =
1758	"== Specify ==
1759	* the reaction related variables for each tray, condenser and reboiler;
1760	* the feed stream of each tray (Inlet);
1761	* the Murphree eficiency for each tray Emv;
1762	* the pump pressure difference;
1763	* the heat supllied in reboiler and condenser;
1764	* the condenser vapor outlet flow (OutletVapour.F);
1765	* the reboiler liquid outlet flow (OutletLiquid.F);
1766	* both splitter outlet flows OR one of the splitter outlet flows and the splitter frac.
1767
1768	== Initial Conditions ==
1769	* the TRAYS temperature (OutletLiquid.T);
1770	* the TRAYS liquid level (Level) OR the TRAYS liquid flow (OutletLiquid.F);
1771	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray;
1772
1773	* the condenser temperature (OutletLiquid.T);
1774	* the condenser liquid level (Level);
1775	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions;
1776
1777	* the reboiler temperature (OutletLiquid.T);
1778	* the reboiler liquid level (Level);
1779	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
1780	";
1781
1782	PARAMETERS
1783	outer PP as Plugin (Type="PP");
1784	outer NComp as Integer (Brief="Number of Components");
1785	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1786	# Section Column Feed Tray - Side Streams Location and Numbering
1787	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1788	NumberOfTrays as Integer (Brief="Number of TRAYS", Default=2);
1789	FeedTrayIndex(NumberOfTrays) as Integer (Brief="Feed Tray Index", Default=0,Hidden=true);
1790	LiqSideTrayIndex(NumberOfTrays) as Integer (Brief="Liquid Side Tray Index", Default=0,Hidden=true);
1791	VapSideTrayIndex(NumberOfTrays) as Integer (Brief="Vapour Side Tray Index", Default=0,Hidden=true);
1792	FeedTrayLocation as Integer (Brief="Feed tray Location", Default=2);
1793	LiquidSideStreamLocation as Integer (Brief="Liquid Side Stream Location", Default=2);
1794	VapourSideStreamLocation as Integer (Brief="Vapour Side Stream Location", Default=2);
1795	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1796	# Section Column Flow Model for Liquid and Vapour
1797	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1798	VapourFlow as Switcher (Valid = ["on", "off"], Default = "off",Hidden=true);
1799	TrayVapourFlow as Switcher (Valid = ["on", "off"], Default = "off",Hidden=true);
1800	TrayLiquidFlow as Switcher (Valid = ["on", "off"], Default = "off",Hidden=true);
1801	alfacond as Real;
1802	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1803	# Section Column Tray Geometry and Auxiliar Parameters
1804	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1805	V as volume (Brief="Total Volume of the tray",Hidden=true);
1806	Q as power (Brief="Rate of heat supply",Hidden=true);
1807	Ap as area (Brief="Plate area = Atray - Adowncomer",Hidden=true);
1808	Ah as area (Brief="Total holes area",Hidden=true);
1809	lw as length (Brief="Weir length",Hidden=true);
1810	g as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
1811	hw as length (Brief="Weir height",Hidden=true);
1812	beta as fraction (Brief="Aeration fraction");
1813	alfa as fraction (Brief="Dry pressure drop coefficient");
1814
1815	VolumeOfTray as volume (Brief="Total Volume of the tray");
1816	HeatSupply as heat_rate (Brief="Rate of heat supply");
1817	PlateArea as area (Brief="Plate area = Atray - Adowncomer");
1818	HolesArea as area (Brief="Total holes area");
1819	WeirLength as length (Brief="Weir length");
1820	WeirHeight as length (Brief="Weir height");
1821
1822	stoic(NComp) as Real(Brief="Stoichiometric matrix");
1823	Hr as energy_mol;
1824
1825	SET
1826
1827	FeedTrayIndex(FeedTrayLocation) =1;
1828	VapSideTrayIndex(FeedTrayLocation) =1;
1829	LiqSideTrayIndex(FeedTrayLocation) =1;
1830
1831	V=VolumeOfTray;
1832	Q=HeatSupply;
1833	Ap=PlateArea;
1834	Ah=HolesArea;
1835	lw=WeirLength;
1836	hw=WeirHeight ;
1837
1838	VARIABLES
1839	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1840	# Tray Initialization For the Whole Column Section
1841	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1842	INITIALIZATION as InitializeSection (Brief = "Column Model Initialization");
1843	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1844	# Column Feed - Side Stream Flow and Murphree Efficiency
1845	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1846	in FeedTray as stream (Brief="Feed stream", PosX=0, PosY=0.55);
1847	VapourDrawOffFlow as flow_mol (Brief = "Stream Molar Flow Rate");
1848	LiquidDrawOffFlow as flow_mol (Brief = "Stream Molar Flow Rate");
1849	MurphreeEff as Real (Brief = "Murphree efficiency");
1850	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1851	# Column Vapour and Liquid Draw Sides
1852	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1853	out VapourDrawOff as vapour_stream (Brief="Vapour Outlet in the section", PosX=1, PosY=0.395,Protected = true);
1854	out LiquidDrawOff as liquid_stream (Brief="Liquid Outlet in the section", PosX=1, PosY=0.44,Protected = true);
1855	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1856	# Heat Ports
1857	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1858	in HeatToReboiler as power (Brief="Heat supplied to Reboiler", PosX=1, PosY=0.89);
1859	in HeatToCondenser as power (Brief="Heat supplied to Condenser", PosX=1, PosY=0.07);
1860	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1861	# Column Devices
1862	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1863	TRAYS(NumberOfTrays) as trayReac;
1864	CONDENSER as condenserReact;
1865	REBOILER as reboilerReact;
1866	SPLITTER as splitter2;
1867	PUMP as pump;
1868	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1869	# Column Outlets
1870	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1871	out VapourDistillate as vapour_stream (Brief="Vapour outlet stream From Condenser", PosX=0.59, PosY=0);
1872	out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Splitter", PosX=1, PosY=0.23);
1873	out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.67, PosY=1);
1874	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1875	# Column Connectors
1876	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1877	in ConnectorCondenserVout as stream (Brief="Connector for Vapour outlet stream From Top Condenser", Hidden=true);
1878	in ConnectorSplitterOut as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true);
1879	in ConnectorReboilerLout as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true);
1880	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1881	# Column Control Ports
1882	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1883	out TRI as control_signal (Brief="Temperature of Reboiler Indicator", PosX=1, PosY=0.83);
1884	out LRI as control_signal (Brief="Level of Reboiler Indicator", PosX=1, PosY=0.855);
1885	out TCI as control_signal (Brief="Temperature of Condenser Indicator", PosX=1, PosY=0.01);
1886	out LCI as control_signal (Brief="Level of Condenser Indicator", PosX=1, PosY=0.04);
1887
1888	CONNECTIONS
1889	#vapor
1890	REBOILER.OutletVapour to TRAYS(NumberOfTrays).InletVapour;
1891	TRAYS(1).OutletVapour to CONDENSER.InletVapour;
1892	TRAYS([2:NumberOfTrays]).OutletVapour to TRAYS([1:NumberOfTrays-1]).InletVapour;
1893
1894	#liquid
1895	CONDENSER.OutletLiquid to SPLITTER.Inlet;
1896	SPLITTER.Outlet2 to PUMP.Inlet;
1897	PUMP.Outlet to TRAYS(1).InletLiquid;
1898	TRAYS([1:NumberOfTrays-1]).OutletLiquid to TRAYS([2:NumberOfTrays]).InletLiquid;
1899	TRAYS(NumberOfTrays).OutletLiquid to REBOILER.InletLiquid;
1900
1901	#Connectors
1902	CONDENSER.OutletVapour to ConnectorCondenserVout;
1903	SPLITTER.Outlet1 to ConnectorSplitterOut;
1904	REBOILER.OutletLiquid to ConnectorReboilerLout;
1905
1906	INITIAL
1907	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1908	# Tray Initialization
1909	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1910	for i in 1:NumberOfTrays do
1911
1912	"The initial temperature of the TRAYS"
1913	TRAYS(i).OutletLiquid.T = INITIALIZATION.TopTemperature+(INITIALIZATION.BottomTemperature-INITIALIZATION.TopTemperature)*((i-1)/(NumberOfTrays-1));
1914
1915	"The initial Level of the TRAYS"
1916	TRAYS(i).Level = INITIALIZATION.LevelFraction*hw;
1917	end
1918
1919	for i in 1:NComp-1 do
1920
1921	for j in 1:NumberOfTrays do
1922
1923	"The initial composition of the TRAYS - Normalized"
1924	TRAYS(j).OutletLiquid.z(i) = INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) +(INITIALIZATION.BottomComposition(i)/sum(INITIALIZATION.BottomComposition)-INITIALIZATION.TopComposition(i)/sum(INITIALIZATION.TopComposition) )*((j-1)/(NumberOfTrays-1));
1925	end
1926
1927	end
1928
1929	EQUATIONS
1930
1931	for i in [1:NumberOfTrays] do
1932
1933	"Murphree Efficiency"
1934	TRAYS(i).OutletVapour.z = MurphreeEff * (TRAYS(i).yideal - TRAYS(i).InletVapour.z) + TRAYS(i).InletVapour.z;
1935
1936	"Energy Holdup"
1937	TRAYS(i).E = TRAYS(i).MLTRAYS(i).OutletLiquid.h + TRAYS(i).MVTRAYS(i).OutletVapour.h - TRAYS(i).OutletLiquid.P*V;
1938
1939	"Component Molar Balance"
1940	diff(TRAYS(i).M)=TRAYS(i).Inlet.FTRAYS(i).Inlet.z + TRAYS(i).InletLiquid.FTRAYS(i).InletLiquid.z + TRAYS(i).InletVapour.FTRAYS(i).InletVapour.z- TRAYS(i).OutletLiquid.FTRAYS(i).OutletLiquid.z - TRAYS(i).OutletVapour.F*TRAYS(i).OutletVapour.z-
1941	TRAYS(i).LiquidSideStream.FTRAYS(i).LiquidSideStream.z-TRAYS(i).VapourSideStream.FTRAYS(i).VapourSideStream.z + stoicTRAYS(i).r3TRAYS(i).ML*TRAYS(i).vL;
1942
1943	"Energy Balance"
1944	diff(TRAYS(i).E) = ( TRAYS(i).Inlet.FTRAYS(i).Inlet.h + TRAYS(i).InletLiquid.FTRAYS(i).InletLiquid.h + TRAYS(i).InletVapour.FTRAYS(i).InletVapour.h- TRAYS(i).OutletLiquid.FTRAYS(i).OutletLiquid.h - TRAYS(i).OutletVapour.F*TRAYS(i).OutletVapour.h
1945	-TRAYS(i).VapourSideStream.FTRAYS(i).VapourSideStream.h - TRAYS(i).LiquidSideStream.FTRAYS(i).LiquidSideStream.h + Q ) + Hr * TRAYS(i).r3 TRAYS(i). vLTRAYS(i).ML;
1946
1947	"Level of clear liquid over the weir"
1948	TRAYS(i).Level = TRAYS(i).ML*TRAYS(i).vL/Ap;
1949
1950	"Geometry Constraint"
1951	V = TRAYS(i).ML* TRAYS(i).vL + TRAYS(i).MV*TRAYS(i).vV;
1952
1953	end
1954
1955	switch VapourFlow
1956	case "on":
1957	"Pressure Drop through the condenser"
1958	CONDENSER.InletVapour.F*TRAYS(1).vV / 'm^2' =
1959	sqrt((TRAYS(1).OutletVapour.P - CONDENSER.OutletLiquid.P + 1e-8 * 'atm')/(TRAYS(1).rhoV*alfacond));
1960	when TRAYS(1).OutletVapour.P < CONDENSER.OutletLiquid.P switchto "off";
1961
1962	case "off":
1963	"Prato selado"
1964	CONDENSER.InletVapour.F = 0.0 * 'mol/s';
1965	when TRAYS(1).OutletVapour.P > CONDENSER.OutletLiquid.P + 1e-3 * 'atm' switchto "on";
1966	end
1967
1968	for i in [1:NumberOfTrays] do
1969
1970	switch TrayLiquidFlow
1971	case "on":
1972	"Francis Equation"
1973	TRAYS(i).OutletLiquid.FTRAYS(i).vL = 1.84'1/s'lw((TRAYS(i).Level-(betahw)+1e-6'm')/(beta))^2;
1974	when TRAYS(i).Level < (beta * hw) switchto "off";
1975
1976	case "off":
1977	"Low level"
1978	TRAYS(i).OutletLiquid.F = 0 * 'mol/h';
1979	when TRAYS(i).Level > (beta * hw) + 1e-6*'m' switchto "on";
1980	end
1981
1982	switch TrayVapourFlow
1983	case "on":
1984	TRAYS(i).InletVapour.FTRAYS(i).vV = sqrt((TRAYS(i).InletVapour.P - TRAYS(i).OutletVapour.P - TRAYS(i).LevelgTRAYS(i).rhoL + 1e-8 'atm')/(TRAYS(i).rhoValfa))Ah;
1985	when TRAYS(i).InletVapour.P < TRAYS(i).OutletVapour.P + TRAYS(i).LevelgTRAYS(i).rhoL switchto "off";
1986
1987	case "off":
1988	TRAYS(i).InletVapour.F = 0 * 'mol/s';
1989	when TRAYS(i).InletVapour.P > TRAYS(i).OutletVapour.P + TRAYS(i).LevelgTRAYS(i).rhoL + 3e-2 * 'atm' switchto "on";
1990	end
1991
1992	end
1993
1994	# Connecting Trays
1995	FeedTray.F*FeedTrayIndex= TRAYS.Inlet.F;
1996	FeedTray.T = TRAYS.Inlet.T;
1997	FeedTray.P = TRAYS.Inlet.P;
1998	FeedTray.z = TRAYS.Inlet.z;
1999	FeedTray.v = TRAYS.Inlet.v;
2000	FeedTray.h = TRAYS.Inlet.h;
2001
2002	HeatToReboiler = REBOILER.InletQ;
2003	HeatToCondenser = CONDENSER.InletQ;
2004
2005	"Reboiler Temperature indicator"
2006	TRI * 'K' = REBOILER.OutletLiquid.T;
2007
2008	"Reboiler Level indicator"
2009	LRIREBOILER.V = REBOILER.LevelREBOILER.Across;
2010
2011	"Condenser Temperature indicator"
2012	TCI * 'K' = CONDENSER.OutletLiquid.T;
2013
2014	"Condenser Level indicator"
2015	LCICONDENSER.V = CONDENSER.LevelCONDENSER.Across;
2016
2017	# Condenser Connector Equations
2018	ConnectorCondenserVout.T = VapourDistillate.T;
2019	ConnectorCondenserVout.P = VapourDistillate.P;
2020	ConnectorCondenserVout.F = VapourDistillate.F;
2021	ConnectorCondenserVout.z = VapourDistillate.z;
2022
2023	# Splitter Connector Equations
2024	ConnectorSplitterOut.T = LiquidDistillate.T;
2025	ConnectorSplitterOut.P = LiquidDistillate.P;
2026	ConnectorSplitterOut.F = LiquidDistillate.F;
2027	ConnectorSplitterOut.z = LiquidDistillate.z;
2028
2029	# Reboiler Connector Equations
2030	ConnectorReboilerLout.T = BottomProduct.T;
2031	ConnectorReboilerLout.P = BottomProduct.P;
2032	ConnectorReboilerLout.F = BottomProduct.F;
2033	ConnectorReboilerLout.z = BottomProduct.z;
2034
2035	VapourDrawOff.F*VapSideTrayIndex= TRAYS.VapourSideStream.F;
2036	VapourDrawOff.T = TRAYS(VapourSideStreamLocation).VapourSideStream.T;
2037	VapourDrawOff.P = TRAYS(VapourSideStreamLocation).VapourSideStream.P;
2038	VapourDrawOff.z = TRAYS(VapourSideStreamLocation).VapourSideStream.z;
2039
2040	LiquidDrawOff.F*LiqSideTrayIndex= TRAYS.LiquidSideStream.F;
2041	LiquidDrawOff.T = TRAYS(LiquidSideStreamLocation).LiquidSideStream.T;
2042	LiquidDrawOff.P = TRAYS(LiquidSideStreamLocation).LiquidSideStream.P;
2043	LiquidDrawOff.z = TRAYS(LiquidSideStreamLocation).LiquidSideStream.z;
2044
2045	VapourDrawOffFlow = VapourDrawOff.F;
2046	LiquidDrawOffFlow = LiquidDrawOff.F;
2047
2048	end
2049
2050	Model Packed_Section_ColumnBasic
2051	ATTRIBUTES
2052	Pallete = false;
2053	Icon = "icon/PackedSectionColumn";
2054	Brief = "Model of a packed column section.";
2055	Info =
2056	"== Model of a packed column section containing ==
2057	* NStages theoretical stages.
2058
2059	== Specify ==
2060	* the feed stream of each tray (Inlet);
2061	* the InletLiquid stream of the top tray;
2062	* the InletVapour stream of the bottom tray;
2063	* the total pressure drop (dP) of the section.
2064
2065	== Initial Conditions ==
2066	* the stages temperature (OutletLiquid.T);
2067	* the stages liquid holdup;
2068	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray.
2069	";
2070
2071	PARAMETERS
2072
2073	outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
2074	outer NComp as Integer;
2075
2076	NumberOfStages as Integer (Brief="Number of Stages", Default=3);
2077	FeedStageIndex(NumberOfStages) as Integer (Brief="Number of Stages", Default=0,Hidden=true);
2078	FeedStageLocation as Integer (Brief="Feed Stage Location", Default=2);
2079
2080	PackingHeight as length (Brief="Height of packing");
2081	HeatSupply as heat_rate (Brief="Rate of heat supply");
2082	ColumnDiameter as length (Brief="Column diameter");
2083	VoidFraction as Real (Brief="Void fraction of packing, m^3/m^3");
2084	ResistanceCoeff as positive (Brief="Resistance coefficient on the liquid load", Default=1);
2085	AreaPerPackingVol as Real (Brief="surface area per packing volume", Unit='m^2/m^3');
2086
2087	V as volume (Brief="Total Volume of the tray",Hidden=true);
2088	g as acceleration (Brief="Gravity Acceleration",Default=9.81,Hidden=true);
2089	hs as length (Brief="Height of the packing stage",Hidden=true);
2090	ColumnArea as area (Brief="Column Sectional Cross Area",Hidden=true);
2091
2092	VARIABLES
2093	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2094	# Stages Initialization For the Whole Column Section
2095	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2096	INITIALIZATION as InitializeStage (Brief = "Column Model Initialization");
2097	in FeedStage as stream (Brief="Feed stream", PosX=0, PosY=0.55);
2098
2099	STAGES(NumberOfStages) as packedStage;
2100	PressureDrop as pressure;
2101
2102	SET
2103	FeedStageIndex(FeedStageLocation) =1;
2104
2105	ColumnArea = ColumnDiameter^23.141590.25;
2106	hs =PackingHeight/NumberOfStages;
2107	V = hs * ColumnArea;
2108
2109	CONNECTIONS
2110
2111	STAGES([2:NumberOfStages]).OutletVapour to STAGES([1:NumberOfStages-1]).InletVapour;
2112	STAGES([1:NumberOfStages-1]).OutletLiquid to STAGES([2:NumberOfStages]).InletLiquid;
2113
2114	EQUATIONS
2115	# Connecting STAGES
2116	FeedStage.F*FeedStageIndex= STAGES.Inlet.F;
2117	FeedStage.T = STAGES.Inlet.T;
2118	FeedStage.P = STAGES.Inlet.P;
2119	FeedStage.z = STAGES.Inlet.z;
2120	FeedStage.v = STAGES.Inlet.v;
2121	FeedStage.h = STAGES.Inlet.h;
2122
2123	STAGES.deltaP = PressureDrop/NumberOfStages;
2124
2125	for i in [1:NumberOfStages] do
2126
2127	"Energy Balance"
2128	diff(STAGES(i).E) = ( STAGES(i).Inlet.FSTAGES(i).Inlet.h + STAGES(i).InletLiquid.FSTAGES(i).InletLiquid.h + STAGES(i).InletVapour.FSTAGES(i).InletVapour.h- STAGES(i).OutletLiquid.FSTAGES(i).OutletLiquid.h
2129	- STAGES(i).OutletVapour.F*STAGES(i).OutletVapour.h + HeatSupply );
2130
2131	"Energy Holdup"
2132	STAGES(i).E = STAGES(i).MLSTAGES(i).OutletLiquid.h + STAGES(i).MVSTAGES(i).OutletVapour.h - STAGES(i).OutletLiquid.P*V;
2133
2134	"Geometry Constraint"
2135	VVoidFraction= STAGES(i).MLSTAGES(i).vL + STAGES(i).MV*STAGES(i).vV;
2136
2137	"Volume flow rate of vapor, m^3/m^2/s"
2138	STAGES(i).uV * ((ColumnArea)VoidFraction - STAGES(i).Al) = STAGES(i).OutletVapour.F STAGES(i).vV;
2139
2140	"Liquid holdup"
2141	STAGES(i).hl = STAGES(i).ML*STAGES(i).vL/V/VoidFraction;
2142
2143	"Particle diameter"
2144	STAGES(i).dp = 6 * (1-VoidFraction)/AreaPerPackingVol;
2145
2146	"Wall Factor"
2147	STAGES(i).invK = (1 + (2STAGES(i).dp/(3ColumnDiameter*(1-VoidFraction))));
2148
2149	"Reynolds number of the vapor stream"
2150	STAGES(i).RevSTAGES(i).invK = STAGES(i).dpSTAGES(i).uVSTAGES(i).rhoV / (STAGES(i).miV(1-VoidFraction));
2151
2152	"Pressure drop and Vapor flow"
2153	STAGES(i).deltaP/hs = ResistanceCoeffAreaPerPackingVolSTAGES(i).uV^2STAGES(i).rhoVSTAGES(i).invK / (2*(VoidFraction-STAGES(i).hl)^3);
2154
2155	"Liquid holdup"
2156	STAGES(i).hl = (12STAGES(i).miLAreaPerPackingVol^2*STAGES(i).uL/STAGES(i).rhoL/g)^1/3;
2157
2158	"Area occupied by the liquid"
2159	STAGES(i).Al = STAGES(i).ML*STAGES(i).vL/hs;
2160
2161	end
2162
2163	INITIAL
2164
2165	for i in 1:NumberOfStages do
2166
2167	"The initial temperature of the STAGES"
2168	STAGES(i).OutletLiquid.T = INITIALIZATION.TopStageTemperature+(INITIALIZATION.BottomStageTemperature-INITIALIZATION.TopStageTemperature)*((i-1)/(NumberOfStages-1));
2169
2170	"The Liquid Holdup of the STAGES"
2171	STAGES(i).ML = INITIALIZATION.LiquidMolarHoldup;
2172
2173	end
2174
2175	for i in 1:NComp-1 do
2176
2177	for j in 1:NumberOfStages do
2178
2179	"The initial composition of the TRAYS"
2180	STAGES(j).OutletLiquid.z(i) = INITIALIZATION.TopStageComposition(i)/sum(INITIALIZATION.TopStageComposition) +(INITIALIZATION.BottomStageComposition(i)/sum(INITIALIZATION.BottomStageComposition)-INITIALIZATION.TopStageComposition(i)/sum(INITIALIZATION.TopStageComposition) )*((j-1)/(NumberOfStages-1));
2181
2182	end
2183
2184	end
2185
2186	end
2187
2188	Model Packed_Section_Column as Packed_Section_ColumnBasic
2189
2190	ATTRIBUTES
2191	Pallete = true;
2192	Icon = "icon/PackedSectionColumn";
2193	Brief = "Model of a packed column section.";
2194	Info =
2195	"== Model of a packed column section containing ==
2196	* NStages theoretical stages.
2197
2198	== Specify ==
2199	* the feed stream of each tray (Inlet);
2200	* the InletLiquid stream of the top tray;
2201	* the InletVapour stream of the bottom tray;
2202	* the total pressure drop (dP) of the section.
2203
2204	== Initial Conditions ==
2205	* the stages temperature (OutletLiquid.T);
2206	* the stages liquid holdup;
2207	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray.
2208	";
2209
2210	VARIABLES
2211
2212	in LiquidInlet as stream (Brief="Liquid Inlet in the section", PosX=0.70, PosY=0);
2213	out VapourOutlet as vapour_stream (Brief="Vapour Outlet in the section", PosX=0.30, PosY=0);
2214
2215	in VapourInlet as stream (Brief="Vapour Inlet in the section", PosX=0.30, PosY=1);
2216	out LiquidOutlet as liquid_stream (Brief="Liquid Outlet in the section", PosX=0.70, PosY=1);
2217
2218	LiquidConnector as stream (Brief="Liquid connection at the middle STAGES", PosX=0.75, PosY=1,Hidden=true);
2219	VapourConnector as stream (Brief="Vapour connection at the middle STAGES", PosX=0.55, PosY=0,Hidden=true);
2220
2221	CONNECTIONS
2222
2223	LiquidConnector to STAGES(1).InletLiquid;
2224	VapourConnector to STAGES(NumberOfStages).InletVapour;
2225
2226	EQUATIONS
2227
2228	LiquidConnector.F= LiquidInlet.F;
2229	LiquidConnector.T = LiquidInlet.T;
2230	LiquidConnector.P = LiquidInlet.P;
2231	LiquidConnector.z = LiquidInlet.z;
2232	LiquidConnector.v = LiquidInlet.v;
2233	LiquidConnector.h = LiquidInlet.h;
2234
2235	VapourConnector.F= VapourInlet.F;
2236	VapourConnector.T = VapourInlet.T;
2237	VapourConnector.P = VapourInlet.P;
2238	VapourConnector.z = VapourInlet.z;
2239	VapourConnector.v = VapourInlet.v;
2240	VapourConnector.h = VapourInlet.h;
2241
2242	LiquidOutlet.F= STAGES(NumberOfStages).OutletLiquid.F;
2243	LiquidOutlet.T = STAGES(NumberOfStages).OutletLiquid.T;
2244	LiquidOutlet.P = STAGES(NumberOfStages).OutletLiquid.P;
2245	LiquidOutlet.z = STAGES(NumberOfStages).OutletLiquid.z;
2246
2247	VapourOutlet.F= STAGES(1).OutletVapour.F;
2248	VapourOutlet.T = STAGES(1).OutletVapour.T;
2249	VapourOutlet.P = STAGES(1).OutletVapour.P;
2250	VapourOutlet.z = STAGES(1).OutletVapour.z;
2251
2252	end
2253
2254	Model PackedDistillation_kettle_cond as Packed_Section_ColumnBasic
2255	ATTRIBUTES
2256	Pallete = true;
2257	Icon = "icon/PackedDistillationKettleCond";
2258	Brief = "Model of a distillation column with dynamic condenser and dynamic reboiler.";
2259	Info =
2260	"== Specify ==
2261	* the feed stream of each tray (Inlet);
2262	* the pump pressure difference;
2263	* the total pressure drop (dP) of the packing;
2264	* the heat supllied in reboiler and condenser;
2265	* the condenser vapor outlet flow (OutletVapour.F);
2266	* the reboiler liquid outlet flow (OutletLiquid.F);
2267	* both splitter outlet flows OR one of the splitter outlet flows and the splitter frac.
2268
2269	== Initial Conditions ==
2270	* the stages temperature (OutletLiquid.T);
2271	* the stages initial molar holdup;
2272	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each stage;
2273
2274	* the condenser temperature (OutletLiquid.T);
2275	* the condenser liquid level (Level);
2276	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions;
2277
2278	* the reboiler temperature (OutletLiquid.T);
2279	* the reboiler liquid level (Level);
2280	* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
2281	";
2282
2283	PARAMETERS
2284
2285	VapourFlow as Switcher (Valid = ["on", "off"], Hidden = true ,Default = "on");
2286
2287	VARIABLES
2288	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2289	# Column Devices
2290	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2291	CONDENSER as condenser (Brief="dymamic condenser");
2292	REBOILER as reboiler (Brief="kettle reboiler");
2293	SPLITTER as splitter2 (Brief="splitter which separate reflux and distillate");
2294	PUMP as pump (Brief="pump in reflux stream");
2295	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2296	# Heat Ports
2297	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2298	in HeatToReboiler as power (Brief="Heat supplied to Reboiler",Protected = true, PosX=1, PosY=0.885);
2299	in HeatToCondenser as power (Brief="Heat supplied to Condenser", Protected = true, PosX=1, PosY=0.034);
2300	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2301	# Column Products
2302	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2303	out VapourDistillate as vapour_stream (Brief="Vapour outlet stream From Top Condenser", PosX=0.66, PosY=0);
2304	out LiquidDistillate as liquid_stream (Brief="Liquid outlet stream From Top Splitter", PosX=1, PosY=0.20);
2305	out BottomProduct as liquid_stream (Brief="Liquid outlet stream From Reboiler", PosX=0.68, PosY=1);
2306	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2307	# Column Connectors
2308	#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2309	in ConnectorCondenserVout as stream (Brief="Connector for Vapour outlet stream From Top Condenser", Hidden=true);
2310	in ConnectorSplitterOut as stream (Brief="Connector for Liquid outlet stream From Top Splitter", Hidden=true);
2311	in ConnectorReboilerLout as stream (Brief="Connector for Liquid outlet stream From Reboiler", Hidden=true);
2312	ConnectorHeatReboiler as power (Brief="Heat supplied to Reboiler", Hidden=true);
2313	ConnectorHeatCondenser as power (Brief="Heat supplied to Condenser", Hidden=true);
2314
2315	CONNECTIONS
2316
2317	#vapor
2318	REBOILER.OutletVapour to STAGES(NumberOfStages).InletVapour;
2319	STAGES(1).OutletVapour to CONDENSER.InletVapour;
2320
2321	#liquid
2322	CONDENSER.OutletLiquid to SPLITTER.Inlet;
2323	SPLITTER.Outlet2 to PUMP.Inlet;
2324	PUMP.Outlet to STAGES(1).InletLiquid;
2325	STAGES(NumberOfStages).OutletLiquid to REBOILER.InletLiquid;
2326
2327	#Connectors
2328	ConnectorHeatReboiler to REBOILER.InletQ;
2329	ConnectorHeatCondenser to CONDENSER.InletQ;
2330
2331	CONDENSER.OutletVapour to ConnectorCondenserVout;
2332	SPLITTER.Outlet1 to ConnectorSplitterOut;
2333	REBOILER.OutletLiquid to ConnectorReboilerLout;
2334
2335	EQUATIONS
2336	HeatToReboiler = ConnectorHeatReboiler;
2337	HeatToCondenser = ConnectorHeatCondenser;
2338
2339	# Condenser Connector Equations
2340	ConnectorCondenserVout.T = VapourDistillate.T;
2341	ConnectorCondenserVout.P = VapourDistillate.P;
2342	ConnectorCondenserVout.F = VapourDistillate.F;
2343	ConnectorCondenserVout.z = VapourDistillate.z;
2344
2345	# Splitter Connector Equations
2346	ConnectorSplitterOut.T = LiquidDistillate.T;
2347	ConnectorSplitterOut.P = LiquidDistillate.P;
2348	ConnectorSplitterOut.F = LiquidDistillate.F;
2349	ConnectorSplitterOut.z = LiquidDistillate.z;
2350
2351	# Reboiler Connector Equations
2352	ConnectorReboilerLout.T = BottomProduct.T;
2353	ConnectorReboilerLout.P = BottomProduct.P;
2354	ConnectorReboilerLout.F = BottomProduct.F;
2355	ConnectorReboilerLout.z = BottomProduct.z;
2356
2357	switch VapourFlow
2358	case "on":
2359	STAGES(NumberOfStages).InletVapour.F*STAGES(NumberOfStages).vV = sqrt((REBOILER.OutletVapour.P - STAGES(NumberOfStages).OutletVapour.P)/
2360	(STAGES(NumberOfStages).rhoVResistanceCoeff20))(ColumnAreaVoidFraction - STAGES(NumberOfStages).Al);
2361	when STAGES(NumberOfStages).InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
2362
2363	case "off":
2364	STAGES(NumberOfStages).InletVapour.F = 0 * 'mol/s';
2365	when REBOILER.OutletVapour.P > STAGES(NumberOfStages).OutletVapour.P + 1e-1 * 'atm' switchto "on";
2366	end
2367
2368	end
2369

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format