Ignore:
Timestamp:
Mar 1, 2007, 3:09:21 PM (15 years ago)
Author:
gerson bicca
Message:

updated cstr reactor model for the new language

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/newlanguage/sample/reactors/sample_cstr.mso

    r82 r169  
    2727
    2828FlowSheet cstr_PropyleneGlycol
    29         PARAMETERS
    30         PP   as CalcObject(Brief="Physical Properties",File="vrpp");
    31         NComp as Integer;
    32         R as Real(Default=8.3145, Unit="J/mol/K");
    3329
    34         DEVICES
    35         s1   as stream_therm;
    36         Reac as cstr;
     30PARAMETERS
     31        PP              as Plugin       (Brief="Physical Properties",File="vrpp");
     32        NComp   as Integer;
     33        R                       as Real         (Default=8.3145, Unit='J/mol/K');
    3734
    38         SET
    39         PP.Components = ["1,2-propylene oxide", "methanol",
    40                 "water", "sulfuric acid", "1,2-propylene glycol" ];
    41         PP.LiquidModel = "PR";
    42         PP.VapourModel = "PR";
     35DEVICES
     36        s1      as source;
     37        Reac    as cstr;
    4338
    44         NComp = PP.NumberOfComponents;
    45         Reac.NReac = 1;
    46         Reac.stoic = [-1, 0, -1, 0, 1];
     39SET
     40
     41        PP.Components   = ["1,2-propylene oxide", "methanol","water","sulfuric acid", "1,2-propylene glycol" ];
     42        PP.LiquidModel  = "PR";
     43        PP.VapourModel  = "PR";
     44
     45        NComp                   = PP.NumberOfComponents;
     46        Reac.NReac      = 1;
     47        Reac.stoic              = [-1, 0, -1, 0, 1];
    4748       
    48         EQUATIONS
     49EQUATIONS
     50
    4951#       Reac.Hr = -36400 * "Btu/lbmol" - 7 * "Btu/lbmol/degR" * (Reac.Fout.T - 528 * "degR");
    5052#       Reac.r = 16.96e12 / "h" * exp(- 32400 * "Btu/lbmol"/R/Reac.Fout.T)
    5153#               * Reac.C(1);
    52         Reac.Hr = 0 * "Btu/lbmol";
    53         Reac.r  = 0 * "mol/l/s";
     54        Reac.Hr = 0 * 'Btu/lbmol';
     55        Reac.r  = 0 * 'mol/l/s';
    5456
    5557        SPECIFY
    56         s1.F = 52126 * "mol/h";
    57         s1.P = 1     * "atm";
    58         s1.T = 534.67* "degR";
    59         s1.z = [0.0468, 0.0783, 0.874, 0.0009, 0];
    60         s1.v = 0.0;
     58        s1.Outlet.F = 52126 * 'mol/h';
     59        s1.Outlet.P = 1     * 'atm';
     60        s1.Outlet.T = 534.67* 'degR';
     61        s1.Outlet.z = [0.0468, 0.0783, 0.874, 0.0009, 0];
     62        s1.Outlet.v = 0.0;
    6163
    62         Reac.Vr    = 300   * "gal";
    63         Reac.Outlet.F= Reac.Inlet.F;
    64         "Adiabatic"
    65         Reac.q     = 0     * "J/s";
    66         Reac.Outlet.P= 1     * "atm";
     64        Reac.Vr                 = 300   * 'gal';
     65        Reac.Outlet.F   = Reac.Inlet.F;
     66
     67"Adiabatic"
     68        Reac.q                  = 0     * 'J/s';
     69        Reac.Outlet.P   = 1     * 'atm';
    6770       
    68         CONNECTIONS
    69         s1  to  Reac.Inlet;
     71CONNECTIONS
     72
     73        s1.Outlet  to  Reac.Inlet;
    7074       
    71         INITIAL
     75INITIAL
     76
    7277        Reac.Outlet.T = Reac.Inlet.T;
    7378        Reac.Outlet.z = Reac.Inlet.z;
    7479
    75         OPTIONS
    76         time = [0:0.1:5]*"h";
     80OPTIONS
     81        TimeStep        = 0.1;
     82        TimeEnd         = 5;
     83
    7784end
Note: See TracChangeset for help on using the changeset viewer.