Changeset 402


Ignore:
Timestamp:
Nov 15, 2007, 11:36:42 PM (15 years ago)
Author:
Rafael de Pelegrini Soares
Message:

More samples on gibbs reactor simple

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/sample/miscellaneous/sample_gibbs_reactor_simple.mso

    r401 r402  
    4747        R  as Real(Brief="Universal gas constant", Unit='J/mol/K', Default=8.314);
    4848        T0 as temperature(Default = 298.15);
    49         P0 as pressure(Default=1);
     49        P0 as Real(Default=1, Unit='bar');
    5050        G0(NComp) as energy_mol         (Brief="Gibbs energy in standard state");
    5151        H0(NComp) as energy_mol         (Brief="Enthalpy in standard state");
     
    5959        P as pressure;
    6060        ni(NComp) as positive(Brief="Initial number of mols", Unit='mol');
    61         n(NComp)  as positive(Brief="Number of mols at equilibrium", Unit='mol', Upper=2);
    62         advance   as Real(Unit='mol', Lower=-2, Upper=2);
     61        n(NComp)  as positive(Brief="Number of mols at equilibrium", Unit='mol', Upper=10);
     62        advance   as Real(Unit='mol', Lower=-10, Upper=10);
    6363
    6464        K  as Real(Brief="Reaction equilibrium constant at T");
     
    113113        Dynamic = false;
    114114end
     115
     116
     117# Ethanol production by ethylene hydratation at 250 degC and 35 bar
     118FlowSheet gibbs_reactor_simple_sample2 as gibbs_reactor_simple
     119        PARAMETERS
     120        PP as Plugin(Brief="Physical Properties",
     121                Type="PP",
     122                Components = ["ethylene", "water", "ethanol"],
     123                LiquidModel = "IdealLiquid",
     124                VapourModel = "PR"
     125        );
     126        NComp as Integer;
     127
     128        SET
     129        NComp = PP.NumberOfComponents;
     130        stoic = [-1, -1, 1];
     131
     132        SPECIFY
     133        T = (250 + 273.15) * 'K';
     134        P = 35 * 'atm';
     135        ni = [1, 5, 0] * 'mol';
     136       
     137        # Expected results:
     138        # advance = 0.21;
     139       
     140        OPTIONS
     141        Dynamic = false;
     142end
     143
     144# Water-gas-shift T = 1100 K, P = 1 bar
     145FlowSheet gibbs_reactor_simple_sample3 as gibbs_reactor_simple
     146        PARAMETERS
     147        PP as Plugin(Brief="Physical Properties",
     148                Type="PP",
     149                Components = ["carbon monoxide", "water", "carbon dioxide", "hydrogen"],
     150                LiquidModel = "IdealLiquid",
     151                VapourModel = "PR"
     152        );
     153        NComp as Integer;
     154
     155        SET
     156        NComp = PP.NumberOfComponents;
     157        stoic = [-1, -1, 1, 1];
     158
     159        SPECIFY
     160        T = 1100 * 'K';
     161        P = 1 * 'bar';
     162        ni = [1, 1, 0, 0] * 'mol';
     163       
     164        # Expected results:
     165        # advance = 0.41;
     166       
     167        OPTIONS
     168        Dynamic = false;
     169end
     170
     171
     172# Water-gas-shift T = 1100 K, P = 10 bar
     173FlowSheet gibbs_reactor_simple_sample4 as gibbs_reactor_simple
     174        PARAMETERS
     175        PP as Plugin(Brief="Physical Properties",
     176                Type="PP",
     177                Components = ["carbon monoxide", "water", "carbon dioxide", "hydrogen"],
     178                LiquidModel = "IdealLiquid",
     179                VapourModel = "PR"
     180        );
     181        NComp as Integer;
     182
     183        SET
     184        NComp = PP.NumberOfComponents;
     185        stoic = [-1, -1, 1, 1];
     186
     187        SPECIFY
     188        T = 1100 * 'K';
     189        P = 1 * 'bar';
     190        ni = [1, 1, 0, 0] * 'mol';
     191       
     192        # Expected results:
     193        # advance = 0.41;
     194       
     195        OPTIONS
     196        Dynamic = false;
     197end
     198
     199
     200# Water-gas-shift T = 1100 K, P = 1 bar excess of water
     201FlowSheet gibbs_reactor_simple_sample5 as gibbs_reactor_simple
     202        PARAMETERS
     203        PP as Plugin(Brief="Physical Properties",
     204                Type="PP",
     205                Components = ["carbon monoxide", "water", "carbon dioxide", "hydrogen"],
     206                LiquidModel = "IdealLiquid",
     207                VapourModel = "PR"
     208        );
     209        NComp as Integer;
     210
     211        SET
     212        NComp = PP.NumberOfComponents;
     213        stoic = [-1, -1, 1, 1];
     214
     215        SPECIFY
     216        T = 1100 * 'K';
     217        P = 1 * 'bar';
     218        ni = [1, 2, 0, 0] * 'mol';
     219       
     220        # Expected results:
     221        # advance = 0.56;
     222       
     223        OPTIONS
     224        Dynamic = false;
     225end
     226
     227# Ammonia synthesis from nitrogen and hydrogen T = 500 degC, P = 1 bar
     228FlowSheet gibbs_reactor_simple_sample6 as gibbs_reactor_simple
     229        PARAMETERS
     230        PP as Plugin(Brief="Physical Properties",
     231                Type="PP",
     232                Components = ["nitrogen", "hydrogen", "ammonia"],
     233                LiquidModel = "IdealLiquid",
     234                VapourModel = "PR"
     235        );
     236        NComp as Integer;
     237
     238        SET
     239        NComp = PP.NumberOfComponents;
     240        stoic = [-1, -1, 2];
     241
     242        SPECIFY
     243        T = (500 + 273.15) * 'K';
     244        P = 1 * 'bar';
     245        ni = [1, 3, 0] * 'mol';
     246       
     247        # Expected results:
     248        # K = 6e-5;
     249       
     250        OPTIONS
     251        Dynamic = false;
     252end
     253
     254# Ammonia synthesis from nitrogen and hydrogen T = 500 degC, P = 300 bar
     255FlowSheet gibbs_reactor_simple_sample7 as gibbs_reactor_simple
     256        PARAMETERS
     257        PP as Plugin(Brief="Physical Properties",
     258                Type="PP",
     259                Components = ["nitrogen", "hydrogen", "ammonia"],
     260                LiquidModel = "IdealLiquid",
     261                VapourModel = "PR"
     262        );
     263        NComp as Integer;
     264
     265        SET
     266        NComp = PP.NumberOfComponents;
     267        stoic = [-1, -3, 2];
     268
     269        SPECIFY
     270        T = (500 + 273.15) * 'K';
     271        P = 300 * 'bar';
     272        ni = [1, 3, 0] * 'mol';
     273       
     274        # Expected results:
     275        # advance = 0.54;
     276       
     277        OPTIONS
     278        Dynamic = false;
     279end
Note: See TracChangeset for help on using the changeset viewer.