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/reboiler.mso

    r721 r735  
    3636* the inlet stream;
    3737* the liquid inlet stream;
    38 * the outlet flows: OutletV.F and OutletL.F;
     38* the outlet flows: OutletVapour.F and OutletLiquid.F;
    3939* the heat supply.
    4040       
    4141== Initial Conditions ==
    4242
    43 * the reboiler temperature (OutletL.T);
     43* the reboiler temperature (OutletLiquid.T);
    4444* the reboiler liquid level (Level);
    45 * (NoComps - 1) OutletL (OR OutletV) compositions.
     45* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
    4646";     
    4747       
    … …  
    5858VARIABLES
    5959
    60 in      InletL                  as stream                               (Brief="Liquid inlet stream", PosX=0.15, PosY=1, Symbol="_{inL}");
    61 out     OutletL         as liquid_stream                (Brief="Liquid outlet stream", PosX=0.40, PosY=1, Symbol="_{outL}");
    62 out     OutletV         as vapour_stream        (Brief="Vapour outlet stream", PosX=0.40, PosY=0, Symbol="_{outV}");
     60in      InletLiquid                     as stream                               (Brief="Liquid inlet stream", PosX=0.15, PosY=1, Symbol="_{inL}");
     61out     OutletLiquid    as liquid_stream                (Brief="Liquid outlet stream", PosX=0.40, PosY=1, Symbol="_{outL}");
     62out     OutletVapour    as vapour_stream        (Brief="Vapour outlet stream", PosX=0.40, PosY=0, Symbol="_{outV}");
    6363in      InletQ                  as power                                (Brief="Heat supplied", PosX=1, PosY=0, Symbol="_{in}");
    6464
    … …  
    7575
    7676        Level                                   = Initial_Level;
    77         OutletL.T                               = Initial_Temperature;
    78         OutletL.z(1:NComp-1)    = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
     77        OutletLiquid.T                          = Initial_Temperature;
     78        OutletLiquid.z(1:NComp-1)       = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
    7979
    8080EQUATIONS
    8181"Component Molar Balance"
    82         diff(M)= InletL.F*InletL.z      - OutletL.F*OutletL.z - OutletV.F*OutletV.z;
     82        diff(M)= InletLiquid.F*InletLiquid.z    - OutletLiquid.F*OutletLiquid.z - OutletVapour.F*OutletVapour.z;
    8383       
    8484"Energy Balance"
    85         diff(E) = InletL.F*InletL.h     - OutletL.F*OutletL.h - OutletV.F*OutletV.h + InletQ;
     85        diff(E) = InletLiquid.F*InletLiquid.h   - OutletLiquid.F*OutletLiquid.h - OutletVapour.F*OutletVapour.h + InletQ;
    8686       
    8787"Molar Holdup"
    88         M = ML*OutletL.z + MV*OutletV.z;
     88        M = ML*OutletLiquid.z + MV*OutletVapour.z;
    8989       
    9090"Energy Holdup"
    91         E = ML*OutletL.h + MV*OutletV.h - OutletL.P*V;
     91        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletLiquid.P*V;
    9292       
    9393"Mol fraction normalisation"
    94         sum(OutletL.z)=1.0;
    95         sum(OutletL.z)=sum(OutletV.z);
     94        sum(OutletLiquid.z)=1.0;
     95        sum(OutletLiquid.z)=sum(OutletVapour.z);
    9696
    9797"Vapour Density"
    98         rhoV = PP.VapourDensity(OutletV.T, OutletV.P, OutletV.z);
     98        rhoV = PP.VapourDensity(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    9999
    100100"Liquid Volume"
    101         vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
     101        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    102102       
    103103"Vapour Volume"
    104         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
     104        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    105105       
    106106"Chemical Equilibrium"
    107         PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z = PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;
     107        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z = PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z;
    108108
    109109"Mechanical Equilibrium"
    110         OutletL.P = OutletV.P;
     110        OutletLiquid.P = OutletVapour.P;
    111111       
    112112"Thermal Equilibrium"
    113         OutletL.T = OutletV.T;
     113        OutletLiquid.T = OutletVapour.T;
    114114       
    115115"Geometry Constraint"
    … …  
    136136       
    137137== Specify ==
    138 * the InletL stream;
    139 * the heat supply OR the outlet temperature (OutletV.T);
     138* the InletLiquid stream;
     139* the heat supply OR the outlet temperature (OutletVapour.T);
    140140";     
    141141
    … …  
    146146
    147147VARIABLES
    148 in      InletL                  as stream                                       (Brief="Liquid inlet stream", PosX=0.3345, PosY=1, Symbol="_{inL}");
    149 out     OutletV         as vapour_stream                (Brief="Vapour outlet stream", PosX=0.3369, PosY=0, Symbol="_{outV}");
     148in      InletLiquid                     as stream                                       (Brief="Liquid inlet stream", PosX=0.3345, PosY=1, Symbol="_{inL}");
     149out     OutletVapour    as vapour_stream                (Brief="Vapour outlet stream", PosX=0.3369, PosY=0, Symbol="_{outV}");
    150150in      InletQ                  as power                                        (Brief="Heat supplied", PosX=1, PosY=0.6111, Symbol="_{in}");
    151151        vV                              as volume_mol                   (Brief="Vapour Molar volume");
    … …  
    155155
    156156"Molar Balance"
    157         InletL.F = OutletV.F;
    158         InletL.z = OutletV.z;
     157        InletLiquid.F = OutletVapour.F;
     158        InletLiquid.z = OutletVapour.z;
    159159       
    160160"Vapour Volume"
    161         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
     161        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    162162       
    163163"Vapour Density"
    164         rhoV = PP.VapourDensity(OutletV.T, OutletV.P, OutletV.z);
     164        rhoV = PP.VapourDensity(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    165165
    166166"Energy Balance"
    167         InletL.F*InletL.h + InletQ = OutletV.F*OutletV.h;
     167        InletLiquid.F*InletLiquid.h + InletQ = OutletVapour.F*OutletVapour.h;
    168168       
    169169"Pressure"
    170         DP = InletL.P - OutletV.P;
     170        DP = InletLiquid.P - OutletVapour.P;
    171171
    172172end
    … …  
    195195       
    196196        VARIABLES
    197 in      InletL as stream(Brief="Liquid inlet stream", PosX=0.3345, PosY=1, Symbol="_{inL}");
    198 out     OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.3369, PosY=0, Symbol="_{outV}");
     197in      InletLiquid as stream(Brief="Liquid inlet stream", PosX=0.3345, PosY=1, Symbol="_{inL}");
     198out     OutletVapour as vapour_stream(Brief="Vapour outlet stream", PosX=0.3369, PosY=0, Symbol="_{outV}");
    199199in      InletQ as power (Brief="Heat supplied", PosX=1, PosY=0.6111, Symbol="_{in}");
    200200
    201201        EQUATIONS
    202202        "Molar Balance"
    203         InletL.F = OutletV.F;
    204         InletL.z = OutletV.z;
     203        InletLiquid.F = OutletVapour.F;
     204        InletLiquid.z = OutletVapour.z;
    205205       
    206206        "Energy Balance"
    207         InletL.F*InletL.h + InletQ = OutletV.F*OutletV.h;
     207        InletLiquid.F*InletLiquid.h + InletQ = OutletVapour.F*OutletVapour.h;
    208208       
    209209        "Pressure"
    210         DP = InletL.P - OutletV.P;
     210        DP = InletLiquid.P - OutletVapour.P;
    211211
    212212        "Fake Vapourisation Fraction"
    213         OutletV.v = 1.0;
     213        OutletVapour.v = 1.0;
    214214       
    215215        "Fake output temperature"
    216         OutletV.T = 300*'K';
     216        OutletVapour.T = 300*'K';
    217217       
    218218        "Pressure Drop through the reboiler"
    219         OutletV.F = k*InletQ;
     219        OutletVapour.F = k*InletQ;
    220220end
    221221
    … …  
    239239* the inlet stream;
    240240* the liquid inlet stream;
    241 * the outlet flows: OutletV.F and OutletL.F;
     241* the outlet flows: OutletVapour.F and OutletLiquid.F;
    242242* the heat supply.
    243243
    244244== Initial Conditions ==
    245 * the reboiler temperature (OutletL.T);
     245* the reboiler temperature (OutletLiquid.T);
    246246* the reboiler liquid level (Level);
    247 * (NoComps - 1) OutletL (OR OutletV) compositions.
     247* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions.
    248248";
    249249       
    … …  
    262262       
    263263VARIABLES
    264 in      InletL  as stream                       (Brief="Liquid inlet stream", PosX=0, PosY=0.5254, Symbol="_{inL}");
    265 out     OutletL as liquid_stream        (Brief="Liquid outlet stream", PosX=0.2413, PosY=1, Symbol="_{outL}");
    266 out     OutletV as vapour_stream        (Brief="Vapour outlet stream", PosX=0.5079, PosY=0, Symbol="_{outV}");
     264in      InletLiquid     as stream                       (Brief="Liquid inlet stream", PosX=0, PosY=0.5254, Symbol="_{inL}");
     265out     OutletLiquid as liquid_stream   (Brief="Liquid outlet stream", PosX=0.2413, PosY=1, Symbol="_{outL}");
     266out     OutletVapour as vapour_stream   (Brief="Vapour outlet stream", PosX=0.5079, PosY=0, Symbol="_{outV}");
    267267        InletQ  as power                        (Brief="Heat supplied", PosX=1, PosY=0.6123, Symbol="_{in}");
    268268
    … …  
    282282
    283283        Level                                   = Initial_Level;
    284         OutletL.T                               = Initial_Temperature;
    285         OutletL.z(1:NComp-1)    = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
     284        OutletLiquid.T                          = Initial_Temperature;
     285        OutletLiquid.z(1:NComp-1)       = Initial_Composition(1:NComp-1)/sum(Initial_Composition);
    286286
    287287EQUATIONS
    288288"Molar Concentration"
    289         OutletL.z = vL * C;
     289        OutletLiquid.z = vL * C;
    290290       
    291291"Reaction"
    292         r3 = exp(-7150*'K'/OutletL.T)*(4.85e4*C(1)*C(2) - 1.23e4*C(3)*C(4)) * 'l/mol/s';
     292        r3 = exp(-7150*'K'/OutletLiquid.T)*(4.85e4*C(1)*C(2) - 1.23e4*C(3)*C(4)) * 'l/mol/s';
    293293
    294294"Component Molar Balance"
    295         diff(M)= InletL.F*InletL.z- OutletL.F*OutletL.z - OutletV.F*OutletV.z + stoic*r3*ML*vL;
     295        diff(M)= InletLiquid.F*InletLiquid.z- OutletLiquid.F*OutletLiquid.z - OutletVapour.F*OutletVapour.z + stoic*r3*ML*vL;
    296296       
    297297"Energy Balance"
    298         diff(E) = InletL.F*InletL.h- OutletL.F*OutletL.h - OutletV.F*OutletV.h + InletQ + Hr * r3 * vL*ML;
     298        diff(E) = InletLiquid.F*InletLiquid.h- OutletLiquid.F*OutletLiquid.h - OutletVapour.F*OutletVapour.h + InletQ + Hr * r3 * vL*ML;
    299299       
    300300"Molar Holdup"
    301         M = ML*OutletL.z + MV*OutletV.z;
     301        M = ML*OutletLiquid.z + MV*OutletVapour.z;
    302302       
    303303"Energy Holdup"
    304         E = ML*OutletL.h + MV*OutletV.h - OutletL.P*V;
     304        E = ML*OutletLiquid.h + MV*OutletVapour.h - OutletLiquid.P*V;
    305305       
    306306"Mol fraction normalisation"
    307         sum(OutletL.z)=1.0;
     307        sum(OutletLiquid.z)=1.0;
    308308       
    309309"Liquid Volume"
    310         vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
     310        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    311311
    312312"Vapour Volume"
    313         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z); 
     313        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);   
    314314
    315315"Vapour Density"
    316         rhoV = PP.VapourDensity(OutletV.T, OutletV.P, OutletV.z);
     316        rhoV = PP.VapourDensity(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    317317       
    318318"Level of liquid phase"
    … …  
    322322       
    323323"Mechanical Equilibrium"
    324         OutletL.P = OutletV.P;
     324        OutletLiquid.P = OutletVapour.P;
    325325       
    326326"Thermal Equilibrium"
    327         OutletL.T = OutletV.T; 
     327        OutletLiquid.T = OutletVapour.T;       
    328328       
    329329"Geometry Constraint"
    … …  
    331331
    332332"Chemical Equilibrium"
    333         PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z =
    334         PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;
    335 
    336         sum(OutletL.z)=sum(OutletV.z);
     333        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z =
     334        PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z;
     335
     336        sum(OutletLiquid.z)=sum(OutletVapour.z);
    337337       
    338338end
Note: See TracChangeset for help on using the changeset viewer.