source: trunk/eml/heat_exchangers/Mheatex.mso @ 840

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

improved documentation of Mheatex model

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 3.3 KB
RevLine 
[78]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.
[1]14*----------------------------------------------------------------------
15* Author: Gerson Balbueno Bicca
16* $Id: Mheatex.mso 544 2008-06-27 15:23:00Z bicca $
17*--------------------------------------------------------------------*#
[45]18
[164]19using "HEX_Engine.mso";
[1]20
21Model Mheatex
[135]22       
23ATTRIBUTES
[307]24        Icon = "icon/Mheatex";
[139]25        Pallete = true;
[284]26        Brief = "Model of a Multistream heat exchangers";
[359]27        Info            =
28"Heat Transfer between multiple hot and cold streams.
29
30== Specify ==
[544]31The user must specify the following variables:
[359]32* The Inlet Hot streams: Any Number
33* The Inlet Cold streams: Any Number
34* The LMTD correction factor
[544]35* The outlet temperature and the outlet pressure for each stream on one side (hot or cold) of the heat exchanger
36* On the other side the user must leave at least one unespecified outlet stream.
[359]37
38== Setting The Mheatex Parameters ==
[544]39* FlowDirection: counter or cocurrent flow
40* Number of hot streams (Nhot)
41* Number of cold streams (Ncold)
[359]42";
43
[1]44PARAMETERS
45
[325]46outer PP                        as Plugin       (Brief="Physical Properties", Type="PP");
47outer NComp             as Integer      (Brief="Number of Components");
48        FlowDirection   as Switcher     (Brief="Flow Direction",Valid=["counter","cocurrent"],Default="cocurrent");
[359]49        Ncold                   as Integer      (Brief="Number of Inlet Cold Streams",Lower=1, Symbol = "N_{cold}");
50        Nhot                    as Integer      (Brief="Number of Inlet Hot Streams",Lower=1, Symbol = "N_{hot}");
[45]51       
[1]52VARIABLES
53
[544]54in  InletHot(Nhot)      as stream       (Brief="Inlet Hot Streams", PosX=0, PosY=0.7156, Symbol="_{hot}^{in}");
55out OutletHot(Nhot)     as streamPH (Brief="Outlet Hot Streams", PosX=1, PosY=0.7156, Symbol="_{hot}^{out}");
56in  InletCold(Ncold)    as stream       (Brief="Inlet Cold Streams", PosX=1, PosY=0.2793, Symbol="_{cold}^{in}");
57out OutletCold(Ncold)   as streamPH     (Brief="Outlet Cold Streams", PosX=0, PosY=0.2793, Symbol="_{cold}^{out}");
[1]58
[359]59        Method  as LMTD_Basic   (Brief="Log Mean Temperature Difference Method", Symbol = " ");
[325]60        Q               as power        (Brief="Heat Transfer", Default=7000, Lower=1e-6, Upper=1e10);
61        UA              as Real     (Brief="UA product",Unit='W/K',Lower=1e-8);
[131]62       
[1]63EQUATIONS
64
65"Hot Flow"
[164]66        OutletHot.F = InletHot.F;       
[1]67
68"Cold Flow"
[164]69        OutletCold.F = InletCold.F;
[1]70       
71"Hot Composition"
[164]72        OutletHot.z = InletHot.z;
[1]73       
74"Cold Composition"
[164]75        OutletCold.z = InletCold.z;
[1]76       
77"Heat Duty Hot Stream"
[164]78        Q =  sum(InletHot.F*(InletHot.h- OutletHot.h));
[1]79       
80"Heat Duty Cold Stream"
[164]81        Q = -sum(InletCold.F*(InletCold.h- OutletCold.h));
[1]82       
83"Heat Duty"
[168]84        Q = UA*Method.LMTD*Method.Fc;
[1]85
[144]86switch FlowDirection
[1]87       
[140]88case    "cocurrent":
89
[131]90"Temperature Difference at Inlet"
[164]91        Method.DT0 = max(InletHot.T) - min(InletCold.T);
[131]92
93"Temperature Difference at Outlet"
[164]94        Method.DTL = min(OutletHot.T) - max(OutletCold.T);
[131]95       
[140]96case "counter":
[1]97       
[131]98"Temperature Difference at Inlet"
[164]99        Method.DT0 = max(InletHot.T) - max(OutletCold.T);
[131]100
101"Temperature Difference at Outlet"
[164]102        Method.DTL = min(OutletHot.T) - min(InletCold.T);
[131]103       
[1]104end
105
106end
Note: See TracBrowser for help on using the repository browser.