initial commit / add all books

15 files changed, 230 insertions, 0 deletions

diff --git a/2510/CH11/EX11.10/Ex11_10.sce b/2510/CH11/EX11.10/Ex11_10.sce
new file mode 100755
index 000000000..234813c06
--- /dev/null
+++ b/2510/CH11/EX11.10/Ex11_10.sce
@@ -0,0 +1,12 @@
+//Variable declaration:
+TH = 140.0+460.0                    //Absolute outside temperature of pipe (ft^2)
+TC = 60.0+460.0                     //Absolute temperature of surrounding atmosphere (ft^2)
+A = 10.0                            //Area of pipe (ft^2)
+E = 0.9                             //Emissivity of pipe
+
+//Calculation:
+Q = E*A*0.173*((TH/100.0)**4-(TC/100.0)**4) //Heat loss due to radiation (Btu/h)
+Q = round(Q*10**-1)/10**-1
+
+//Result:
+printf("The heat loss due to radiation is : %f Btu/h.",Q)
diff --git a/2510/CH11/EX11.11/Ex11_11.sce b/2510/CH11/EX11.11/Ex11_11.sce
new file mode 100755
index 000000000..ded3aa878
--- /dev/null
+++ b/2510/CH11/EX11.11/Ex11_11.sce
@@ -0,0 +1,12 @@
+//Variable declaration:
+//Froma example 11.10:
+Q = 880.0                   //Heat loss due to radiation (Btu/h)
+A = 10.0                    //Area of pipe (ft^2)
+TH = 140.0                  //Absolute outside temperature of pipe (°F)
+TC = 60.0                   //Absolute temperature of surrounding atmosphere (°F)
+
+//Calculation:
+hr = Q/(A*(TH-TC))          //Radiation heat transfer coefficient (Btu/h.ft^2.°F)
+
+//Result:
+printf("The radiation heat transfer coefficient is : %.1f Btu/h.ft^2.°F.",hr)
diff --git a/2510/CH11/EX11.12/Ex11_12.sce b/2510/CH11/EX11.12/Ex11_12.sce
new file mode 100755
index 000000000..d808e5af5
--- /dev/null
+++ b/2510/CH11/EX11.12/Ex11_12.sce
@@ -0,0 +1,35 @@
+//Variable declaration:
+D = 0.0833                          //Diameter of tube (ft)
+L = 2.0                             //Length of tube (ft)
+h = 2.8                             //Heat transfer coefficient (Btu/h.ft^2.°F)
+Ta1 = 1500.0+460.0                  //Temperature of hot air in furnace (°R)
+Ta2 = 1350.0+460.0                  //Temperature of hot air in the furnace brick walls (°R)
+Tt = 600.0+460.0                    //Surface temperature of tube (°R)
+E = 0.6                             //Surface emissivity of tube
+s = 0.1713*10**-8                   //Stefan-Boltzmann constant
+pi = %pi
+
+//Calculation:
+//Case 1:
+A = pi*D*L                          //Area of tube (ft^2)
+Qc = round(h*A*(Ta1-Tt)*10**-1)/10**-1         //Convection heat transfer from air to tube (Btu/h)
+Qr = round(E*s*A*(Ta2**4-Tt**4)*10**-2)/10**-2 //Radiation feat transfer from wall to tube (Btu/h)
+Q = Qr+Qc                           //Total heat transfer (Btu/h)
+//Case 2:
+Qp = Qr/Q*100                       //Radiation percent 
+//Case 3:
+hr = Qr/(A*(Ta2-Tt))                //Radiation heat transfer coefficient (Btu/h.ft^2.°F)
+//Case 4:
+T = Ta2-Tt                          //Temperature difference (°F)
+
+//Result:
+printf("1. The convective heat transferred to the metal tube is : %f Btu/h.",Qc)
+printf("   The radiative heat transferred to the metal tube is : %f Btu/h.",Qr)
+printf("   The total heat transferred to the metal tube is : %f Btu/h .",Q)
+printf("2. The percent of total heat transferred by radiation is : %.1f %%.",Qp)
+printf("3. The radiation heat transfer coefficient is : %.1f Btu/h.ft^2.°F.",hr)
+if (T > 200) then
+    printf("4. The use of the approximation Equation (11.30), hr = 4EsTav^3, is not appropriate.")
+elseif (T < 200) then
+    printf("4. The use of the approximation Equation (11.30), hr = 4EsTav^3, is appropriate.")
+end
diff --git a/2510/CH11/EX11.13/Ex11_13.sce b/2510/CH11/EX11.13/Ex11_13.sce
new file mode 100755
index 000000000..813920638
--- /dev/null
+++ b/2510/CH11/EX11.13/Ex11_13.sce
@@ -0,0 +1,12 @@
+//Variable declaration:
+Q = 5.0                         //Radiation heat transfer (W)
+E = 1.0                         //Emissivity of filament
+s = 5.669*10**-8                //Stefan-Boltzmann constant
+T1 = 900.0+273.0                //Light bulb temperature (K)
+T2 = 150.0+273.0                //Glass bulb temperature (K)
+
+//Calculation:
+A = Q/(E*s*(T1**4-T2**4))       //Surface area of the filament (m^2)
+
+//Result:
+printf("The surface area of the filament is : %.2f cm^2",A*10**4)
diff --git a/2510/CH11/EX11.14/Ex11_14.sce b/2510/CH11/EX11.14/Ex11_14.sce
new file mode 100755
index 000000000..91ab7fd9a
--- /dev/null
+++ b/2510/CH11/EX11.14/Ex11_14.sce
@@ -0,0 +1,25 @@
+//Variable declaration:
+T1 = 127.0+273.0                    //Surface temperature (K)
+T2 = 20.0+273.0                     //Wall temperature (K)
+T3 = 22.0+273.0                     //Air temperature (K)
+s = 5.669*10**-8                    //Stefan-Boltzmann constant
+e = 0.76                            //Surface emissivity of anodized aluminium
+D = 0.06                            //Diameter of %pipe (m)
+L = 100.0                           //Length of %pipe (m)
+h = 15.0                            //%pipe convective heat transfer coefficient (W/m^2.K)
+
+//Calculation:
+Eb = s*T1**4                        //Emissive energy of %pipe (W/m^2)
+E = e*Eb                            //Emissive power from surface of %pipe (W/m^2)
+A = %pi*D*L                          //Surface area of %pipe (m^2)
+Qc = h*A*(T1-T3)                    //Convection heat transfer to air (W)
+Qr = e*s*A*(T1**4-T2**4)            //Radiation heat transfer rate (W)
+Q = Qc+Qr                           //Total heat transfer rate (Btu/h)
+Tav = (T1+T2)/2.0                   //Average temperature (K)
+hr = 4*e*s*Tav**3                   //Radiation heat transfer coefficient (W/m^2.K)
+
+//Result:
+printf("The emissive power from surface of %%pipe is : %.0f W/m^2.",E)
+printf("The convection heat transfer to air is : %.1f kW.",Qc/10**3)
+printf("The radiation heat transfer rate is : %.1f kW",Qr/10**3)
+printf("The radiation heat transfer coefficient is : %.1f W/m^2.K.",hr)
diff --git a/2510/CH11/EX11.15/Ex11_15.sce b/2510/CH11/EX11.15/Ex11_15.sce
new file mode 100755
index 000000000..6d61ab3f2
--- /dev/null
+++ b/2510/CH11/EX11.15/Ex11_15.sce
@@ -0,0 +1,10 @@
+//Variable declaration:
+//From example 11.14:
+Qc = 15.0                       //Convection heat transfer coefficient (W/m^2.K)
+hr = 7.2                        //Radiation heat transfer coefficient (W/m^2.K)
+
+//Calculation:
+X = hr/(Qc+hr)*100.0            //Percent heat transfer by radiation (%)
+
+//Result:
+printf("The percent heat transfer by radiation is : %.1f %%.",X)
diff --git a/2510/CH11/EX11.16/Ex11_16.sce b/2510/CH11/EX11.16/Ex11_16.sce
new file mode 100755
index 000000000..9b7dbbfec
--- /dev/null
+++ b/2510/CH11/EX11.16/Ex11_16.sce
@@ -0,0 +1,15 @@
+//Variable declaration:
+FV = 1.0                            //Correction factor
+//From example 11.9:
+FE = 0.358                          //Emissivity correction factor
+TH = 300.0+460.0                    //Absolute temperature of external surface (°R)
+TC = 75.0+460.0                     //Absolute temperature of duct (°R)
+AH = 0.622                          //Area of pipe (ft^2)
+s = 0.173*10**-8                    //Stefan-Boltzmann constant
+
+//Calculation:
+Q = FV*FE*AH*s*(TH**4-TC**4)        //Heat transfer rate (Btu/h.ft)
+
+//Result:
+printf("The heat transfer rate is : %.2f Btu/h.ft",Q)
+printf("Since, Q obtained in (11.9) is 96.96 Btu/h.ft, the solution does not match with book.")
diff --git a/2510/CH11/EX11.17/Ex11_17.sce b/2510/CH11/EX11.17/Ex11_17.sce
new file mode 100755
index 000000000..46e7f8b55
--- /dev/null
+++ b/2510/CH11/EX11.17/Ex11_17.sce
@@ -0,0 +1,24 @@
+//Variable declaration:
+//From figure 11.2:
+L = 1.0                                 //Space between plates (m)
+X = 0.5                                 //Length of plate (m)
+Y = 2.0                                 //Width of plate (m)
+s = 5.669*10**-8                        //Stefan-Boltzmann constant
+TH = 2000.0+273.0                       //Temperature of hotter plate (K)
+TC = 1000.0+273.0                       //Temperature of colder plate (K)
+Btu = 0.2934*10**-3                     //Btu/h in a KW
+
+//Calculation:
+A = X*Y                                 //Area of plate (m^2)
+Z1 = Y/L                                //Ratio of width with space
+Z2 = X/L                                //Ratio of length with space
+//From figure 11.2:
+FV = 0.18                               //Correction factor
+FE = 1.0                                //Emissivity correction factor
+Q1 = FV*FE*s*A*(TH**4-TC**4)            //Net radiant heat exchange between plates (kW)
+Q2 = Q1/Btu                             //Net radiant heat exchange between plates in Btu/h (Btu/h)
+Q1 = round(Q1*10**-2)/10**-2
+
+//Result:
+printf("The net radiant heat exchange between plates is : %f kW.",Q1)
+printf("The net radiant heat exchange between plates in Btu/h is : %.2f  x 10^8 Btu/h.",Q2/10**8)
diff --git a/2510/CH11/EX11.3/Ex11_3.sce b/2510/CH11/EX11.3/Ex11_3.sce
new file mode 100755
index 000000000..0af8b179b
--- /dev/null
+++ b/2510/CH11/EX11.3/Ex11_3.sce
@@ -0,0 +1,9 @@
+//Variable declaration:
+syms l                                 //Wavelength (mu.m)
+I = 40*exp(-l**2)                       //Intensity of radiation (Btu/h.ft^2.mu.m)
+
+//Calculation:
+E = eval(integrate(I, l,0,%inf))                 //Total emissive power (Btu/h.ft^2)
+
+//Result:
+printf("The total emissive power is : %.1f Btu/h.ft^2.",E)
diff --git a/2510/CH11/EX11.4/Ex11_4.sce b/2510/CH11/EX11.4/Ex11_4.sce
new file mode 100755
index 000000000..a31d594c6
--- /dev/null
+++ b/2510/CH11/EX11.4/Ex11_4.sce
@@ -0,0 +1,14 @@
+//Variable declaration:
+l = 0.25                            //Wavelength (mu.m)
+//From equation 11.4:
+lT = 2884                           //Product of wavelength and absolute temperature (mu.m.°R)
+
+//Calculation:
+T = lT/l                            //Sun's temperature (°R)
+T1 = round(T * 10**-2)/10**-2
+T = T - 460
+T460 = round(T * 10**-3)/10**-3
+
+//Result:
+printf("The Sun s temperature is : %f  °R.",T1)
+printf("The Sun s temperature in fahrenheit scale is : %f °F.",T460)
diff --git a/2510/CH11/EX11.5/Ex11_5.sce b/2510/CH11/EX11.5/Ex11_5.sce
new file mode 100755
index 000000000..d2c3cfd5b
--- /dev/null
+++ b/2510/CH11/EX11.5/Ex11_5.sce
@@ -0,0 +1,11 @@
+//Variable declaration:
+T1 = 1500.0+460.0                   //Absolute temperature 1 (°R)
+T2 = 1000.0+460.0                   //Absolute temperature 2 (°R)
+
+//Calculation:
+X = T1**4/T2**4                     //Ratio of quantity of heat transferred
+x = 100*(T1**4-T2**4)/T2**4         //Percentage increase in heat transfer (%)
+
+//Result:
+printf("The ratio of the quantity/rate of heat transferred is : %.2f .",X)
+printf("The percentage increase in heat transfer is : %.0f %%",x)
diff --git a/2510/CH11/EX11.6/Ex11_6.sce b/2510/CH11/EX11.6/Ex11_6.sce
new file mode 100755
index 000000000..46239f3ac
--- /dev/null
+++ b/2510/CH11/EX11.6/Ex11_6.sce
@@ -0,0 +1,10 @@
+//Variable declaration:
+T1 = 1200.0+460.0                       //Absolute temperature of wall 1 (°R)
+T2 = 800.0+460.0                        //Absolute temperature of wall 2 (°R)
+
+//Calculation:
+//From equation 11.23:
+X = 0.173*((T1/100.0)**4-(T2/100.0)**4) //Heat removed from colder wall (Btu/h.ft^2)
+
+//Result:
+printf("The heat removed from the colder wall to maintain a steady-state is : %.0f Btu/h.ft^2.",X)
diff --git a/2510/CH11/EX11.7/Ex11_7.sce b/2510/CH11/EX11.7/Ex11_7.sce
new file mode 100755
index 000000000..afaca937d
--- /dev/null
+++ b/2510/CH11/EX11.7/Ex11_7.sce
@@ -0,0 +1,13 @@
+//Variable declaration:
+s = 0.173                       //Stefan-Boltzmann constant (Btu/h.ft^2.°R)
+EH = 0.5                        //Energy transferred from hotter body (Btu/h.ft^2)
+EC = 0.75                       //Energy transferred to colder body (Btu/h.ft^2)
+TH = 1660.0                     //Absolute temperature of hotter body (°R)
+TC = 1260.0                     //Absolute temperature of colder body (°R)
+
+//Calculation:
+E = s*((TH/100.0)**4-(TC/100.0)**4)/((1.0/EH)+(1.0/EC)-1.0) //Net energy exchange per unit area (Btu/h.ft^2)
+E = round(E*10**-1)/10**-1
+
+//Result:
+printf("The net energy exchange per unit area is : %f Btu/h.ft^2.",E)
diff --git a/2510/CH11/EX11.8/Ex11_8.sce b/2510/CH11/EX11.8/Ex11_8.sce
new file mode 100755
index 000000000..936f56aaf
--- /dev/null
+++ b/2510/CH11/EX11.8/Ex11_8.sce
@@ -0,0 +1,10 @@
+//Variable declaration:
+//From example 11.6-11.7:
+E1 = 8776.0                     //Energy exchange between black bodies (Btu/h.ft^2)
+E2 = 3760.0                     //Energy exchange between non-black bodies (Btu/h.ft^2)
+
+//Calculation:
+D = (E1-E2)/E1*100              //Percent difference in energy (%)
+
+//Result:
+printf("The percent difference relative to the black body is: %.1f %%.",D)
diff --git a/2510/CH11/EX11.9/Ex11_9.sce b/2510/CH11/EX11.9/Ex11_9.sce
new file mode 100755
index 000000000..504fb2449
--- /dev/null
+++ b/2510/CH11/EX11.9/Ex11_9.sce
@@ -0,0 +1,18 @@
+//Variable declaration:
+s = 0.173*10**-8                        //Stefan-Boltzmann constant (Btu/h.ft^2.°R)
+TH = 300.0+460.0                        //Absolute temperature of external surface (°R)
+TC = 75.0+460.0                         //Absolute temperature of duct (°R)
+//From Table 6.2:
+AH = 0.622                              //External surface area of pipe (ft^2)
+//From Table 11.2:
+EH = 0.44                               //Emissivity of oxidized steel
+AC = 4.0*1.0*1.0                        //External surface area of duct (ft^2)
+EC = 0.23                               //Emissivity of galvanized zinc
+
+//Calculation:
+FE = 1.0/(1.0/EH+((AH/AC)*(1.0/EC-1.0))) //Emissivity correction factor
+Q = FE*AH*s*(TH**4-TC**4)                //Net radiation heat transfer (Btu/h.ft)
+
+//Result:
+printf("The net radiation heat transfer is : %.2f Btu/h.ft^2.",Q)
+printf("There is a calculation error in book.")