source: trunk/eml/costs/column_cost.mso @ 987

Last change on this file since 987 was 732, checked in by gerson bicca, 14 years ago

costs models are obsoletes

File size: 18.3 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* Author: Núbia do Carmo Ferreira
17* $Id: column.mso 210 2007-03-15 12:52:28Z arge $
18*---------------------------------------------------------------------#
19
20using "stage_separators/column";
21
22
23
24
25                                Needs to be Updated !!!!!!!
26
27
28#----------------------------------------------------------------------
29* Model of cost of a distillation column containing:
30*       - NTrays like tray;
31*       - a kettle reboiler;
32*       - dymamic condenser;
33*       - a splitter which separate reflux and distillate;
34*       - a pump in reflux stream;
35*
36*       - This model is valid if:
37*          0.91m < Di < 7.32m
38*          17.57m < Lt < 51.82m
39*          0.6m < D < 4.8m
40*
41*---------------------------------------------------------------------#
42
43Model Distillation_kettle_cond_cost as Distillation_kettle_cond
44        ATTRIBUTES
45        Pallete         = true;
46        Icon            = "icon/DistillationKettleCond";
47        Brief           = "Model of a distillation column with dynamic condenser and dynamic reboiler.";
48        Info            =
49"== Specify ==
50* the feed stream of each tray (Inlet);
51* the Murphree eficiency for each tray Emv;
52* the pump pressure difference;
53* the heat supllied in reboiler and condenser;
54* the condenser vapor outlet flow (OutletV.F);
55* the reboiler liquid outlet flow (OutletL.F);
56* both splitter outlet flows OR one of the splitter outlet flows and the splitter frac.
57* all necessary dimensions and materials for cost evaluation
58       
59== Initial Conditions ==
60* the trays temperature (OutletL.T);
61* the trays liquid level (Level) OR the trays liquid flow (OutletL.F);
62* (NoComps - 1) OutletL (OR OutletV) compositions for each tray;
63       
64* the condenser temperature (OutletL.T);
65* the condenser liquid level (Level);
66* (NoComps - 1) OutletL (OR OutletV) compositions;
67       
68* the reboiler temperature (OutletL.T);
69* the reboiler liquid level (Level);
70* (NoComps - 1) OutletL (OR OutletV) compositions.
71";
72
73        PARAMETERS
74        Material as Switcher (Valid = ["Stainless steel 316", "Stainless steel 304", "Carpenter 20 CB_3", "Nickel 200", "Monel 400", "Inconel 600", "Incoloy 825", "Titanium"],
75                        Default = "Stainless steel 304");
76        Tray_Type as Switcher (Valid = ["Valve", "Grid", "Bubble cap", "Sieve"], Default= "Valve");
77        Tray_Material as Switcher (Valid = ["Stainless steel 304", "Stainless steel 316", "Carpenter 20 CB_3", "Monel"], Default= "Stainless steel 304");
78        Cost(12,4) as Real;
79        Di as length    (Brief="Internal Diameter for the hoof of the towers");
80        Lt as length    (Brief="Hoof Length");
81        Tb as length    (Brief="Bottom thickness");
82        Tp as length    (Brief="Wall thickness");       
83        D  as length    (Brief="Tower diameter");
84        dens_mass_material  as dens_mass (Brief = "Mass Density of the Material");
85        Pi      as constant    (Brief="Pi Number",Default=3.14159265);
86       
87        VARIABLES
88        Cs as currency  (Brief="Total Cost of the hoof of the towers");
89        Cb as currency  (Brief="Base Cost for the hoof of the distillation tower");
90        Fm as positive  (Brief="Cost Factor based on the construction material");
91        Cpl as currency (Brief="Cost for stairs and platform");
92        Cbt as currency (Brief="Base Cost for the trays of the tower");
93        Ftm as positive (Brief="Cost Factor based on the construction material of the trays");
94        Fnt as positive (Brief="Cost Factor based on the number of trays in the tower");
95        Ftt as positive (Brief="Cost Factor based on the type of the tray");
96        Ct as currency  (Brief="Total Cost");
97        Ws as mass              (Brief="Equipment Weight");
98       
99        EQUATIONS
100
101        "Total Cost of the hoof of the towers"
102        Cs = Cb*Fm;
103
104        "Base Cost for the hoof of the distillation tower"
105        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Ws/'kg') + Cost(1,3)*(ln(Ws/'kg'))^2 + Cost(1,4)*(Lt/Di)*(ln(Tb/Tp)));
106
107        "Cost for stairs and platform"
108        Cpl = 'US$'*Cost(2,1)*(Di^0.63316)/'m^0.63316'*(Lt^0.80161)/'m^0.80161';
109       
110        "Base Cost for the trays of the tower"
111        Cbt = 'US$'*Cost(6,1)*exp(Cost(6,2)*D/'m');
112
113        switch Tray_Material
114                case "Stainless steel 304":
115
116        "Cost Factor based on the construction material of the trays"
117        Ftm = Cost(7,1) + Cost(7,2)*D/'m';
118       
119                case "Stainless steel 316":
120               
121        "Cost Factor based on the construction material of the trays"
122        Ftm = Cost(8,1) + Cost(8,2)*D/'m';             
123
124                case "Carpenter 20 CB_3":
125       
126        "Cost Factor based on the construction material of the trays"
127        Ftm = Cost(9,1) + Cost(9,2)*D/'m';
128
129                case "Monel":
130
131        "Cost Factor based on the construction material of the trays"
132        Ftm = Cost(10,1) + Cost(10,2)*D/'m';
133
134        end
135
136        "Cost Factor based on the number of trays in the tower"
137        Fnt = Cost(11,1)/(Cost(11,2))^NTrays;
138
139        "Total Cost"
140        Ct = Cb*Fm + NTrays*Cbt*Ftm*Ftt*Fnt + Cpl;
141
142        "Cost Factor based on the type of the tray"
143        Ftt = Cost(12,1);
144
145        "Cost Factor based on the construction material"
146        Fm = Cost(5,1);
147
148        "Equipment Weight"
149        Ws = dens_mass_material*Di*(Lt + 0.8116*Di)*Tp*Pi;
150end
151
152#* -------------------------------------------------------------------
153* Distillation Column cost model with:
154*
155*       - NTrays like tray;
156*       - a vessel in the bottom of column;
157*       - a splitter who separate the bottom product and the stream to reboiler;
158*       - steady state reboiler (thermosyphon);
159*       - a steady state condenser with subcooling;
160*       - a vessel drum (layed cilinder);
161*       - a splitter which separate reflux and distillate;
162*       - a pump in reflux stream.
163*
164* ------------------------------------------------------------------#
165Model Distillation_thermosyphon_subcooling_cost as Distillation_thermosyphon_subcooling
166        ATTRIBUTES
167        Pallete         = true;
168        Icon            = "icon/DistillationThermosyphonSubcooling";
169        Brief           = "Model of a distillation column with steady condenser and steady reboiler.";
170        Info            =
171"== Specify ==
172* the feed stream of each tray (Inlet);
173* the Murphree eficiency for each tray Emv;
174* the pump head;
175* the condenser pressure drop;
176* the heat supllied in top and bottom tanks;
177* the heat supllied in condenser and reboiler;
178* the Outlet1 flow in the bottom splitter (spbottom.Outlet1.F) that corresponds to the bottom product;
179* both  top splitter outlet flows OR one of the splitter outlet flows and the splitter frac.
180* all necessary dimensions and materials for cost evaluation
181
182== Initial Conditions ==
183* the trays temperature (OutletL.T);
184* the trays liquid level (Level) OR the trays liquid flow (OutletL.F);
185* (NoComps - 1) OutletL (OR OutletV) compositions for each tray;
186       
187* the top tank temperature (OutletL.T);
188* the top tank liquid level (Level);
189* (NoComps - 1) OutletL (OR OutletV) compositions;
190       
191* the bottom tank temperature (OutletL.T);
192* the bottom tank liquid level (Level);
193* (NoComps - 1) OutletL (OR OutletV) compositions.
194";
195
196        PARAMETERS
197        Material as Switcher (Valid = ["Stainless steel 316", "Stainless steel 304", "Carpenter 20 CB_3", "Nickel 200", "Monel 400", "Inconel 600", "Incoloy 825", "Titanium"],
198                        Default = "Stainless steel 304");
199        Tray_Type as Switcher (Valid = ["Valve", "Grid", "Bubble cap", "Sieve"], Default= "Valve");
200        Tray_Material as Switcher (Valid = ["Stainless steel 304", "Stainless steel 316", "Carpenter 20 CB_3", "Monel"], Default= "Stainless steel 304");
201        Cost(12,4) as Real;
202
203        Di as length    (Brief="Internal Diameter for the hoof of the towers");
204        Lt as length    (Brief="Hoof Length");
205        Tb as length    (Brief="Bottom thickness");
206        Tp as length    (Brief="Wall thickness");       
207        D  as length    (Brief="Tower diameter");
208        dens_mass_material  as dens_mass (Brief = "Mass Density of the Material");
209        Pi      as constant    (Brief="Pi Number",Default=3.14159265);
210       
211        VARIABLES
212        Cs as currency  (Brief="Total Cost of the hoof of the towers");
213        Cb as currency  (Brief="Base Cost for the hoof of the distillation tower");
214        Fm as positive  (Brief="Cost Factor based on the construction material");
215        Cpl as currency (Brief="Cost for stairs and platform");
216        Cbt as currency (Brief="Base Cost for the trays of the tower");
217        Ftm as positive (Brief="Cost Factor based on the construction material of the trays");
218        Fnt as positive (Brief="Cost Factor based on the number of trays in the tower");
219        Ftt as positive (Brief="Cost Factor based on the type of the tray");
220        Ct as currency  (Brief="Total Cost");
221        Ws as mass              (Brief="Equipment Weight");
222
223        EQUATIONS
224
225        "Total Cost of the hoof of the towers"
226        Cs = Cb*Fm;
227
228        "Base Cost for the hoof of the distillation tower"
229        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Ws/'kg') + Cost(1,3)*(ln(Ws/'kg'))^2 + Cost(1,4)*(Lt/Di)*(ln(Tb/Tp)));
230
231        "Cost for stairs and platform"
232        Cpl = 'US$'*Cost(2,1)*(Di^0.63316)/'m^0.63316'*(Lt^0.80161)/'m^0.80161';
233       
234        "Base Cost for the trays of the tower"
235        Cbt = 'US$'*Cost(6,1)*exp(Cost(6,2)*D/'m');
236
237        switch Tray_Material
238                case "Stainless steel 304":
239
240        "Cost Factor based on the construction material of the trays"
241        Ftm = Cost(7,1) + Cost(7,2)*D/'m';
242       
243                case "Stainless steel 316":
244               
245        "Cost Factor based on the construction material of the trays"
246        Ftm = Cost(8,1) + Cost(8,2)*D/'m';             
247
248                case "Carpenter 20 CB_3":
249       
250        "Cost Factor based on the construction material of the trays"
251        Ftm = Cost(9,1) + Cost(9,2)*D/'m';
252
253                case "Monel":
254
255        "Cost Factor based on the construction material of the trays"
256        Ftm = Cost(10,1) + Cost(10,2)*D/'m';
257
258        end
259
260        "Cost Factor based on the number of trays in the tower"
261        Fnt = Cost(11,1)/(Cost(11,2))^NTrays;
262
263        "Total Cost"
264        Ct = Cb*Fm + NTrays*Cbt*Ftm*Ftt*Fnt + Cpl;
265
266        "Cost Factor based on the type of the tray"
267        Ftt = Cost(12,1);
268
269        "Cost Factor based on the construction material"
270        Fm = Cost(5,1);
271
272        "Equipment Weight"
273        Ws = dens_mass_material*Di*(Lt + 0.8116*Di)*Tp*Pi;
274end
275
276#* -------------------------------------------------------------------
277* Distillation Column model with:
278*
279*       - NTrays like tray;
280*       - a vessel in the bottom of column;
281*       - a splitter who separate the bottom product and the stream to reboiler;
282*       - steady state reboiler (thermosyphon);
283*       - a dynamic condenser without subcooling;
284*       - a splitter which separate reflux and distillate;
285*       - a pump in reflux stream.
286*
287* ------------------------------------------------------------------#
288Model Distillation_thermosyphon_cond_cost as Distillation_thermosyphon_cond
289        ATTRIBUTES
290        Pallete         = true;
291        Icon            = "icon/DistillationThermosyphonCond";
292        Brief           = "Model of a distillation column with dynamic condenser and steady reboiler.";
293        Info            =
294"== Specify ==
295* the feed stream of each tray (Inlet);
296* the Murphree eficiency for each tray Emv;
297* the pump head;
298* the condenser vapor outlet flow (OutletV.F);
299* the heat supllied in bottom tank;
300* the heat supllied in condenser and reboiler;
301* the Outlet1 flow in the bottom splitter (spbottom.Outlet1.F) that corresponds to the bottom product;
302* all necessary dimensions and materials for cost evaluation
303
304== Initial Conditions ==
305* the trays temperature (OutletL.T);
306* the trays liquid level (Level) OR the trays liquid flow (OutletL.F);
307* (NoComps - 1) OutletL (OR OutletV) compositions for each tray;
308       
309* the condenser temperature (OutletL.T);
310* the condenser liquid level (Level);
311* (NoComps - 1) OutletL (OR OutletV) compositions;
312       
313* the bottom tank temperature (OutletL.T);
314* the bottom tank liquid level (Level);
315* (NoComps - 1) OutletL (OR OutletV) compositions.
316";
317
318        PARAMETERS
319        Material as Switcher (Valid = ["Stainless steel 316", "Stainless steel 304", "Carpenter 20 CB_3", "Nickel 200", "Monel 400", "Inconel 600", "Incoloy 825", "Titanium"],
320                        Default = "Stainless steel 304");
321        Tray_Type as Switcher (Valid = ["Valve", "Grid", "Bubble cap", "Sieve"], Default= "Valve");
322        Tray_Material as Switcher (Valid = ["Stainless steel 304", "Stainless steel 316", "Carpenter 20 CB_3", "Monel"], Default= "Stainless steel 304");
323        Cost(12,4) as Real;
324        Di as length    (Brief="Internal Diameter for the hoof of the towers");
325        Lt as length    (Brief="Hoof Length");
326        Tb as length    (Brief="Bottom thickness");
327        Tp as length    (Brief="Wall thickness");       
328        D  as length    (Brief="Tower diameter");
329        dens_mass_material  as dens_mass (Brief = "Mass Density of the Material");
330        Pi      as constant    (Brief="Pi Number",Default=3.14159265);
331       
332        VARIABLES
333        Cs as currency  (Brief="Total Cost of the hoof of the towers");
334        Cb as currency  (Brief="Base Cost for the hoof of the distillation tower");
335        Fm as positive  (Brief="Cost Factor based on the construction material");
336        Cpl as currency (Brief="Cost for stairs and platform");
337        Cbt as currency (Brief="Base Cost for the trays of the tower");
338        Ftm as positive (Brief="Cost Factor based on the construction material of the trays");
339        Fnt as positive (Brief="Cost Factor based on the number of trays in the tower");
340        Ftt as positive (Brief="Cost Factor based on the type of the tray");
341        Ct as currency  (Brief="Total Cost");
342        Ws as mass              (Brief="Equipment Weight");
343       
344        EQUATIONS
345
346        "Total Cost of the hoof of the towers"
347        Cs = Cb*Fm;
348
349        "Base Cost for the hoof of the distillation tower"
350        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Ws/'kg') + Cost(1,3)*(ln(Ws/'kg'))^2 + Cost(1,4)*(Lt/Di)*(ln(Tb/Tp)));
351
352        "Cost for stairs and platform"
353        Cpl = 'US$'*Cost(2,1)*(Di^0.63316)/'m^0.63316'*(Lt^0.80161)/'m^0.80161';
354       
355        "Base Cost for the trays of the tower"
356        Cbt = 'US$'*Cost(6,1)*exp(Cost(6,2)*D/'m');
357
358        switch Tray_Material
359                case "Stainless steel 304":
360
361        "Cost Factor based on the construction material of the trays"
362        Ftm = Cost(7,1) + Cost(7,2)*D/'m';
363       
364                case "Stainless steel 316":
365               
366        "Cost Factor based on the construction material of the trays"
367        Ftm = Cost(8,1) + Cost(8,2)*D/'m';             
368
369                case "Carpenter 20 CB_3":
370       
371        "Cost Factor based on the construction material of the trays"
372        Ftm = Cost(9,1) + Cost(9,2)*D/'m';
373
374                case "Monel":
375
376        "Cost Factor based on the construction material of the trays"
377        Ftm = Cost(10,1) + Cost(10,2)*D/'m';
378
379        end
380
381        "Cost Factor based on the number of trays in the tower"
382        Fnt = Cost(11,1)/(Cost(11,2))^NTrays;
383
384        "Total Cost"
385        Ct = Cb*Fm + NTrays*Cbt*Ftm*Ftt*Fnt + Cpl;
386
387        "Cost Factor based on the type of the tray"
388        Ftt = Cost(12,1);
389
390        "Cost Factor based on the construction material"
391        Fm = Cost(5,1);
392
393        "Equipment Weight"
394        Ws = dens_mass_material*Di*(Lt + 0.8116*Di)*Tp*Pi;
395end
396
397#* -------------------------------------------------------------------
398* Distillation Column model with:
399*
400*       - NTrays like tray;
401*       - a kettle reboiler;
402*       - a steady state condenser with subcooling;
403*       - a vessel drum (layed cilinder);
404*       - a splitter which separate reflux and distillate;
405*       - a pump in reflux stream.
406*
407* ------------------------------------------------------------------#
408Model Distillation_kettle_subcooling_cost as Distillation_kettle_subcooling
409        ATTRIBUTES
410        Pallete         = true;
411        Icon            = "icon/DistillationKettleSubcooling";
412        Brief           = "Model of a distillation column with steady condenser and dynamic reboiler.";
413        Info            =
414"== Specify ==
415* the feed stream of each tray (Inlet);
416* the Murphree eficiency for each tray (Emv);
417* the pump pressure difference;
418* the heat supllied in reboiler and condenser;
419* the heat supllied in the top tank;
420* the condenser pressure drop;
421* the reboiler liquid outlet flow (OutletL.F);
422* both splitter outlet flows OR one of the splitter outlet flows and the splitter frac.
423* all necessary dimensions and materials for cost evaluation
424
425== Initial Conditions ==
426* the trays temperature (OutletL.T);
427* the trays liquid level (Level) OR the trays liquid flow (OutletL.F);
428* (NoComps - 1) OutletL (OR OutletV) compositions for each tray;
429       
430* the top tank temperature (OutletL.T);
431* the top tank liquid level (Level);
432* (NoComps - 1) OutletL (OR OutletV) compositions;
433       
434* the reboiler temperature (OutletL.T);
435* the reboiler liquid level (Level);
436* (NoComps - 1) OutletL (OR OutletV) compositions.
437";
438
439        PARAMETERS
440        Material as Switcher (Valid = ["Stainless steel 316", "Stainless steel 304", "Carpenter 20 CB_3", "Nickel 200", "Monel 400", "Inconel 600", "Incoloy 825", "Titanium"],
441                        Default = "Stainless steel 304");
442        Tray_Type as Switcher (Valid = ["Valve", "Grid", "Bubble cap", "Sieve"], Default= "Valve");
443        Tray_Material as Switcher (Valid = ["Stainless steel 304", "Stainless steel 316", "Carpenter 20 CB_3", "Monel"], Default= "Stainless steel 304");
444        Cost(12,4) as Real;
445        Di as length    (Brief="Internal Diameter for the hoof of the towers");
446        Lt as length    (Brief="Hoof Length");
447        Tb as length    (Brief="Bottom thickness");
448        Tp as length    (Brief="Wall thickness");       
449        D  as length    (Brief="Tower diameter");
450        dens_mass_material  as dens_mass (Brief = "Mass Density of the Material");
451        Pi      as constant    (Brief="Pi Number",Default=3.14159265);
452       
453        VARIABLES
454        Cs as currency  (Brief="Total Cost of the hoof of the towers");
455        Cb as currency  (Brief="Base Cost for the hoof of the distillation tower");
456        Fm as positive  (Brief="Cost Factor based on the construction material");
457        Cpl as currency (Brief="Cost for stairs and platform");
458        Cbt as currency (Brief="Base Cost for the trays of the tower");
459        Ftm as positive (Brief="Cost Factor based on the construction material of the trays");
460        Fnt as positive (Brief="Cost Factor based on the number of trays in the tower");
461        Ftt as positive (Brief="Cost Factor based on the type of the tray");
462        Ct as currency  (Brief="Total Cost");
463        Ws as mass              (Brief="Equipment Weight");
464       
465        EQUATIONS
466
467        "Total Cost of the hoof of the towers"
468        Cs = Cb*Fm;
469
470        "Base Cost for the hoof of the distillation tower"
471        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Ws/'kg') + Cost(1,3)*(ln(Ws/'kg'))^2 + Cost(1,4)*(Lt/Di)*(ln(Tb/Tp)));
472
473        "Cost for stairs and platform"
474        Cpl = 'US$'*Cost(2,1)*(Di^0.63316)/'m^0.63316'*(Lt^0.80161)/'m^0.80161';
475       
476        "Base Cost for the trays of the tower"
477        Cbt = 'US$'*Cost(6,1)*exp(Cost(6,2)*D/'m');
478
479        switch Tray_Material
480                case "Stainless steel 304":
481
482        "Cost Factor based on the construction material of the trays"
483        Ftm = Cost(7,1) + Cost(7,2)*D/'m';
484       
485                case "Stainless steel 316":
486               
487        "Cost Factor based on the construction material of the trays"
488        Ftm = Cost(8,1) + Cost(8,2)*D/'m';             
489
490                case "Carpenter 20 CB_3":
491       
492        "Cost Factor based on the construction material of the trays"
493        Ftm = Cost(9,1) + Cost(9,2)*D/'m';
494
495                case "Monel":
496
497        "Cost Factor based on the construction material of the trays"
498        Ftm = Cost(10,1) + Cost(10,2)*D/'m';
499
500        end
501
502        "Cost Factor based on the number of trays in the tower"
503        Fnt = Cost(11,1)/(Cost(11,2))^NTrays;
504
505        "Total Cost"
506        Ct = Cb*Fm + NTrays*Cbt*Ftm*Ftt*Fnt + Cpl;
507
508        "Cost Factor based on the type of the tray"
509        Ftt = Cost(12,1);
510
511        "Cost Factor based on the construction material"
512        Fm = Cost(5,1);
513
514        "Equipment Weight"
515        Ws = dens_mass_material*Di*(Lt + 0.8116*Di)*Tp*Pi;
516end
517
518*#
Note: See TracBrowser for help on using the repository browser.