source: branches/gui/eml/costs/reboiler_cost.mso @ 918

Last change on this file since 918 was 918, checked in by Rafael de Pelegrini Soares, 12 years ago

Fixed some flash and cost models

File size: 3.1 KB
RevLine 
[918]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* Model of costs for a reboiler
17*--------------------------------------------------------------------
18*
19*       Streams:
20*               * a liquid inlet stream
21*               * a liquid outlet stream
22*               * a vapour outlet stream
23*               * a feed stream
24*
25*       Assumptions:
26*               * perfect mixing of both phases
27*               * thermodynamics equilibrium
28*               * no liquid entrainment in the vapour stream
29*
30*       Specify:
31*               * the Feed stream
32*               * the Liquid inlet stream
33*               * the outlet flows: OutletV.F and OutletL.F
34*
35*       Initial:
36*               * the reboiler temperature (OutletL.T)
37*               * the reboiler liquid level (Ll)
38*               * (NoComps - 1) OutletL (OR OutletV) compositions
39*
40*       - This Model is valid if 14m^2 < Across < 1100 m^2
41*
42*----------------------------------------------------------------------
43* Author: Núbia do Carmo Ferreira
44* $Id: reboiler.mso 210 2007-03-15 12:52:28Z arge $
45*--------------------------------------------------------------------*#
46
47using "stage_separators/reboiler";
48
49Model reboiler_cost as reboiler
50        ATTRIBUTES
51        Pallete         = true;
52        Icon            = "icon/Reboiler";
53
54        PARAMETERS
55        Material as Switcher (Valid = ["Stainless steel 316", "Stainless steel 304", "Stainless steel 347", "Nickel 200", "Monel 400", "Inconel 600", "Incoloy 825", "Titanium", "Hastelloy"],
56                        Default = "Stainless steel 316");
57        Cost(8,3) as Real;
58       
59        VARIABLES
60        Ce as currency  (Brief="Capital Cost");
61        Cb as currency  (Brief="Basic Cost");
62        Fd as positive  (Brief="Cost Factor based on the type of the heat exchanger");
63        Fp as positive  (Brief="Cost Factor based on the project pressure");
64        Fm as positive  (Brief="Cost Factor based on the construction material");
65        P  as pressure  (Brief="Average pressure");
66
67        EQUATIONS
68        "Average pressure"
69        P = 0.5*(InletLiquid.P + OutletLiquid.P);
70       
71        "Capital Cost"
72        Ce = Cb*Fd*Fp*Fm;
73
74        "Basic Cost"
75        Cb = 'US$'*exp(Cost(1,1) + Cost(1,2)*ln(Across/'m^2') + Cost(1,3)*(ln(Across/'m^2'))^2);
76
77        "Cost Factor based on the type of the heat exchanger"
78        Fd = exp(Cost(3,1) + Cost(3,2)*ln(Across/'m^2') + Cost(3,3)*ln(Across/'m^2'));
79
80        if P <= 700 * 'kPa' then
81                "Cost Factor based on the project pressure" # verificar
82                Fp = Cost(5,1) + Cost(5,2)*ln(Across/'m^2');
83        else if 700 * 'kPa' < P and P < 2100 * 'kPa' then
84                Fp = Cost(5,1) + Cost(5,2)*ln(Across/'m^2');
85        else if 2100 * 'kPa' < P and P < 4200 * 'kPa' then
86                Fp = Cost(6,1) + Cost(6,2)*ln(Across/'m^2');
87        else if 4200 * 'kPa' < P and P < 6200 * 'kPa' then
88                Fp = Cost(7,1) + Cost(7,2)*ln(Across/'m^2');
89        else
90                Fp = Cost(7,1) + Cost(7,2)*ln(Across/'m^2'); # verificar
91        end
92        end
93        end
94        end
95
96        "Cost Factor based on the construction material"
97        Fm = Cost(8,1) + Cost(8,2)*ln(Across/'m^2');
98       
99end
Note: See TracBrowser for help on using the repository browser.