Changeset 706 for trunk/eml


Ignore:
Timestamp:
Jan 26, 2009, 8:08:28 PM (14 years ago)
Author:
Argimiro Resende Secchi
Message:

Removing temporary variable from tray model.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/eml/stage_separators/tray.mso

    r522 r706  
    144144        beta as fraction (Brief="Aeration fraction");
    145145        alfa as fraction (Brief="Dry pressure drop coefficient");
    146         w as Real (Brief="Feehery's correlation coefficient", Unit='1/m^4', Default=1);
    147         btray as Real (Brief="Elgue's correlation coefficient", Unit='kg/m/mol^2', Default=1);
    148         fw as Real      (Brief = "Olsen's correlation coefficient", Default=1);
     146        w as Real (Brief="Feehery correlation coefficient", Unit='1/m^4', Default=1);
     147        btray as Real (Brief="Elgue correlation coefficient", Unit='kg/m/mol^2', Default=1);
     148        fw as Real      (Brief = "Olsen correlation coefficient", Default=1);
    149149        Np as Real      (Brief = "Number of liquid passes in the tray", Default=1);
    150150        Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
     
    161161        rhoL as dens_mass;
    162162        rhoV as dens_mass;
    163 
    164         btemp as Real (Brief="Temporary variable of Roffel's liquid flow equation");
    165163       
    166164        EQUATIONS
     
    187185                       
    188186                                case "Roffel_Fl":
    189                                 OutletL.F = 2/3*sqrt(2*g)*rhoL/sum(Mw*OutletL.z)*lw*(2*btemp-1)*(ML*sum(Mw*OutletL.z)/(Ap*1.3)/rhoL/(2*btemp-1))^1.5;
     187                                OutletL.F = 2/3*rhoL/sum(Mw*OutletL.z)*lw*(ML*sum(Mw*OutletL.z)/(Ap*1.3)/rhoL)^1.5*sqrt(2*g/
     188                                                        (2*(1 - 0.3593/'Pa^0.0888545'*(OutletV.F*sum(Mw*OutletV.z)/(Ap*1.3)/sqrt(rhoV))^0.177709)-1)); 
    190189                        end
    191190                when Level < (beta * hw) switchto "off";
     
    197196        end
    198197       
    199         btemp = 1 - 0.3593/'Pa^0.0888545'*(OutletV.F*sum(Mw*OutletV.z)/(Ap*1.3)/sqrt(rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709';
    200 
    201198        switch VapourFlow
    202199                case "on":
Note: See TracChangeset for help on using the changeset viewer.