source: mso/eml/stage_separators/condenser.mso @ 37

Last change on this file since 37 was 1, checked in by Rafael de Pelegrini Soares, 16 years ago

Initial import of the library

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 3.1 KB
Line 
1#*-------------------------------------------------------------------
2* Model of a dynamic condenser
3*--------------------------------------------------------------------
4*
5*       Streams:
6*               * a vapour inlet stream
7*               * a liquid outlet stream
8*
9*       Assumptions:
10*               * perfect mixing of both phases
11*               * thermodynamics equilibrium
12*
13*       Specify:
14*               * the Inlet stream
15*               * the Outlet flows
16*
17*       Initial:
18*               * the condenser temperature (OutletL.T)
19*               * the condenser level (Ll)
20*               * (NoComps - 1) Outlet compositions
21*
22*----------------------------------------------------------------------
23* Author: Paula B. Staudt
24* $Id: condenser.mso 1 2006-06-20 17:33:53Z rafael $
25*--------------------------------------------------------------------*#
26
27using "streams";
28
29Model condenser
30        PARAMETERS
31ext PP as CalcObject;
32ext NComp as Integer;
33        V as volume (Brief="Condenser total volume");
34        Across as area (Brief="Cross Section Area of reboiler");
35
36        VARIABLES
37in      InletV as stream; #(Brief="Vapour inlet stream");
38out     OutletL as stream_therm; #(Brief="Liquid outlet stream");
39out     OutletV as stream_therm; #(Brief="Vapour outlet stream");
40in      Q as heat_rate (Brief="Heat supplied");
41
42        M(NComp) as mol (Brief="Molar Holdup in the tray");
43        ML as mol (Brief="Molar liquid holdup");
44        MV as mol (Brief="Molar vapour holdup");
45        E as energy (Brief="Total Energy Holdup on tray");
46        vL as volume_mol (Brief="Liquid Molar Volume");
47        vV as volume_mol (Brief="Vapour Molar volume");
48        Level as length (Brief="Level of liquid phase");
49
50        EQUATIONS
51        "Component Molar Balance"
52        diff(M) = InletV.F*InletV.z - OutletL.F*OutletL.z
53                                - OutletV.F*OutletV.z;
54
55        "Energy Balance"
56        diff(E) = InletV.F*InletV.h - OutletL.F*OutletL.h
57                                - OutletV.F*OutletV.h + Q;
58
59        "Molar Holdup"
60        M = ML*OutletL.z + MV*OutletV.z;
61       
62        "Energy Holdup"
63        E = ML*OutletL.h + MV*OutletV.h - OutletV.P*V;
64       
65        "Mol fraction normalisation"
66        sum(OutletL.z)=1.0;
67        sum(OutletL.z)=sum(OutletV.z);
68
69        "Liquid Volume"
70        vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
71        "Vapour Volume"
72        vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
73
74        "Chemical Equilibrium"
75        PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z =
76                PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;
77
78        "Thermal Equilibrium"
79        OutletL.T = OutletV.T;
80
81        "Mechanical Equilibrium"
82        OutletV.P = OutletL.P;
83
84        "Geometry Constraint"
85        V = ML*vL + MV*vV;
86
87        "Level of liquid phase"
88        Level = ML*vL/Across;
89       
90        "Vapourisation Fraction"
91        OutletL.v = 0.0;
92        OutletV.v = 1.0;
93end
94
95
96#*----------------------------------------------------------------------
97* Model of a  Steady State condenser with no thermodynamics equilibrium
98*---------------------------------------------------------------------*#
99Model condenserSteady
100        PARAMETERS
101ext PP as CalcObject;
102ext NComp as Integer;
103        DP as press_delta (Brief="Pressure Drop in the condenser");
104
105        VARIABLES
106in      InletV as stream; #(Brief="Vapour inlet stream");
107out     OutletL as stream_therm; #(Brief="Liquid outlet stream");
108in      Q as heat_rate (Brief="Heat supplied");
109
110        EQUATIONS
111        "Molar Balance"
112        InletV.F = OutletL.F;
113        InletV.z = OutletL.z;
114               
115        "Energy Balance"
116        InletV.F*InletV.h = OutletL.F*OutletL.h + Q;
117       
118        "Pressure"
119        DP = InletV.P - OutletL.P;
120
121        "Vapourisation Fraction"
122        OutletL.v = 0.0;
123end
Note: See TracBrowser for help on using the repository browser.