Changeset 100


Ignore:
Timestamp:
Jan 9, 2007, 12:15:56 PM (15 years ago)
Author:
gerson bicca
Message:

updated

Location:
mso
Files:
2 added
21 edited

Legend:

Unmodified
Added
Removed
  • mso/eml/heat_exchangers/DoublePipe.mso

    r78 r100  
    1919
    2020using "HEX_Engine";
    21 #=====================================================================
     21#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    2222#       Basic Models for Double Pipe Heat Exchangers
    23 #=====================================================================
     23#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    2424
    2525Model DoublePipe_Basic
     
    7070        Properties.Cold.Average.Mw = sum(M*Inlet.Cold.z);
    7171
    72 
    73 
    7472if Inlet.Cold.v equal 0
    7573        then   
     
    208206end
    209207
    210 #=====================================================================
     208#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    211209#       Thermal Details
    212 #=====================================================================
     210#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    213211"Hot Stream Heat Capacity"
    214212        Details.Ch =Inlet.Hot.F*Properties.Hot.Average.Cp;
     
    225223"Heat Capacity Ratio"   
    226224        Details.Cr*Details.Cmax   = Details.Cmin;
    227 #=====================================================================
     225#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    228226#       Energy Balance
    229 #=====================================================================
     227#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    230228"Energy Balance Hot Stream"
    231229        Details.Q = Inlet.Hot.F*(Inlet.Hot.h-Outlet.Hot.h);
     
    234232        Details.Q = Inlet.Cold.F*(Outlet.Cold.h - Inlet.Cold.h);
    235233
    236 #=====================================================================
     234#--------------------------------------------------------------------
    237235#       Material Balance
    238 #=====================================================================
     236#--------------------------------------------------------------------
    239237"Flow Mass Inlet Cold Stream"
    240238        Properties.Cold.Inlet.Fw        =  sum(M*Inlet.Cold.z)*Inlet.Cold.F;
     
    255253        Inlet.Cold.F = Outlet.Cold.F;
    256254
    257 #======================================
     255#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    258256#       Constraints
    259 #======================================
     257#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    260258"Hot Stream Molar Fraction Constraint"
    261259        Outlet.Hot.z=Inlet.Hot.z;
     
    294292
    295293end
    296 
    297 "Overall Heat Transfer Coefficient"
    298 #       Details.U*(Resistances.Rtube+Resistances.Rwall+Resistances.Rshell)=1;
    299         Details.U=1/(Resistances.Rtube+Resistances.Rwall+Resistances.Rshell);
    300294
    301295end
     
    355349"Outer Pipe Film Coefficient"
    356350        Outer.HeatTransfer.hcoeff= HE.PipeFilmCoeff(Outer.HeatTransfer.Re,Outer.HeatTransfer.PR,Properties.Hot.Average.K,Outer.HeatTransfer.Dh,Lpipe)*Outer.HeatTransfer.Phi;
    357 #       Outer.HeatTransfer.hcoeff= (0.027*Outer.HeatTransfer.Re^(4/5)*Outer.HeatTransfer.PR^(1/3)*Properties.Hot.Average.K/Outer.HeatTransfer.Dh)*Outer.HeatTransfer.Phi;
    358351
    359352"Inner Pipe Film Coefficient"
    360353        Inner.HeatTransfer.hcoeff= HE.PipeFilmCoeff(Inner.HeatTransfer.Re,Inner.HeatTransfer.PR,Properties.Cold.Average.K,DiInner,Lpipe)*Inner.HeatTransfer.Phi;
    361 #       Inner.HeatTransfer.hcoeff= (0.027*Inner.HeatTransfer.Re^(4/5)*Inner.HeatTransfer.PR^(1/3)*Properties.Cold.Average.K/DiInner)*Inner.HeatTransfer.Phi;
    362354
    363355"Outer Pipe Pressure Drop"
     
    407399"Inner Pipe Film Coefficient"
    408400        Inner.HeatTransfer.hcoeff= HE.PipeFilmCoeff(Inner.HeatTransfer.Re,Inner.HeatTransfer.PR,Properties.Hot.Average.K,DiInner,Lpipe)*Inner.HeatTransfer.Phi;
    409 #       Inner.HeatTransfer.hcoeff= (0.027*Inner.HeatTransfer.Re^(4/5)*Inner.HeatTransfer.PR^(1/3)*Properties.Hot.Average.K/DiInner)*Inner.HeatTransfer.Phi;
    410401
    411402"Outer Pipe Film Coefficient"
    412403        Outer.HeatTransfer.hcoeff= HE.PipeFilmCoeff(Outer.HeatTransfer.Re,Outer.HeatTransfer.PR,Properties.Cold.Average.K,Outer.HeatTransfer.Dh,Lpipe)*Outer.HeatTransfer.Phi;
    413 #       Outer.HeatTransfer.hcoeff= (0.027*Outer.HeatTransfer.Re^(4/5)*Outer.HeatTransfer.PR^(1/3)*Properties.Cold.Average.K/Outer.HeatTransfer.Dh)*Outer.HeatTransfer.Phi;
    414404
    415405"Outer Pipe Pressure Drop"
     
    446436        Outer.HeatTransfer.Vmean*(Outer.HeatTransfer.As*Properties.Cold.Average.rho)= Properties.Cold.Inlet.Fw;
    447437       
    448        
    449438"Inner Pipe Velocity"
    450439        Inner.HeatTransfer.Vmean*(Inner.HeatTransfer.As*Properties.Hot.Average.rho)     = Properties.Hot.Inlet.Fw;
     
    461450        Resistances.Rshell*(Outer.HeatTransfer.hcoeff)=1;
    462451
    463 
     452"Overall Heat Transfer Coefficient Clean"
     453        Details.Uc*(Resistances.Rtube+Resistances.Rwall+Resistances.Rshell)=1;
     454
     455"Overall Heat Transfer Coefficient Dirty"
     456        Details.Ud*(Resistances.Rfi*(DoInner/DiInner) + Resistances.Rfo + Resistances.Rtube + Resistances.Rwall + Resistances.Rshell)=1;
     457       
    464458end
    465459
    466460Model DoublePipe_Basic_NTU                      as DoublePipe
    467 #=====================================================================
     461#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    468462#       Basic Model Double Pipe Heat Exchanger - NTU Method
    469 #=====================================================================
     463#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    470464VARIABLES
    471465
     
    481475
    482476Model DoublePipe_Basic_LMTD                     as DoublePipe
    483 #=====================================================================
     477#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    484478#       Basic Model for Double Pipe Heat Exchanger- LMTD Method
    485 #=====================================================================
     479#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    486480VARIABLES
    487481
     
    518512
    519513"Exchange Surface Area"
    520         Details.Q = Details.U*Pi*DoInner*Lpipe*LMTD;
     514        Details.Q = Details.Ud*Pi*DoInner*Lpipe*LMTD;
    521515
    522516end
     
    552546       
    553547        then   
    554 "Effectiveness in Cocurrent Flow"
     548"Effectiveness"
    555549        Eft = 1-exp(-Details.NTU);
    556550       
     
    599593        Kwall           as conductivity (Brief="Tube Wall Material Thermal Conductivity",Default=1.0);
    600594
     595VARIABLES
     596
     597Unity(Npipe)  as DoublePipe_Basic;
     598
    601599SET
    602600        Pi      = 3.14159265;
    603601        Hside   = HE.FluidAlocation();
    604 
    605 VARIABLES
    606 
    607 Unity(Npipe)  as DoublePipe_Basic;
     602       
     603#"Inner Pipe Cross Sectional Area for Flow"
     604        Unity.Inner.HeatTransfer.As=Pi*DiInner*DiInner/4;
     605       
     606#"Outer Pipe Cross Sectional Area for Flow"
     607        Unity.Outer.HeatTransfer.As=Pi*(DiOuter*DiOuter-DoInner*DoInner)/4;
     608       
     609#"Inner Pipe Hydraulic Diameter for Heat Transfer"
     610        Unity.Inner.HeatTransfer.Dh=DiInner;
     611       
     612#"Outer Pipe Hydraulic Diameter for Heat Transfer"
     613        Unity.Outer.HeatTransfer.Dh=(DiOuter*DiOuter-DoInner*DoInner)/DoInner;
     614
     615#"Inner Pipe Hydraulic Diameter for Pressure Drop"
     616        Unity.Inner.PressureDrop.Dh=DiInner;
     617       
     618#"Outer Pipe Hydraulic Diameter for Pressure Drop"
     619        Unity.Outer.PressureDrop.Dh=DiOuter-DoInner;
    608620
    609621EQUATIONS
    610622
    611623for i in [1:Npipe]
     624
     625"Overall Heat Transfer Coefficient Clean"
     626        Unity(i).Details.Uc*(Unity(i).Resistances.Rtube+Unity(i).Resistances.Rwall+Unity(i).Resistances.Rshell)=1;
     627
     628"Overall Heat Transfer Coefficient Dirty"
     629        Unity(i).Details.Ud*(Unity(i).Resistances.Rfi*(DoInner/DiInner) + Unity(i).Resistances.Rfo + Unity(i).Resistances.Rtube + Unity(i).Resistances.Rwall + Unity(i).Resistances.Rshell)=1;
    612630
    613631"Exchange Surface Area"
    614632        Unity(i).Details.A=Pi*DoInner*Lpipe;
    615633       
    616 "Inner Pipe Cross Sectional Area for Flow"
    617         Unity(i).Inner.HeatTransfer.As=Pi*DiInner*DiInner/4;
    618        
    619 "Outer Pipe Cross Sectional Area for Flow"
    620         Unity(i).Outer.HeatTransfer.As=Pi*(DiOuter*DiOuter-DoInner*DoInner)/4;
    621        
    622 "Inner Pipe Hydraulic Diameter for Heat Transfer"
    623         Unity(i).Inner.HeatTransfer.Dh=DiInner;
    624        
    625 "Outer Pipe Hydraulic Diameter for Heat Transfer"
    626         Unity(i).Outer.HeatTransfer.Dh=(DiOuter*DiOuter-DoInner*DoInner)/DoInner;
    627 
    628 "Inner Pipe Hydraulic Diameter for Pressure Drop"
    629         Unity(i).Inner.PressureDrop.Dh=DiInner;
    630        
    631 "Outer Pipe Hydraulic Diameter for Pressure Drop"
    632         Unity(i).Outer.PressureDrop.Dh=DiOuter-DoInner;
    633 
    634634if Hside equal 1
    635635       
     
    693693
    694694"Outer Pipe Phi correction"
    695         Unity(i).Outer.HeatTransfer.Phi                 = HE.PhiCorrection(Unity(i).Properties.Cold.Average.Mu,Unity(i).Properties.Cold.Wall.Mu);
     695        Unity(i).Outer.HeatTransfer.Phi         = HE.PhiCorrection(Unity(i).Properties.Cold.Average.Mu,Unity(i).Properties.Cold.Wall.Mu);
    696696       
    697697"Inner Pipe Phi correction"
     
    739739
    740740Model Multitubular_Basic_LMTD           as Multitubular_Basic
    741 #=====================================================================
     741#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    742742#       Basic Model for Double Pipe Heat Exchanger- LMTD Method
    743 #=====================================================================
     743#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
    744744VARIABLES
    745745
     
    777777
    778778"Exchange Surface Area"
    779         Unity(i).Details.Q = Unity(i).Details.U*Unity(i).Details.A*LMTD(i);
     779        Unity(i).Details.Q = Unity(i).Details.Ud*Unity(i).Details.A*LMTD(i);
    780780
    781781end
  • mso/eml/heat_exchangers/HEX_Engine.mso

    r78 r100  
    146146#===================================================================== 
    147147VARIABLES
    148 
    149148PdTube                  as press_delta  (Brief="Tube Pressure Drop",Default=0.01, Lower=1e-10);
    150149Pdtotal                 as press_delta  (Brief="Total Pressure Drop",Default=0.01, Lower=1e-10);
     
    186185Vnozzle_in      as velocity             (Brief="Inlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
    187186Vnozzle_out     as velocity             (Brief="Outlet Nozzle Velocity",Default=1, Upper=1e5, Lower=0);
    188 RVsquare_in     as flux_mass    (Brief="Inlet Nozzle rho-V^2");
    189 RVsquare_out    as flux_mass    (Brief="Outlet Nozzle rho-V^2");
     187RVsquare_out    as positive (Brief = "Outlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = "kg/s^2/m");
     188RVsquare_in     as positive (Brief = "Inlet Nozzle rho-V^2", Default=1, Upper=1e6, Unit = "kg/s^2/m");
    190189
    191190EQUATIONS
    192 
    193191"Shell Side Total Pressure Drop"
    194192        Pdtotal = PdCross + PdEndZones + Pdnozzle_in + Pdnozzle_out + Pdwindow;
     
    233231#===================================================================== 
    234232VARIABLES
    235 
    236233Rtube   as positive     (Brief="Tube Resistance",Unit="m^2*K/kW",Lower=1e-6);
    237234Rwall   as positive     (Brief="Wall Resistance",Unit="m^2*K/kW",Lower=1e-6);
    238235Rshell  as positive     (Brief="Shell Resistance",Unit="m^2*K/kW",Lower=1e-6);
    239 
     236Rfi             as positive     (Brief="Inside Fouling Resistance",Unit="m^2*K/kW",Default=1e-6,Lower=0);
     237Rfo             as positive     (Brief="Outside Fouling Resistance",Unit="m^2*K/kW",Default=1e-6,Lower=0);
    240238end
    241239
     
    247245A               as area                         (Brief="Exchange Surface Area");
    248246Q               as power                        (Brief="Heat Transfer", Default=7000, Lower=1e-6, Upper=1e10);
    249 U               as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient",Default=1,Lower=1e-6,Upper=1e10);
     247Uc              as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Clean",Default=1,Lower=1e-6,Upper=1e10);
     248Ud              as heat_trans_coeff (Brief="Overall Heat Transfer Coefficient Dirty",Default=1,Lower=1e-6,Upper=1e10);
    250249Ch      as positive                     (Brief="Hot Stream Heat Capacity",Lower=1e-3,Default=1e3,Unit="W/K");
    251250Cc      as positive                     (Brief="Cold Stream Heat Capacity",Lower=1e-3,Default=1e3,Unit="W/K");
     
    257256EQUATIONS
    258257"Number of Units Transference"
    259         NTU*Cmin = U*A;
     258        NTU*Cmin = Ud*A;
    260259end
    261260
     
    298297Model DoublePipe_HeatTransfer
    299298#=====================================================================
    300 #       Heat Transfer
     299#       Double Pipe Heat Transfer Block
    301300#===================================================================== 
    302301PARAMETERS
     
    314313Model DoublePipe_PressureDrop
    315314#=====================================================================
    316 #       Pressure Drop
     315#       #       Double Pipe Pressure Drop Block
    317316#===================================================================== 
    318317PARAMETERS
     
    326325
    327326Model Main_DoublePipe
     327#=====================================================================
     328#       Double Pipe Main Variables
     329#===================================================================== 
    328330VARIABLES
    329331HeatTransfer as DoublePipe_HeatTransfer;
  • mso/eml/heat_exchangers/HeatExchangerDetailed.mso

    r78 r100  
    640640        Shell.HeatTransfer.Sm = HE.CrossFlowArea(Baffles.Ls);
    641641       
    642 "Overall Heat Transfer Coefficient"
    643         Details.U=1/(Dotube/(Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Shell.HeatTransfer.hshell)));
     642"Overall Heat Transfer Coefficient Dirty"
     643        Details.Ud=1/(Dotube/(Tubes.HeatTransfer.htube*Ditube)+Resistances.Rfo+Resistances.Rfi*(Dotube/Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Shell.HeatTransfer.hshell)));
     644
     645"Overall Heat Transfer Coefficient Clean"
     646        (1/Details.Ud)=(1/Details.Uc)+Resistances.Rfo+Resistances.Rfi*(Dotube/Ditube);
    644647
    645648"Exchange Surface Area"
     
    700703
    701704"Exchange Surface Area"
    702         Details.Q   = Details.U*Details.A*MTD;
     705        Details.Q   = Details.Ud*Details.A*MTD;
    703706       
    704707"Mean Temperature Difference"   
     
    10841087        Unity(i).Resistances.Rshell*(Unity(i).Shell.HeatTransfer.hshell)=1;
    10851088       
    1086 "Overall Heat Transfer Coefficient"
    1087         Unity(i).Details.U=1/(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)));
     1089"Overall Heat Transfer Coefficient Clean"
     1090        Unity(i).Details.Uc=1/(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)));
     1091
     1092"Overall Heat Transfer Coefficient Dirty"
     1093        (1/Unity(i).Details.Ud)=(1/Unity(i).Details.Uc)+Unity(i).Resistances.Rfo+Unity(i).Resistances.Rfi*(Dotube/Ditube);
    10881094
    10891095"Exchange Surface Area"
     
    15601566        Unity(i).Resistances.Rshell*(Unity(i).Shell.HeatTransfer.hshell)=1;
    15611567       
    1562 "Overall Heat Transfer Coefficient"
    1563         Unity(i).Details.U*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;
     1568"Overall Heat Transfer Coefficient Clean"
     1569        Unity(i).Details.Uc*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;
     1570
     1571"Overall Heat Transfer Coefficient Dirty"
     1572        (1/Unity(i).Details.Ud)=(1/Unity(i).Details.Uc)+Unity(i).Resistances.Rfo+Unity(i).Resistances.Rfi*(Dotube/Ditube);
    15641573
    15651574"Exchange Surface Area"
     
    16191628
    16201629"Exchange Surface Area"
    1621 #       Unity(i).Details.Q   = Unity(i).Details.U*Unity(i).Details.A*Fc(i)*LMTD(i);
    1622         Unity(i).Details.Q   = Unity(i).Details.U*Pi*Dotube*Ntt*Ltube*Fc(i)*LMTD(i);
     1630        Unity(i).Details.Q   = Unity(i).Details.Ud*Pi*Dotube*Ntt*Ltube*Fc(i)*LMTD(i);
    16231631       
    16241632"Mean Temperature Difference"   
  • mso/eml/heat_exchangers/HeatExchangerDiscretized.mso

    r78 r100  
    365365EQUATIONS
    366366"Exchange Surface Area"
    367         Details.Q   = Details.U*Details.A*MTD;
     367        Details.Q   = Details.Ud*Details.A*MTD;
    368368       
    369369"Mean Temperature Difference"   
     
    424424PdropShellNozzle        as pressure                     (Brief="Total Shell Side Nozzles Pressure Drop");
    425425PdropShell                      as pressure                     (Brief="Total Shell Side Pressure Drop");
    426 Uaverage                        as heat_trans_coeff (Brief="Average Overall Heat Transfer Coefficient",Default=1,Lower=1e-6,Upper=1e10);
     426Udaverage                       as heat_trans_coeff (Brief="Average Overall Heat Transfer Coefficient Dirty",Default=1,Lower=1e-6,Upper=1e10);
     427Ucaverage                       as heat_trans_coeff (Brief="Average Overall Heat Transfer Coefficient Clean",Default=1,Lower=1e-6,Upper=1e10);
    427428hshellaverage           as heat_trans_coeff     (Brief="Average Shell Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
    428429htubeaverage            as heat_trans_coeff     (Brief="Average Tube Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
     
    506507        Sumary.htubeaverage             = sum(Unity.Tubes.HeatTransfer.htube)/Sumary.Zones;
    507508
    508 "Average Overall Heat Transfer Coefficient"
    509         Sumary.Uaverage                 = sum(Unity.Details.U)/Sumary.Zones;
     509"Average Overall Heat Transfer Coefficient Dirty"
     510        Sumary.Udaverage                = sum(Unity.Details.Ud)/Sumary.Zones;
     511       
     512"Average Overall Heat Transfer Coefficient Clean"
     513        Sumary.Ucaverage                = sum(Unity.Details.Uc)/Sumary.Zones;
    510514
    511515"Area Total"
     
    822826        Unity(i).Resistances.Rshell*(Unity(i).Shell.HeatTransfer.hshell)        =       1;
    823827       
    824 "Overall Heat Transfer Coefficient"
    825         Unity(i).Details.U*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;
     828"Overall Heat Transfer Coefficient Dirty"
     829        Unity(i).Details.Ud*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;
     830
     831"Overall Heat Transfer Coefficient Clean"
     832        (1/Unity(i).Details.Ud)=(1/Unity(i).Details.Uc)+Unity(i).Resistances.Rfo+Unity(i).Resistances.Rfi*(Dotube/Ditube);
    826833
    827834"Exchange Surface Area"
     
    831838        Ltube = Unity(i).Baffles.Lsi+Unity(i).Baffles.Lso+Unity(i).Baffles.Ls*(Nb-1);
    832839       
    833 
    834840"Js Factor"     
    835841        Unity(i).Shell.HeatTransfer.Js                  =       1;
     
    853859        Unity(i).Shell.HeatTransfer.Jtotal              =       Unity(i).Shell.HeatTransfer.Jc*Unity(i).Shell.HeatTransfer.Jl*Unity(i).Shell.HeatTransfer.Jb*Unity(i).Shell.HeatTransfer.Jr*Unity(i).Shell.HeatTransfer.Js;
    854860end
    855 
    856861
    857862"Velocity Tube Side Inlet Nozzle"
     
    991996        Sumary.htubeaverage             = sum(Unity.Tubes.HeatTransfer.htube)/Sumary.Zones;
    992997
    993 "Average Overall Heat Transfer Coefficient"
    994         Sumary.Uaverage                 = sum(Unity.Details.U)/Sumary.Zones;
     998"Average Overall Heat Transfer Coefficient Dirty"
     999        Sumary.Udaverage                = sum(Unity.Details.Ud)/Sumary.Zones;
     1000       
     1001"Average Overall Heat Transfer Coefficient Clean"
     1002        Sumary.Ucaverage                = sum(Unity.Details.Uc)/Sumary.Zones;   
    9951003
    9961004"Area Total"
     
    10321040"Shell Side Cross Flow Area"
    10331041        Unity(Nb+1).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(Nb+1).Baffles.Lso);
    1034 
    1035 #"LMTD Correction Factor"
    1036 #       Unity.Fc = HE.EshellCorrectionFactor(Unity(1).Inlet.Hot.T,Unity(Nb+1).Outlet.Hot.T,Unity(Nb+1).Inlet.Cold.T,Unity(1).Outlet.Cold.T);
    10371042
    10381043if side equal 1
     
    13111316        Unity(i).Resistances.Rshell*(Unity(i).Shell.HeatTransfer.hshell)        =       1;
    13121317       
    1313 "Overall Heat Transfer Coefficient"
    1314         Unity(i).Details.U*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;
     1318"Overall Heat Transfer Coefficient Dirty"
     1319        Unity(i).Details.Ud*(Dotube/(Unity(i).Tubes.HeatTransfer.htube*Ditube)+(Dotube*ln(Dotube/Ditube)/(2*Kwall))+(1/(Unity(i).Shell.HeatTransfer.hshell)))=1;
     1320
     1321"Overall Heat Transfer Coefficient Clean"
     1322        (1/Unity(i).Details.Ud)=(1/Unity(i).Details.Uc)+Unity(i).Resistances.Rfo+Unity(i).Resistances.Rfi*(Dotube/Ditube);
    13151323
    13161324"Exchange Surface Area"
  • mso/eml/heat_exchangers/HeatExchangerSimplified.mso

    r78 r100  
    349349        PressureDrop.Cold.Pdrop  = Inlet.Cold.P*PressureDrop.Cold.FPdrop;
    350350       
    351        
    352351end
    353352
     
    407406
    408407"Exchange Surface Area"
    409         Details.Q = Details.U*Details.A*MTD;   
     408        Details.Q = Details.Ud*Details.A*MTD;   
    410409       
    411410"Mean Temperature Difference"   
  • mso/eml/heat_exchangers/Mheatex.mso

    r78 r100  
    8787SET
    8888
    89 Side = HE.FlowDir();
     89# Flow Direction
     90        Side = HE.FlowDir();
    9091
    91 Inlet.Ncold  = Ncold;
    92 Outlet.Ncold = Ncold;
     92# Inlet Ncold Parameters
     93        Inlet.Ncold  = Ncold;
     94       
     95# Outlet Ncold Parameters
     96        Outlet.Ncold = Ncold;
    9397
    94 Inlet.Nhot  = Nhot ;
    95 Outlet.Nhot = Nhot ;
    96 
     98# Inlet Nhot Parameters
     99        Inlet.Nhot  = Nhot;
     100       
     101# Outlet Nhot Parameters
     102        Outlet.Nhot = Nhot ;
    97103
    98104EQUATIONS
  • mso/sample/heat_exchangers/Eshell_Detailed_LMTD.mso

    r89 r100  
    5353        NComp                           = PP.NumberOfComponents;
    5454       
    55         exchanger.HE.HotSide            = "Shell";
     55        exchanger.HE.HotSide  = "Shell";
    5656
    5757#   LMTD Correction Factor
     
    8585exchanger.Ditube                = 0.013395  *"m";
    8686exchanger.Dotube                = 0.015875      *"m";
    87 exchanger.Kwall                         = 0.57          *"kW/m/K";
     87exchanger.Kwall                         = 0.057         *"kW/m/K";
    8888exchanger.Donozzle_Tube     = 0.203     *"m";
    8989exchanger.Dinozzle_Tube         = 0.203         *"m";
     
    107107#   Cold Stream
    108108#============================================   
    109         streamcold_in.F  = 75                   * "mol/s";
    110         streamcold_in.T  = 333   * "K";
     109        streamcold_in.F  = 75   * "mol/s";
     110        streamcold_in.T  = 333  * "K";
    111111        streamcold_in.z  = [1];
    112112        streamcold_in.P  = 2210*"kPa";
     
    116116        exchanger.Baffles.Ls    = 0.622  *"m";
    117117        exchanger.Baffles.Lsi   = 0.807  *"m";
     118#=====================================================================
     119#       Fouling
     120#=====================================================================
     121exchanger.Resistances.Rfi = 0*"m^2*K/kW";
     122exchanger.Resistances.Rfo = 0*"m^2*K/kW";
    118123
    119124OPTIONS
  • mso/sample/heat_exchangers/Eshell_Detailed_NTU.mso

    r89 r100  
    2525
    2626using "heat_exchangers/HeatExchangerDetailed.mso";
    27 #===============================================================
    28 #   erro ao trocar de modelo termodinamico
    29 #===============================================================
     27
    3028FlowSheet Exchanger_E_shell_Detailed_NTU
    3129       
     
    105103#       Baffle Spacing
    106104#=====================================================================
    107         exchanger.Baffles.Ls    = 0.622         *"m";
     105        exchanger.Baffles.Ls    = 0.622  *"m";
    108106        exchanger.Baffles.Lsi   = 0.807  *"m";
     107#=====================================================================
     108#       Fouling
     109#=====================================================================
     110exchanger.Resistances.Rfi = 0.0018*"m^2*K/kW";
     111exchanger.Resistances.Rfo = 0.0021*"m^2*K/kW";
    109112
    110113OPTIONS
  • mso/sample/heat_exchangers/Eshell_Discretized_LMTD.mso

    r89 r100  
    118118        exchanger.Unity.Baffles.Ls      = 0.622         *"m";
    119119        exchanger.Unity.Baffles.Lsi     = 0.807  *"m";
    120 
     120#=====================================================================
     121#       Fouling
     122#=====================================================================
     123exchanger.Unity.Resistances.Rfi = 0*"m^2*K/kW";
     124exchanger.Unity.Resistances.Rfo = 0*"m^2*K/kW";
    121125OPTIONS
    122126#============================================
  • mso/sample/heat_exchangers/Eshell_Discretized_NTU.mso

    r89 r100  
    114114        exchanger.Unity.Baffles.Ls      = 0.622         *"m";
    115115        exchanger.Unity.Baffles.Lsi     = 0.807  *"m";
     116#=====================================================================
     117#       Fouling
     118#=====================================================================
     119exchanger.Unity.Resistances.Rfi = 0*"m^2*K/kW";
     120exchanger.Unity.Resistances.Rfo = 0*"m^2*K/kW";
    116121
    117122OPTIONS
  • mso/sample/heat_exchangers/MheaterTeste.mso

    r89 r100  
    4949SET
    5050
    51         PP.LiquidModel  = "RK";
    52         PP.VapourModel  = "RK";
     51        PP.LiquidModel  = "PR";
     52        PP.VapourModel  = "PR";
    5353        PP.Components   = ["water","methanol","benzene"];
    5454        NComp                   = PP.NumberOfComponents;
    5555       
    5656       
    57         Mheater.HE.FlowDirection  = "Cocurrent";
     57        Mheater.HE.FlowDirection  = "Counter";
    5858        Mheater.Ncold = 2;
    5959        Mheater.Nhot  = 2;
  • mso/sample/heat_exchangers/Mheater_project.mso

    r89 r100  
    2929* $Id$
    3030*--------------------------------------------------------------------*#
    31 
    3231using "heat_exchangers/Mheatex";
    3332
     
    5251
    5352        PP                      as CalcObject   (File="vrpp");
    54 #       HE                      as CalcObject   (File="heatex");
    5553        NComp           as Integer              (Brief="Number Components");
    5654        Ncold           as Integer              (Brief="Number of Inlet Cold Streams");
     
    6563        Mheater.Nhot            = 2;
    6664    Mheater.Ncold               = 2;
    67         Mheater.HE.FlowDirection      = "Cocurrent";
    68 #       HE.FlowDirection      = "Counter";
    69 
     65        Mheater.HE.FlowDirection      = "Counter";
    7066       
    7167SPECIFY
  • mso/sample/heat_exchangers/Multipass_Detailed.mso

    r89 r100  
    115115        exchanger.Unity.Baffles.Lsi     = 0.807  *"m";
    116116
     117#=====================================================================
     118#       Fouling
     119#=====================================================================
     120exchanger.Unity.Resistances.Rfi = 0*"m^2*K/kW";
     121exchanger.Unity.Resistances.Rfo = 0*"m^2*K/kW";
    117122
    118123OPTIONS
  • mso/sample/heat_exchangers/Sample_DoublePipe.mso

    r89 r100  
    7373#               Hot Stream
    7474#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    75 
    7675        streamhot.F = 10                        * "kmol/h";
    7776        streamhot.T = (68+273.15)   * "K";
     
    8079       
    8180#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    82 #    Cold Stream
     81#       Cold Stream
    8382#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
    84 
    8583        streamcold.F = 5                        * "kmol/h";
    8684        streamcold.P = 5.1              * "bar";
    8785        streamcold.T = (23+273.15)      * "K";
    8886    streamcold.z = [0,1];
     87#=====================================================================
     88#       Fouling
     89#=====================================================================
     90exchanger.Resistances.Rfi = 0.0018*"m^2*K/kW";
     91exchanger.Resistances.Rfo = 0.0021*"m^2*K/kW";
    8992
    9093OPTIONS
    91 mode            = "steady";
     94
     95mode = "steady";
    9296
    9397end
  • mso/sample/heat_exchangers/Sample_Simplified.mso

    r89 r100  
    5454
    5555        exchanger.Details.A                             = 65.031*"m^2";
    56         exchanger.Details.U                             = 0.75*"W/(m^2*K)";
     56        exchanger.Details.Ud                            = 0.75*"W/(m^2*K)";
     57        exchanger.Details.Uc                            = 0.75*"W/(m^2*K)";
    5758        exchanger.PressureDrop.Hot.Pdrop        = 0*"atm";
    5859        exchanger.PressureDrop.Cold.Pdrop   = 0*"atm";
     
    105106
    106107        exchanger.Details.A                             = 65.031*"m^2";
    107         exchanger.Details.U                             = 0.75*"W/(m^2*K)";
     108        exchanger.Details.Ud                            = 0.75*"W/(m^2*K)";
     109        exchanger.Details.Uc                            = 0.75*"W/(m^2*K)";
    108110        exchanger.PressureDrop.Hot.Pdrop        = 0*"atm";
    109111        exchanger.PressureDrop.Cold.Pdrop   = 0*"atm";
  • mso/sample/heat_exchangers/sampleEshell.mso

    r89 r100  
    5151SPECIFY
    5252
    53 exchanger.Details.U                             = 210*"W/(m^2*K)";
     53exchanger.Details.Ud                            = 210*"W/(m^2*K)";
     54exchanger.Details.Uc                            = 210*"W/(m^2*K)";
    5455exchanger.PressureDrop.Hot.Pdrop        = 0.1*"kPa";
    5556exchanger.PressureDrop.Cold.Pdrop   = 0.2*"kPa";
  • mso/sample/heat_exchangers/sampleEshell_LMTD.mso

    r89 r100  
    5656SPECIFY
    5757
    58 exchanger.Details.U                             = 210*"W/(m^2*K)";
     58exchanger.Details.Ud                            = 210*"W/(m^2*K)";
     59exchanger.Details.Uc                            = 210*"W/(m^2*K)";
    5960exchanger.PressureDrop.Hot.Pdrop        = 0.1*"kPa";
    6061exchanger.PressureDrop.Cold.Pdrop   = 0.2*"kPa";
    6162
    6263streamhot.F             = 20    * "kmol/h";
    63 streamhot.T             = 450     * "K";
    64 streamhot.P             = 120    * "kPa";
     64streamhot.T             = 450   * "K";
     65streamhot.P             = 120   * "kPa";
    6566streamhot.z             = [1,0,0,0];
    6667
  • mso/sample/heat_exchangers/sampleLMTD.mso

    r89 r100  
    5656       
    5757"Overall Heat Transfer Coefficient"     
    58         exchanger.Details.U = 210*"W/(m^2*K)";
     58        exchanger.Details.Ud = 210*"W/(m^2*K)";
     59        exchanger.Details.Uc = 210*"W/(m^2*K)";
    5960
    6061exchanger.PressureDrop.Hot.Pdrop        = 0.1*"kPa";
  • mso/sample/heat_exchangers/sampleNTU.mso

    r89 r100  
    5454SPECIFY
    5555
    56 exchanger.Details.U                             = 210*"W/(m^2*K)";
     56exchanger.Details.Ud                            = 210*"W/(m^2*K)";
     57exchanger.Details.Uc                            = 210*"W/(m^2*K)";
    5758exchanger.PressureDrop.Hot.Pdrop        = 0.2*"kPa";
    5859exchanger.PressureDrop.Cold.Pdrop   = 0.2*"kPa";
  • mso/sample/heat_exchangers/samples1.mso

    r89 r100  
    7777#============================================
    7878        exchanger.Details.A                     = 65    *"m^2";
    79         exchanger.Details.U                 = 0.74  *"W/m^2/K";
     79        exchanger.Details.Ud                = 0.74  *"W/m^2/K";
     80        exchanger.Details.Uc                = 0.74  *"W/m^2/K";
    8081        exchanger.PressureDrop.Hot.Pdrop        = 0*"atm";
    8182        exchanger.PressureDrop.Cold.Pdrop   = 0*"atm";
  • mso/sample/heat_exchangers/samples2.mso

    r89 r100  
    7979#============================================
    8080        exchanger.Details.A                     = 65    *"m^2";
    81         exchanger.Details.U                 = 0.74  *"W/m^2/K";
     81        exchanger.Details.Ud                = 0.74  *"W/m^2/K";
     82        exchanger.Details.Uc                = 0.74  *"W/m^2/K";
    8283        exchanger.PressureDrop.Hot.Pdrop        = 0   * "atm";
    8384        exchanger.PressureDrop.Cold.Pdrop       = 0   * "atm";
Note: See TracChangeset for help on using the changeset viewer.