 Timestamp:
 Dec 5, 2006, 2:46:19 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

mso/eml/reactors/pfr.mso
r47 r69 45 45 q(NDisc) as heat_rate; 46 46 M(NComp, NDisc) as mol (Brief = "Molar holdup"); 47 C(NComp, NDisc) as conc_mol(Brief="Components concentration"); 47 Mt(NDisc) as mol (Brief = "Molar holdup"); 48 C(NComp, NDisc) as conc_mol(Brief="Components concentration", Lower=1e6); 48 49 E(NDisc) as energy (Brief="Total Energy Holdup on element"); 49 50 … … 70 71 71 72 for z in [1:NDisc] 72 for c in [1:NComp ]73 for c in [1:NComp1] 73 74 "Component Molar Balance" 74 75 diff(M(c,z)) = (str(z).F*str(z).z(c)  str(z+1).F*str(z+1).z(c)) … … 81 82 82 83 "Energy Holdup" 83 E(z) = sum(M(:,z))*str(z+1).h  str(z+1).P*Across*L/NDisc;84 E(z) = Mt(z)*str(z+1).h  str(z+1).P*Across*L/NDisc; 84 85 85 86 "mass flow is considered constant" 86 87 str(z+1).F*vol(z+1) = str(z).F*vol(z); # FIXME: is this correct? No (constant velocity: only for equimolar) 87 #rho(z+1)*vel(z+1) = rho(z)*vel(z); # FIXME: this is correct! But does not converge.88 #rho(z+1)*vel(z+1) = rho(z)*vel(z); # FIXME: this is correct! But does not converge. 88 89 89 90 "Molar concentration" 90 91 C(:,z) * Across*L/NDisc = M(:,z); 92 93 "Sum of M" 94 Mt(z) = sum(M(:,z)); 95 96 "Geometrical constraint" 97 Across*L/NDisc = Mt(z) * vol(z); 91 98 92 99 "Molar fraction" 93 str(z+1).z = C(:,z) * vol(z+1);94 end 100 str(z+1).z * Mt(z) = M(:,z); 101 end 95 102 96 103 for z in [1:NDisc+1] … … 99 106 100 107 "Specific Mass" 101 rho(z) = sum(str(z).z*Mw)/vol(z);108 rho(z) = PP.VapourDensity(str(z).T, str(z).P, str(z).z); 102 109 103 110 "Velocity"
Note: See TracChangeset
for help on using the changeset viewer.