Changeset 44 for mso/eml


Ignore:
Timestamp:
Nov 5, 2006, 3:58:23 PM (16 years ago)
Author:
Argimiro Resende Secchi
Message:

Fix part of PFR.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • mso/eml/reactors/pfr.mso

    r36 r44  
    1313
    1414#*
    15 * Generic PFR model
     15* Generic PFR model with constant mass holdup
    1616*
    1717* Requires the information of:
     
    8383                E(z) = sum(M(:,z))*str(z+1).h - str(z+1).P*Across*L/NDisc;
    8484
    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       
    8889                "Molar concentration"
    8990                C(:,z) * Across*L/NDisc = M(:,z);
     91
    9092                "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);
    9294        end     
    93        
     95
    9496        for z in [1:NDisc+1]
    9597                "Specific Volume"
    9698                vol(z) = PP.VapourVolume(str(z).T, str(z).P, str(z).z);
     99               
    97100                "Specific Mass"
    98101                rho(z) = sum(str(z).z*Mw)/vol(z);
Note: See TracChangeset for help on using the changeset viewer.