Changeset 44 for mso/eml/reactors

Timestamp:

Nov 5, 2006, 3:58:23 PM (17 years ago)

Author:

Argimiro Resende Secchi

Message:

Fix part of PFR.

File:

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

mso/eml/reactors/pfr.mso

@@ -13 +13 @@
 
 #*
-* Generic PFR model
+* Generic PFR model with constant mass holdup
 *
 * Requires the information of:
@@ -83 +83 @@
                 E(z) = sum(M(:,z))*str(z+1).h - str(z+1).P*Across*L/NDisc;
 
-                "Volume flow is considered constant"
-                str(z+1).F*vol(z+1) = str(z).F*vol(z); # FIXME: is this correct?
-                
+                "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.
+        
                 "Molar concentration"
                 C(:,z) * Across*L/NDisc = M(:,z);
+
                 "Molar fraction"
-                str(z+1).z = C(:,z) * vol(z); # FIXME: is this correct?
+                str(z+1).z = C(:,z) * vol(z+1);
         end     
-        
+
         for z in [1:NDisc+1]
                 "Specific Volume"
                 vol(z) = PP.VapourVolume(str(z).T, str(z).P, str(z).z);
+                
                 "Specific Mass"
                 rho(z) = sum(str(z).z*Mw)/vol(z);