Changeset 516 for branches/rate


Ignore:
Timestamp:
May 15, 2008, 6:37:10 PM (14 years ago)
Author:
Argimiro Resende Secchi
Message:

Updating tray rate.

File:
1 edited

Legend:

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

    r507 r516  
    494494       
    495495        VARIABLES
    496         NL(NComp) as flux_mol   (Brief = "Stream Molar Flux Rate on Liquid Phase");
    497         NV(NComp) as flux_mol   (Brief = "Stream Molar Flux Rate on Vapour Phase");
     496        NL(NComp) as flow_mol   (Brief = "Stream Molar Rate on Liquid Phase");
     497        NV(NComp) as flow_mol   (Brief = "Stream Molar Rate on Vapour Phase");
    498498        T as temperature                (Brief = "Stream Temperature");
    499499        P as pressure                   (Brief = "Stream Pressure");
     
    503503        htL as heat_trans_coeff (Brief = "Heat Transference Coefficient on Liquid Phase");
    504504        htV as heat_trans_coeff (Brief = "Heat Transference Coefficient on Vapour Phase");     
    505         E_liq as heat_flux      (Brief = "Liquid Energy on interface");
    506     E_vap as heat_flux      (Brief = "Vapour Energy on interface");     
     505        E_liq as heat_rate      (Brief = "Liquid Energy Rate at interface");
     506    E_vap as heat_rate      (Brief = "Vapour Energy Rate at interface");       
    507507        hL as enth_mol          (Brief = "Liquid Molar Enthalpy");
    508508        hV as enth_mol          (Brief = "Vapour Molar Enthalpy");
     
    565565        "Component Molar Balance"
    566566        diff(M_liq)=Inlet.F*Inlet.z + InletL.F*InletL.z
    567         - OutletL.F*OutletL.z + interf.a*interf.NL;
     567        - OutletL.F*OutletL.z + interf.NL;
    568568       
    569569        diff(M_vap)=InletFV.F*InletFV.z + InletV.F*InletV.z
    570         - OutletV.F*OutletV.z - interf.a*interf.NV;
     570        - OutletV.F*OutletV.z - interf.NV;
    571571       
    572572        "Energy Balance"
    573573        diff(E_liq) = Inlet.F*Inlet.h + InletL.F*InletL.h
    574                 - OutletL.F*OutletL.h  + Q + interf.a*interf.E_liq;
     574                - OutletL.F*OutletL.h  + Q + interf.E_liq;
    575575       
    576576        diff(E_vap) = InletFV.F*InletFV.h + InletV.F*InletV.h
    577                 - OutletV.F*OutletV.h  - interf.a*interf.E_vap;
     577                - OutletV.F*OutletV.h  - interf.E_vap;
    578578       
    579579        "Molar Holdup"
     
    587587        E_vap = MV*(OutletV.h - OutletV.P*vV);
    588588       
    589         "Energy on interface"
    590         interf.E_liq = interf.htL*(interf.T-OutletL.T)+sum(interf.NL)*interf.hL;       
    591        
    592         interf.E_vap = interf.htV*(OutletV.T-interf.T)+sum(interf.NV)*interf.hV;
     589        "Energy Rate through the interface"
     590        interf.E_liq = interf.htL*interf.a*(interf.T-OutletL.T)+sum(interf.NL)*interf.hL;       
     591       
     592        interf.E_vap = interf.htV*interf.a*(OutletV.T-interf.T)+sum(interf.NV)*interf.hV;
    593593       
    594594        "Mass Conservation"
    595         interf.NL=interf.NV;
     595        interf.NL = interf.NV;
    596596       
    597597        "Energy Conservation"
     
    619619        Level = ML*vL/Ap;
    620620
    621         "Total Mass Transfer Fluxes"
    622         interf.NL(1:NC1)=sumt(interf.kL*(interf.x(1:NC1)-OutletL.z(1:NC1)))/vL+
     621        "Total Mass Transfer Rates"
     622        interf.NL(1:NC1)=interf.a*sumt(interf.kL*(interf.x(1:NC1)-OutletL.z(1:NC1)))/vL+
    623623                OutletL.z(1:NC1)*sum(interf.NL);
    624624       
    625         interf.NV(1:NC1)=sumt(interf.kV*(OutletV.z(1:NC1)-interf.y(1:NC1)))/vV+
     625        interf.NV(1:NC1)=interf.a*sumt(interf.kV*(OutletV.z(1:NC1)-interf.y(1:NC1)))/vV+
    626626                OutletV.z(1:NC1)*sum(interf.NV);
    627627
Note: See TracChangeset for help on using the changeset viewer.