Changeset 207


Ignore:
Timestamp:
Mar 15, 2007, 2:17:35 AM (15 years ago)
Author:
Argimiro Resende Secchi
Message:

Fix initial conditions of sample_cstr.

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/reactors/cstr.mso

    r77 r207  
    4141
    4242        EQUATIONS
     43        "Molar concentration"
     44        Outlet.z * M = C * Vr;
     45
    4346        "Molar fraction"
    44         Outlet.z * M = C * Vr;
     47        sum(Outlet.z) = 1;
     48
    4549        "Vapourisation Fraction"
    4650        Outlet.v = Inlet.v;
     
    6569        EQUATIONS
    6670        "Component Molar Balance"
    67         diff(Outlet.z*M) = (Inlet.F*Inlet.z - Outlet.F*Outlet.z) + sumt(stoic*r)*Vr;
     71        diff(Outlet.z*M) = (Inlet.F*Inlet.z - Outlet.F*Outlet.z) + sumt(stoic*r)*Vr; #FIXME sumt()
    6872
    6973        "Reactor Energy Balance"
    7074        diff(M*Outlet.h) = Inlet.F*Inlet.h - Outlet.F*Outlet.h +
    71                 sum(Hr*sum(stoic*r))*Vr - q;
     75                sum(Hr*sum(stoic*r))*Vr - q; # FIXME sum()
    7276end
  • trunk/sample/reactors/sample_cstr.mso

    r82 r207  
    3838        SET
    3939        PP.Components = ["1,2-propylene oxide", "methanol",
    40                 "water", "sulfuric acid", "1,2-propylene glycol" ];
     40                "water", "1,2-propylene glycol" ];
    4141        PP.LiquidModel = "PR";
    4242        PP.VapourModel = "PR";
     
    4444        NComp = PP.NumberOfComponents;
    4545        Reac.NReac = 1;
    46         Reac.stoic = [-1, 0, -1, 0, 1];
     46        Reac.stoic = [-1, 0, -1, 1];
    4747       
    4848        EQUATIONS
    49 #       Reac.Hr = -36400 * "Btu/lbmol" - 7 * "Btu/lbmol/degR" * (Reac.Fout.T - 528 * "degR");
    50 #       Reac.r = 16.96e12 / "h" * exp(- 32400 * "Btu/lbmol"/R/Reac.Fout.T)
    51 #               * Reac.C(1);
    52         Reac.Hr = 0 * "Btu/lbmol";
    53         Reac.r  = 0 * "mol/l/s";
     49        Reac.Hr = -36400 * "Btu/lbmol" - 7 * "Btu/lbmol/degR" * (Reac.Outlet.T - 528 * "degR");
     50        Reac.r = 16.96e12 / "h" * exp(- 32400 * "Btu/lbmol"/R/Reac.Outlet.T)
     51                * Reac.C(1);
    5452
    5553        SPECIFY
     
    5755        s1.P = 1     * "atm";
    5856        s1.T = 534.67* "degR";
    59         s1.z = [0.0468, 0.0783, 0.874, 0.0009, 0];
     57        s1.z = [0.0468, 0.0783, 0.874, 0];
    6058        s1.v = 0.0;
    6159
     
    7068       
    7169        INITIAL
    72         Reac.Outlet.T = Reac.Inlet.T;
    73         Reac.Outlet.z = Reac.Inlet.z;
     70        Reac.Outlet.T = 500*"degR";
     71        Reac.M = 100 * "kmol";
     72        Reac.Outlet.z(1:3) = [0.0468, 0.0783, 0.874];
    7473
    7574        OPTIONS
Note: See TracChangeset for help on using the changeset viewer.