Ignore:
Timestamp:
Mar 20, 2009, 4:59:07 PM (14 years ago)
Author:
gerson bicca
Message:

PFR reactor updates

Location:
branches/gui/sample/reactors
Files:
1 added
1 edited

Legend:

Unmodified
Added
Removed
  • branches/gui/sample/reactors/sample_pfr.mso

    r585 r745  
    4949DEVICES
    5050
    51         s1      as source;
     51        s1      as simple_source;
    5252        Reac as pfr;
    5353
     
    5656       
    5757        s1.ValidPhases = "Vapour-Liquid";
    58         s1.CompositionBasis = "Molar";
     58        #s1.CompositionBasis = "Molar";
    5959       
    6060        Reac.NDisc      = 10;
    61         Reac.Across     = 0.7 * 'in^2';
    62         Reac.L                  = 2.28 * 'm';
     61        Reac.Dpipe      = 3 * 'in';
     62        Reac.L                  = 15 * 'm';
    6363        Reac.NReac      = 1;
     64        Reac.Roughness  = 4.572E-5*'m';
    6465
    6566        # Reaction 1: A -> B + C
     
    7677        Reac.Hr(1,z) = -80.77 * 'kJ/mol';
    7778
    78 "Pressure Drop (no pressure drop)"
    79         Reac.str(z+1).P = Reac.str(z).P;
    80 
    8179end
    8280
     
    8583        s1.T = 1035 * 'K';
    8684        s1.P = 1.6  * 'atm';
    87         s1.Composition = [1, 0, 0];
     85        s1.MolarComposition = [1, 0, 0];
    8886
    8987        Reac.q     = 0     * 'J/s';
     
    9896       
    9997        Reac.str(z).T = Reac.Inlet.T;
     98        Reac.vel(z) = Reac.vel(1);
     99        #Reac.str(z).F = Reac.str(1).F;
    100100        Reac.str(z).z(1:NComp-1) = Reac.Inlet.z(1:NComp-1);
    101101
    102102end
    103        
     103
     104VARIABLES
     105
     106SOMA as Real;
     107
     108EQUATIONS
     109
     110SOMA = sum(Reac.Outlet.z);
     111
    104112OPTIONS
    105113        # This model can be solved in both steady or dynamic
     
    107115        TimeStep        = 0.05;
    108116        TimeEnd         = 3;
     117        DAESolver(File="sundials", RelativeAccuracy=1e-2);
     118        NLASolver(RelativeAccuracy=1e-4);
     119        #InitialFile = "PFR_AceticAnhydride";
    109120        #Dynamic = false;
    110121        #DAESolver = "dassl";
Note: See TracChangeset for help on using the changeset viewer.