Ignore:
Timestamp:
Mar 2, 2007, 10:06:53 AM (16 years ago)
Author:
gerson bicca
Message:

some modifications in the models for the new language (reactors model)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/newlanguage/sample/reactors/fogler/chap4/molarflow_pfr.mso

    r82 r171  
    4848        PARAMETERS
    4949        NComp   as Integer;
    50         ka              as Real                 (Brief="Specific rate of reaction", Unit="1/min");
     50        ka              as Real                 (Brief="Specific rate of reaction", Unit='1/min');
    5151        Kc              as conc_mol     (Brief="Equilibrium constant");
    5252       
    5353        VARIABLES
    54     F(NComp)as flow_mol         (Brief="Molar flow", Unit="mol/min");
    55         C(NComp)as conc_mol     (Brief="Molar concentration", Unit="mol/l");
    56         r(NComp)as reaction_mol (Brief="Reaction rate", Unit="mol/min/l");
    57         Fo(NComp) as flow_mol   (Brief="Input molar flow of A", Unit="mol/min");
    58     Ft          as flow_mol     (Brief="Total molar flow", Unit="mol/min");
    59         Cto     as conc_mol     (Brief="Initial concentration", Unit="mol/l");
    60         V               as volume               (Brief="Reactor volume", Unit="l");
     54    F(NComp)as flow_mol         (Brief="Molar flow", Unit='mol/min');
     55        C(NComp)as conc_mol     (Brief="Molar concentration", Unit='mol/l');
     56        r(NComp)as reaction_mol (Brief="Reaction rate", Unit='mol/min/l');
     57        Fo(NComp) as flow_mol   (Brief="Input molar flow of A", Unit='mol/min');
     58    Ft          as flow_mol     (Brief="Total molar flow", Unit='mol/min');
     59        Cto     as conc_mol     (Brief="Initial concentration", Unit='mol/l');
     60        V               as volume               (Brief="Reactor volume", Unit='l');
    6161       
    6262        EQUATIONS
    6363        "Change time in V"
    64         V = time*"l/s";
     64        V = time*'l/s';
    6565
    6666        "Molar balance"
    67         diff(F(1))= r(1)*"l/s";
     67        diff(F(1))= r(1)*'l/s';
    6868       
    6969        "Reaction rate of A"
     
    8484        SET
    8585        NComp = 2; # components A and B
    86         ka = 2.7*"1/min";
    87         Kc = 1.2*"mol/l";
     86        ka = 2.7*'1/min';
     87        Kc = 1.2*'mol/l';
    8888       
    8989    SPECIFY
    9090        "Input total molar concentration"
    91         Cto = 0.1*"mol/l";
     91        Cto = 0.1*'mol/l';
    9292       
    9393        "Input molar flow"
    94         Fo = [10, 0]*"mol/min";
     94        Fo = [10, 0]*'mol/min';
    9595       
    9696        INITIAL
    9797        "Molar flow of A"
    98         F(1) = 10*"mol/min";
     98        F(1) = 10*'mol/min';
    9999       
    100100        OPTIONS
    101         time=[0:1:100];
     101        TimeStep=1;
     102        TimeEnd=100;
    102103        DAESolver = "dassl";
    103104end
Note: See TracChangeset for help on using the changeset viewer.