initial commit / add all books

1 files changed, 30 insertions, 0 deletions

diff --git a/2510/CH16/EX16.10/Ex16_10.sce b/2510/CH16/EX16.10/Ex16_10.sce
new file mode 100755
index 000000000..c338ed1f9
--- /dev/null
+++ b/2510/CH16/EX16.10/Ex16_10.sce
@@ -0,0 +1,30 @@
+//Variable declaration:
+t2 = 84.0                       //Temperature of water leaving the tube ( C)
+t1 = 16.0                       //Temperature of water entering the tube ( C)
+m1 = 10000.0/3600.0             //Mass flowrate of water (kg/s)
+T2 = 94.0                       //Temperature of oil leaving the shell ( C)
+T1 = 160.0                      //Temperature of oil entering the shell ( C)
+
+//Calculation:
+Tw = (t1+t2)/2.0                //Average bulk temperature of water ( C)
+To = (T1+T2)/2.0                //Average bulk temperature of oil ( C)
+//From table 16.1:
+p1 = 987.0                      //Density of water (kg/m^3)
+cp1 = 4176.0                    //Heat capacity of water (J/kg. C)
+p2 = 822.0                      //Density of oil (kg/m^3)
+Q = m1*cp1*(t2-t1)              //Heat load (W)
+cp2 = 4820.0                    //Heat capacity of oil (J/kg. C)
+m2 = Q/(cp2*(T1-T2))            //Mass flowrate of oil (kg/s)
+DT1 = T2-t1                     //Temperature driving force 1 ( C)
+DT2 = T1-t2                     //Temperature driving force 2 ( C)
+DTlm1 = ((DT1-DT2)/log(DT1/DT2))  //Log mean temperature driving force for ideal countercurrent heat exchanger ( C)
+P = (t2-t1)/(T1 - t1)           //Dimensionless ratio P
+R = (T1-T2)/(t2-t1)             //Dimensionless ratio R
+//From figure 16.7:
+F = 0.965                       //Correction factor
+DTlm2 = F*DTlm1                 //Log mean temperature driving force for 1-4 shell and tube exchanger ( C)
+
+//Result:
+printf("1. The heat load is : %.3f  MW .",Q/10**6)
+printf("2. The countercurrent flow log mean temperature difference is : %.0f C .",DTlm1)
+printf("3. The F correction factor and the corrected log mean temperature difference is : %.1f C .",DTlm2)