source: branches/gui/sample/costs/sample_distillation_kettle_cond_cost.mso @ 941

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

costs models are obsoletes

File size: 3.9 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* Sample file for column cost model
17*--------------------------------------------------------------------
18*
19* This sample file needs VRTherm DEMO(www.vrtech.com.br) to run
20* SectionColumn_Test and needs VRTherm full to run the distillation
21* column flowsheet.
22*
23*----------------------------------------------------------------------
24* Author: Núbia do Carmo Ferreira
25* $Id: sample_column.mso 247 2007-04-24 13:44:18Z rafael $
26*--------------------------------------------------------------------#
27
28
29
30Needs to be Updated
31
32
33using "costs/column_cost";
34
35FlowSheet sample_Distillation_kettle_cond_cost
36       
37        PARAMETERS
38        PP      as Plugin(Brief="Physical Properties",
39                Type="PP",
40                Components = [ "isobutane", "n-pentane", "propylene",
41                "benzene", "isobutene" ],
42                LiquidModel = "PR",
43                VapourModel = "PR"
44        );
45        NComp   as Integer;
46       
47       
48        SET
49        NComp = PP.NumberOfComponents;
50
51        DEVICES
52        col as Distillation_kettle_cond_cost;
53        feed as source;
54        zero as stream;
55        Qc as energy_source;
56        Qr as energy_source;
57       
58        CONNECTIONS
59        feed.Outlet to col.trays(5).Inlet;
60        zero            to col.reb.Inlet;
61        zero            to col.trays([1:4]).Inlet;
62        zero            to col.trays([6:col.NTrays]).Inlet;
63        Qc.OutletQ      to col.cond.InletQ;
64        Qr.OutletQ      to col.reb.InletQ;
65       
66        SPECIFY
67        feed.F = 113.4 * 'kmol/h';
68        feed.T = 291 * 'K';
69        feed.P = 168.3 * 'kPa';
70        feed.Composition = 1/NComp;
71       
72        zero.F = 0 * 'kmol/h';
73        zero.T = 300 * 'K';
74        zero.P = 1 * 'atm';
75        zero.z = 1/NComp;
76        zero.v = 0;
77        zero.h = 0 * 'J/mol';
78       
79        col.sptop.Outlet2.F = 85 * 'kmol/h';
80        col.reb.OutletL.F = 28.4 * 'kmol/h';
81        col.sptop.frac = 0.444445;
82        col.cond.OutletV.F = 0 * 'kmol/h';
83        Qr.OutletQ = 3.7743e6 * 'kJ/h';
84        Qc.OutletQ = -3.71e6 * 'kJ/h';
85        col.pump1.dP = 16 * 'kPa';
86        col.trays.Emv = 1;
87        col.alfaTopo = 2;
88       
89        SET
90        col.NTrays = 8;
91        col.cond.V = 2 * 'm^3';
92        col.cond.Across = 1 * 'm^2';
93        col.trays.V = 4 * 'ft^3';
94        col.trays.Ah = 0.394 * 'ft^2';
95        col.trays.lw = 20.94 * 'in';
96        col.trays.hw = 0.125 * 'ft';
97        col.trays.Q = 0 * 'kW';
98        col.trays.beta = 0.6;
99        col.trays.alfa = 4;
100        col.trays.Ap = 3.94 * 'ft^2';
101        col.reb.V = 2 * 'm^3';
102        col.reb.Across = 1 * 'm^2';
103       
104        #cost
105        col.Material = "Stainless steel 316";
106        col.Tray_Material = "Stainless steel 304";
107        col.Tray_Type = "Valve";
108        col.Di = 0.9 * 'm';
109        col.Lt = 17.5 * 'm';
110        col.Tb = 0.01250 * 'm';
111        col.Tp = 0.01349 * 'm';
112        col.D  = 2.4 * 'm';
113        col.dens_mass_material = 8000 * 'kg/m^3';
114
115        col.Cost(1,:) = [6.950,0.1808,0.02468,0.01580];
116        col.Cost(2,:) = [834.86,0,0,0];
117        col.Cost(3,:) = [6.448,0.21887,0.02297,0];
118        col.Cost(4,:) = [1017.0,0,0,0];
119        col.Cost(5,:) = [1.7,0,0,0]; #For "Stainless steel 304" - tray material
120        col.Cost(6,:) = [278.38,0.1739,0,0];
121        col.Cost(7,:) = [1.189,0.1894,0,0];
122        col.Cost(8,:) = [1.401,0.2376,0,0];
123        col.Cost(9,:) = [1.525,0.2585,0,0];
124        col.Cost(10,:) = [2.306,0.3674,0,0];
125        col.Cost(11,:) = [2.25,1.0414,0,0];
126        col.Cost(12,:) = [1,0,0,0]; #For valve trays
127       
128        INITIAL
129        # condenser
130        col.cond.OutletL.T = 260 *'K';
131        col.cond.Level = 1 * 'm';
132        col.cond.OutletL.z([1:4]) = [0.65, 0.05, 0.01, 0.01];
133
134        # reboiler
135        col.reb.OutletL.T = 300 *'K';
136        col.reb.Level = 1 * 'm';
137        col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01];
138
139        # column trays
140        col.trays.OutletL.T = [290:(300-290)/(col.NTrays-1):300] * 'K';
141        col.trays.Level = 1.2 * col.trays.hw;
142        col.trays.OutletL.z([1:4]) = [0.5, 0.05, 0.01, 0.01];
143
144        OPTIONS
145        TimeStep = 0.1;
146        TimeEnd = 50;
147        #Dynamic = false;       
148end
149*#
Note: See TracBrowser for help on using the repository browser.