Changeset 1005 for trunk/eml/heat_exchangers/HeatExchangerDetailed.mso

Timestamp:

Jun 12, 2017, 10:29:43 AM (5 years ago)

Author:

Argimiro Resende Secchi

Message:

Fixed convergence problem in detailed heat exchanger.

File:

• trunk/eml/heat_exchangers/HeatExchangerDetailed.mso (modified) (5 diffs)

trunk/eml/heat_exchangers/HeatExchangerDetailed.mso

@@ -107 +107 @@
 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++#
 Pi                                                      as constant   (Brief="Pi Number", Hidden = true, Default=3.14159265, Symbol = "\pi");
+pi                                                      as constant   (Brief="Pi Number", Unit = 'rad', Default=3.14159265, Symbol = "\pi");
 Aonozzle_Shell          as area                 (Brief="Shell Outlet Nozzle Area", Hidden = true, Lower=1E-6 ,  Symbol ="A_{nozzle\_out }^{shell}");
 Ainozzle_Shell          as area                 (Brief="Shell Inlet Nozzle Area", Hidden = true, Lower=1E-6 ,  Symbol ="A_{nozzle\_in }^{shell}");
@@ -130 +131 @@
 #"Pi Number"
         Pi      = 3.14159265;
+        pi = Pi * 'rad';
 
 #"Baffle cut angle in degrees"
-        Ods     = (360/Pi)*acos(1-0.02*Baffles.BaffleCut);
+        Ods     = (360/pi)*acos(1-0.02*Baffles.BaffleCut);
 
 #"Baffle cut angle relative to the centerline"
-        Octl    = (360/Pi)*acos((Shell.ShellID/(Shell.ShellID - Clearances.BundleToShell - Tubes.TubeOD))*(1-0.02*Baffles.BaffleCut));
+        Octl    = (360/pi)*acos((Shell.ShellID/(Shell.ShellID - Clearances.BundleToShell - Tubes.TubeOD))*(1-0.02*Baffles.BaffleCut));
 
 #"Fraction of number of tubes in baffle window"
-        Ftw     = (Octl/360)-sin(Octl*Pi/180)/(2*Pi);
+        Ftw     = (Octl/360)-sin(Octl*pi/180)/(2*Pi);
 
 #"Shell to baffle leakage area"
@@ -150 +152 @@
 
 #"comments"
-        Dw      = (4*abs((Pi*Shell.ShellID*Shell.ShellID*((Ods/360)-sin(Ods*Pi/180)/(2*Pi))/4)-(Tubes.NumberOfTubes*Pi*Tubes.TubeOD*Tubes.TubeOD*Ftw/4)))/(Pi*Tubes.TubeOD*Tubes.NumberOfTubes*Ftw+ Pi*Shell.ShellID*Ods/360);
+        Dw      = (4*abs((Pi*Shell.ShellID*Shell.ShellID*((Ods/360)-sin(Ods*pi/180)/(2*Pi))/4)-(Tubes.NumberOfTubes*Pi*Tubes.TubeOD*Tubes.TubeOD*Ftw/4)))/(Pi*Tubes.TubeOD*Tubes.NumberOfTubes*Ftw+ Pi*Shell.ShellID*Ods/360);
 
 #"Tube Side Inlet Nozzle Area"
@@ -180 +182 @@
 "Shell Stream Average Temperature"
         Shell.Properties.Average.T = 0.5*InletShell.T + 0.5*OutletShell.T;
+#       Shell.Properties.Average.T = InletShell.T;
         
 "Tube Stream Average Temperature"
         Tubes.Properties.Average.T = 0.5*InletTube.T + 0.5*OutletTube.T;
+#       Tubes.Properties.Average.T = InletTube.T;
         
 "Shell Stream Average Pressure"
-        Shell.Properties.Average.P = 0.5*InletShell.P+0.5*OutletShell.P;
+#       Shell.Properties.Average.P = 0.5*InletShell.P+0.5*OutletShell.P;
+        Shell.Properties.Average.P = InletShell.P;
         
 "Tube Stream Average Pressure"
-        Tubes.Properties.Average.P = 0.5*InletTube.P+0.5*OutletTube.P;
+#       Tubes.Properties.Average.P = 0.5*InletTube.P+0.5*OutletTube.P;
+        Tubes.Properties.Average.P = InletTube.P;
         
 "Shell Stream Average Molecular Weight"
@@ -355 +361 @@
 
 "Mass Velocity in Window Zone"
-        mw      = Shell.Properties.Inlet.Fw/sqrt(abs(Sm*abs((Pi*Shell.ShellID*Shell.ShellID*((Ods/360)-sin(Ods*Pi/180)/(2*Pi))/4)-(Tubes.NumberOfTubes*Pi*Tubes.TubeOD*Tubes.TubeOD*Ftw/4))));
+        mw      = Shell.Properties.Inlet.Fw/sqrt(abs(Sm*abs((Pi*Shell.ShellID*Shell.ShellID*((Ods/360)-sin(Ods*pi/180)/(2*Pi))/4)-(Tubes.NumberOfTubes*Pi*Tubes.TubeOD*Tubes.TubeOD*Ftw/4))));
 
 switch TubeFlowRegime