Ignore:
Timestamp:
Mar 1, 2007, 3:35:36 PM (16 years ago)
Author:
gerson bicca
Message:

updated pfr reactor model for the new language

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/newlanguage/eml/reactors/pfr.mso

    r77 r170  
    3030
    3131Model pfr
    32         PARAMETERS
    33 ext PP   as CalcObject (Brief = "External Physical Properties");
    34 ext     NComp as Integer(Brief="Number of components");
    35         NReac as Integer(Brief="Number of reactions");
    36         stoic(NComp, NReac) as Real (Brief = "Stoichiometric Matrix");
    37         NDisc as Integer(Brief="Number of points of discretization", Default=10);
    38         Mw(NComp)  as molweight (Brief="Component Mol Weight");
     32
     33ATTRIBUTES
     34        Pallete         = true;
     35        Brief           = "Brief";
     36        Info            =
     37        "write some information";
     38
     39PARAMETERS
     40
     41outer PP                                                        as Plugin               (Brief = "External Physical Properties");
     42outer   NComp                                           as Integer              (Brief="Number of components");
     43                NReac                                           as Integer              (Brief="Number of reactions");
     44                stoic(NComp, NReac)     as Real                         (Brief = "Stoichiometric Matrix");
     45                NDisc                                           as Integer              (Brief="Number of points of discretization", Default=10);
     46                Mw(NComp)                               as molweight    (Brief="Component Mol Weight");
     47                L                                                               as length               (Brief="Reactor Length");
     48                Across                                          as area                 (Brief="Cross section area");
     49
     50SET
    3951       
    40         L as length(Brief="Reactor Length");
    41         Across as area(Brief="Cross section area");
    42 
    43         SET
    4452        Mw   = PP.MolecularWeight();
    4553
    46         VARIABLES
    47 in      Inlet  as stream;
    48 out     Outlet as stream;
    49         str(NDisc+1) as stream_therm;
     54VARIABLES
     55
     56in              Inlet   as stream       (Brief = "Inlet Stream");
     57out     Outlet  as stream       (Brief = "Outlet Stream");
     58
     59        str(NDisc+1) as streamPH;
    5060        vel(NDisc+1) as velocity;
    5161        vol(NDisc+1) as vol_mol;
    5262        rho(NDisc+1) as dens_mass;
    5363       
    54         q(NDisc)    as heat_rate;
    55         M(NComp, NDisc)    as mol (Brief = "Molar holdup");
    56         Mt(NDisc)    as mol (Brief = "Molar holdup");
    57         C(NComp, NDisc)  as conc_mol(Brief="Components concentration", Lower=-1e-6);
    58         E(NDisc) as energy (Brief="Total Energy Holdup on element");
    59        
    60         r(NReac, NDisc) as reaction_mol;
    61         Hr(NReac, NDisc) as heat_reaction;
     64        q(NDisc)                                as heat_rate;
     65        M(NComp, NDisc)         as mol                  (Brief = "Molar holdup");
     66        Mt(NDisc)                       as mol                  (Brief = "Molar holdup");
     67        C(NComp, NDisc)         as conc_mol     (Brief="Components concentration", Lower=-1e-6);
     68        E(NDisc)                                as energy               (Brief="Total Energy Holdup on element");
     69        r(NReac, NDisc)                 as reaction_mol;
     70        Hr(NReac, NDisc)        as heat_reaction;
    6271
    63         EQUATIONS
    64         "Vapourisation Fraction"
     72EQUATIONS
     73
     74"Vapourisation Fraction"
    6575        str.v = Inlet.v;
    6676
    67         "Inlet boundary"
     77"Inlet boundary"
    6878        str(1).F = Inlet.F;
    6979        str(1).T = Inlet.T;
     
    7181        str(1).z = Inlet.z;
    7282
    73         "Outlet boundary"
     83"Outlet boundary"
    7484        Outlet.F = str(NDisc+1).F;
    7585        Outlet.T = str(NDisc+1).T;
Note: See TracChangeset for help on using the changeset viewer.