Changeset 69 for mso/eml

Timestamp:

Dec 5, 2006, 2:46:19 PM (17 years ago)

Author:

Rafael de Pelegrini Soares

Message:

Updated PFR model

File:

• mso/eml/reactors/pfr.mso (modified) (4 diffs)

mso/eml/reactors/pfr.mso

@@ -45 +45 @@
         q(NDisc)    as heat_rate;
         M(NComp, NDisc)    as mol (Brief = "Molar holdup");
-        C(NComp, NDisc)  as conc_mol(Brief="Components concentration");
+        Mt(NDisc)    as mol (Brief = "Molar holdup");
+        C(NComp, NDisc)  as conc_mol(Brief="Components concentration", Lower=-1e-6);
         E(NDisc) as energy (Brief="Total Energy Holdup on element");
         
@@ -70 +71 @@
         
         for z in [1:NDisc]
-                for c in [1:NComp]
+                for c in [1:NComp-1]
                         "Component Molar Balance"
                         diff(M(c,z)) = (str(z).F*str(z).z(c) - str(z+1).F*str(z+1).z(c))
@@ -81 +82 @@
 
                 "Energy Holdup"
-                E(z) = sum(M(:,z))*str(z+1).h - str(z+1).P*Across*L/NDisc;
+                E(z) = Mt(z)*str(z+1).h - str(z+1).P*Across*L/NDisc;
 
                 "mass flow is considered constant"
                 str(z+1).F*vol(z+1) = str(z).F*vol(z); # FIXME: is this correct? No (constant velocity: only for equimolar)
-#               rho(z+1)*vel(z+1) = rho(z)*vel(z);  # FIXME: this is correct! But does not converge.
+                #rho(z+1)*vel(z+1) = rho(z)*vel(z);  # FIXME: this is correct! But does not converge.
         
                 "Molar concentration"
                 C(:,z) * Across*L/NDisc = M(:,z);
+                
+                "Sum of M"
+                Mt(z) = sum(M(:,z));
+                
+                "Geometrical constraint"
+                Across*L/NDisc = Mt(z) * vol(z);
 
                 "Molar fraction"
-                str(z+1).z = C(:,z) * vol(z+1);
-        end     
+                str(z+1).z * Mt(z) = M(:,z);
+        end
 
         for z in [1:NDisc+1]
@@ -99 +106 @@
                 
                 "Specific Mass"
-                rho(z) = sum(str(z).z*Mw)/vol(z);
+                rho(z) = PP.VapourDensity(str(z).T, str(z).P, str(z).z);
 
                 "Velocity"