Changeset 45

Timestamp:

Nov 7, 2006, 2:11:10 PM (17 years ago)

Author:

gerson bicca

Message:

updated heat exchanger model

Location:

mso

Files:

• eml/heat_exchangers/HEX_Engine.mso (modified) (1 diff)
• eml/heat_exchangers/HeatExchangerDetailed.mso (modified) (36 diffs)
• eml/heat_exchangers/HeatExchangerDiscretized.mso (modified) (16 diffs)
• eml/heat_exchangers/HeatExchangerSimplified.mso (modified) (13 diffs)
• eml/heat_exchangers/Mheatex.mso (modified) (3 diffs)
• eml/heat_exchangers/heater.mso (modified) (4 diffs)
• eml/heat_exchangers/heatex.dll (modified) (previous)
• sample/heat_exchangers/Eshell_Detailed_LMTD.mso (modified) (3 diffs)
• sample/heat_exchangers/Eshell_Detailed_NTU.mso (modified) (3 diffs)
• sample/heat_exchangers/Eshell_Discretized_LMTD.mso (modified) (5 diffs)
• sample/heat_exchangers/Eshell_Discretized_NTU.mso (modified) (3 diffs)
• sample/heat_exchangers/MheaterTeste.mso (modified) (3 diffs)
• sample/heat_exchangers/Mheater_project.mso (modified) (3 diffs)
• sample/heat_exchangers/Multipass_Detailed.mso (modified) (3 diffs)
• sample/heat_exchangers/NTU_Method.mso (modified) (1 diff)
• sample/heat_exchangers/Sample_Simplified.mso (modified) (5 diffs)
• sample/heat_exchangers/heater_project.mso (modified) (6 diffs)
• sample/heat_exchangers/sampleEshell.mso (modified) (2 diffs)
• sample/heat_exchangers/sampleEshell_LMTD.mso (modified) (3 diffs)
• sample/heat_exchangers/sampleLMTD.mso (modified) (2 diffs)
• sample/heat_exchangers/sampleNTU.mso (modified) (3 diffs)
• sample/heat_exchangers/sample_cooler.mso (modified) (4 diffs)
• sample/heat_exchangers/sample_heater.mso (modified) (7 diffs)
• sample/heat_exchangers/samples1.mso (modified) (3 diffs)
• sample/heat_exchangers/samples2.mso (modified) (2 diffs)

mso/eml/heat_exchangers/HEX_Engine.mso

@@ -147 +147 @@
 #=====================================================================  
 VARIABLES
-Re                      as positive                     (Brief="Tube Side Reynolds Number",Default=1000,Lower=1);
-htube           as heat_trans_coeff     (Brief="Tube Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
-PR                      as positive                     (Brief="Tube Side Prandtl Number",Default=0.5,Lower=1e-8);
-PRw             as positive                     (Brief="Tube Side Prandtl Number at Wall Temperature",Default=0.5,Lower=1e-8);
-Phi             as positive                     (Brief="Phi Correction",Default=1,Lower=1e-3);
-Vtube           as velocity                     (Brief="Tube Side Velocity",Lower=1e-8);
+Re        as positive             (Brief="Tube Side Reynolds Number",Default=1000,Lower=1);
+htube as heat_trans_coeff (Brief="Tube Side Film Coefficient",Default=1,Lower=1e-12, Upper=1e6);
+PR        as positive             (Brief="Tube Side Prandtl Number",Default=0.5,Lower=1e-8);
+PRw   as positive                 (Brief="Tube Side Prandtl Number at Wall Temperature",Default=0.5,Lower=1e-8);
+Phi   as positive                 (Brief="Phi Correction",Default=1,Lower=1e-3);
+Vtube as velocity                 (Brief="Tube Side Velocity",Lower=1e-8);
 end
 

mso/eml/heat_exchangers/HeatExchangerDetailed.mso

@@ -8 +8 @@
 
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
-#       Basic Model for Detailed Heat Exchangers
+#       Basic Model for Detailed Shell and Tubes Heat Exchangers
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
 
@@ -16 +16 @@
 
 ext PP          as CalcObject   (Brief="External Physical Properties");
-ext HE          as CalcObject   (Brief="STHE Calculations",     File="heatex");
 ext NComp       as Integer      (Brief="Number of Components");
-  M(NComp)      as molweight    (Brief="Component Mol Weight");
+        HE              as CalcObject   (Brief="STHE Calculations",File="heatex");
+        M(NComp)        as molweight    (Brief="Component Mol Weight");
         
 VARIABLES
+
 in  Inlet               as Inlet_Main_Stream;   
 out Outlet              as Outlet_Main_Stream;  
@@ -39 +40 @@
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
 
-"HotStream Stream Average Temperature"
+"Hot Stream Average Temperature"
         Properties.Hot.Average.T = 0.5*Inlet.Hot.T + 0.5*Outlet.Hot.T;
         
@@ -45 +46 @@
         Properties.Cold.Average.T = 0.5*Inlet.Cold.T + 0.5*Outlet.Cold.T;
         
-"HotStream Stream Average Pressure"
+"Hot Stream Average Pressure"
         Properties.Hot.Average.P = 0.5*Inlet.Hot.P+0.5*Outlet.Hot.P;
         
@@ -62 +63 @@
         then    
         
-"Heat Capacity Cold Stream"
+"Cold Stream Average Heat Capacity"
         Properties.Cold.Average.Cp      =       PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+
+"Cold Stream Inlet Heat Capacity"
         Properties.Cold.Inlet.Cp        =       PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+
+"Cold Stream Outlet Heat Capacity"
         Properties.Cold.Outlet.Cp       =       PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
-"Mass Density Cold Stream"
+"Cold Stream Average Mass Density"
         Properties.Cold.Average.rho =   PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+
+"Cold Stream Inlet Mass Density"
         Properties.Cold.Inlet.rho       =       PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+
+"Cold Stream Outlet Mass Density"
         Properties.Cold.Outlet.rho      =       PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
-"Viscosity Cold Stream"
+"Cold Stream Average Viscosity"
         Properties.Cold.Average.Mu      =       PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+
+"Cold Stream inlet Viscosity"
         Properties.Cold.Inlet.Mu        =       PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+        
+"Cold Stream Outlet Viscosity"
         Properties.Cold.Outlet.Mu       =       PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
-"Conductivity Cold Stream"
+"Cold Stream Average Conductivity"
         Properties.Cold.Average.K       =       PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+
+"Cold Stream Inlet Conductivity"        
         Properties.Cold.Inlet.K         =       PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+
+"Cold Stream Outlet Conductivity"
         Properties.Cold.Outlet.K        =       PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
-"Heat Capacity Cold Stream"
+"Cold Stream Heat Capacity at Wall Temperature"
         Properties.Cold.Wall.Cp         =       PP.LiquidCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
         
-"Viscosity Cold Stream"
+"Cold Stream Viscosity at Wall Temperature"
         Properties.Cold.Wall.Mu         =       PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
 
-"Conductivity Cold Stream"
+"Cold Stream Conductivity at Wall Temperature"
         Properties.Cold.Wall.K          =       PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
 
@@ -94 +111 @@
         else
 
-"Heat Capacity Cold Stream"
+"Cold Stream Average Heat Capacity"
         Properties.Cold.Average.Cp      =       PP.VapourCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+
+"Cold Stream Inlet Heat Capacity"       
         Properties.Cold.Inlet.Cp        =       PP.VapourCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+
+"Cold Stream Outlet Heat Capacity"      
         Properties.Cold.Outlet.Cp       =       PP.VapourCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
-"Mass Density Cold Stream"
+"Cold Stream Average Mass Density"
         Properties.Cold.Average.rho =   PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+
+"Cold Stream Inlet Mass Density"
         Properties.Cold.Inlet.rho       =       PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+
+"Cold Stream Outlet Mass Density"       
         Properties.Cold.Outlet.rho      =       PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
-"Viscosity Cold Stream"
+"Cold Stream Average Viscosity "
         Properties.Cold.Average.Mu      =       PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+
+"Cold Stream Inlet Viscosity "  
         Properties.Cold.Inlet.Mu        =       PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+
+"Cold Stream Outlet Viscosity " 
         Properties.Cold.Outlet.Mu       =       PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
-"Conductivity Cold Stream"
+"Cold Stream Average Conductivity "
         Properties.Cold.Average.K       =       PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+
+"Cold Stream Inlet Conductivity "
         Properties.Cold.Inlet.K         =       PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+
+"Cold Stream Outlet Conductivity "
         Properties.Cold.Outlet.K        =       PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
         
-"Heat Capacity Cold Stream"
+"Cold Stream Heat Capacity at Wall Temperature"
         Properties.Cold.Wall.Cp         =       PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
 
 
-"Viscosity Cold Stream"
+"Cold Stream Viscosity at Wall Temperature"
         Properties.Cold.Wall.Mu         =       PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
 
-"Conductivity Cold Stream"
+"Cold Stream Conductivity at Wall Temperature"
         Properties.Cold.Wall.K          =       PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
         
@@ -271 +304 @@
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
 ext PP          as CalcObject   (Brief="External Physical Properties");
-ext     HE              as CalcObject   (Brief="STHE Calculations",File="heatex");
 side        as Integer          (Brief="Fluid Alocation Flag",Lower=0,Upper=1);
 Pi                      as constant     (Brief="Pi Number",Default=3.14159265);
-
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
 #                               Shell Geometrical Parameters                                            #
@@ -417 +448 @@
         then
         
-"Wall Temperature"
-#       Properties.Hot.Wall.Twall = HE.WallTemperature(Properties.Hot.Average.T,Properties.Cold.Average.T,Tubes.HeatTransfer.htube,Shell.HeatTransfer.hshell);
-
-#       Properties.Cold.Wall.Twall = HE.WallTemperature(Properties.Cold.Average.T,Properties.Hot.Average.T,Tubes.HeatTransfer.htube,Shell.HeatTransfer.hshell);
-
-Properties.Hot.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
-
-Properties.Cold.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
+"Hot Wall Temperature"
+        Properties.Hot.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
+
+"ColdWall Temperature"
+        Properties.Cold.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
 
 "Tube Side Velocity"
-        Tubes.HeatTransfer.Vtube  = HE.TubeVelocity(Properties.Cold.Inlet.Fw,Properties.Cold.Average.rho);
+        Tubes.HeatTransfer.Vtube  = Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Properties.Cold.Average.rho*Ntt);
 
 "Tube Side Reynolds Number"
-        Tubes.HeatTransfer.Re = HE.TubeReynoldsNumber(Properties.Cold.Average.rho,Tubes.HeatTransfer.Vtube,Properties.Cold.Average.Mu);
+        Tubes.HeatTransfer.Re     =     (Properties.Cold.Average.rho*Tubes.HeatTransfer.Vtube*Ditube)/Properties.Cold.Average.Mu;
         
 "Tube Side Prandtl Number"
-        Tubes.HeatTransfer.PR = HE.PrandtlNumber(Properties.Cold.Average.K,Properties.Cold.Average.Cp,Properties.Cold.Average.Mw,Properties.Cold.Average.Mu);
+        Tubes.HeatTransfer.PR = ((Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Average.Mu)/Properties.Cold.Average.K;
 
 "Tube Side Prandtl Number at Wall Temperature"
-        Tubes.HeatTransfer.PRw = HE.PrandtlNumber(Properties.Cold.Wall.K,Properties.Cold.Wall.Cp,Properties.Cold.Average.Mw,Properties.Cold.Wall.Mu);
+        Tubes.HeatTransfer.PRw = ((Properties.Cold.Wall.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Wall.Mu)/Properties.Cold.Wall.K;
 
 "Tube Side Film Coefficient"
-        Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Cold.Average.K,Tubes.HeatTransfer.Phi);
-        
+        Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Cold.Average.K)*Tubes.HeatTransfer.Phi;
+
 "Shell Side Prandtl Number"
-        Shell.HeatTransfer.PR=HE.PrandtlNumber(Properties.Hot.Average.K,Properties.Hot.Average.Cp,Properties.Hot.Average.Mw,Properties.Hot.Average.Mu); 
+        Shell.HeatTransfer.PR = ((Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Average.Mu)/Properties.Hot.Average.K;
 
 "Shell Side Prandtl Number at Wall Temperature"
-        Shell.HeatTransfer.PRw = HE.PrandtlNumber(Properties.Hot.Wall.K,Properties.Hot.Wall.Cp,Properties.Hot.Average.Mw,Properties.Hot.Wall.Mu);       
+        Shell.HeatTransfer.PRw = ((Properties.Hot.Wall.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Wall.Mu)/Properties.Hot.Wall.K;
 
 "Tube Side Pressure Drop"
@@ -483 +511 @@
         else
 
-"Wall Temperature"
-#       Properties.Cold.Wall.Twall = HE.WallTemperature(Properties.Cold.Average.T,Properties.Hot.Average.T,Tubes.HeatTransfer.htube,Shell.HeatTransfer.hshell);
-
-#       Properties.Hot.Wall.Twall = HE.WallTemperature(Properties.Hot.Average.T,Properties.Cold.Average.T,Tubes.HeatTransfer.htube,Shell.HeatTransfer.hshell);
-
-
-Properties.Hot.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
-
-Properties.Cold.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
+"Hot Wall Temperature"
+        Properties.Hot.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
+
+"Cold Wall Temperature"
+        Properties.Cold.Wall.Twall  = (Properties.Hot.Average.T+Properties.Cold.Average.T)/2;
 
 "Tube Side Velocity"
-        Tubes.HeatTransfer.Vtube=HE.TubeVelocity(Properties.Hot.Inlet.Fw,Properties.Hot.Average.rho);
+        Tubes.HeatTransfer.Vtube  = Properties.Hot.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Properties.Hot.Average.rho*Ntt);
 
 "Tube Side Reynolds Number"
-        Tubes.HeatTransfer.Re=HE.TubeReynoldsNumber(Properties.Hot.Average.rho,Tubes.HeatTransfer.Vtube,Properties.Hot.Average.Mu);
-        
+        Tubes.HeatTransfer.Re=(Properties.Hot.Average.rho*Tubes.HeatTransfer.Vtube*Ditube)/Properties.Hot.Average.Mu;
+
 "Tube Side Prandtl Number"
-        Tubes.HeatTransfer.PR=HE.PrandtlNumber(Properties.Hot.Average.K,Properties.Hot.Average.Cp,Properties.Hot.Average.Mw,Properties.Hot.Average.Mu);
+        Tubes.HeatTransfer.PR = ((Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Average.Mu)/Properties.Hot.Average.K;
 
 "Tube Side Prandtl Number at Wall"
-        Tubes.HeatTransfer.PRw=HE.PrandtlNumber(Properties.Hot.Wall.K,Properties.Hot.Wall.Cp,Properties.Hot.Average.Mw,Properties.Hot.Wall.Mu);
+        Tubes.HeatTransfer.PRw = ((Properties.Hot.Wall.Cp/Properties.Hot.Average.Mw)*Properties.Hot.Wall.Mu)/Properties.Hot.Wall.K;
 
 "Tube Side Film Coefficient"
-        Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Hot.Average.K,Tubes.HeatTransfer.Phi);
+        Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Tubes.HeatTransfer.Re,Tubes.HeatTransfer.PR,Properties.Hot.Average.K)*Tubes.HeatTransfer.Phi;
         
 "Shell Side Prandtl Number"
-        Shell.HeatTransfer.PR=HE.PrandtlNumber(Properties.Cold.Average.K,Properties.Cold.Average.Cp,Properties.Cold.Average.Mw,Properties.Cold.Average.Mu);
+        Shell.HeatTransfer.PR = ((Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Average.Mu)/Properties.Cold.Average.K;
 
 "Shell Side Prandtl Number at Wall"
-        Shell.HeatTransfer.PRw=HE.PrandtlNumber(Properties.Cold.Wall.K,Properties.Cold.Wall.Cp,Properties.Cold.Average.Mw,Properties.Cold.Wall.Mu);
+        Shell.HeatTransfer.PRw=((Properties.Cold.Wall.Cp/Properties.Cold.Average.Mw)*Properties.Cold.Wall.Mu)/Properties.Cold.Wall.K;
 
 "Tube Side Pressure Drop"
@@ -546 +570 @@
 "Pressure Drop Cold Stream"
         Outlet.Cold.P  = Inlet.Cold.P - Shell.PressureDrop.Pdtotal;
-
 
 end
@@ -598 +621 @@
         
 VARIABLES
-LMTD            as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=5);
-Fc                      as positive             (Brief="LMTD Correction Factor",Lower=0.4);
-MTD                     as temp_delta   (Brief="Mean Temperature Difference",Lower=5);
+DT0     as temp_delta   (Brief="Temperature Difference at Inlet",Lower=1);
+DTL             as temp_delta   (Brief="Temperature Difference at Outlet",Lower=1);
+LMTD    as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=5);
+Fc              as positive             (Brief="LMTD Correction Factor",Lower=0.4);
+MTD             as temp_delta   (Brief="Mean Temperature Difference",Lower=5);
 
 EQUATIONS
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+#                       Log Mean Temperature Difference
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+if abs(DT0 - DTL) > 0.05*max(abs([DT0,DTL])) 
+        
+        then
+"Log Mean Temperature Difference"
+        LMTD= (DT0-DTL)/ln(DT0/DTL);
+
+        else
+        
+if DT0*DTL equal 0
+        
+        then
+"Log Mean Temperature Difference"
+        LMTD = 0.5*(DT0+DTL);
+        
+        else
+"Log Mean Temperature Difference"
+        LMTD = 0.5*(DT0+DTL)*(1-(DT0-DTL)^2/(DT0*DTL)*(1+(DT0-DTL)^2/(DT0*DTL)/2)/12);
+        
+end
+        
+end
 
 "Exchange Surface Area"
@@ -608 +657 @@
         
 "Mean Temperature Difference"   
-        MTD   = Fc*LMTD;
+        MTD = Fc*LMTD;
         
 end
@@ -630 +679 @@
         
 "Shell Side Reynolds Number"
-        Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Hot.Average.Mu);
+        Shell.HeatTransfer.Re=(Dotube*Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Hot.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Shell.HeatTransfer.Ji,Properties.Hot.Average.Cp,Properties.Hot.Average.Mw,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Shell.HeatTransfer.PR,Shell.HeatTransfer.Jtotal,Shell.HeatTransfer.Phi);
+        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*(Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -648 +697 @@
 
 "Shell Side Reynolds Number"
-        Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Cold.Average.Mu);
+        Shell.HeatTransfer.Re=(Dotube*Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Cold.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Shell.HeatTransfer.Ji,Properties.Cold.Average.Cp,Properties.Cold.Average.Mw,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Shell.HeatTransfer.PR,Shell.HeatTransfer.Jtotal,Shell.HeatTransfer.Phi);
+        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*(Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -691 +740 @@
         
 "Shell Side Reynolds Number"
-        Shell.HeatTransfer.Re= HE.ShellReynoldsNumber(Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Hot.Average.Mu);
+        Shell.HeatTransfer.Re=(Dotube*Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Hot.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Shell.HeatTransfer.hshell= HE.ShellFilmCoeff(Shell.HeatTransfer.Ji,Properties.Hot.Average.Cp,Properties.Hot.Average.Mw,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Shell.HeatTransfer.PR,Shell.HeatTransfer.Jtotal,Shell.HeatTransfer.Phi);
+        Shell.HeatTransfer.hshell= Shell.HeatTransfer.Ji*(Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*(Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -709 +758 @@
 
 "Shell Side Reynolds Number"
-        Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Cold.Average.Mu);
+        Shell.HeatTransfer.Re=(Dotube*Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Cold.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Shell.HeatTransfer.Ji,Properties.Cold.Average.Cp,Properties.Cold.Average.Mw,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Shell.HeatTransfer.PR,Shell.HeatTransfer.Jtotal,Shell.HeatTransfer.Phi);
+        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*(Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -735 +784 @@
 PARAMETERS
 
-Nshell          as Integer              (Brief="N Shell in Series",Default=2);
-ext     HE              as CalcObject   (Brief="STHE Calculations",File="heatex.dll");
-ext PP          as CalcObject   (Brief="External Physical Properties");
-side        as Integer          (Brief="Fluid Alocation",Lower=0,Upper=1);
-Pi                      as constant     (Brief="Pi Number",Default=3.14159265);
+        Nshell          as Integer              (Brief="N Shell in Series",Default=2);
+ext PP                  as CalcObject   (Brief="External Physical Properties");
+        HE                      as CalcObject   (Brief="heatex Calculations", File="heatex");
+        side        as Integer          (Brief="Fluid Alocation",Lower=0,Upper=1);
+        Pi                      as constant     (Brief="Pi Number",Default=3.14159265);
 
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
@@ -780 +829 @@
 
 Unity([1:Nshell-1]).Outlet.Hot  to Unity([2:Nshell]).Inlet.Hot;
-Unity([2:Nshell]).Outlet.Cold to Unity([1:Nshell-1]).Inlet.Cold;
+Unity([2:Nshell]).Outlet.Cold   to Unity([1:Nshell-1]).Inlet.Cold;
 
 EQUATIONS
@@ -790 +839 @@
 
 "Ji Factor"
-        Unity(i).Shell.HeatTransfer.Ji = HE.JiFactor(Unity(i).Shell.HeatTransfer.Re);
+        Unity(i).Shell.HeatTransfer.Ji =HE.JiFactor(Unity(i).Shell.HeatTransfer.Re);
 
 "Jc Factor"
@@ -841 +890 @@
 
 "Tube Side Velocity"
-        Unity(i).Tubes.HeatTransfer.Vtube  = HE.TubeVelocity(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Properties.Cold.Average.rho);
+        Unity(i).Tubes.HeatTransfer.Vtube  = Unity(i).Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Cold.Average.rho*Ntt);
 
 "Tube Side Reynolds Number"
-        Unity(i).Tubes.HeatTransfer.Re = HE.TubeReynoldsNumber(Unity(i).Properties.Cold.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Properties.Cold.Average.Mu);
-        
+        Unity(i).Tubes.HeatTransfer.Re=(Unity(i).Properties.Cold.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Cold.Average.Mu;
+
 "Tube Side Prandtl Number"
-        Unity(i).Tubes.HeatTransfer.PR = HE.PrandtlNumber(Unity(i).Properties.Cold.Average.K,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.PR = ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
 
 "Tube Side Prandtl Number at Wall Temperature"
-        Unity(i).Tubes.HeatTransfer.PRw = HE.PrandtlNumber(Unity(i).Properties.Cold.Wall.K,Unity(i).Properties.Cold.Wall.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Wall.Mu);
+        Unity(i).Tubes.HeatTransfer.PRw = ((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
 
 "Tube Side Film Coefficient"
-        Unity(i).Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K,Unity(i).Tubes.HeatTransfer.Phi);
+        Unity(i).Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K)*Unity(i).Tubes.HeatTransfer.Phi;
         
 "Shell Side Prandtl Number"
-        Unity(i).Shell.HeatTransfer.PR=HE.PrandtlNumber(Unity(i).Properties.Hot.Average.K,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Average.Mu);    
+        Unity(i).Shell.HeatTransfer.PR=((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;  
 
 "Shell Side Prandtl Number at Wall Temperature"
-        Unity(i).Shell.HeatTransfer.PRw = HE.PrandtlNumber(Unity(i).Properties.Hot.Wall.K,Unity(i).Properties.Hot.Wall.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Wall.Mu);  
+        Unity(i).Shell.HeatTransfer.PRw = ((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;        
 
 "Tube Side Pressure Drop"
@@ -904 +953 @@
 
 "Tube Side Velocity"
-        Unity(i).Tubes.HeatTransfer.Vtube       =       HE.TubeVelocity(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Properties.Hot.Average.rho);
+        Unity(i).Tubes.HeatTransfer.Vtube  = Unity(i).Properties.Hot.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Hot.Average.rho*Ntt);
 
 "Tube Side Reynolds Number"
-        Unity(i).Tubes.HeatTransfer.Re          =       HE.TubeReynoldsNumber(Unity(i).Properties.Hot.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.Re          =       (Unity(i).Properties.Hot.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Hot.Average.Mu;
         
 "Tube Side Prandtl Number"
-        Unity(i).Tubes.HeatTransfer.PR          =       HE.PrandtlNumber(Unity(i).Properties.Hot.Average.K,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.PR          =       ((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;
         
 "Tube Side Prandtl Number at Wall"
-        Unity(i).Tubes.HeatTransfer.PRw     =  HE.PrandtlNumber(Unity(i).Properties.Hot.Wall.K,Unity(i).Properties.Hot.Wall.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Wall.Mu);
+        Unity(i).Tubes.HeatTransfer.PRw     =  ((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;
 
 "Tube Side Film Coefficient"
-        Unity(i).Tubes.HeatTransfer.htube       =       HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K,Unity(i).Tubes.HeatTransfer.Phi);
+        Unity(i).Tubes.HeatTransfer.htube       =       HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K)*Unity(i).Tubes.HeatTransfer.Phi;
         
 "Shell Side Prandtl Number"
-        Unity(i).Shell.HeatTransfer.PR          =       HE.PrandtlNumber(Unity(i).Properties.Cold.Average.K,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Shell.HeatTransfer.PR          =       ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
         
 "Shell Side Prandtl Number at Wall"
-        Unity(i).Shell.HeatTransfer.PRw     =   HE.PrandtlNumber(Unity(i).Properties.Cold.Wall.K,Unity(i).Properties.Cold.Wall.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Wall.Mu);
+        Unity(i).Shell.HeatTransfer.PRw     =   ((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
 
 "Tube Side Pressure Drop"
@@ -1009 +1058 @@
         
 "Shell Side Reynolds Number"
-        Unity(i).Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Hot.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Unity(i).Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Unity(i).Shell.HeatTransfer.Ji,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Shell.HeatTransfer.PR,Unity(i).Shell.HeatTransfer.Jtotal,Unity(i).Shell.HeatTransfer.Phi);
+        Unity(i).Shell.HeatTransfer.hshell=Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*(Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -1027 +1076 @@
 
 "Shell Side Reynolds Number"
-        Unity(i).Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Cold.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Unity(i).Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Unity(i).Shell.HeatTransfer.Ji,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Shell.HeatTransfer.PR,Unity(i).Shell.HeatTransfer.Jtotal,Unity(i).Shell.HeatTransfer.Phi);
+        Unity(i).Shell.HeatTransfer.hshell=Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*(Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -1087 +1136 @@
         Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
 
-"Counter Flow LMTD"
-        LMTD = HE.CounterLMTD(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
-
+"Temperature Difference at Inlet"
+        DT0 = Inlet.Hot.T - Outlet.Cold.T;
+
+"Temperature Difference at Outlet"
+        DTL = Outlet.Hot.T - Inlet.Cold.T;
+        
 "Js Factor"     
         Shell.HeatTransfer.Js = HE.JsFactor(Shell.HeatTransfer.Re,Baffles.Lsi,Baffles.Lso,Baffles.Ls);
@@ -1098 +1150 @@
         
 "Shell Side Reynolds Number"
-        Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Hot.Average.Mu);
+        Shell.HeatTransfer.Re=(Dotube*Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Hot.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Shell.HeatTransfer.Ji,Properties.Hot.Average.Cp,Properties.Hot.Average.Mw,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Shell.HeatTransfer.PR,Shell.HeatTransfer.Jtotal,Shell.HeatTransfer.Phi);
+        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*(Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -1116 +1168 @@
 
 "Shell Side Reynolds Number"
-        Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Cold.Average.Mu);
+        Shell.HeatTransfer.Re=(Dotube*Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Cold.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Shell.HeatTransfer.Ji,Properties.Cold.Average.Cp,Properties.Cold.Average.Mw,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Shell.HeatTransfer.PR,Shell.HeatTransfer.Jtotal,Shell.HeatTransfer.Phi);
+        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*(Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -1147 +1199 @@
         Fc = HE.FshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
 
-"Counter Flow LMTD"
-        LMTD = HE.CounterLMTD(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
-
+"Temperature Difference at Inlet"
+        DT0 = Inlet.Hot.T - Outlet.Cold.T;
+
+"Temperature Difference at Outlet"
+        DTL = Outlet.Hot.T - Inlet.Cold.T;
+        
 "Js Factor"     
         Shell.HeatTransfer.Js = HE.JsFactor(Shell.HeatTransfer.Re,Baffles.Lsi,Baffles.Lso,Baffles.Ls);
@@ -1158 +1213 @@
         
 "Shell Side Reynolds Number"
-        Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Hot.Average.Mu);
+        Shell.HeatTransfer.Re=(Dotube*Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Hot.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Shell.HeatTransfer.Ji,Properties.Hot.Average.Cp,Properties.Hot.Average.Mw,Properties.Hot.Inlet.Fw,Shell.HeatTransfer.Sm,Shell.HeatTransfer.PR,Shell.HeatTransfer.Jtotal,Shell.HeatTransfer.Phi);
+        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Hot.Average.Cp/Properties.Hot.Average.Mw)*(Properties.Hot.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -1176 +1231 @@
 
 "Shell Side Reynolds Number"
-        Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Properties.Cold.Average.Mu);
+        Shell.HeatTransfer.Re=(Dotube*Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)/Properties.Cold.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Shell.HeatTransfer.Ji,Properties.Cold.Average.Cp,Properties.Cold.Average.Mw,Properties.Cold.Inlet.Fw,Shell.HeatTransfer.Sm,Shell.HeatTransfer.PR,Shell.HeatTransfer.Jtotal,Shell.HeatTransfer.Phi);
+        Shell.HeatTransfer.hshell=Shell.HeatTransfer.Ji*(Properties.Cold.Average.Cp/Properties.Cold.Average.Mw)*(Properties.Cold.Inlet.Fw/Shell.HeatTransfer.Sm)*(Shell.HeatTransfer.PR^(-2/3))*Shell.HeatTransfer.Jtotal*Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -1193 +1248 @@
 end
 
-
-
 end
 
@@ -1204 +1257 @@
 
 Nshell          as Integer              (Brief="N Shell in Series",Default=2);
-ext     HE              as CalcObject   (Brief="STHE Calculations",File="heatex.dll");
+        HE              as CalcObject   (Brief="STHE Calculations",File="heatex");
 ext PP          as CalcObject   (Brief="External Physical Properties");
 side        as Integer          (Brief="Fluid Alocation",Lower=0,Upper=1);
@@ -1313 +1366 @@
 
 "Tube Side Velocity"
-        Unity(i).Tubes.HeatTransfer.Vtube  = HE.TubeVelocity(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Properties.Cold.Average.rho);
+        Unity(i).Tubes.HeatTransfer.Vtube  = Unity(i).Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Cold.Average.rho*Ntt);
 
 "Tube Side Reynolds Number"
-        Unity(i).Tubes.HeatTransfer.Re = HE.TubeReynoldsNumber(Unity(i).Properties.Cold.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Properties.Cold.Average.Mu);
-        
+        Unity(i).Tubes.HeatTransfer.Re=(Unity(i).Properties.Cold.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Cold.Average.Mu;
+
 "Tube Side Prandtl Number"
-        Unity(i).Tubes.HeatTransfer.PR = HE.PrandtlNumber(Unity(i).Properties.Cold.Average.K,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.PR = ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
 
 "Tube Side Prandtl Number at Wall"
-        Unity(i).Tubes.HeatTransfer.PRw = HE.PrandtlNumber(Unity(i).Properties.Cold.Wall.K,Unity(i).Properties.Cold.Wall.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Wall.Mu);
+        Unity(i).Tubes.HeatTransfer.PRw = ((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
 
 "Tube Side Film Coefficient"
-        Unity(i).Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K,Unity(i).Tubes.HeatTransfer.Phi);
+        Unity(i).Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K)*Unity(i).Tubes.HeatTransfer.Phi;
         
 "Shell Side Prandtl Number"
-        Unity(i).Shell.HeatTransfer.PR=HE.PrandtlNumber(Unity(i).Properties.Hot.Average.K,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Average.Mu);    
+        Unity(i).Shell.HeatTransfer.PR=((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;  
 
 "Shell Side Prandtl Number at Wall"
-        Unity(i).Shell.HeatTransfer.PRw=HE.PrandtlNumber(Unity(i).Properties.Hot.Wall.K,Unity(i).Properties.Hot.Wall.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Wall.Mu);    
+        Unity(i).Shell.HeatTransfer.PRw=((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;  
 
 "Tube Side Pressure Drop"
@@ -1376 +1429 @@
 
 "Tube Side Velocity"
-        Unity(i).Tubes.HeatTransfer.Vtube=HE.TubeVelocity(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Properties.Hot.Average.rho);
+        Unity(i).Tubes.HeatTransfer.Vtube  = Unity(i).Properties.Hot.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Hot.Average.rho*Ntt);
 
 "Tube Side Reynolds Number"
-        Unity(i).Tubes.HeatTransfer.Re=HE.TubeReynoldsNumber(Unity(i).Properties.Hot.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.Re=(Unity(i).Properties.Hot.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Hot.Average.Mu;
         
 "Tube Side Prandtl Number"
-        Unity(i).Tubes.HeatTransfer.PR=HE.PrandtlNumber(Unity(i).Properties.Hot.Average.K,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.PR=((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;
 
 "Tube Side Prandtl Number at Wall"
-        Unity(i).Tubes.HeatTransfer.PRw=HE.PrandtlNumber(Unity(i).Properties.Hot.Wall.K,Unity(i).Properties.Hot.Wall.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Wall.Mu);
+        Unity(i).Tubes.HeatTransfer.PRw=((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;
 
 "Tube Side Film Coefficient"
-        Unity(i).Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K,Unity(i).Tubes.HeatTransfer.Phi);
+        Unity(i).Tubes.HeatTransfer.htube= HE.TubeFilmCoeff(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K)*Unity(i).Tubes.HeatTransfer.Phi;
         
 "Shell Side Prandtl Number"
-        Unity(i).Shell.HeatTransfer.PR=HE.PrandtlNumber(Unity(i).Properties.Cold.Average.K,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Shell.HeatTransfer.PR=((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
 
 "Shell Side Prandtl Number at Wall"
-        Unity(i).Shell.HeatTransfer.PRw=HE.PrandtlNumber(Unity(i).Properties.Cold.Wall.K,Unity(i).Properties.Cold.Wall.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Wall.Mu);
+        Unity(i).Shell.HeatTransfer.PRw=((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
 
 "Tube Side Pressure Drop"
@@ -1432 +1485 @@
 end
 
+if side equal 1 
+        
+        then 
+
+"Shell Side inlet Nozzle rho-V^2"
+        Unity(i).Shell.PressureDrop.RVsquare_in = Unity(i).Properties.Hot.Inlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_in)^2;
+
+"Shell Side Outlet Nozzle rho-V^2"
+        Unity(i).Shell.PressureDrop.RVsquare_out = Unity(i).Properties.Hot.Outlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_out)^2;
+        
+        else
+        
+"Shell Side inlet Nozzle rho-V^2"
+        Unity(i).Shell.PressureDrop.RVsquare_in = Unity(i).Properties.Cold.Inlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_in)^2;
+
+"Shell Side Outlet Nozzle rho-V^2"
+        Unity(i).Shell.PressureDrop.RVsquare_out = Unity(i).Properties.Cold.Outlet.rho*(Unity(i).Shell.PressureDrop.Vnozzle_out)^2;
+
+end
+
 "Tube Resistance"       
         Unity(i).Resistances.Rtube*(Unity(i).Tubes.HeatTransfer.htube*Ditube) = Dotube;
@@ -1464 +1537 @@
         
 "Shell Side Reynolds Number"
-        Unity(i).Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Hot.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Unity(i).Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Unity(i).Shell.HeatTransfer.Ji,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Shell.HeatTransfer.PR,Unity(i).Shell.HeatTransfer.Jtotal,Unity(i).Shell.HeatTransfer.Phi);
+        Unity(i).Shell.HeatTransfer.hshell=Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*(Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"
@@ -1482 +1555 @@
 
 "Shell Side Reynolds Number"
-        Unity(i).Shell.HeatTransfer.Re=HE.ShellReynoldsNumber(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Cold.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Unity(i).Shell.HeatTransfer.hshell=HE.ShellFilmCoeff(Unity(i).Shell.HeatTransfer.Ji,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Shell.HeatTransfer.PR,Unity(i).Shell.HeatTransfer.Jtotal,Unity(i).Shell.HeatTransfer.Phi);
+        Unity(i).Shell.HeatTransfer.hshell=Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*(Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
 
 "Shell Pressure Drop Cross Flow"

mso/eml/heat_exchangers/HeatExchangerDiscretized.mso

@@ -16 +16 @@
 
 ext PP      as CalcObject(Brief="External Physical Properties");
-ext HE      as CalcObject(Brief="STHE Calculations",File="heatex");
+        HE          as CalcObject(Brief="STHE Calculations",File="heatex");
 ext NComp   as Integer   (Brief="Number of Components");
   M(NComp)  as molweight (Brief="Component Mol Weight");
@@ -279 +279 @@
         
 VARIABLES
-LMTD            as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=5);
-Fc                      as positive             (Brief="LMTD Correction Factor",Lower=0.5);
-MTD                     as temp_delta   (Brief="Mean Temperature Difference",Lower=1);
+DT0     as temp_delta   (Brief="Temperature Difference at Inlet",Lower=1);
+DTL             as temp_delta   (Brief="Temperature Difference at Outlet",Lower=1);
+LMTD    as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=5);
+Fc              as positive             (Brief="LMTD Correction Factor",Lower=0.5);
+MTD             as temp_delta   (Brief="Mean Temperature Difference",Lower=1);
 
 EQUATIONS
@@ -290 +292 @@
         MTD   = Fc*LMTD;
         
-"LMTD Correction Factor"
-        Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
-
-"LMTD"
-        LMTD = HE.CounterLMTD(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
-
+#"LMTD Correction Factor"
+#       Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
+
+"Temperature Difference at Inlet"
+        DT0 = Inlet.Hot.T - Outlet.Cold.T;
+
+"Temperature Difference at Outlet"
+        DTL = Outlet.Hot.T - Inlet.Cold.T;
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+#                       Log Mean Temperature Difference
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+if abs(DT0 - DTL) > 0.05*max(abs([DT0,DTL])) 
+        
+        then
+"Log Mean Temperature Difference"
+        LMTD= (DT0-DTL)/ln(DT0/DTL);
+
+        else
+        
+if DT0*DTL equal 0
+        
+        then
+"Log Mean Temperature Difference"
+        LMTD = 0.5*(DT0+DTL);
+        
+        else
+"Log Mean Temperature Difference"
+        LMTD = 0.5*(DT0+DTL)*(1-(DT0-DTL)^2/(DT0*DTL)*(1+(DT0-DTL)^2/(DT0*DTL)/2)/12);
+        
+end
+        
+end
+
+        
 end
 
@@ -332 +362 @@
 PARAMETERS
 
-ext     HE      as CalcObject   (Brief="STHE Calculations",File="heatex.dll");
+        HE      as CalcObject   (Brief="STHE Calculations",File="heatex");
 ext PP  as CalcObject   (Brief="External Physical Properties");
 side    as Integer              (Brief="Fluid Alocation",Lower=0,Upper=1);
@@ -601 +631 @@
         
 "Shell Side Reynolds Number"
-        Unity(i).Shell.HeatTransfer.Re                          =       HE.ShellReynoldsNumber(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Hot.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Unity(i).Shell.HeatTransfer.hshell                      =       HE.ShellFilmCoeff(Unity(i).Shell.HeatTransfer.Ji,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Shell.HeatTransfer.PR,Unity(i).Shell.HeatTransfer.Jtotal,Unity(i).Shell.HeatTransfer.Phi);
+        Unity(i).Shell.HeatTransfer.hshell =Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*(Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
 
 
@@ -617 +647 @@
 
 "Tube Side Velocity"
-        Unity(i).Tubes.HeatTransfer.Vtube       =       HE.TubeVelocity(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Properties.Cold.Average.rho);
+        Unity(i).Tubes.HeatTransfer.Vtube  = Unity(i).Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Cold.Average.rho*Ntt);
 
 "Tube Side Reynolds Number"
-        Unity(i).Tubes.HeatTransfer.Re          =       HE.TubeReynoldsNumber(Unity(i).Properties.Cold.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.Re          =       (Unity(i).Properties.Cold.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Cold.Average.Mu;
         
 "Tube Side Prandtl Number"
-        Unity(i).Tubes.HeatTransfer.PR          =       HE.PrandtlNumber(Unity(i).Properties.Cold.Average.K,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.PR          = ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
 
 "Tube Side Prandtl Number at Wall"
-        Unity(i).Tubes.HeatTransfer.PRw         =       HE.PrandtlNumber(Unity(i).Properties.Cold.Wall.K,Unity(i).Properties.Cold.Wall.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Wall.Mu);
+        Unity(i).Tubes.HeatTransfer.PRw         =       ((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
 
 "Tube Side Film Coefficient"
-        Unity(i).Tubes.HeatTransfer.htube       =       HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K,
-        Unity(i).Tubes.HeatTransfer.Phi,Sumary.Lz(i));
+        Unity(i).Tubes.HeatTransfer.htube       =       HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K,Sumary.Lz(i))*Unity(i).Tubes.HeatTransfer.Phi;
         
 "Shell Side Prandtl Number"
-        Unity(i).Shell.HeatTransfer.PR          =       HE.PrandtlNumber(Unity(i).Properties.Hot.Average.K,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Average.Mu);   
+        Unity(i).Shell.HeatTransfer.PR          =       ((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K; 
 
 "Shell Side Prandtl Number at Wall"
-        Unity(i).Shell.HeatTransfer.PRw         =       HE.PrandtlNumber(Unity(i).Properties.Hot.Wall.K,Unity(i).Properties.Hot.Wall.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Wall.Mu);    
+        Unity(i).Shell.HeatTransfer.PRw         =       ((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;  
 
 "Tube Side Pressure Drop"
@@ -657 +686 @@
         
 "Shell Side Reynolds Number"
-        Unity(i).Shell.HeatTransfer.Re                  =       HE.ShellReynoldsNumber(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Cold.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Unity(i).Shell.HeatTransfer.hshell              =       HE.ShellFilmCoeff(Unity(i).Shell.HeatTransfer.Ji,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,
-        Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Shell.HeatTransfer.PR,Unity(i).Shell.HeatTransfer.Jtotal,Unity(i).Shell.HeatTransfer.Phi);
+        Unity(i).Shell.HeatTransfer.hshell      =Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*
+        (Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
 
 
@@ -667 +696 @@
         Unity(i).Shell.PressureDrop.Pdwindow    =       HE.DeltaPwindowIncremental(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Cold.Average.rho,Unity(i).Properties.Cold.Average.Mu,Unity(i).Baffles.Ls);
 
-
 "Hot Wall Temperature"
         Unity(i).Properties.Hot.Wall.Twall      =       (Unity(i).Properties.Hot.Average.T+Unity(i).Properties.Cold.Average.T)/2;
@@ -675 +703 @@
 
 "Tube Side Velocity"
-        Unity(i).Tubes.HeatTransfer.Vtube               =       HE.TubeVelocity(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Properties.Hot.Average.rho);
+        Unity(i).Tubes.HeatTransfer.Vtube               = Unity(i).Properties.Hot.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Hot.Average.rho*Ntt);
 
 "Tube Side Reynolds Number"
-        Unity(i).Tubes.HeatTransfer.Re                  =       HE.TubeReynoldsNumber(Unity(i).Properties.Hot.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.Re                  =       (Unity(i).Properties.Hot.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Hot.Average.Mu;
         
 "Tube Side Prandtl Number"
-        Unity(i).Tubes.HeatTransfer.PR                  =       HE.PrandtlNumber(Unity(i).Properties.Hot.Average.K,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.PR                  =       ((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;
         
 "Tube Side Prandtl Number at Wall"
-        Unity(i).Tubes.HeatTransfer.PRw                 =       HE.PrandtlNumber(Unity(i).Properties.Hot.Wall.K,Unity(i).Properties.Hot.Wall.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Wall.Mu);
-        
+        Unity(i).Tubes.HeatTransfer.PRw                 =       ((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;
         
 "Tube Side Film Coefficient"
-        Unity(i).Tubes.HeatTransfer.htube               =       HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K,
-        Unity(i).Tubes.HeatTransfer.Phi,Sumary.Lz(i));
+        Unity(i).Tubes.HeatTransfer.htube               =       HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K,Sumary.Lz(i))*Unity(i).Tubes.HeatTransfer.Phi;
         
 "Shell Side Prandtl Number"
-        Unity(i).Shell.HeatTransfer.PR                  =       HE.PrandtlNumber(Unity(i).Properties.Cold.Average.K,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Shell.HeatTransfer.PR                  =       ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
 
 "Shell Side Prandtl Number at Wall"
-        Unity(i).Shell.HeatTransfer.PRw                 =       HE.PrandtlNumber(Unity(i).Properties.Cold.Wall.K,Unity(i).Properties.Cold.Wall.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Wall.Mu);
+        Unity(i).Shell.HeatTransfer.PRw                 =       ((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
 
 "Tube Side Pressure Drop"
@@ -821 +847 @@
 PARAMETERS
 
-ext     HE      as CalcObject   (Brief="STHE Calculations",File="heatex.dll");
+        HE      as CalcObject   (Brief="STHE Calculations",File="heatex");
 ext PP  as CalcObject   (Brief="External Physical Properties");
 side    as Integer              (Brief="Fluid Alocation",Lower=0,Upper=1);
@@ -929 +955 @@
         Unity(Nb+1).Shell.HeatTransfer.Sm = HE.CrossFlowArea(Unity(Nb+1).Baffles.Lso);
 
+"LMTD Correction Factor"
+        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);
 
 if side equal 1
@@ -1079 +1107 @@
 
 for i in [1:Nb+1]
-        
+
+
 if side equal 1 
         
@@ -1090 +1119 @@
         
 "Shell Side Reynolds Number"
-        Unity(i).Shell.HeatTransfer.Re                          =       HE.ShellReynoldsNumber(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Hot.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Unity(i).Shell.HeatTransfer.hshell                      =       HE.ShellFilmCoeff(Unity(i).Shell.HeatTransfer.Ji,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Shell.HeatTransfer.PR,Unity(i).Shell.HeatTransfer.Jtotal,Unity(i).Shell.HeatTransfer.Phi);
-
+        Unity(i).Shell.HeatTransfer.hshell =Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*(Unity(i).Properties.Hot.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
 
 "Shell Pressure Baffle Window"
@@ -1106 +1134 @@
 
 "Tube Side Velocity"
-        Unity(i).Tubes.HeatTransfer.Vtube       =       HE.TubeVelocity(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Properties.Cold.Average.rho);
+        Unity(i).Tubes.HeatTransfer.Vtube       =       Unity(i).Properties.Cold.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Cold.Average.rho*Ntt);
 
 "Tube Side Reynolds Number"
-        Unity(i).Tubes.HeatTransfer.Re          =       HE.TubeReynoldsNumber(Unity(i).Properties.Cold.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.Re          =       (Unity(i).Properties.Cold.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Cold.Average.Mu;
         
 "Tube Side Prandtl Number"
-        Unity(i).Tubes.HeatTransfer.PR          =       HE.PrandtlNumber(Unity(i).Properties.Cold.Average.K,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.PR          =       ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
 
 "Tube Side Prandtl Number at Wall"
-        Unity(i).Tubes.HeatTransfer.PRw         =       HE.PrandtlNumber(Unity(i).Properties.Cold.Wall.K,Unity(i).Properties.Cold.Wall.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Wall.Mu);
+        Unity(i).Tubes.HeatTransfer.PRw         =       ((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
 
 "Tube Side Film Coefficient"
-        Unity(i).Tubes.HeatTransfer.htube       =       HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K,
-        Unity(i).Tubes.HeatTransfer.Phi,Sumary.Lz(i));
+        Unity(i).Tubes.HeatTransfer.htube       =       HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Cold.Average.K,Sumary.Lz(i))*Unity(i).Tubes.HeatTransfer.Phi;
         
 "Shell Side Prandtl Number"
-        Unity(i).Shell.HeatTransfer.PR          =       HE.PrandtlNumber(Unity(i).Properties.Hot.Average.K,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Average.Mu);   
+        Unity(i).Shell.HeatTransfer.PR          =       ((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K; 
 
 "Shell Side Prandtl Number at Wall"
-        Unity(i).Shell.HeatTransfer.PRw         =       HE.PrandtlNumber(Unity(i).Properties.Hot.Wall.K,Unity(i).Properties.Hot.Wall.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Wall.Mu);    
+        Unity(i).Shell.HeatTransfer.PRw         =       ((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;  
 
 "Tube Side Pressure Drop"
@@ -1146 +1173 @@
         
 "Shell Side Reynolds Number"
-        Unity(i).Shell.HeatTransfer.Re                  =       HE.ShellReynoldsNumber(Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Shell.HeatTransfer.Re=(Dotube*Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)/Unity(i).Properties.Cold.Average.Mu;
 
 "Shell Heat Transfer Coefficient"
-        Unity(i).Shell.HeatTransfer.hshell              =       HE.ShellFilmCoeff(Unity(i).Shell.HeatTransfer.Ji,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,
-        Unity(i).Properties.Cold.Inlet.Fw,Unity(i).Shell.HeatTransfer.Sm,Unity(i).Shell.HeatTransfer.PR,Unity(i).Shell.HeatTransfer.Jtotal,Unity(i).Shell.HeatTransfer.Phi);
+        Unity(i).Shell.HeatTransfer.hshell      =Unity(i).Shell.HeatTransfer.Ji*(Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*
+        (Unity(i).Properties.Cold.Inlet.Fw/Unity(i).Shell.HeatTransfer.Sm)*(Unity(i).Shell.HeatTransfer.PR^(-2/3))*Unity(i).Shell.HeatTransfer.Jtotal*Unity(i).Shell.HeatTransfer.Phi;
 
 
@@ -1164 +1191 @@
 
 "Tube Side Velocity"
-        Unity(i).Tubes.HeatTransfer.Vtube               =       HE.TubeVelocity(Unity(i).Properties.Hot.Inlet.Fw,Unity(i).Properties.Hot.Average.rho);
+        Unity(i).Tubes.HeatTransfer.Vtube               = Unity(i).Properties.Hot.Inlet.Fw*Tpass/((Pi*Ditube*Ditube/4)*Unity(i).Properties.Hot.Average.rho*Ntt);
 
 "Tube Side Reynolds Number"
-        Unity(i).Tubes.HeatTransfer.Re                  =       HE.TubeReynoldsNumber(Unity(i).Properties.Hot.Average.rho,Unity(i).Tubes.HeatTransfer.Vtube,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.Re                  =       (Unity(i).Properties.Hot.Average.rho*Unity(i).Tubes.HeatTransfer.Vtube*Ditube)/Unity(i).Properties.Hot.Average.Mu;
         
 "Tube Side Prandtl Number"
-        Unity(i).Tubes.HeatTransfer.PR                  =       HE.PrandtlNumber(Unity(i).Properties.Hot.Average.K,Unity(i).Properties.Hot.Average.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Average.Mu);
+        Unity(i).Tubes.HeatTransfer.PR                  =((Unity(i).Properties.Hot.Average.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Average.Mu)/Unity(i).Properties.Hot.Average.K;
         
 "Tube Side Prandtl Number at Wall"
-        Unity(i).Tubes.HeatTransfer.PRw                 =       HE.PrandtlNumber(Unity(i).Properties.Hot.Wall.K,Unity(i).Properties.Hot.Wall.Cp,Unity(i).Properties.Hot.Average.Mw,Unity(i).Properties.Hot.Wall.Mu);
+        Unity(i).Tubes.HeatTransfer.PRw                 =       ((Unity(i).Properties.Hot.Wall.Cp/Unity(i).Properties.Hot.Average.Mw)*Unity(i).Properties.Hot.Wall.Mu)/Unity(i).Properties.Hot.Wall.K;
         
         
 "Tube Side Film Coefficient"
-        Unity(i).Tubes.HeatTransfer.htube               =       HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K,
-        Unity(i).Tubes.HeatTransfer.Phi,Sumary.Lz(i));
+        Unity(i).Tubes.HeatTransfer.htube               =       HE.TubeFilmCoeffIncremental(Unity(i).Tubes.HeatTransfer.Re,Unity(i).Tubes.HeatTransfer.PR,Unity(i).Properties.Hot.Average.K,Sumary.Lz(i))*Unity(i).Tubes.HeatTransfer.Phi;
         
 "Shell Side Prandtl Number"
-        Unity(i).Shell.HeatTransfer.PR                  =       HE.PrandtlNumber(Unity(i).Properties.Cold.Average.K,Unity(i).Properties.Cold.Average.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Average.Mu);
+        Unity(i).Shell.HeatTransfer.PR                  =       ((Unity(i).Properties.Cold.Average.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Average.Mu)/Unity(i).Properties.Cold.Average.K;
 
 "Shell Side Prandtl Number at Wall"
-        Unity(i).Shell.HeatTransfer.PRw                 =       HE.PrandtlNumber(Unity(i).Properties.Cold.Wall.K,Unity(i).Properties.Cold.Wall.Cp,Unity(i).Properties.Cold.Average.Mw,Unity(i).Properties.Cold.Wall.Mu);
+        Unity(i).Shell.HeatTransfer.PRw                 =       ((Unity(i).Properties.Cold.Wall.Cp/Unity(i).Properties.Cold.Average.Mw)*Unity(i).Properties.Cold.Wall.Mu)/Unity(i).Properties.Cold.Wall.K;
 
 "Tube Side Pressure Drop"

mso/eml/heat_exchangers/HeatExchangerSimplified.mso

@@ -1 @@
- 
+#--------------------------------------------------------------------
+# Author: Gerson Balbueno Bicca 
+# $Id$
+#-------------------------------------------------------------------- 
 using "HEX_Engine";
 #=====================================================================
@@ -7 +10 @@
 PARAMETERS
 ext PP          as CalcObject   (Brief="External Physical Properties");
-ext HE          as CalcObject   (Brief="STHE Calculations",File="heatex.dll");
+        HE              as CalcObject   (Brief="STHE Calculations",File="heatex");
 ext NComp       as Integer      (Brief="Number of Components");
         M(NComp)        as molweight    (Brief="Component Mol Weight");
@@ -54 +57 @@
         then    
 "Heat Capacity Cold Stream"
-        Properties.Cold.Average.Cp      =       PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+        Properties.Cold.Average.Cp              =       PP.LiquidCp(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
         Properties.Cold.Inlet.Cp                =       PP.LiquidCp(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
         Properties.Cold.Outlet.Cp               =       PP.LiquidCp(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
@@ -61 +64 @@
         Properties.Cold.Average.rho     =       PP.LiquidDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
         Properties.Cold.Inlet.rho               =       PP.LiquidDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
-        Properties.Cold.Outlet.rho      =       PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
+        Properties.Cold.Outlet.rho              =       PP.LiquidDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
 "Viscosity Cold Stream"
-        Properties.Cold.Average.Mu =    PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
-        Properties.Cold.Inlet.Mu                =               PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
-        Properties.Cold.Outlet.Mu       =       PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
+        Properties.Cold.Average.Mu              =       PP.LiquidViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+        Properties.Cold.Inlet.Mu                =       PP.LiquidViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+        Properties.Cold.Outlet.Mu               =       PP.LiquidViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
 "Conductivity Cold Stream"
-        Properties.Cold.Average.K        =      PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
-        Properties.Cold.Inlet.K                         =               PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
-        Properties.Cold.Outlet.K                =               PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
+        Properties.Cold.Average.K               =       PP.LiquidThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+        Properties.Cold.Inlet.K                 =       PP.LiquidThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+        Properties.Cold.Outlet.K                =       PP.LiquidThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
 "Heat Capacity Cold Stream"
@@ -77 +80 @@
         
 "Viscosity Cold Stream"
-        Properties.Cold.Wall.Mu =       PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
+        Properties.Cold.Wall.Mu                 =       PP.LiquidViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
 
 "Conductivity Cold Stream"
-        Properties.Cold.Wall.K =        PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
+        Properties.Cold.Wall.K                  =       PP.LiquidThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
 
 
@@ -92 +95 @@
 "Mass Density Cold Stream"
         Properties.Cold.Average.rho     =       PP.VapourDensity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
-        Properties.Cold.Inlet.rho       =       PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
-        Properties.Cold.Outlet.rho      =       PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
+        Properties.Cold.Inlet.rho               =       PP.VapourDensity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+        Properties.Cold.Outlet.rho              =       PP.VapourDensity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
 "Viscosity Cold Stream"
-        Properties.Cold.Average.Mu =    PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
-        Properties.Cold.Inlet.Mu =              PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
-        Properties.Cold.Outlet.Mu =             PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
+        Properties.Cold.Average.Mu              =       PP.VapourViscosity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+        Properties.Cold.Inlet.Mu                =       PP.VapourViscosity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+        Properties.Cold.Outlet.Mu               =       PP.VapourViscosity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
 
 "Conductivity Cold Stream"
-        Properties.Cold.Average.K =     PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
-        Properties.Cold.Inlet.K =               PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
-        Properties.Cold.Outlet.K =              PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
+        Properties.Cold.Average.K               =       PP.VapourThermalConductivity(Properties.Cold.Average.T,Properties.Cold.Average.P,Inlet.Cold.z);
+        Properties.Cold.Inlet.K                 =       PP.VapourThermalConductivity(Inlet.Cold.T,Inlet.Cold.P,Inlet.Cold.z);
+        Properties.Cold.Outlet.K                =       PP.VapourThermalConductivity(Outlet.Cold.T,Outlet.Cold.P,Outlet.Cold.z);
         
 "Heat Capacity Cold Stream"
-        Properties.Cold.Wall.Cp         =       PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
+        Properties.Cold.Wall.Cp                 =       PP.VapourCp(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
 
 
 "Viscosity Cold Stream"
-        Properties.Cold.Wall.Mu =       PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
+        Properties.Cold.Wall.Mu                 =       PP.VapourViscosity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
 
 "Conductivity Cold Stream"
-        Properties.Cold.Wall.K =        PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
+        Properties.Cold.Wall.K                  =       PP.VapourThermalConductivity(Properties.Cold.Wall.Twall,Properties.Cold.Average.P,Inlet.Cold.z);
         
         
@@ -140 +143 @@
 "Conductivity Hot Stream"
         Properties.Hot.Average.K        =               PP.LiquidThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);    
-        Properties.Hot.Inlet.K  =               PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);      
+        Properties.Hot.Inlet.K          =               PP.LiquidThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);      
         Properties.Hot.Outlet.K         =               PP.LiquidThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
 
 "Heat Capacity Hot Stream"
-        Properties.Hot.Wall.Cp  =               PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
+        Properties.Hot.Wall.Cp          =               PP.LiquidCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
 
 "Viscosity Hot Stream"
-        Properties.Hot.Wall.Mu  =               PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);     
+        Properties.Hot.Wall.Mu          =               PP.LiquidViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);     
 
 "Conductivity Hot Stream"
-        Properties.Hot.Wall.K   =               PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
+        Properties.Hot.Wall.K           =               PP.LiquidThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
         
 
@@ -172 +175 @@
 "Conductivity Hot Stream"
         Properties.Hot.Average.K        =               PP.VapourThermalConductivity(Properties.Hot.Average.T,Properties.Hot.Average.P,Inlet.Hot.z);    
-        Properties.Hot.Inlet.K  =               PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);      
+        Properties.Hot.Inlet.K          =               PP.VapourThermalConductivity(Inlet.Hot.T,Inlet.Hot.P,Inlet.Hot.z);      
         Properties.Hot.Outlet.K         =               PP.VapourThermalConductivity(Outlet.Hot.T,Outlet.Hot.P,Outlet.Hot.z);   
 
 "Heat Capacity Hot Stream"
-        Properties.Hot.Wall.Cp  =               PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
+        Properties.Hot.Wall.Cp          =               PP.VapourCp(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
 
 "Viscosity Hot Stream"
-        Properties.Hot.Wall.Mu  =               PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
+        Properties.Hot.Wall.Mu          =               PP.VapourViscosity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);
 
 "Conductivity Hot Stream"
-        Properties.Hot.Wall.K   =               PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
+        Properties.Hot.Wall.K           =               PP.VapourThermalConductivity(Properties.Hot.Wall.Twall,Properties.Hot.Average.P,Inlet.Hot.z);   
 
 
@@ -195 +198 @@
 "Cold Stream Heat Capacity"
         Details.Cc =Inlet.Cold.F*Properties.Cold.Average.Cp;
-        
-"Heat Capacity Ratio"
-
-Details.Cmin  = min([Details.Ch,Details.Cc]);
-Details.Cmax  = max([Details.Ch,Details.Cc]);
-Details.Cr    = Details.Cmin/Details.Cmax;
+
+"Minimum Heat Capacity"
+        Details.Cmin  = min([Details.Ch,Details.Cc]);
+
+"Maximum Heat Capacity"
+        Details.Cmax  = max([Details.Ch,Details.Cc]);
+
+"Heat Capacity Ratio"   
+        Details.Cr    = Details.Cmin/Details.Cmax;
 #=====================================================================
 #       Energy Balance
@@ -285 +291 @@
 #       Basic Model for Heat Exchangers - LMTD Method
 #=====================================================================
-
 VARIABLES
 
+DT0     as temp_delta   (Brief="Temperature Difference at Inlet",Lower=1);
+DTL             as temp_delta   (Brief="Temperature Difference at Outlet",Lower=1);
 LMTD    as temp_delta   (Brief="Logarithmic Mean Temperature Difference",Lower=1);
 Fc              as positive             (Brief="LMTD Correction Factor",Lower=0.5);
-MTD             as temp_delta   (Brief="Mean Temperature Difference");
-
-EQUATIONS
+MTD             as temp_delta   (Brief="Mean Temperature Difference",Lower=1);
+
+EQUATIONS
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+#                       Log Mean Temperature Difference
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+if abs(DT0 - DTL) > 0.05*max(abs([DT0,DTL])) 
+        
+        then
+"Log Mean Temperature Difference"
+        LMTD= (DT0-DTL)/ln(DT0/DTL);
+
+        else
+        
+if DT0*DTL equal 0
+        
+        then
+"Log Mean Temperature Difference"
+        LMTD = 0.5*(DT0+DTL);
+        
+        else
+"Log Mean Temperature Difference"
+        LMTD = 0.5*(DT0+DTL)*(1-(DT0-DTL)^2/(DT0*DTL)*(1+(DT0-DTL)^2/(DT0*DTL)/2)/12);
+        
+end
+        
+end
 
 "Exchange Surface Area"
@@ -313 +344 @@
 Model HeatExchanger_LMTD        as Heatex_Basic_LMTD
 
-EQUATIONS
-
-"Cocurrent Flow LMTD"
-        LMTD = HE.LogMeanTemperature(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
+PARAMETERS
+
+        Side as Integer         (Brief="Flow Direction",Lower=0,Upper=1);
+
+SET
+
+Side = HE.FlowDir(); # Return Flow Direction
+
+EQUATIONS
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+#       Flow Direction 
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+if Side equal 0
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+#       Cocurrent Flow
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#   
+        then
+"Temperature Difference at Inlet"
+        DT0 = Inlet.Hot.T - Inlet.Cold.T;
+
+"Temperature Difference at Outlet"
+        DTL = Outlet.Hot.T - Outlet.Cold.T;
+
+        else
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+#       Counter Flow
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#           
+"Temperature Difference at Inlet"
+        DT0 = Inlet.Hot.T - Outlet.Cold.T;
+
+"Temperature Difference at Outlet"
+        DTL = Outlet.Hot.T - Inlet.Cold.T;
+end
 
 end
@@ -325 +385 @@
 #=====================================================================  
 EQUATIONS
-
-"Counter Flow LMTD"
-        LMTD = HE.CounterLMTD(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
-
+"Temperature Difference at Inlet"
+        DT0 = Inlet.Hot.T - Outlet.Cold.T;
+
+"Temperature Difference at Outlet"
+        DTL = Outlet.Hot.T - Inlet.Cold.T;
+        
 "LMTD Correction Factor"
         Fc = HE.EshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
@@ -339 +401 @@
 #=====================================================================
 EQUATIONS
-
-"Counter Flow LMTD"
-        LMTD = HE.CounterLMTD(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);
-
+"Temperature Difference at Inlet"
+        DT0 = Inlet.Hot.T - Outlet.Cold.T;
+
+"Temperature Difference at Outlet"
+        DTL = Outlet.Hot.T - Inlet.Cold.T;
+        
 "LMTD Correction Factor"
         Fc = HE.FshellCorrectionFactor(Inlet.Hot.T,Outlet.Hot.T,Inlet.Cold.T,Outlet.Cold.T);

mso/eml/heat_exchangers/Mheatex.mso

@@ -13 +13 @@
 * $Id$
 *--------------------------------------------------------------------*#
+
 using "streams.mso";
 
-Model Parameters
-#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
-#       Heat Exchanger External Parameters 
-#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
-PARAMETERS
-        
-ext PP                  as CalcObject   (Brief="Physical Properties");
-ext     NComp           as Integer              (Brief="Number of Components"); 
-ext     Ncold           as Integer              (Brief="Number of Inlet Cold Streams",Lower=1);
-ext     Nhot            as Integer              (Brief="Number of Inlet Hot Streams",Lower=1);
-        
-end     
-
-Model Inlet_Main_Stream         as Parameters
+Model Inlet_Main_Stream         
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
 #       Inlet Streams
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+PARAMETERS
+
+        Ncold   as Integer      (Brief="Number of Inlet Cold Streams",Lower=1);
+        Nhot    as Integer      (Brief="Number of Inlet Hot Streams",Lower=1);
+        
 VARIABLES
 
-Hot(Nhot)               as stream;# Inlet Hot Streams
-Cold(Ncold)     as stream;# Inlet Cold Streams
+        Hot  (Nhot)     as stream;# Inlet Hot Streams
+        Cold (Ncold)    as stream;# Inlet Cold Streams
         
 end
 
-Model Outlet_Main_Stream        as Parameters
+Model Outlet_Main_Stream        
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
 #       Outlet Streams
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+PARAMETERS
+
+        Ncold   as Integer      (Brief="Number of Inlet Cold Streams",Lower=1);
+        Nhot    as Integer      (Brief="Number of Inlet Hot Streams",Lower=1);
+        
 VARIABLES
 
-Hot(Nhot)               as stream_therm;# Outlet Hot Streams
-Cold(Ncold)     as stream_therm;# Outlet Cold Streams
+        Hot  (Nhot)     as stream_therm;# Outlet Hot Streams
+        Cold (Ncold)    as stream_therm;# Outlet Cold Streams
         
 end
@@ -56 +54 @@
 PARAMETERS
 
-ext     HE              as CalcObject   (Brief="Cold Box Calculations",File="Heatex.dll");
+ext PP          as CalcObject   (Brief="Physical Properties");
+ext     NComp   as Integer              (Brief="Number of Components"); 
+        HE              as CalcObject   (Brief="Cold Box Calculations",File="heatex");
         Side    as Integer              (Brief="Flow Direction",Lower=0,Upper=1);
-
-SET
-
-Side = HE.FlowDir();
-
+        Ncold   as Integer              (Brief="Number of Inlet Cold Streams",Lower=1);
+        Nhot    as Integer              (Brief="Number of Inlet Hot Streams",Lower=1);
+        
 VARIABLES
 
@@ -71 +69 @@
         LMTD    as temp_delta   (Brief="Logarithmic Mean Temperature Difference");
         UA      as positive     (Unit="W/K");
-        
+
+SET
+
+Side = HE.FlowDir();
+
+Inlet.Ncold  = Ncold;
+Outlet.Ncold = Ncold;
+
+Inlet.Nhot  = Nhot ;
+Outlet.Nhot = Nhot ;
+
+
 EQUATIONS
 

mso/eml/heat_exchangers/heater.mso

@@ -22 +22 @@
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#   
 PARAMETERS
-        ext PP                  as CalcObject   (Brief="Physical Properties");
-        ext     NComp           as Integer              (Brief="Number of Components"); 
-        ext     Ninlet      as Integer          (Brief="Number of Inlet Streams",Lower=1);
+        ext PP          as CalcObject(Brief="Physical Properties");
+        ext     NComp   as Integer       (Brief="Number of Components");        
+                Ninlet  as Integer       (Brief="Number of Inlet Streams",Lower=1);
         
 VARIABLES
         in  Inlet(Ninlet)       as stream;              #(Brief="Inlet Streams")
         out Outlet              as stream_therm;#(Brief="Outlet Stream") 
-        Q                               as power                (Brief="Heat Transfer");
+        Q                               as power                (Brief = "Heat Transfer");
         Vfrac                           as fraction     (Brief = "Vapor fraction Outlet Stream");
         Lfrac                           as fraction     (Brief = "Liquid fraction Outlet Stream");
@@ -45 +45 @@
 end
 
-
 "Outlet Vapourisation Fraction"
         Outlet.v = PP.VapourFraction(Outlet.T,Outlet.P,Outlet.z);
@@ -58 +57 @@
 
 Model Heater            as Heater_Cooler_Basic
-
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+#       Heater 
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#   
 EQUATIONS
 
@@ -67 +68 @@
 
 Model Cooler            as Heater_Cooler_Basic
-
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
+#       Cooler 
+#+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#   
 EQUATIONS
 

mso/sample/heat_exchangers/Eshell_Detailed_LMTD.mso

@@ +1 @@
+
+#===============================================================
+#   Heat Exchanger TEMA E Shell - LMTD Method
+#  converge após várias tentativas
+#===============================================================
 
 using "Heat_exchangers/HeatExchangerDetailed.mso";
 
-#===============================================================
-#   Heat Exchanger TEMA E Shell - LMTD Method
-#===============================================================        
 FlowSheet Exchanger_E_shell_Detailed_LMTD
         
 DEVICES
+
         exchanger               as E_Shell_LMTD_Det;
         streamhot_in    as streamTP;
@@ -13 +16 @@
         
 CONNECTIONS
+
         streamhot_in    to exchanger.Inlet.Hot;
         streamcold_in   to exchanger.Inlet.Cold;
         
 PARAMETERS
+
         PP                      as CalcObject   (File="vrpp");
         NComp           as Integer;
-        HE                      as CalcObject   (File="Heatex");
 
 SET
-
+#===============================================================
+#   Heat Exchanger Options
+#===============================================================        
         PP.Components           = ["water"];
-        HE.HotSide              = "Shell";
-        HE.LMTDcorrection       = "Bowmann";
-#       HE.LMTDcorrection       = "Fakeri";
         PP.LiquidModel          = "RK";
         PP.VapourModel          = "RK";
         NComp                           = PP.NumberOfComponents;
+        
+        exchanger.HE.HotSide            = "Shell";
+
+#   LMTD Correction Factor
+        exchanger.HE.LMTDcorrection     = "Fakeri";
+#       exchanger.HE.LMTDcorrection     = "Bowmann";
+
+#   Heat Transfer Correlation
+#       exchanger.HE.LaminarFlow                = "Schlunder";
+#       exchanger.HE.TurbulentFlow              = "SiederTate";
+#       exchanger.HE.TransitionFlow     = "GnielinskiII";
+
 #=====================================================================
 #       Shell Geometrical Parameters
 #=====================================================================       
-exchanger.Tpass                         = 2;
+exchanger.Tpass                         = 4;
 exchanger.Dishell               = 0.75  *"m";
 exchanger.Lcf                   = 0.043  *"m";
@@ -80 +95 @@
 #       Baffle Spacing
 #=====================================================================
-        exchanger.Baffles.Ls    = 0.622         *"m";
+        exchanger.Baffles.Ls    = 0.622  *"m";
         exchanger.Baffles.Lsi   = 0.807  *"m";
 
 OPTIONS
- mode                           = "steady";
- relativeAccuracy       = 1e-6;
+
+ mode = "steady";
 
 end

mso/sample/heat_exchangers/Eshell_Detailed_NTU.mso

@@ -1 +1 @@
 
 using "Heat_exchangers/HeatExchangerDetailed.mso";
-
 #===============================================================
 #   Heat Exchanger TEMA E Shell - NTU Method
+#   erro ao trocar de modelo termodinamico
 #===============================================================
 FlowSheet Exchanger_E_shell_Detailed_NTU
@@ -19 +19 @@
         PP                      as CalcObject   (File="vrpp");
         NComp           as Integer;
-        HE                      as CalcObject   (File="Heatex");
+
+SET
+        PP.Components                   = ["benzene","toluene"];
+        exchanger.HE.HotSide    = "Shell";
+        PP.LiquidModel                  = "PR";
+        PP.VapourModel                  = "PR";
+        NComp                                   = PP.NumberOfComponents;
         
-SET
-        PP.Components           = ["benzene","toluene"];
-        HE.HotSide              = "Shell";
-        PP.LiquidModel          = "PR";
-        PP.VapourModel          = "PR";
-        NComp                           = PP.NumberOfComponents;
+#   Heat Transfer Correlation
+        exchanger.HE.LaminarFlow                = "Schlunder";
+        exchanger.HE.TurbulentFlow              = "SiederTate";
+        exchanger.HE.TransitionFlow     = "GnielinskiII";
+
 #=====================================================================
 #       Shell Geometrical Parameters
 #=====================================================================    
-exchanger.Tpass                         = 4;
+exchanger.Tpass                         = 2;
 exchanger.Dishell               = 0.75  *"m";
 exchanger.Lcf                   = 0.043 *"m";
@@ -81 +86 @@
 
 OPTIONS
- mode                           = "steady";
- relativeAccuracy       = 1e-6;
+
+ mode   = "steady";
 
 end

mso/sample/heat_exchangers/Eshell_Discretized_LMTD.mso

@@ -24 +24 @@
 #    Parâmetros
 #===============================================================
-        PP                      as CalcObject   (File="vrpp.dll");
+        PP                      as CalcObject   (File="vrpp");
         NComp           as Integer;
-        HE                      as CalcObject   (File="Heatex.dll");
         
 SET
@@ -33 +32 @@
 #===============================================================
         PP.Components           = ["benzene","toluene"];
-        HE.HotSide              = "Shell";
-        HE.LMTDcorrection       = "Fakeri";
+        exchanger.HE.HotSide            = "Shell";
         PP.LiquidModel          = "PR";
         PP.VapourModel          = "PR";
         NComp                           = PP.NumberOfComponents;
+        
+        exchanger.HE.LMTDcorrection     = "Fakeri";
+#       exchanger.HE.LaminarFlow        = "Schlunder";
+#       exchanger.HE.TurbulentFlow      = "SiederTate";
+#       exchanger.HE.TransitionFlow     = "GnielinskiII";
 #=====================================================================
 #       Parâmetros geométricos do casco
@@ -76 +79 @@
         streamhot_in.F  = 40    * "mol/s";
         streamhot_in.T  = 393.15* "K";
-        streamhot_in.v  = 0    ;
         streamhot_in.z  = [1,0]     ;
+        streamhot_in.v  = 0     ;
         streamhot_in.P  = 740     * "kPa";
 #============================================
@@ -84 +87 @@
         streamcold_in.F  = 122                  * "mol/s";
         streamcold_in.T  = 333.15   * "K";
+        streamcold_in.z  = [0,1];
         streamcold_in.v  = 0;
-        streamcold_in.z  = [0,1];
         streamcold_in.P  = 2210*"kPa";
 #=====================================================================
@@ -98 +101 @@
 #============================================
  mode                   = "steady";
- relativeAccuracy = 1e-6;
 
 end

mso/sample/heat_exchangers/Eshell_Discretized_NTU.mso

@@ -24 +24 @@
 #    Parâmetros
 #===============================================================
-        PP                      as CalcObject   (File="vrpp.dll");
+        PP                      as CalcObject   (File="vrpp");
         NComp           as Integer;
-        HE                      as CalcObject   (File="Heatex.dll");
         
 SET
@@ -33 +32 @@
 #===============================================================
         PP.Components           = ["benzene","toluene"];
-        HE.HotSide              = "Shell";
-#       HE.TurbulentFlow        = "SiederTate";
-        PP.LiquidModel          = "PR";
-        PP.VapourModel          = "PR";
+        exchanger.HE.HotSide            = "Shell";
+        exchanger.HE.TurbulentFlow      = "SiederTate";
+        PP.LiquidModel          = "RK";
+        PP.VapourModel          = "RK";
         NComp                           = PP.NumberOfComponents;
 #=====================================================================
@@ -97 +96 @@
 #   Simulation Options
 #============================================
- mode                   = "steady";
- relativeAccuracy = 1e-6;
+ mode = "steady";
 
 end

mso/sample/heat_exchangers/MheaterTeste.mso

@@ -29 +29 @@
         streamcold2             as stream_therm;
         
+PARAMETERS
+
+        PP                      as CalcObject   (File="vrpp");
+        NComp           as Integer              (Brief="Number Components");
+        
+SET
+
+        PP.LiquidModel  = "RK";
+        PP.VapourModel  = "RK";
+        PP.Components   = ["water","methanol","benzene"];
+        NComp                   = PP.NumberOfComponents;
+        
+        
+        Mheater.HE.FlowDirection  = "Cocurrent";
+        Mheater.Ncold = 2;
+        Mheater.Nhot  = 2;
+        
 CONNECTIONS
 
@@ -36 +53 @@
         streamcold2     to Mheater.Inlet.Cold(2);
         
-PARAMETERS
 
-        PP                      as CalcObject   (File="vrpp.dll");
-        HE                      as CalcObject   (File="heatex.dll");
-        NComp           as Integer              (Brief="Number Components");
-        Ncold           as Integer              (Brief="Number of Inlet Cold Streams");
-        Nhot            as Integer              (Brief="Number of Inlet Hot Streams");
-        
-SET
-
-        PP.LiquidModel  = "RK";
-        PP.VapourModel  = "RK";
-        PP.Components   = ["water","methanol","benzene"];
-        NComp                   = PP.NumberOfComponents;
-        Nhot            = 2;
-    Ncold               = 2;
-#       HE.FlowDirection      = "Cocurrent";
-        HE.FlowDirection      = "Counter";
 
         
@@ -104 +104 @@
 OPTIONS
   
-        mode                    = "steady";
-        relativeAccuracy        = 1e-8;
+mode    = "steady";
         
 end

mso/sample/heat_exchangers/Mheater_project.mso

@@ -38 +38 @@
 PARAMETERS
 
-        PP                      as CalcObject   (File="vrpp.dll");
-        HE                      as CalcObject   (File="heatex.dll");
+        PP                      as CalcObject   (File="vrpp");
+#       HE                      as CalcObject   (File="heatex");
         NComp           as Integer              (Brief="Number Components");
         Ncold           as Integer              (Brief="Number of Inlet Cold Streams");
@@ -50 +50 @@
         PP.Components   = ["water","methanol","benzene"];
         NComp                   = PP.NumberOfComponents;
-        Nhot            = 2;
-    Ncold               = 2;
-        HE.FlowDirection      = "Cocurrent";
+        Mheater.Nhot            = 2;
+    Mheater.Ncold               = 2;
+        Mheater.HE.FlowDirection      = "Cocurrent";
 #       HE.FlowDirection      = "Counter";
 
@@ -105 +105 @@
   
         mode                    = "steady";
-        outputLevel             = "all";
         relativeAccuracy        = 1e-8;
         

mso/sample/heat_exchangers/Multipass_Detailed.mso

@@ -24 +24 @@
 #    Parâmetros
 #===============================================================
-        PP                      as CalcObject   (File="vrpp.dll");
+        PP                      as CalcObject   (File="vrpp");
         NComp           as Integer;
-        HE                      as CalcObject   (File="Heatex.dll");
         
 SET
@@ -33 +32 @@
 #===============================================================
         PP.Components           = ["benzene","toluene"];
-        HE.HotSide              = "Shell";
-#       HE.TurbulentFlow        = "SiederTate";
+        exchanger.HE.HotSide            = "Shell";
+        exchanger.HE.TurbulentFlow      = "SiederTate";
         PP.LiquidModel          = "PR";
         PP.VapourModel          = "PR";
@@ -100 +99 @@
 #============================================
  mode                   = "steady";
- relativeAccuracy = 1e-6;
 
 end

mso/sample/heat_exchangers/NTU_Method.mso

@@ -20 +20 @@
 PARAMETERS
 
-        PP                      as CalcObject   (File="vrpp");
-        HE                      as CalcObject   (File="heatex");
+        PP                      as CalcObject   (File="vrpp.dll");
+        HE                      as CalcObject   (File="heatex.dll");
         NComp           as Integer;
         

mso/sample/heat_exchangers/Sample_Simplified.mso

@@ -16 +16 @@
 PARAMETERS
         PP                      as CalcObject   (File="vrpp");
-        HE              as CalcObject   (File="Heatex");
         NComp           as Integer;
 SET
@@ -25 +24 @@
         NComp                           = PP.NumberOfComponents;
         
-#       HE.FlowDirection      = "Cocurrent";
-        HE.FlowDirection      = "Counter";
+#       exchanger.HE.FlowDirection      = "Cocurrent";
+        exchanger.HE.FlowDirection      = "Counter";
 
 SPECIFY
@@ -68 +67 @@
 PARAMETERS
         PP                      as CalcObject   (File="vrpp");
-        HE              as CalcObject   (File="Heatex");
         NComp           as Integer;
 SET
@@ -77 +75 @@
         NComp                           = PP.NumberOfComponents;
         
-#       HE.FlowDirection      = "Cocurrent";
-        HE.FlowDirection      = "Counter";
+        exchanger.HE.FlowDirection      = "Cocurrent";
+#       exchanger.HE.FlowDirection      = "Counter";
 
 SPECIFY
@@ -100 +98 @@
         
 OPTIONS
- mode                     = "steady";
- relativeAccuracy = 1e-6;
-
+ mode  = "steady";
 end
 

mso/sample/heat_exchangers/heater_project.mso

@@ -25 +25 @@
 PARAMETERS
 
-        PP                      as CalcObject (File="vrpp.dll");
+        PP                      as CalcObject (File="vrpp");
         NComp           as Integer;
-        Ninlet      as Integer     (Brief="Number of Inlet Streams");
         
 SET
@@ -36 +35 @@
         PP.Components   = ["water"];
         NComp                   = PP.NumberOfComponents;
-        Ninlet          = 1;
+        heater.Ninlet          = 1;
 
 SPECIFY
@@ -55 +54 @@
   
         mode                    = "steady";
-        relativeAccuracy        = 1e-8;
-        
         
 end
@@ -77 +74 @@
         PP                      as CalcObject (File="vrpp.dll");
         NComp   as Integer;
-        Ninlet      as Integer     (Brief="Number of Inlet Streams");
         
 SET
@@ -85 +81 @@
         PP.Components   = ["water"];
         NComp                   = PP.NumberOfComponents;
-        Ninlet          = 1;
+        cooler.Ninlet          = 1;
 
 SPECIFY
@@ -103 +99 @@
 OPTIONS
   
-        mode                    = "steady";
-        relativeAccuracy        = 1e-8;
+        mode    = "steady";
         
         

mso/sample/heat_exchangers/sampleEshell.mso

@@ -7 +7 @@
 
 PP               as CalcObject (Brief="Physical Properties",File="vrpp");
-HE           as CalcObject (Brief="Heat Exchangers Calculations",File="Heatex");
-NComp as Integer;
+NComp    as Integer;
         
 SET
@@ -40 +39 @@
 streamcold.T            = 300 * "K";    
 streamcold.z            = [0,0.5, 0.1, 0.4];
-exchanger.Outlet.Cold.T = 350*"K";
+exchanger.Outlet.Cold.T = 347*"K";
         
 OPTIONS

mso/sample/heat_exchangers/sampleEshell_LMTD.mso

@@ -6 +6 @@
 PARAMETERS
 
-PP               as CalcObject (Brief="Physical Properties",File="vrpp");
-HE           as CalcObject (Brief="Heat Exchangers Calculations",File="Heatex");
-NComp as Integer;
+PP              as CalcObject (Brief="Physical Properties",File="vrpp");
+NComp   as Integer;
         
-SET
-
-PP.LiquidModel          = "PR";
-PP.VapourModel          = "PR";
-PP.Components   = ["water","n-butane", "benzene", "n-octane" ];
-NComp                                           = PP.NumberOfComponents;
 
 DEVICES
@@ -21 +14 @@
 streamhot       as streamTP;
 streamcold  as streamTP;
+
+SET
+
+PP.LiquidModel  = "PR";
+PP.VapourModel  = "PR";
+PP.Components   = ["water","n-butane", "benzene", "n-octane" ];
+NComp                   = PP.NumberOfComponents;
+
+        exchanger.HE.LMTDcorrection     = "Bowmann";
+#       HE.LMTDcorrection       = "Fakeri";
         
+
 CONNECTIONS
 
-streamhot               to exchanger.Inlet.Hot;
+streamhot       to exchanger.Inlet.Hot;
 streamcold      to exchanger.Inlet.Cold;
         
@@ -32 +36 @@
 exchanger.PressureDrop.Hot.Pdrop        = 0.1*"kPa";
 exchanger.PressureDrop.Cold.Pdrop   = 0.2*"kPa";
+
 streamhot.F             = 20    * "kmol/h";
 streamhot.T             = 450     * "K";
 streamhot.P             = 120    * "kPa";
 streamhot.z             = [1,0,0,0];
+
 streamcold.F            = 10 * "kmol/h";
 streamcold.P            = 120 * "kPa";
 streamcold.T            = 300 * "K";    
 streamcold.z            = [0,0.5, 0.1, 0.4];
-exchanger.Outlet.Cold.T = 350*"K";
+
+exchanger.Outlet.Cold.T = 330*"K";
         
 OPTIONS
- mode                     = "steady";
+ mode   = "steady";
  end

mso/sample/heat_exchangers/sampleLMTD.mso

@@ -7 +7 @@
 
 PP               as CalcObject (Brief="Physical Properties",File="vrpp");
-HE           as CalcObject (Brief="Heat Exchangers Calculations",File="Heatex");
-NComp as Integer;
+NComp   as Integer;
         
-SET
 
-PP.LiquidModel          = "PR";
-PP.VapourModel          = "PR";
-PP.Components   = ["water","n-butane", "benzene", "n-octane" ];
-NComp                                           = PP.NumberOfComponents;
-HE.FlowDirection    = "Cocurrent";
         
 DEVICES
@@ -22 +15 @@
 streamhot       as streamTP;
 streamcold  as streamTP;
-        
+SET
+
+PP.LiquidModel          = "PR";
+PP.VapourModel          = "PR";
+PP.Components           = ["water","n-butane", "benzene", "n-octane" ];
+NComp                           = PP.NumberOfComponents;
+exchanger.HE.FlowDirection    = "Cocurrent";    
 CONNECTIONS
 
-streamhot               to exchanger.Inlet.Hot;
+streamhot       to exchanger.Inlet.Hot;
 streamcold      to exchanger.Inlet.Cold;
         
 SPECIFY
 
-exchanger.Fc                                                                            = 1;
-exchanger.Details.U                                                             = 210*"W/(m^2*K)";
+"LMTD Correction Factor"
+        exchanger.Fc = 1;
+        
+"Overall Heat Transfer Coefficient"     
+        exchanger.Details.U = 210*"W/(m^2*K)";
+
 exchanger.PressureDrop.Hot.Pdrop        = 0.1*"kPa";
 exchanger.PressureDrop.Cold.Pdrop   = 0.2*"kPa";
-streamhot.F                                                                                     = 20    * "kmol/h";
-streamhot.T                                                                                     = 450     * "K";
+
+streamhot.F             = 20    * "kmol/h";
+streamhot.T             = 450     * "K";
 streamhot.P             = 120    * "kPa";
 streamhot.z             = [1,0,0,0];
+
 streamcold.F            = 10 * "kmol/h";
 streamcold.P            = 120 * "kPa";
 streamcold.T            = 300 * "K";    
 streamcold.z            = [0,0.5, 0.1, 0.4];
-exchanger.Outlet.Cold.T = 350*"K";
+
+exchanger.Outlet.Cold.T = 340*"K";
         
 OPTIONS

mso/sample/heat_exchangers/sampleNTU.mso

@@ -6 +6 @@
 PARAMETERS
 
-PP               as CalcObject (Brief="Physical Properties",File="vrpp");
-HE           as CalcObject (Brief="Heat Exchangers Calculations",File="Heatex");
+PP        as CalcObject (Brief="Physical Properties",File="vrpp");
 NComp as Integer;
         
+
+        
+DEVICES
+exchanger       as HeatExchanger_NTU;
+streamhot       as streamTP;
+streamcold  as streamTP;
+
 SET
 
@@ -16 +22 @@
 PP.Components   = ["water","n-butane", "benzene", "n-octane" ];
 NComp                                           = PP.NumberOfComponents;
-HE.FlowDirection    = "Cocurrent";
-        
-DEVICES
-exchanger       as HeatExchanger_NTU;
-streamhot       as streamTP;
-streamcold  as streamTP;
-        
+exchanger.HE.FlowDirection    = "Cocurrent";
+
 CONNECTIONS
 
 streamhot               to exchanger.Inlet.Hot;
-streamcold      to exchanger.Inlet.Cold;
+streamcold              to exchanger.Inlet.Cold;
         
 SPECIFY
@@ -41 +42 @@
 streamcold.T            = 300 * "K";    
 streamcold.z            = [0,0.5, 0.1, 0.4];
-exchanger.Outlet.Cold.T = 350*"K";
+exchanger.Outlet.Cold.T = 340*"K";
         
 OPTIONS
- mode                     = "steady";
+ mode   = "steady";
  end

mso/sample/heat_exchangers/sample_cooler.mso

@@ -14 +14 @@
 DEVICES
 
+
         cooler                  as Cooler;
         streamcold1     as stream_therm;
@@ -27 +28 @@
         PP                      as CalcObject (File="vrpp");
         NComp           as Integer;
-        Ninlet      as Integer     (Brief="Number of Inlet Streams");
         
 SET
@@ -36 +36 @@
         PP.Components   = ["water","methanol"];
         NComp                   = PP.NumberOfComponents;
-        Ninlet          = 2;
+        cooler.Ninlet   = 2;
 
 SPECIFY
@@ -55 +55 @@
         cooler.Outlet.P         = 1 * "atm";
 
-        
-        
 OPTIONS
-  
-        mode                    = "steady";
+        mode                    = "steady";
         relativeAccuracy        = 1e-8;
         

mso/sample/heat_exchangers/sample_heater.mso

@@ -25 +25 @@
 PARAMETERS
 
-        PP                      as CalcObject (File="vrpp");
-        NComp           as Integer;
-        Ninlet      as Integer     (Brief="Number of Inlet Streams");
+        PP              as CalcObject (File="vrpp");
+        NComp   as Integer;
         
 SET
@@ -36 +35 @@
         PP.Components   = ["water"];
         NComp                   = PP.NumberOfComponents;
-        Ninlet          = 1;
+        heater.Ninlet   = 1;
 
 SPECIFY
@@ -50 +49 @@
         heater.Outlet.T         = 363*"K";
         heater.Outlet.P         = 1 * "atm";
+
         
         
@@ -55 +55 @@
   
         mode                    = "steady";
-        outputLevel             = "all";
-        relativeAccuracy        = 1e-8;
-        
         
 end
@@ -78 +75 @@
         PP                      as CalcObject (File="vrpp");
         NComp   as Integer;
-        Ninlet      as Integer     (Brief="Number of Inlet Streams");
         
 SET
@@ -86 +82 @@
         PP.Components   = ["water"];
         NComp                   = PP.NumberOfComponents;
-        Ninlet          = 1;
+        cooler.Ninlet   = 1;
 
 SPECIFY
@@ -105 +101 @@
   
         mode                    = "steady";
-        outputLevel             = "all";
-        relativeAccuracy        = 1e-8;
-        
         
 end

mso/sample/heat_exchangers/samples1.mso

@@ -20 +20 @@
 PARAMETERS
 
-        PP                      as CalcObject   (File="vrpp.dll");
+        PP                      as CalcObject   (File="vrpp");
         NComp           as Integer;
-        HE              as CalcObject   (File="Heatex.dll");
         
 SET
@@ -30 +29 @@
         PP.Components           = ["water"];
         NComp                           = PP.NumberOfComponents;
-        HE.FlowDirection      = "Counter";
+        exchanger.HE.FlowDirection      = "Counter";
         
 SPECIFY
@@ -64 +63 @@
 #============================================
  mode                   = "steady";
- outputLevel    = "high";
- relativeAccuracy = 1e-8;
 
 end

mso/sample/heat_exchangers/samples2.mso

@@ -20 +20 @@
 PARAMETERS
 
-        PP                      as CalcObject   (File="vrpp.dll");
-        HE                      as CalcObject   (File="heatex.dll");
+        PP              as CalcObject   (File="vrpp");
         NComp   as Integer;
         
@@ -65 +64 @@
 #============================================
  mode                   = "steady";
- outputLevel    = "all";
- relativeAccuracy = 1e-8;
 
 end