source: branches/gui/eml/costs/column_cost.mso @ 554

Last change on this file since 554 was 354, checked in by Argimiro Resende Secchi, 15 years ago

Fixing more wiki notation.

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