 Timestamp:
 Nov 5, 2006, 3:58:23 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

mso/eml/reactors/pfr.mso
r36 r44 13 13 14 14 #* 15 * Generic PFR model 15 * Generic PFR model with constant mass holdup 16 16 * 17 17 * Requires the information of: … … 83 83 E(z) = sum(M(:,z))*str(z+1).h  str(z+1).P*Across*L/NDisc; 84 84 85 "Volume flow is considered constant" 86 str(z+1).F*vol(z+1) = str(z).F*vol(z); # FIXME: is this correct? 87 85 "mass flow is considered constant" 86 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 88 89 "Molar concentration" 89 90 C(:,z) * Across*L/NDisc = M(:,z); 91 90 92 "Molar fraction" 91 str(z+1).z = C(:,z) * vol(z ); # FIXME: is this correct?93 str(z+1).z = C(:,z) * vol(z+1); 92 94 end 93 95 94 96 for z in [1:NDisc+1] 95 97 "Specific Volume" 96 98 vol(z) = PP.VapourVolume(str(z).T, str(z).P, str(z).z); 99 97 100 "Specific Mass" 98 101 rho(z) = sum(str(z).z*Mw)/vol(z);
Note: See TracChangeset
for help on using the changeset viewer.