Changeset 507 for branches


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

Updating tray rate model.

Location:
branches/rate
Files:
2 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
  • branches/rate/sample/stage_separators/sample_tray.mso

    r506 r507  
    271271        t1.interf.kL=0.001*'m/s';
    272272        t1.interf.kV=0.001*'m/s';
    273         t1.ML=3*'mol';
    274         t1.MV=3*'mol';
    275273
    276274        SET
     
    285283        t1.Q = 0 * 'kW';
    286284       
    287        
    288285        INITIAL
    289286        t1.OutletL.T = 290 *'K';
    290287        t1.OutletV.T = 290 *'K';
    291         #t1.Level = 0.9*'m' ;
     288        t1.Level = 0.9*'m' ;
    292289        t1.OutletL.z(1) = 0.5;
    293290        t1.OutletV.z(1) = 0.5;
    294         #t1.MV = 0.01*'mol' ;
    295        
    296         OPTIONS
     291        t1.MV = 0.01*'mol' ;
     292       
     293        OPTIONS
     294        TimeStep = 0.1;
    297295        TimeEnd = 100;
    298 end
     296        TimeUnit = 's';
     297#       GuessFile = "tray_Test.rlt";
     298end
Note: See TracChangeset for help on using the changeset viewer.