Changeset 516
 Timestamp:
 May 15, 2008, 6:37:10 PM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/rate/eml/stage_separators/tray.mso
r507 r516 494 494 495 495 VARIABLES 496 NL(NComp) as fl ux_mol (Brief = "Stream Molar FluxRate on Liquid Phase");497 NV(NComp) as fl ux_mol (Brief = "Stream Molar FluxRate 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"); 498 498 T as temperature (Brief = "Stream Temperature"); 499 499 P as pressure (Brief = "Stream Pressure"); … … 503 503 htL as heat_trans_coeff (Brief = "Heat Transference Coefficient on Liquid Phase"); 504 504 htV as heat_trans_coeff (Brief = "Heat Transference Coefficient on Vapour Phase"); 505 E_liq as heat_ flux (Brief = "Liquid Energy oninterface");506 E_vap as heat_ flux (Brief = "Vapour Energy oninterface");505 E_liq as heat_rate (Brief = "Liquid Energy Rate at interface"); 506 E_vap as heat_rate (Brief = "Vapour Energy Rate at interface"); 507 507 hL as enth_mol (Brief = "Liquid Molar Enthalpy"); 508 508 hV as enth_mol (Brief = "Vapour Molar Enthalpy"); … … 565 565 "Component Molar Balance" 566 566 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; 568 568 569 569 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; 571 571 572 572 "Energy Balance" 573 573 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; 575 575 576 576 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; 578 578 579 579 "Molar Holdup" … … 587 587 E_vap = MV*(OutletV.h  OutletV.P*vV); 588 588 589 "Energy oninterface"590 interf.E_liq = interf.htL* (interf.TOutletL.T)+sum(interf.NL)*interf.hL;591 592 interf.E_vap = interf.htV* (OutletV.Tinterf.T)+sum(interf.NV)*interf.hV;589 "Energy Rate through the interface" 590 interf.E_liq = interf.htL*interf.a*(interf.TOutletL.T)+sum(interf.NL)*interf.hL; 591 592 interf.E_vap = interf.htV*interf.a*(OutletV.Tinterf.T)+sum(interf.NV)*interf.hV; 593 593 594 594 "Mass Conservation" 595 interf.NL =interf.NV;595 interf.NL = interf.NV; 596 596 597 597 "Energy Conservation" … … 619 619 Level = ML*vL/Ap; 620 620 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+ 623 623 OutletL.z(1:NC1)*sum(interf.NL); 624 624 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+ 626 626 OutletV.z(1:NC1)*sum(interf.NV); 627 627
Note: See TracChangeset
for help on using the changeset viewer.