Changeset 410 for trunk

Timestamp:

Nov 20, 2007, 1:30:10 PM (16 years ago)

Author:

Rafael de Pelegrini Soares

Message:

Added symbols for the simple gibbs reactor

File:

• trunk/sample/miscellaneous/sample_gibbs_reactor_simple.mso (modified) (1 diff)

trunk/sample/miscellaneous/sample_gibbs_reactor_simple.mso

@@ -43 +43 @@
 outer NComp     as Integer (Brief="Number of components", Default=1);
 
-        nu(NComp) as Real;
+        nu(NComp) as Real(Symbol="\nu_i");
 
         R  as Real(Brief="Universal gas constant", Unit='J/mol/K', Default=8.314);
-        T0 as temperature(Default = 298.15);
-        P0 as Real(Default=1, Unit='bar');
-        G0(NComp) as energy_mol         (Brief="Gibbs energy in standard state");
-        H0(NComp) as energy_mol         (Brief="Enthalpy in standard state");
-        
-        SET
-        G0 = PP.IdealGasGibbsOfFormationAt25C();
-        H0 = PP.IdealGasEnthalpyOfFormationAt25C();
+        T0 as temperature(Default = 298.15, Symbol="T_{298}");
+        P0 as Real(Default=1, Unit='bar', Symbol="P^0");
+        g0(NComp) as energy_mol         (Brief="Gibbs energy in standard state", Symbol="g^0_{298,i}");
+        h0(NComp) as energy_mol         (Brief="Enthalpy in standard state", Symbol="h^0_{298,i}");
+        
+        SET
+        g0 = PP.IdealGasGibbsOfFormationAt25C();
+        h0 = PP.IdealGasEnthalpyOfFormationAt25C();
 
         VARIABLES
         T as temperature;
         P as pressure;
-        n0(NComp) as positive(Brief="Initial number of mols", Unit='mol');
-        n(NComp)  as positive(Brief="Number of mols at equilibrium", Unit='mol');
-        advance   as Real(Unit='mol');
+        n0(NComp) as positive(Brief="Initial number of mols", Unit='mol', Symbol="n_{i,0}");
+        n(NComp)  as positive(Brief="Number of mols at equilibrium", Unit='mol', Symbol="n_i");
+        advance   as Real(Unit='mol', Symbol="\epsilon");
 
         K  as positive(Brief="Reaction equilibrium constant at T");
-        K0 as positive(Brief="Reaction equilibrium constant at T0");
-        
-        phi(NComp)      as fugacity(Brief="Fugacity coefficient", Default=1);
+        K0 as positive(Brief="Reaction equilibrium constant at T0", Symbol="K_{298}");
+        
+        phi(NComp)      as fugacity(Brief="Fugacity coefficient", Default=1, Symbol="\phi_i");
         
         EQUATIONS
         "Equilibrium constant at 298 K"
-        K0 = exp(-sum(nu*G0)/(R*T0));
+        K0 = exp(-sum(nu*g0)/(R*T0));
         "Equilibrium constant at T"
-        K = K0 * exp(-sum(nu*H0)/R*(1/T - 1/T0));
+        K = K0 * exp(-sum(nu*h0)/R*(1/T - 1/T0));
 
         "Equilibrium rule"