Ignore:
Timestamp:
Feb 26, 2009, 11:00:39 PM (13 years ago)
Author:
gerson bicca
Message:

updates (some samples are obsoletes)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/gui/eml/stage_separators/condenser.mso

    r721 r735  
    3232== Specify ==
    3333* the inlet stream;
    34 * the outlet flows: OutletV.F and OutletL.F;
     34* the outlet flows: OutletVapour.F and OutletLiquid.F;
    3535* the heat supply.
    3636       
    3737== Initial Conditions ==
    38 * the condenser temperature (OutletL.T);
     38* the condenser temperature (OutletLiquid.T);
    3939* the condenser liquid level (Level);
    40 * (NoComps - 1) OutletL (OR OutletV) compositions.
     40* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
    4141";     
    4242       
     
    5353       
    5454VARIABLES
    55 in              InletV          as stream                               (Brief="Vapour inlet stream", PosX=0.15, PosY=0, Symbol="_{inV}");
    56 out     OutletL         as liquid_stream                (Brief="Liquid outlet stream", PosX=0.4513, PosY=1, Symbol="_{outL}");
    57 out     OutletV         as vapour_stream        (Brief="Vapour outlet stream", PosX=0.4723, PosY=0, Symbol="_{outV}");
     55in              InletVapour             as stream                               (Brief="Vapour inlet stream", PosX=0.15, PosY=0, Symbol="_{inV}");
     56out     OutletLiquid    as liquid_stream                (Brief="Liquid outlet stream", PosX=0.4513, PosY=1, Symbol="_{outL}");
     57out     OutletVapour    as vapour_stream        (Brief="Vapour outlet stream", PosX=0.4723, PosY=0, Symbol="_{outV}");
    5858in              InletQ          as power                                (Brief="Cold supplied", PosX=1, PosY=0, Symbol="_{in}");
    5959
     
    6969
    7070        Level                                   = Initial_Level;
    71         OutletL.T                               = Initial_Temperature;
    72         OutletL.z(1:NComp-1)    = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
     71        OutletLiquid.T                          = Initial_Temperature;
     72        OutletLiquid.z(1:NComp-1)       = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
    7373       
    7474EQUATIONS
    7575"Component Molar Balance"
    76         diff(M) = InletV.F*InletV.z - OutletL.F*OutletL.z- OutletV.F*OutletV.z;
     76        diff(M) = InletVapour.F*InletVapour.z - OutletLiquid.F*OutletLiquid.z- OutletVapour.F*OutletVapour.z;
    7777
    7878"Energy Balance"
    79         diff(E) = InletV.F*InletV.h - OutletL.F*OutletL.h- OutletV.F*OutletV.h + InletQ;
     79        diff(E) = InletVapour.F*InletVapour.h - OutletLiquid.F*OutletLiquid.h- OutletVapour.F*OutletVapour.h + InletQ;
    8080
    8181"Molar Holdup"
    82         M = ML*OutletL.z + MV*OutletV.z;
     82        M = ML*OutletLiquid.z + MV*OutletVapour.z;
    8383       
    8484"Energy Holdup"
    85         E = ML*OutletL.h + MV*OutletV.h - OutletV.P*V;
     85        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletVapour.P*V;
    8686       
    8787"Mol fraction normalisation"
    88         sum(OutletL.z)=1.0;
    89         sum(OutletL.z)=sum(OutletV.z);
     88        sum(OutletLiquid.z)=1.0;
     89        sum(OutletLiquid.z)=sum(OutletVapour.z);
    9090
    9191"Liquid Volume"
    92         vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
     92        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    9393       
    9494"Vapour Volume"
    95         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
     95        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    9696
    9797"Chemical Equilibrium"
    98         PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z =
    99                 PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;
     98        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z =
     99                PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z;
    100100
    101101"Thermal Equilibrium"
    102         OutletL.T = OutletV.T;
     102        OutletLiquid.T = OutletVapour.T;
    103103
    104104"Mechanical Equilibrium"
    105         OutletV.P = OutletL.P;
     105        OutletVapour.P = OutletLiquid.P;
    106106
    107107"Geometry Constraint"
     
    138138
    139139VARIABLES
    140 in      InletV          as stream                               (Brief="Vapour inlet stream", PosX=0.3431, PosY=0, Symbol="_{inV}");
    141 out     OutletL as liquid_stream                (Brief="Liquid outlet stream", PosX=0.34375, PosY=1, Symbol="_{outL}");
     140in      InletVapour             as stream                               (Brief="Vapour inlet stream", PosX=0.3431, PosY=0, Symbol="_{inV}");
     141out     OutletLiquid as liquid_stream           (Brief="Liquid outlet stream", PosX=0.34375, PosY=1, Symbol="_{outL}");
    142142in      InletQ          as power                                (Brief="Cold supplied", PosX=1, PosY=0.5974, Symbol="_{in}");
    143143        DP                      as press_delta          (Brief="Pressure Drop in the condenser",Default=0);
     
    146146
    147147"Molar Balance"
    148         InletV.F = OutletL.F;
    149         InletV.z = OutletL.z;
     148        InletVapour.F = OutletLiquid.F;
     149        InletVapour.z = OutletLiquid.z;
    150150
    151151"Energy Balance"
    152         InletV.F*InletV.h = OutletL.F*OutletL.h + InletQ;
     152        InletVapour.F*InletVapour.h = OutletLiquid.F*OutletLiquid.h + InletQ;
    153153
    154154"Pressure"
    155         DP = InletV.P - OutletL.P;
     155        DP = InletVapour.P - OutletLiquid.P;
    156156
    157157end
     
    174174* the reaction related variables;
    175175* the inlet stream;
    176 * the outlet flows: OutletV.F and OutletL.F;
     176* the outlet flows: OutletVapour.F and OutletLiquid.F;
    177177* the heat supply.
    178178
    179179== Initial Conditions ==
    180 * the condenser temperature (OutletL.T);
     180* the condenser temperature (OutletLiquid.T);
    181181* the condenser liquid level (Level);
    182 * (NoComps - 1) OutletL (OR OutletV) compositions.
     182* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
    183183";
    184184       
     
    198198VARIABLES
    199199
    200 in      InletV          as stream                       (Brief="Vapour inlet stream", PosX=0.1164, PosY=0, Symbol="_{inV}");
    201 out     OutletL         as liquid_stream        (Brief="Liquid outlet stream", PosX=0.4513, PosY=1, Symbol="_{outL}");
    202 out     OutletV         as vapour_stream        (Brief="Vapour outlet stream", PosX=0.4723, PosY=0, Symbol="_{outV}");
     200in      InletVapour             as stream                       (Brief="Vapour inlet stream", PosX=0.1164, PosY=0, Symbol="_{inV}");
     201out     OutletLiquid    as liquid_stream        (Brief="Liquid outlet stream", PosX=0.4513, PosY=1, Symbol="_{outL}");
     202out     OutletVapour    as vapour_stream        (Brief="Vapour outlet stream", PosX=0.4723, PosY=0, Symbol="_{outV}");
    203203        InletQ          as power                        (Brief="Cold supplied", PosX=1, PosY=0.6311, Symbol="_{in}");
    204204
     
    217217
    218218        Level                                   = Initial_Level;
    219         OutletL.T                               = Initial_Temperature;
    220         OutletL.z(1:NComp-1)    = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
     219        OutletLiquid.T                          = Initial_Temperature;
     220        OutletLiquid.z(1:NComp-1)       = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
    221221
    222222EQUATIONS
    223223"Molar Concentration"
    224         OutletL.z = vL * C;
     224        OutletLiquid.z = vL * C;
    225225       
    226226"Reaction"
    227         r3 = exp(-7150*'K'/OutletL.T)*(4.85e4*C(1)*C(2) - 1.23e4*C(3)*C(4)) * 'l/mol/s';
     227        r3 = exp(-7150*'K'/OutletLiquid.T)*(4.85e4*C(1)*C(2) - 1.23e4*C(3)*C(4)) * 'l/mol/s';
    228228       
    229229"Component Molar Balance"
    230         diff(M) = InletV.F*InletV.z - OutletL.F*OutletL.z - OutletV.F*OutletV.z + stoic*r3*ML*vL;
     230        diff(M) = InletVapour.F*InletVapour.z - OutletLiquid.F*OutletLiquid.z - OutletVapour.F*OutletVapour.z + stoic*r3*ML*vL;
    231231
    232232"Energy Balance"
    233         diff(E) = InletV.F*InletV.h - OutletL.F*OutletL.h- OutletV.F*OutletV.h + InletQ + Hr * r3 * ML*vL;
     233        diff(E) = InletVapour.F*InletVapour.h - OutletLiquid.F*OutletLiquid.h- OutletVapour.F*OutletVapour.h + InletQ + Hr * r3 * ML*vL;
    234234
    235235"Molar Holdup"
    236         M = ML*OutletL.z + MV*OutletV.z;
     236        M = ML*OutletLiquid.z + MV*OutletVapour.z;
    237237       
    238238"Energy Holdup"
    239         E = ML*OutletL.h + MV*OutletV.h - OutletV.P*V;
     239        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletVapour.P*V;
    240240       
    241241"Mol fraction normalisation"
    242         sum(OutletL.z)=1.0;
     242        sum(OutletLiquid.z)=1.0;
    243243
    244244"Liquid Volume"
    245         vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
     245        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    246246
    247247"Vapour Volume"
    248         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
     248        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    249249
    250250"Thermal Equilibrium"
    251         OutletL.T = OutletV.T;
     251        OutletLiquid.T = OutletVapour.T;
    252252
    253253"Mechanical Equilibrium"
    254         OutletV.P = OutletL.P;
     254        OutletVapour.P = OutletLiquid.P;
    255255
    256256"Geometry Constraint"
     
    263263       
    264264"Chemical Equilibrium"
    265         PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z =
    266         PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;
    267 
    268         sum(OutletL.z)=sum(OutletV.z);
     265        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z =
     266        PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z;
     267
     268        sum(OutletLiquid.z)=sum(OutletVapour.z);
    269269
    270270end
Note: See TracChangeset for help on using the changeset viewer.