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

    r709 r735  
    4747        VapourSideStream        as vapour_stream        (Brief="vapour Sidestream", Hidden=true, Symbol="_{outV}");
    4848
    49 in      InletL                  as stream                               (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
    50 in      InletV                  as stream                               (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
    51 out     OutletL         as liquid_stream        (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");
    52 out     OutletV         as vapour_stream        (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");
     49in      InletLiquid                     as stream                               (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
     50in      InletVapour                     as stream                               (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
     51out     OutletLiquid    as liquid_stream        (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");
     52out     OutletVapour    as vapour_stream        (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");
    5353
    5454
     
    6565EQUATIONS
    6666"Component Molar Balance"
    67         diff(M)=Inlet.F*Inlet.z + InletL.F*InletL.z + InletV.F*InletV.z- OutletL.F*OutletL.z - OutletV.F*OutletV.z-
     67        diff(M)=Inlet.F*Inlet.z + InletLiquid.F*InletLiquid.z + InletVapour.F*InletVapour.z- OutletLiquid.F*OutletLiquid.z - OutletVapour.F*OutletVapour.z-
    6868        LiquidSideStream.F*LiquidSideStream.z-VapourSideStream.F*VapourSideStream.z;
    6969
    7070"Molar Holdup"
    71         M = ML*OutletL.z + MV*OutletV.z;
     71        M = ML*OutletLiquid.z + MV*OutletVapour.z;
    7272
    7373"Mol fraction normalisation"
    74         sum(OutletL.z)= 1.0;
    75         sum(OutletL.z)= sum(OutletV.z);
     74        sum(OutletLiquid.z)= 1.0;
     75        sum(OutletLiquid.z)= sum(OutletVapour.z);
    7676
    7777"Liquid Volume"
    78         vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
     78        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    7979
    8080"Vapour Volume"
    81         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
     81        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    8282
    8383"Chemical Equilibrium"
    84         PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z = PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, yideal)*yideal;
     84        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z = PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, yideal)*yideal;
    8585
    8686"Thermal Equilibrium"
    87         OutletV.T = OutletL.T;
     87        OutletVapour.T = OutletLiquid.T;
    8888
    8989"Mechanical Equilibrium"
    90         OutletV.P = OutletL.P;
     90        OutletVapour.P = OutletLiquid.P;
    9191       
    9292"Thermal Equilibrium Vapour Side Stream"
    93         OutletV.T = VapourSideStream.T;
     93        OutletVapour.T = VapourSideStream.T;
    9494
    9595"Thermal Equilibrium Liquid Side Stream"
    96         OutletL.T = LiquidSideStream.T;
     96        OutletLiquid.T = LiquidSideStream.T;
    9797
    9898"Mechanical Equilibrium Vapour Side Stream"
    99         OutletV.P= VapourSideStream.P;
     99        OutletVapour.P= VapourSideStream.P;
    100100
    101101"Mechanical Equilibrium Liquid Side Stream"
    102         OutletL.P = LiquidSideStream.P;
     102        OutletLiquid.P = LiquidSideStream.P;
    103103
    104104"Composition Liquid Side Stream"
    105         OutletL.z= LiquidSideStream.z;
     105        OutletLiquid.z= LiquidSideStream.z;
    106106       
    107107"Composition Vapour Side Stream"
    108         OutletV.z= VapourSideStream.z;
     108        OutletVapour.z= VapourSideStream.z;
    109109
    110110end
     
    120120* the Liquid inlet stream
    121121* the Vapour inlet stream
    122 * the Vapour outlet flow (OutletV.F)
     122* the Vapour outlet flow (OutletVapour.F)
    123123       
    124124== Initial ==
    125 * the plate temperature (OutletL.T)
    126 * the liquid height (Level) OR the liquid flow OutletL.F
    127 * (NoComps - 1) OutletL compositions
     125* the plate temperature (OutletLiquid.T)
     126* the liquid height (Level) OR the liquid flow OutletLiquid.F
     127* (NoComps - 1) OutletLiquid compositions
    128128
    129129== Options ==
     
    149149
    150150"Liquid Density"
    151         rhoL = PP.LiquidDensity(OutletL.T, OutletL.P, OutletL.z);
     151        rhoL = PP.LiquidDensity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    152152       
    153153"Vapour Density"
    154         rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z);
     154        rhoV = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z);
    155155
    156156end
     
    176176* the Liquid inlet stream;
    177177* the Vapour inlet stream;
    178 * the Vapour outlet flow (OutletV.F);
     178* the Vapour outlet flow (OutletVapour.F);
    179179* the reaction related variables.
    180180       
    181181== Initial ==
    182 * the plate temperature (OutletL.T)
    183 * the liquid height (Level) OR the liquid flow OutletL.F
    184 * (NoComps - 1) OutletL compositions
     182* the plate temperature (OutletLiquid.T)
     183* the liquid height (Level) OR the liquid flow OutletLiquid.F
     184* (NoComps - 1) OutletLiquid compositions
    185185";
    186186
     
    197197       
    198198
    199 in              InletL          as      stream                          (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
    200 in              InletV          as      stream                          (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
    201 out     OutletL         as      liquid_stream           (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");
    202 out     OutletV         as      vapour_stream   (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");
     199in              InletLiquid             as      stream                          (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
     200in              InletVapour             as      stream                          (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
     201out     OutletLiquid    as      liquid_stream           (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");
     202out     OutletVapour    as      vapour_stream   (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");
    203203
    204204        yideal(NComp)   as fraction;
     
    221221
    222222"Molar Concentration"
    223         OutletL.z = vL * C;
     223        OutletLiquid.z = vL * C;
    224224       
    225225"Reaction"
    226         r3 = exp(-7150*'K'/OutletL.T)*(4.85e4*C(1)*C(2) - 1.23e4*C(3)*C(4))*'l/mol/s';
     226        r3 = exp(-7150*'K'/OutletLiquid.T)*(4.85e4*C(1)*C(2) - 1.23e4*C(3)*C(4))*'l/mol/s';
    227227       
    228228"Molar Holdup"
    229         M = ML*OutletL.z + MV*OutletV.z;
     229        M = ML*OutletLiquid.z + MV*OutletVapour.z;
    230230
    231231"Thermal Equilibrium Vapour Side Stream"
    232         OutletV.T = VapourSideStream.T;
     232        OutletVapour.T = VapourSideStream.T;
    233233
    234234"Thermal Equilibrium Liquid Side Stream"
    235         OutletL.T = LiquidSideStream.T;
     235        OutletLiquid.T = LiquidSideStream.T;
    236236
    237237"Mechanical Equilibrium Vapour Side Stream"
    238         OutletV.P= VapourSideStream.P;
     238        OutletVapour.P= VapourSideStream.P;
    239239
    240240"Mechanical Equilibrium Liquid Side Stream"
    241         OutletL.P = LiquidSideStream.P;
     241        OutletLiquid.P = LiquidSideStream.P;
    242242
    243243"Composition Liquid Side Stream"
    244         OutletL.z= LiquidSideStream.z;
     244        OutletLiquid.z= LiquidSideStream.z;
    245245       
    246246"Composition Vapour Side Stream"
    247         OutletV.z= VapourSideStream.z;
     247        OutletVapour.z= VapourSideStream.z;
    248248
    249249"Mol fraction normalisation"
    250         sum(OutletL.z)= 1.0;
     250        sum(OutletLiquid.z)= 1.0;
    251251
    252252"Liquid Volume"
    253         vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
     253        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    254254
    255255"Vapour Volume"
    256         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
     256        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    257257
    258258"Thermal Equilibrium"
    259         OutletV.T = OutletL.T;
     259        OutletVapour.T = OutletLiquid.T;
    260260
    261261"Mechanical Equilibrium"
    262         OutletV.P = OutletL.P;
     262        OutletVapour.P = OutletLiquid.P;
    263263
    264264        Vol = ML*vL;
    265265       
    266266"Liquid Density"
    267         rhoL = PP.LiquidDensity(OutletL.T, OutletL.P, OutletL.z);
     267        rhoL = PP.LiquidDensity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    268268
    269269"Vapour Density"
    270         rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z);
     270        rhoV = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z);
    271271
    272272"Chemical Equilibrium"
    273         PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z =       PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, yideal)*yideal;
    274 
    275         sum(OutletL.z)= sum(OutletV.z);
     273        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z =   PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, yideal)*yideal;
     274
     275        sum(OutletLiquid.z)= sum(OutletVapour.z);
    276276
    277277end
     
    293293       
    294294== Initial ==
    295 * the plate temperature (OutletL.T)
     295* the plate temperature (OutletLiquid.T)
    296296* the liquid molar holdup ML
    297 * (NoComps - 1) OutletL compositions
     297* (NoComps - 1) OutletLiquid compositions
    298298";     
    299299       
     
    313313
    314314        Inlet                   as stream                               (Brief="Feed stream", PosX=0, PosY=0.4932, Symbol="_{in}");
    315 in      InletL                  as stream                               (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
    316 in      InletV                  as stream                               (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
    317 out     OutletL         as liquid_stream        (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");
    318 out     OutletV         as vapour_stream        (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");
     315in      InletLiquid                     as stream                               (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
     316in      InletVapour                     as stream                               (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
     317out     OutletLiquid    as liquid_stream        (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");
     318out     OutletVapour    as vapour_stream        (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");
    319319
    320320        M(NComp) as mol (Brief="Molar Holdup in the tray", Default=0.01, Lower=0, Upper=100);
     
    346346
    347347"Component Molar Balance"
    348         diff(M)=Inlet.F*Inlet.z + InletL.F*InletL.z + InletV.F*InletV.z- OutletL.F*OutletL.z - OutletV.F*OutletV.z;
     348        diff(M)=Inlet.F*Inlet.z + InletLiquid.F*InletLiquid.z + InletVapour.F*InletVapour.z- OutletLiquid.F*OutletLiquid.z - OutletVapour.F*OutletVapour.z;
    349349
    350350"Molar Holdup"
    351         M = ML*OutletL.z + MV*OutletV.z;
     351        M = ML*OutletLiquid.z + MV*OutletVapour.z;
    352352       
    353353"Mol fraction normalisation"
    354         sum(OutletL.z)= 1.0;
     354        sum(OutletLiquid.z)= 1.0;
    355355       
    356356"Liquid Volume"
    357         vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
     357        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    358358
    359359"Vapour Volume"
    360         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
     360        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    361361       
    362362"Chemical Equilibrium"
    363         PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z =       PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z;
     363        PP.LiquidFugacityCoefficient(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z)*OutletLiquid.z =   PP.VapourFugacityCoefficient(OutletVapour.T, OutletVapour.P, OutletVapour.z)*OutletVapour.z;
    364364       
    365365"Thermal Equilibrium"
    366         OutletV.T = OutletL.T;
     366        OutletVapour.T = OutletLiquid.T;
    367367       
    368368"Mechanical Equilibrium"
    369         OutletL.P = OutletV.P;
     369        OutletLiquid.P = OutletVapour.P;
    370370       
    371371"Liquid Density"
    372         rhoL = PP.LiquidDensity(OutletL.T, OutletL.P, OutletL.z);
     372        rhoL = PP.LiquidDensity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    373373
    374374"Vapour Density"
    375         rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z);
     375        rhoV = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z);
    376376
    377377"Liquid viscosity"
    378         miL = PP.LiquidViscosity(OutletL.T, OutletL.P, OutletL.z);
     378        miL = PP.LiquidViscosity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    379379
    380380"Vapour viscosity"
    381         miV = PP.VapourViscosity(InletV.T, InletV.P, InletV.z);
     381        miV = PP.VapourViscosity(InletVapour.T, InletVapour.P, InletVapour.z);
    382382
    383383"Volume flow rate of liquid, m^3/m^2/s"
    384         uL * Al = OutletL.F * vL;
    385        
    386         deltaP = InletV.P - OutletV.P;
     384        uL * Al = OutletLiquid.F * vL;
     385       
     386        deltaP = InletVapour.P - OutletVapour.P;
    387387       
    388388end
     
    458458in      Inlet as stream (Brief="Feed stream", PosX=0, PosY=0.4932, Symbol="_{in}");
    459459in      InletFV as stream (Brief="Feed stream", PosX=0, PosY=0.4932, Symbol="_{in}");
    460 in      InletL as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
    461 in      InletV as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
    462 out     OutletL as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");
    463 out     OutletV as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");
     460in      InletLiquid as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0, Symbol="_{inL}");
     461in      InletVapour as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1, Symbol="_{inV}");
     462out     OutletLiquid as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1, Symbol="_{outL}");
     463out     OutletVapour as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0, Symbol="_{outV}");
    464464
    465465        M_liq(NComp) as mol (Brief="Liquid Molar Holdup in the tray");
     
    479479        EQUATIONS
    480480        "Component Molar Balance"
    481         diff(M_liq)=Inlet.F*Inlet.z + InletL.F*InletL.z
    482         - OutletL.F*OutletL.z + interf.NL;
    483        
    484         diff(M_vap)=InletFV.F*InletFV.z + InletV.F*InletV.z
    485         - OutletV.F*OutletV.z - interf.NV;
     481        diff(M_liq)=Inlet.F*Inlet.z + InletLiquid.F*InletLiquid.z
     482        - OutletLiquid.F*OutletLiquid.z + interf.NL;
     483       
     484        diff(M_vap)=InletFV.F*InletFV.z + InletVapour.F*InletVapour.z
     485        - OutletVapour.F*OutletVapour.z - interf.NV;
    486486       
    487487        "Energy Balance"
    488         diff(E_liq) = Inlet.F*Inlet.h + InletL.F*InletL.h
    489                 - OutletL.F*OutletL.h  + Q + interf.E_liq;
    490        
    491         diff(E_vap) = InletFV.F*InletFV.h + InletV.F*InletV.h
    492                 - OutletV.F*OutletV.h  - interf.E_vap;
     488        diff(E_liq) = Inlet.F*Inlet.h + InletLiquid.F*InletLiquid.h
     489                - OutletLiquid.F*OutletLiquid.h  + Q + interf.E_liq;
     490       
     491        diff(E_vap) = InletFV.F*InletFV.h + InletVapour.F*InletVapour.h
     492                - OutletVapour.F*OutletVapour.h  - interf.E_vap;
    493493       
    494494        "Molar Holdup"
    495         M_liq = ML*OutletL.z;
    496        
    497         M_vap = MV*OutletV.z;
     495        M_liq = ML*OutletLiquid.z;
     496       
     497        M_vap = MV*OutletVapour.z;
    498498       
    499499        "Energy Holdup"
    500         E_liq = ML*(OutletL.h - OutletL.P*vL);
    501        
    502         E_vap = MV*(OutletV.h - OutletV.P*vV);
     500        E_liq = ML*(OutletLiquid.h - OutletLiquid.P*vL);
     501       
     502        E_vap = MV*(OutletVapour.h - OutletVapour.P*vV);
    503503       
    504504        "Energy Rate through the interface"
    505         interf.E_liq = interf.htL*interf.a*(interf.T-OutletL.T)+sum(interf.NL)*interf.hL;       
    506        
    507         interf.E_vap = interf.htV*interf.a*(OutletV.T-interf.T)+sum(interf.NV)*interf.hV;
     505        interf.E_liq = interf.htL*interf.a*(interf.T-OutletLiquid.T)+sum(interf.NL)*interf.hL; 
     506       
     507        interf.E_vap = interf.htV*interf.a*(OutletVapour.T-interf.T)+sum(interf.NV)*interf.hV;
    508508       
    509509        "Mass Conservation"
     
    514514       
    515515        "Mol fraction normalisation"
    516         sum(OutletL.z)= 1.0;
    517         sum(OutletL.z)= sum(OutletV.z);
     516        sum(OutletLiquid.z)= 1.0;
     517        sum(OutletLiquid.z)= sum(OutletVapour.z);
    518518        sum(interf.x)=1.0;
    519519        sum(interf.x)=sum(interf.y);
    520520       
    521521        "Liquid Volume"
    522         vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z);
     522        vL = PP.LiquidVolume(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    523523        "Vapour Volume"
    524         vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z);
     524        vV = PP.VapourVolume(OutletVapour.T, OutletVapour.P, OutletVapour.z);
    525525       
    526526        "Chemical Equilibrium"
     
    535535
    536536        "Total Mass Transfer Rates"
    537         interf.NL(1:NC1)=interf.a*sumt(interf.kL*(interf.x(1:NC1)-OutletL.z(1:NC1)))/vL+
    538                 OutletL.z(1:NC1)*sum(interf.NL);
     537        interf.NL(1:NC1)=interf.a*sumt(interf.kL*(interf.x(1:NC1)-OutletLiquid.z(1:NC1)))/vL+
     538                OutletLiquid.z(1:NC1)*sum(interf.NL);
    539539
    540540#       interf.NL(1:NC1)=0.01*'kmol/s';
    541541       
    542         interf.NV(1:NC1)=interf.a*sumt(interf.kV*(OutletV.z(1:NC1)-interf.y(1:NC1)))/vV+
    543                 OutletV.z(1:NC1)*sum(interf.NV);
     542        interf.NV(1:NC1)=interf.a*sumt(interf.kV*(OutletVapour.z(1:NC1)-interf.y(1:NC1)))/vV+
     543                OutletVapour.z(1:NC1)*sum(interf.NV);
    544544
    545545        "Mechanical Equilibrium"
    546         OutletV.P = OutletL.P;
    547         interf.P=OutletL.P;
     546        OutletVapour.P = OutletLiquid.P;
     547        interf.P=OutletLiquid.P;
    548548end
    549549
     
    558558* the Liquid inlet stream
    559559* the Vapour inlet stream
    560 * the Vapour outlet flow (OutletV.F)
     560* the Vapour outlet flow (OutletVapour.F)
    561561       
    562562== Initial ==
    563 * the plate temperature of both phases (OutletL.T and OutletV.T)
     563* the plate temperature of both phases (OutletLiquid.T and OutletVapour.T)
    564564* the liquid height (Level) OR the liquid flow holdup (ML)
    565565* the vapor holdup (MV)
    566 * (NoComps - 1) OutletL compositions
     566* (NoComps - 1) OutletLiquid compositions
    567567";
    568568
     
    584584        EQUATIONS
    585585        "Liquid Density"
    586         rhoL = PP.LiquidDensity(OutletL.T, OutletL.P, OutletL.z);
     586        rhoL = PP.LiquidDensity(OutletLiquid.T, OutletLiquid.P, OutletLiquid.z);
    587587        "Vapour Density"
    588         rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z);
     588        rhoV = PP.VapourDensity(InletVapour.T, InletVapour.P, InletVapour.z);
    589589
    590590        switch LiquidFlow
    591591                case "on":
    592592                "Francis Equation"
    593 #               OutletL.F*vL = 1.84*'m^0.5/s'*lw*((Level-(beta*hw))/(beta))^1.5;
    594                 OutletL.F*vL = 1.84*'1/s'*lw*((Level-(beta*hw))/(beta))^2;
     593#               OutletLiquid.F*vL = 1.84*'m^0.5/s'*lw*((Level-(beta*hw))/(beta))^1.5;
     594                OutletLiquid.F*vL = 1.84*'1/s'*lw*((Level-(beta*hw))/(beta))^2;
    595595                when Level < (beta * hw) switchto "off";
    596596               
    597597                case "off":
    598598                "Low level"
    599                 OutletL.F = 0 * 'mol/h';
     599                OutletLiquid.F = 0 * 'mol/h';
    600600                when Level > (beta * hw) + 1e-6*'m' switchto "on";
    601601        end
     
    603603        switch VapourFlow
    604604                case "on":
    605                 InletV.F*vV = sqrt((InletV.P - OutletV.P)/(rhoV*alfa))*Ah;
    606                 when InletV.F < 1e-6 * 'kmol/h' switchto "off";
     605                InletVapour.F*vV = sqrt((InletVapour.P - OutletVapour.P)/(rhoV*alfa))*Ah;
     606                when InletVapour.F < 1e-6 * 'kmol/h' switchto "off";
    607607               
    608608                case "off":
    609                 InletV.F = 0 * 'mol/s';
    610                 when InletV.P > OutletV.P + Level*g*rhoL + 1e-1 * 'atm' switchto "on";
     609                InletVapour.F = 0 * 'mol/s';
     610                when InletVapour.P > OutletVapour.P + Level*g*rhoL + 1e-1 * 'atm' switchto "on";
    611611        end     
    612612end
Note: See TracChangeset for help on using the changeset viewer.