Ignore:
Timestamp:
Mar 7, 2007, 10:17:29 AM (17 years ago)
Author:
Rafael de Pelegrini Soares
Message:

Updated distillation column samples to run on the new language

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/newlanguage/sample/stage_separators/sample_column.mso

    r177 r187  
    5555       
    5656        SPECIFY
    57         feed.F = 113.4 * "kmol/h";
    58         feed.T = 291 * "K";
    59         feed.P = 168.3 * "kPa";
     57        feed.F = 113.4 * 'kmol/h';
     58        feed.T = 291 * 'K';
     59        feed.P = 168.3 * 'kPa';
    6060        feed.z = [0.5, 0.5];
    6161
    62         zero.F = 0 * "kmol/h";
    63         zero.T = 300 * "K";
    64         zero.P = 1 * "atm";
     62        zero.F = 0 * 'kmol/h';
     63        zero.T = 300 * 'K';
     64        zero.P = 1 * 'atm';
    6565        zero.z = [0.5, 0.5];
    6666        zero.v = 0;
    67         zero.h = 0 * "J/mol";
    68        
    69         cond.F = 68 * "kmol/h";
    70         cond.P = 150 * "kPa";
    71         cond.T = 281.75 * "K";
     67        zero.h = 0 * 'J/mol';
     68       
     69        cond.F = 68 * 'kmol/h';
     70        cond.P = 150 * 'kPa';
     71        cond.T = 281.75 * 'K';
    7272        cond.z = [0.6664, 0.3336];
    7373
    74         reb.F = 153 * "kmol/h";
    75         reb.P = 185 * "kPa";
    76         reb.T = 328.12 * "K";
     74        reb.F = 153 * 'kmol/h';
     75        reb.P = 185 * 'kPa';
     76        reb.T = 328.12 * 'K';
    7777        reb.z = [0.001848, 0.9982];
    7878       
    … …  
    8282        sec.NTrays = 2;
    8383        #COLUMN
    84         sec.trays.V = 4 * "ft^3";
    85         sec.trays.Ah = 0.394 * "ft^2";
    86         sec.trays.lw = 20.94 * "in";
    87         sec.trays.hw = 0.125 * "ft";
    88         sec.trays.Q = 0 * "kW";
     84        sec.trays.V = 4 * 'ft^3';
     85        sec.trays.Ah = 0.394 * 'ft^2';
     86        sec.trays.lw = 20.94 * 'in';
     87        sec.trays.hw = 0.125 * 'ft';
     88        sec.trays.Q = 0 * 'kW';
    8989        sec.trays.beta = 0.6;
    9090        sec.trays.alfa = 4;
    91         sec.trays.Ap = 3.94 * "ft^2";
     91        sec.trays.Ap = 3.94 * 'ft^2';
    9292
    9393        INITIAL
    94         sec.trays.OutletL.T = 290 *"K";
     94        sec.trays.OutletL.T = 290 *'K';
    9595        sec.trays.Level = 0.9 * sec.trays.hw;
    9696        sec.trays.OutletL.z(1) = 0.5;
    … …  
    139139       
    140140        SPECIFY
    141         feed.F = 113.4 * "kmol/h";
    142         feed.T = 291 * "K";
    143         feed.P = 168.3 * "kPa";
     141        feed.F = 113.4 * 'kmol/h';
     142        feed.T = 291 * 'K';
     143        feed.P = 168.3 * 'kPa';
    144144        feed.z = [0.5, 0.5];
    145145
    146         zero.F = 0 * "kmol/h";
    147         zero.T = 300 * "K";
    148         zero.P = 1 * "atm";
     146        zero.F = 0 * 'kmol/h';
     147        zero.T = 300 * 'K';
     148        zero.P = 1 * 'atm';
    149149        zero.z = [0.5, 0.5];
    150150        zero.v = 0;
    151         zero.h = 0 * "J/mol";
    152        
    153         cond.F = 68 * "kmol/h";
    154         cond.P = 150 * "kPa";
    155         cond.T = 281.75 * "K";
     151        zero.h = 0 * 'J/mol';
     152       
     153        cond.F = 68 * 'kmol/h';
     154        cond.P = 150 * 'kPa';
     155        cond.T = 281.75 * 'K';
    156156        cond.z = [0.6664, 0.3336];
    157157
    158         reb.F = 153 * "kmol/h";
    159         reb.P = 185 * "kPa";
    160         reb.T = 328.12 * "K";
     158        reb.F = 153 * 'kmol/h';
     159        reb.P = 185 * 'kPa';
     160        reb.T = 328.12 * 'K';
    161161        reb.z = [0.001848, 0.9982];
    162162       
    … …  
    166166        sec.NTrays = 8;
    167167        #COLUMN
    168         sec.trays.V = 4 * "ft^3";
    169         sec.trays.Ah = 0.394 * "ft^2";
    170         sec.trays.lw = 20.94 * "in";
    171         sec.trays.hw = 0.125 * "ft";
    172         sec.trays.Q = 0 * "kW";
     168        sec.trays.V = 4 * 'ft^3';
     169        sec.trays.Ah = 0.394 * 'ft^2';
     170        sec.trays.lw = 20.94 * 'in';
     171        sec.trays.hw = 0.125 * 'ft';
     172        sec.trays.Q = 0 * 'kW';
    173173        sec.trays.beta = 0.6;
    174174        sec.trays.alfa = 4;
    175         sec.trays.Ap = 3.94 * "ft^2";
     175        sec.trays.Ap = 3.94 * 'ft^2';
    176176
    177177        INITIAL
    178         sec.trays.OutletL.T = [290:(330-290)/(sec.NTrays-1):330] *"K";
     178        sec.trays.OutletL.T = [290:(330-290)/(sec.NTrays-1):330] *'K';
    179179        sec.trays.Level = 0.3 * sec.trays.hw;
    180180        sec.trays.OutletL.z(1) = 0.5;
    … …  
    212212       
    213213        CONNECTIONS
    214         feed to col.trays(5).Inlet;
     214        feed.Outlet to col.trays(5).Inlet;
    215215        zero to col.reb.Inlet;
    216216        zero to col.trays([1:4]).Inlet;
    … …  
    220220       
    221221        SPECIFY
    222         feed.Outlet.F = 113.4 * "kmol/h";
    223         feed.Outlet.T = 291 * "K";
    224         feed.Outlet.P = 168.3 * "kPa";
     222        feed.Outlet.F = 113.4 * 'kmol/h';
     223        feed.Outlet.T = 291 * 'K';
     224        feed.Outlet.P = 168.3 * 'kPa';
    225225        feed.Outlet.z = 1/NComp;
    226226       
    227         zero.F = 0 * "kmol/h";
    228         zero.T = 300 * "K";
    229         zero.P = 1 * "atm";
     227        zero.F = 0 * 'kmol/h';
     228        zero.T = 300 * 'K';
     229        zero.P = 1 * 'atm';
    230230        zero.z = 1/NComp;
    231231        zero.v = 0;
    232         zero.h = 0 * "J/mol";
    233        
    234         col.sptop.Outlet2.F = 85 * "kmol/h";
    235         col.reb.OutletL.F = 28.4 * "kmol/h";
     232        zero.h = 0 * 'J/mol';
     233       
     234        col.sptop.Outlet2.F = 85 * 'kmol/h';
     235        col.reb.OutletL.F = 28.4 * 'kmol/h';
    236236        col.sptop.frac = 0.444445;
    237         col.cond.OutletV.F = 0 * "kmol/h";
    238         Qr = 3.7743e6 * "kJ/h";
    239         Qc = -3.71e6 * "kJ/h";
    240         col.pump1.dP = 16 * "kPa";
     237        col.cond.OutletV.F = 0 * 'kmol/h';
     238        Qr = 3.7743e6 * 'kJ/h';
     239        Qc = -3.71e6 * 'kJ/h';
     240        col.pump1.dP = 16 * 'kPa';
    241241        col.trays.Emv = 1;
    242242       
    243243        SET
    244244        col.NTrays = 8;
    245         col.cond.V = 2 * "m^3";
    246         col.cond.Across = 1 * "m^2";
    247         col.trays.V = 4 * "ft^3";
    248         col.trays.Ah = 0.394 * "ft^2";
    249         col.trays.lw = 20.94 * "in";
    250         col.trays.hw = 0.125 * "ft";
    251         col.trays.Q = 0 * "kW";
     245        col.cond.V = 2 * 'm^3';
     246        col.cond.Across = 1 * 'm^2';
     247        col.trays.V = 4 * 'ft^3';
     248        col.trays.Ah = 0.394 * 'ft^2';
     249        col.trays.lw = 20.94 * 'in';
     250        col.trays.hw = 0.125 * 'ft';
     251        col.trays.Q = 0 * 'kW';
    252252        col.trays.beta = 0.6;
    253253        col.trays.alfa = 4;
    254         col.trays.Ap = 3.94 * "ft^2";
    255         col.reb.V = 2 * "m^3";
    256         col.reb.Across = 1 * "m^2";
     254        col.trays.Ap = 3.94 * 'ft^2';
     255        col.reb.V = 2 * 'm^3';
     256        col.reb.Across = 1 * 'm^2';
    257257       
    258258        INITIAL
    259259        # condenser
    260         col.cond.OutletL.T = 260 *"K";
    261         col.cond.Level = 1 * "m";
     260        col.cond.OutletL.T = 260 *'K';
     261        col.cond.Level = 1 * 'm';
    262262        col.cond.OutletL.z([1:4]) = [0.65, 0.05, 0.01, 0.01];
    263263
    264264        # reboiler
    265         col.reb.OutletL.T = 330 *"K";
    266         col.reb.Level = 1 * "m";
     265        col.reb.OutletL.T = 330 *'K';
     266        col.reb.Level = 1 * 'm';
    267267        col.reb.OutletL.z([1:4]) = [0.1, 0.7, 0.01, 0.01];
    268268
    269269        # column trays
    270         col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * "K";
     270        col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K';
    271271        col.trays.Level = 1.2 * col.trays.hw;
    272272        col.trays.OutletL.z([1:4]) = [0.5, 0.05, 0.01, 0.01];
    … …  
    274274        OPTIONS
    275275        RelativeAccuracy = 1e-3;
    276         TimeStep = 0.01;
     276        TimeStep = 0.1;
    277277        TimeEnd = 50;
    278278        #time = [0:0.01:1, 2:50];
Note: See TracChangeset for help on using the changeset viewer.