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/eml/reactors/cstr.mso

    r77 r169  
    1616* Model of cstr reactor
    1717*--------------------------------------------------------------------
    18 *
    19 *
    20 *
    21 *----------------------------------------------------------------------
    2218* Author: Paula B. Staudt
    2319* $Id$
     
    2723
    2824Model cstr_basic
    29         PARAMETERS
    30 ext     NComp as Integer;
    31         NReac as Integer;
    32         stoic(NComp, NReac) as Real (Brief = "Stoichiometric Matrix");
    3325
    34         VARIABLES
    35 in      Inlet  as stream;
    36 out     Outlet as stream_therm;
    37         q    as heat_rate;
    38         Vr   as volume (Brief="Reacting Volume");
    39         M    as mol (Brief = "Molar total amount");
    40         C(NComp)  as conc_mol(Brief="Components concentration");
     26ATTRIBUTES
     27        Pallete         = false;
     28        Brief           = "Brief";
     29        Info            =
     30        "write some information";
     31       
     32PARAMETERS
    4133
    42         EQUATIONS
    43         "Molar fraction"
     34outer   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
     38VARIABLES
     39
     40in              Inlet   as stream               (Brief="Inlet Stream");
     41out     Outlet  as streamPH     (Brief="Outlet Stream");
     42
     43        q                       as heat_rate    (Brief="Heat");
     44        Vr                      as volume               (Brief="Reacting Volume");
     45        M                       as mol                  (Brief = "Molar total amount");
     46        C(NComp)        as conc_mol     (Brief="Components concentration");
     47
     48EQUATIONS
     49
     50"Molar fraction"
    4451        Outlet.z * M = C * Vr;
    45         "Vapourisation Fraction"
     52       
     53"Vapourisation Fraction"
    4654        Outlet.v = Inlet.v;
     55
    4756end
    4857
    4958Model cstr_yield as cstr_basic
    50         PARAMETERS
    51         yield(NReac) as fraction;
     59
     60ATTRIBUTES
     61        Pallete         = true;
     62        Brief           = "Brief";
     63        Info            =
     64        "write some information";
    5265       
    53         EQUATIONS
     66PARAMETERS
     67       
     68        yield(NReac) as fraction (Brief="Yield");
     69       
     70EQUATIONS
     71
    5472        Outlet.z*Outlet.F = Inlet.z*Inlet.F * sumt(stoic*(1-yield));
    55         "Reactor Energy Balance"
     73
     74"Reactor Energy Balance"
    5675        diff(M*Outlet.h) = Inlet.F*Inlet.h - Outlet.F*Outlet.h - q;
     76
    5777end
    5878
    5979
    6080Model cstr as cstr_basic
    61         VARIABLES
    62         r(NReac) as reaction_mol;
    63         Hr(NReac) as heat_reaction;
    6481       
    65         EQUATIONS
    66         "Component Molar Balance"
     82ATTRIBUTES
     83        Pallete         = true;
     84        Brief           = "Brief";
     85        Info            =
     86        "write some information";
     87       
     88VARIABLES
     89       
     90        r(NReac)        as reaction_mol         (Brief="Molar Reaction Rate");
     91        Hr(NReac)       as heat_reaction        (Brief="Heat Reaction");
     92       
     93EQUATIONS
     94
     95"Component Molar Balance"
    6796        diff(Outlet.z*M) = (Inlet.F*Inlet.z - Outlet.F*Outlet.z) + sumt(stoic*r)*Vr;
    6897
    69         "Reactor Energy Balance"
    70         diff(M*Outlet.h) = Inlet.F*Inlet.h - Outlet.F*Outlet.h +
    71                 sum(Hr*sum(stoic*r))*Vr - q;
     98"Reactor Energy Balance"
     99        diff(M*Outlet.h) = Inlet.F*Inlet.h - Outlet.F*Outlet.h +sum(Hr*sum(stoic*r))*Vr - q;
     100
    72101end
Note: See TracChangeset for help on using the changeset viewer.