Changeset 507 for branches/rate/eml


Ignore:
Timestamp:
Apr 28, 2008, 4:17:48 PM (15 years ago)
Author:
Argimiro Resende Secchi
Message:

Updating tray rate model.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/rate/eml/stage_separators/tray.mso

    r505 r507  
    507507        hL as enth_mol          (Brief = "Liquid Molar Enthalpy");
    508508        hV as enth_mol          (Brief = "Vapour Molar Enthalpy");
    509         ConcL as conc_mol       (Brief = "Liquid Molar Concentration");
    510         ConcV as conc_mol       (Brief = "Vapour Molar Concentration");
    511509        kL(NC1,NC1) as velocity (Brief = "Mass Transfer Coefficients");
    512510        kV(NC1,NC1) as velocity (Brief = "Mass Transfer Coefficients");
     
    600598        interf.E_liq = interf.E_vap;
    601599       
    602        
    603600        "Mol fraction normalisation"
    604601        sum(OutletL.z)= 1.0;
     
    621618        "Level of clear liquid over the weir"
    622619        Level = ML*vL/Ap;
    623        
     620
    624621        "Total Mass Transfer Fluxes"
    625         interf.NL(1:NC1)=interf.ConcL*sumt(interf.kL*(interf.x(1:NC1)-OutletL.z(1:NC1)))+
     622        interf.NL(1:NC1)=sumt(interf.kL*(interf.x(1:NC1)-OutletL.z(1:NC1)))/vL+
    626623                OutletL.z(1:NC1)*sum(interf.NL);
    627624       
    628         interf.NV(1:NC1)=interf.ConcV*sumt(interf.kV*(OutletV.z(1:NC1)-interf.y(1:NC1)))+
     625        interf.NV(1:NC1)=sumt(interf.kV*(OutletV.z(1:NC1)-interf.y(1:NC1)))/vV+
    629626                OutletV.z(1:NC1)*sum(interf.NV);
    630627
Note: See TracChangeset for help on using the changeset viewer.