source: trunk/sample/heat_exchangers/Sample_hairpin.mso @ 623

Last change on this file since 623 was 484, checked in by gerson bicca, 15 years ago

Hairpin: added option to solve convergence problems from a simplified model
File size: 5.4 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* Hairpin Heat Exchanger - NTU and LMTD Method
17*--------------------------------------------------------------------
18*
19* This sample file needs VRTherm (www.vrtech.com.br) to run.
20*
21*----------------------------------------------------------------------
22* Author: Gerson B. Bicca
23* $Id: Sample_hairpin.mso  $
24*--------------------------------------------------------------------*#
25
26using "heat_exchangers/Hairpin.mso";
27
28FlowSheet LMTD_Method
29
30PARAMETERS
31        PP                      as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
32                                                                        Components = ["n-hexane","water"]);
33        NComp           as Integer;
34       
35DEVICES
36
37        exchanger       as Hairpin_LMTD;
38        InletHot        as simple_source;
39        InletCold       as simple_source;
40        OutletHot       as sink;
41        OutletCold      as sink;
42       
43       
44CONNECTIONS
45
46        InletHot.Outlet                 to exchanger.InletInner;
47        InletCold.Outlet                        to exchanger.InletOuter;
48        exchanger.OutletOuter   to OutletCold.Inlet;
49        exchanger.OutletInner   to OutletHot.Inlet;
50       
51SET
52
53NComp   = PP.NumberOfComponents;
54
55#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
56#   Options
57#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
58
59exchanger.FlowDirection                                         = "counter";
60exchanger.InnerTurbulentCorrelation     = "SiederTate";
61
62exchanger.OuterTurbulentCorrelation     = "SiederTate";
63exchanger.OuterLaminarCorrelation       = "Schlunder";
64
65#exchanger.CalculationApproach = "Simplified";
66exchanger.CalculationApproach = "Full";
67exchanger.Qestimated                                    = 4*'kW';
68#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
69# Hairpin Geometrical Parameters and Alocation
70#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
71
72exchanger.DoInner               =       0.04826*'m';
73exchanger.DiInner               =       0.04089*'m';
74exchanger.DiOuter               =       0.07793*'m';
75exchanger.Kwall                 =       0.057 *'kW/m/K';
76exchanger.Lpipe                 =       3*'m';
77
78#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
79#       Fouling
80#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
81exchanger.Rfi = 0.0018*'m^2*K/kW';
82exchanger.Rfo = 0.0021*'m^2*K/kW';
83
84SPECIFY
85#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
86#               Hot Stream
87#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
88        InletHot.Outlet.F = 10                          * 'kmol/h';
89        InletHot.Outlet.T = (68+273.15)   * 'K';
90        InletHot.Outlet.P = 5.4                         * 'bar';
91        InletHot.Outlet.z = [1,0];
92       
93#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
94#       Cold Stream
95#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
96        InletCold.Outlet.F = 5                          * 'kmol/h';
97        InletCold.Outlet.P = 5.1                * 'bar';
98        InletCold.Outlet.T = (23+273.15)        * 'K';
99    InletCold.Outlet.z = [0,1];
100       
101OPTIONS
102
103Dynamic  = false;
104GuessFile = "GuessSampleHairpin";
105
106end
107
108FlowSheet NTU_Method
109
110PARAMETERS
111
112        PP                      as Plugin (Type ="PP", LiquidModel = "PR", VapourModel = "PR",
113                                                                        Components = ["n-hexane","water"]);
114        NComp           as Integer;
115       
116DEVICES
117
118        exchanger       as Hairpin_NTU;
119        InletHot        as simple_source;
120        InletCold       as simple_source;
121        OutletHot       as sink;
122        OutletCold      as sink;
123       
124       
125CONNECTIONS
126
127        InletHot.Outlet                 to exchanger.InletInner;
128        InletCold.Outlet                        to exchanger.InletOuter;
129        exchanger.OutletOuter   to OutletCold.Inlet;
130        exchanger.OutletInner   to OutletHot.Inlet;
131       
132SET
133
134NComp   = PP.NumberOfComponents;
135
136#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
137#   Options
138#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
139
140exchanger.FlowDirection  = "counter";
141exchanger.InnerTurbulentCorrelation  = "SiederTate";
142
143exchanger.OuterTurbulentCorrelation  = "SiederTate";
144exchanger.OuterLaminarCorrelation     = "Schlunder";
145
146#exchanger.CalculationApproach = "Simplified";
147exchanger.CalculationApproach = "Full";
148
149exchanger.Qestimated                                    = 2*'kW';
150exchanger.Eftestimated                          = 0.4;
151#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
152# Hairpin Geometrical Parameters and Alocation
153#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
154
155exchanger.DoInner               =       0.04826*'m';
156exchanger.DiInner               =       0.04089*'m';
157exchanger.DiOuter               =       0.07793*'m';
158exchanger.Kwall                 =       0.057 *'kW/m/K';
159exchanger.Lpipe                 =       3*'m';
160
161#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
162#       Fouling
163#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
164exchanger.Rfi = 0.0018*'m^2*K/kW';
165exchanger.Rfo = 0.0021*'m^2*K/kW';
166
167SPECIFY
168#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
169#               Hot Stream
170#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
171        InletHot.Outlet.F = 10                          * 'kmol/h';
172        InletHot.Outlet.T = (68+273.15)   * 'K';
173        InletHot.Outlet.P = 5.4                         * 'bar';
174        InletHot.Outlet.z = [1,0];
175       
176#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
177#       Cold Stream
178#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
179        InletCold.Outlet.F = 5                          * 'kmol/h';
180        InletCold.Outlet.P = 5.1                * 'bar';
181        InletCold.Outlet.T = (23+273.15)        * 'K';
182    InletCold.Outlet.z = [0,1];
183
184OPTIONS
185
186Dynamic = false;
187
188GuessFile = "GuessPipe";
189
190end
Note: See TracBrowser for help on using the repository browser.