initial commit / add all books

13 files changed, 219 insertions, 0 deletions

diff --git a/2510/CH9/EX9.1/Ex9_1.sce b/2510/CH9/EX9.1/Ex9_1.sce
new file mode 100755
index 000000000..7f44e52ad
--- /dev/null
+++ b/2510/CH9/EX9.1/Ex9_1.sce
@@ -0,0 +1,15 @@
+//Variable declaration:
+D = 1.0                     //Diamete of vessel (ft)
+L = 1.5                     //Length of vessel (ft)
+T1 = 390.0                  //Surface temperature of vessel (°F)
+T2 = 50.0                   //Surrounding temperature of vessel (°F)
+h = 4.0                     //Convective heat transfer coefficient (Btu/h.ft.°F)
+pi = %pi
+
+//Calculation:
+A = pi*D*L+2*pi*(D/2)**2    //Total heat transfer area (ft^2)
+Q = h*A*(T1-T2)             //Rate of heat transfer (Btu/h)
+R = 1/(h*A)                 //Thermal resistance (°F.h/Btu)
+
+//Result:
+printf("The thermal resistance of vessel wal is : %.4f  °F.h/Btu.",R)
diff --git a/2510/CH9/EX9.10/Ex9_10.sce b/2510/CH9/EX9.10/Ex9_10.sce
new file mode 100755
index 000000000..07324c8e2
--- /dev/null
+++ b/2510/CH9/EX9.10/Ex9_10.sce
@@ -0,0 +1,19 @@
+//Variable declaration:
+D = 0.902/12.0                  //Inside diameter of tube (ft)
+T_in = 60.0                     //Temperature water entering the tube (°F)
+T_out = 70.0                    //Temperature water leaving the tube (°F)
+V = 7.0                         //Average wave velocity water (ft/s)
+p = 62.3                        //Density of water (lb/ft^3)
+mu = 2.51/3600.0                //Dynamic viscosity of water (lb/ft.s)
+Cp = 1.0                        //Viscosity of centipoise (Btu/lb.°F)
+k = 0.34                        //Thermal conductivity of water (Btu/h.ft.°F)
+
+//Calculation:
+Re = D*V*p/mu                   //Reynolds Number
+Pr = Cp*mu/k*3600               //Prandtl number
+//From equation 9.26:
+Nu = 0.023*(Re**0.8)*(Pr**0.4)  //Nusselt number
+h = (k/D)*Nu                    //Average film heat transfer coefficient (Btu/h.ft^2.°F)
+
+//Result:
+printf("The required average film heat transfer coefficient is : %.0f  Btu/h.ft^2.°F.",h)
diff --git a/2510/CH9/EX9.11/Ex9_11.sce b/2510/CH9/EX9.11/Ex9_11.sce
new file mode 100755
index 000000000..b37c3251d
--- /dev/null
+++ b/2510/CH9/EX9.11/Ex9_11.sce
@@ -0,0 +1,21 @@
+//Variable declaration:
+P = 1.0132 * 10**5              //Air pressure (Pa)
+T = 300.0+273.0                 //Air temperature (K)
+V = 5.0                         //Air flow velocity (m/s)
+D = 2.54/100.0                  //Diameter of tube (m)
+R = 287.0                       //Gas constant (m^2/s^2.K)
+//From Appendix:
+Pr = 0.713                      //Prandtl number of nitrogen
+mu = 1.784*10**(-5)             //Dynamic viscosity of nitrogen (kg/m.s)
+k = 0.0262                      //Thermal conductivity of nitrogen (W/m.K)
+Cp = 1.041                      //Heat capacity of nitrogen (kJ/kg.K)
+
+//Calculation:
+p = P/(R*T)                     //Density of air
+Re = D*V*p/mu                   //Reynolds number
+//From table 9.5:
+Nu = 0.023*(Re**0.8)*(Pr**0.3)  //Nusselt number
+h = (k/D)*Nu                    //Heat transfer coefficient (W/m^2.K)
+
+//Result:
+printf("The required Heat transfer coefficient is : %.2f  W/m^2.K.",h)
diff --git a/2510/CH9/EX9.12/Ex9_12.sce b/2510/CH9/EX9.12/Ex9_12.sce
new file mode 100755
index 000000000..b553e9162
--- /dev/null
+++ b/2510/CH9/EX9.12/Ex9_12.sce
@@ -0,0 +1,24 @@
+//Variable declaration:
+T1 = 15.0                   //Water entering temperature (°C)
+T2 = 60.0                   //Water leaving temperature (°C)
+D = 0.022                   //Inside diameter of tube (m)
+V = 0.355                   //Average water flow velocity (m/s)
+TC = 150.0                  //Outside wall temperature (°C)
+//From Appendix:
+p = 993.0                   //Density of water (kg/m^3)
+mu = 0.000683               //Dynamic viscosity of water (kg/m.s)
+Cp = 4.17*10**3             //Heat capacity of water (J/kg.K)
+k = 0.63                    //Thermal conductivity of water (W/m.K)
+
+//Calculation:
+Tav1 = (T1+T2)/2.0          //Average bulk temperature of water (°C)
+Re = D*V*p/mu               //Reynolds number
+Pr = Cp*mu/k                //Prandtl number
+Tav2 = (Tav1+TC)/2.0        //Fluid's average wall temperature (°C)
+//From Appendix:
+mu_w = 0.000306             //Dynamic viscosity of fluid at wall (kg/m.s)
+//From Table 9.5:
+h = (k/D)*0.027*Re**0.8*Pr**0.33*(mu/mu_w)**0.14    //Heat transfer coefficient for water (W/m^2.K)
+
+//Result:
+printf("The heat transfer coefficient for water is : %.1f  W/m^2.K.",h)
diff --git a/2510/CH9/EX9.13/Ex9_13.sce b/2510/CH9/EX9.13/Ex9_13.sce
new file mode 100755
index 000000000..ec9f8f882
--- /dev/null
+++ b/2510/CH9/EX9.13/Ex9_13.sce
@@ -0,0 +1,11 @@
+//Variable declaration:
+//From example 9.7:
+h = 38.7                    //Average heat transfer coefficient (W/m^2.K)
+L = 1.2                     //Length of plate (m)
+k = 0.025                   //Thermal conductivity of air (W/m)
+
+//Calculation:
+Bi = h*L/k                  //Average Biot number
+
+//Result:
+printf("The average Biot number is : %.0f ",Bi)
diff --git a/2510/CH9/EX9.14/Ex9_14.sce b/2510/CH9/EX9.14/Ex9_14.sce
new file mode 100755
index 000000000..b66fe45a7
--- /dev/null
+++ b/2510/CH9/EX9.14/Ex9_14.sce
@@ -0,0 +1,26 @@
+//Variable declaration:
+k = 60.0                    //Thermal conductivity of rod (W/m.K)
+p = 7850.0                  //Density of rod (kg/m^3)
+Cp = 434.0                  //Heat capacity of rod (J/kg.K)
+h = 140.0                   //Convection heat transfer coefficient (W/m^2.K)
+D = 0.01                    //Diameter of rod (m)
+kf = 0.6                    //Thermal conductivity of fluid (W/m.K)
+L = 2.5                     //Length of rod (m)
+Ts = 250.0                  //Surface temperature of rod (°C)
+Tf = 25.0                   //Fluid temperature (°C)
+
+//Calculation:
+//Case 1:
+a = k/(p*Cp)                //Thermal diffusivity of bare rod (m^2/s)
+//Case 2:
+Nu = h*D/kf                 //Nusselt number
+//Case 3:
+Bi = h*D/k                  //Biot number of bare rod
+//Case 4:
+Q = h*(%pi*D*L)*(Ts-Tf)      //Heat transferred from rod to fluid (W)
+
+//Result:
+printf("1. The thermal diffusivity of the bare rod is : %.2f  x 10^-5 m^2/s.",a/10**-5)
+printf("2. The nusselt number is : %.2f .",Nu)
+printf("3. The Biot number is : %.4f .",Bi)
+printf("4. The heat transferred from the rod to the fluid is : %.0f  W.",Q)
diff --git a/2510/CH9/EX9.2/Ex9_2.sce b/2510/CH9/EX9.2/Ex9_2.sce
new file mode 100755
index 000000000..77b5d1425
--- /dev/null
+++ b/2510/CH9/EX9.2/Ex9_2.sce
@@ -0,0 +1,14 @@
+//Variable declaration:
+//From example 9.1:
+R = 0.0398                      //Theral resistance (°F.h/Btu)
+Btu = 3.412                     //Btu/h in a watt
+C = 1.8                         //Change in degree fahrenheit for a degree change in celsius
+K = 1                           //Change in degree celsius for a unit change in Kelvin
+
+//Calculation:
+Rc = R*Btu/C                    //Thermal resistance in degree cesius per watt (°C/W)
+Rk = Rc/K                       //Thermal resistance in Kelvin per watt (K/W)
+
+//Result:
+printf("The thermal resistance in °C/W is : %.3f  °C/W.",Rc)
+printf("The thermal resistance in K/W is : %.3f  K/W.",Rk)
diff --git a/2510/CH9/EX9.3/Ex9_3.sce b/2510/CH9/EX9.3/Ex9_3.sce
new file mode 100755
index 000000000..9ce1307c9
--- /dev/null
+++ b/2510/CH9/EX9.3/Ex9_3.sce
@@ -0,0 +1,14 @@
+//Variable declaration:
+h = 48.0                    //Convective heat transfer coefficient (Btu/h.ft.°F)
+A = 2*1.5                   //Total heat transfer area (ft^2)
+Ts = 530.0                  //Surface temperature of plate (°F)
+Tm = 105.0                  //Maintained temperature of opposite side of plate (°F)
+kW = 3.4123*10**3           //Units kW in a Btu/h
+
+//Calculation:
+Q = h*A*(Ts-Tm)             //Heat transfer rate in Btu/h (Btu/h)
+Q1 = Q/kW                   //Heat transfer rate in kW (kW)
+
+//Result:
+printf("The heat transfer rate in Btu/h is : %f  Btu/h.",Q)
+printf("The heat transfer rate in kW is : %.2f  kW.",Q1)
diff --git a/2510/CH9/EX9.4/Ex9_4.sce b/2510/CH9/EX9.4/Ex9_4.sce
new file mode 100755
index 000000000..69a540fab
--- /dev/null
+++ b/2510/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,12 @@
+//Variable declaration:
+TS = 10+273                     //Outer surface temperature of wall (K)
+Q = 3000.0                      //Heat transfer rate (W)
+h = 100.0                       //Convection coefficient of air (W/m^2)
+A = 3.0                         //Area of glass window (m^2)
+
+//Calculation:
+TM = TS-Q/(h*A)                 //Bulk temperature of fluid (K)
+
+//Result:
+printf("The bulk temperature of fluid is : %f  K.",TM)
+printf("The bulk temperature of fluid is : %f  °C.",TM-273)
diff --git a/2510/CH9/EX9.5/Ex9_5.sce b/2510/CH9/EX9.5/Ex9_5.sce
new file mode 100755
index 000000000..5103569f0
--- /dev/null
+++ b/2510/CH9/EX9.5/Ex9_5.sce
@@ -0,0 +1,12 @@
+//Variable declaration:
+h = 24.0                    //Plant operating hour per day (h/day)
+d = 350.0                   //Plant operating day per year (day/yr)
+
+//Calculation:
+N = h*d                     //Operating hours per year (h/yr)
+//From example 9.1:
+Q = 8545.0                  //Rate of energy loss (Btu/h)
+Qy = Q*N                    //Steady-state energy loss yearly (Btu/yr)
+
+//Result:
+printf("The yearly steady-state energy loss is : %.2f  x 10^7 Btu/yr.",Qy/10**7)
diff --git a/2510/CH9/EX9.7/Ex9_7.sce b/2510/CH9/EX9.7/Ex9_7.sce
new file mode 100755
index 000000000..5364b82bb
--- /dev/null
+++ b/2510/CH9/EX9.7/Ex9_7.sce
@@ -0,0 +1,21 @@
+
+//Variable declaration:
+x = 0.3                         //Length from the leading age of the plate (m)
+L = 1.2                         //Length of plate (m)
+TS = 58.0                       //Surface temperature of plate (°C)
+Ta = 21.0                       //Temperature of flowing air (°C)
+
+//Calculation:
+hx = 25/x**0.4                  //Local heat transfer coefficient at 0.3m (W/m^2.K) (Part 1)
+syms y                          //Length
+hy = 25/y**0.4                  //hx at the end of the plate (W/m^2.K)
+h = integrate(hy, y,0,L)/L    //Average heat transfer coefficient (W/m^2.K)
+Q = hx*(TS-Ta)                  //Heat flux at 0.3m from leading edge of plate (W/m^2)
+hL = 25/L**0.4                  //Local heat transfer coefficient at plate end (W/m^2.K) (Part 2) 
+r = h/hL                        //Ratio h/hL at the end of the plate
+
+//Result:
+printf("1. The heat flux at 0.3 m from the leading edge of the plate is : %.0f  W/m^2.",Q)
+printf("2. The local heat transfer coefficient at the end of the plate is : %.1f  W/m^2.K.",hL)
+disp("3. The ratio h/hL at the end of plate is : ")
+disp(r)
diff --git a/2510/CH9/EX9.8/Ex9_8.sce b/2510/CH9/EX9.8/Ex9_8.sce
new file mode 100755
index 000000000..b34740b1f
--- /dev/null
+++ b/2510/CH9/EX9.8/Ex9_8.sce
@@ -0,0 +1,15 @@
+//Variable declaration:
+//From example 9.7:
+b = 1.0                     //Width of plate (m)
+L = 1.2                     //Length of plate (m)
+TS = 58.0                   //Surface temperture of plate (°C)
+Ta = 21.0                   //Air flow temperature (°C)
+h = 38.7                    //Average heat transfer coefficient (W/m^2.K)
+
+//Calculation:
+A = b*L                     //Area for heat transfer for the entire plate (m^2)
+Q = h*A*(TS-Ta)             //Rate of heat transfer over the whole length of the plate (W)
+Q = round(Q*10**-1)/10**-1
+
+//Result:
+printf("The rate of heat transfer over the whole length of the plate is : %.1f W.",Q)
diff --git a/2510/CH9/EX9.9/Ex9_9.sce b/2510/CH9/EX9.9/Ex9_9.sce
new file mode 100755
index 000000000..b8de3f493
--- /dev/null
+++ b/2510/CH9/EX9.9/Ex9_9.sce
@@ -0,0 +1,15 @@
+//Variable declaration:
+m = 0.075                           //Mass rate of air flow (kg/s)
+D = 0.225                           //Diameter of tube (m)
+mu = 208*10**-7                     //Dynamic viscosity of fluid (N)
+Pr = 0.71                           //Prandtl number
+k = 0.030                           //Thermal conductivity of air (W/m.K)
+
+//Calculation:
+Re = 4*m/(%pi*D*mu)                  //Reynolds number
+//From equation 9.26:
+Nu = 0.023*(Re**0.8)*(Pr**0.3)      //Nusselt number
+h = (k/D)*Nu                        //Heat transfer coefficient of air (W/m^2.K)
+
+//Result:
+printf("The Heat transfer coefficient of air is : %.2f  W/m^2.K.",h)