Changeset 859


Ignore:
Timestamp:
Oct 9, 2009, 5:15:53 PM (13 years ago)
Author:
gerson bicca
Message:

updated samples

Location:
trunk/sample
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/sample/processes/Sample_Process.mso

    r853 r859  
    126126       
    127127        Reac.str(z).T = Reac.Inlet.T;
    128         Reac.str(z).z(1:NComp-1) = Reac.Inlet.z(1:NComp-1);
     128        Reac.str(z).z(1:NComp) = Reac.Inlet.z(1:NComp);
    129129
    130130end
  • trunk/sample/reactors/sample_pfr.mso

    r634 r859  
    3939PARAMETERS
    4040
    41         PP              as Plugin (Brief="Physical Properties",
     41
     42
     43        PP as Plugin (Brief="Physical Properties",
     44
    4245                Type = "PP",
    43                 Components = ["acetone", "acetic anhydride", "methane" ],
     46
     47                Components = ["acetone", "ketene", "methane" ],
     48
    4449                LiquidModel = "PR",
     50
    4551                VapourModel = "PR"
     52
    4653        );
     54
    4755        NComp   as Integer;
     56
     57SET
     58
     59        NComp = PP.NumberOfComponents;
     60
     61
    4862
    4963DEVICES
    5064
    5165        s1      as source;
    52         Reac as pfr;
     66
     67        Reac    as pfr;
     68
     69
     70
     71CONNECTIONS
     72
     73        s1.Outlet  to  Reac.Inlet;
     74
     75       
    5376
    5477SET
    55         NComp = PP.NumberOfComponents;
    56        
    57         Reac.NDisc      = 10;
    58         Reac.Across     = 0.7 * 'in^2';
    59         Reac.L                  = 2.28 * 'm';
     78
     79        Reac.NDisc      = 15;
     80
     81        Reac.L          = 2.28 * 'm';
     82
     83        # Area was chosen to match the original example data
     84
     85        Reac.Across = 1.27*'m^3'/Reac.L;
     86
     87
     88
     89        # Only one reaction
     90
    6091        Reac.NReac      = 1;
    6192
    62         # Reaction 1: A -> B + C
     93        # Reaction 1: acetone -> ketene + methane
     94
    6395        Reac.stoic(:,1) = [-1, 1, 1];
     96
     97
     98
     99SPECIFY
     100
     101        # Feed specification
     102
     103        s1.Fw = 8000 * 'kg/h';
     104
     105        s1.T = 1035 * 'K';
     106
     107        s1.P = 1.6  * 'atm';
     108
     109        s1.Composition = [0.98, 0.01, 0.01];
     110
     111
     112
     113        # Adiabatic
     114
     115        Reac.q     = 0     * 'J/s';
     116
    64117       
     118
    65119EQUATIONS
    66120
    67 for z in [1:Reac.NDisc]
     121
     122
     123for z in [1:Reac.NDisc]
     124
    68125       
     126
    69127"Reaction Rate = k*C(1)"
     128
    70129        Reac.r(1,z)  = exp(34.34 - (34222 * 'K') / Reac.str(z).T)*'1/s' * Reac.C(1,z);
    71130
     131
     132
    72133"Heat of reaction"
     134
    73135        Reac.Hr(1,z) = -80.77 * 'kJ/mol';
    74136
     137
     138
    75139"Pressure Drop (no pressure drop)"
     140
    76141        Reac.str(z+1).P = Reac.str(z).P;
     142
     143
    77144
    78145end
    79146
    80 SPECIFY
    81         s1.F = 138/1000 * 'kmol/h';
    82         s1.T = 1035 * 'K';
    83         s1.P = 1.6  * 'atm';
    84         s1.Composition = [1, 0, 0];
    85147
    86 "Adiabatic"
    87         Reac.q     = 0     * 'J/s';
    88        
    89 CONNECTIONS
    90148
    91         s1.Outlet  to  Reac.Inlet;
    92        
    93149INITIAL
    94150
    95 for z in [2:Reac.NDisc+1]
    96        
     151
     152
     153for z in [2:Reac.NDisc+1]
     154
    97155        Reac.str(z).T = Reac.Inlet.T;
    98         Reac.str(z).z(1:NComp-1) = Reac.Inlet.z(1:NComp-1);
     156
     157        Reac.str(z).z(1:NComp) = Reac.Inlet.z(1:NComp);
    99158
    100159end
     160
    101161       
     162
    102163OPTIONS
    103         # This model can be solved in both steady or dynamic
     164
    104165        Dynamic         = true;
     166
    105167        TimeStep        = 0.05;
    106         TimeEnd         = 3;
    107         #Dynamic = false;
    108         #DAESolver = "dassl";
     168
     169        TimeEnd         = 10;
     170
     171        GuessFile = "PFR_AceticAnhydride";
     172
     173        Dynamic = false; # Requires a GuessFile from a dynamic simulation
    109174end
Note: See TracChangeset for help on using the changeset viewer.