summaryrefslogtreecommitdiff
path: root/2510/CH9
diff options
context:
space:
mode:
Diffstat (limited to '2510/CH9')
-rwxr-xr-x2510/CH9/EX9.1/Ex9_1.sce15
-rwxr-xr-x2510/CH9/EX9.10/Ex9_10.sce19
-rwxr-xr-x2510/CH9/EX9.11/Ex9_11.sce21
-rwxr-xr-x2510/CH9/EX9.12/Ex9_12.sce24
-rwxr-xr-x2510/CH9/EX9.13/Ex9_13.sce11
-rwxr-xr-x2510/CH9/EX9.14/Ex9_14.sce26
-rwxr-xr-x2510/CH9/EX9.2/Ex9_2.sce14
-rwxr-xr-x2510/CH9/EX9.3/Ex9_3.sce14
-rwxr-xr-x2510/CH9/EX9.4/Ex9_4.sce12
-rwxr-xr-x2510/CH9/EX9.5/Ex9_5.sce12
-rwxr-xr-x2510/CH9/EX9.7/Ex9_7.sce21
-rwxr-xr-x2510/CH9/EX9.8/Ex9_8.sce15
-rwxr-xr-x2510/CH9/EX9.9/Ex9_9.sce15
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)