initial commit / add all books

14 files changed, 486 insertions, 0 deletions

diff --git a/1244/CH1/EX1.1/Example11.sce b/1244/CH1/EX1.1/Example11.sce
new file mode 100755
index 000000000..336809ca4
--- /dev/null
+++ b/1244/CH1/EX1.1/Example11.sce
@@ -0,0 +1,37 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.1 ")

+

+//Temperature Inside in F

+Ti = 55;

+//Temperature outside in F

+To = 45;

+//Thickness of the wall in ft

+t = 1;

+//Heat loss through the wall in Btu/h-ft2

+q = 3.4;

+

+//Converting Btu/h-ft2 to W/m2

+disp("Heat loss through the wall in W/m2 is")

+//Heat loss through the wall in W/m2 

+qdash = (q*0.2931)/0.0929

+

+//Heat loss for a 100ft2 surface over a 24-h period

+disp("Heat loss for a 100ft2 surface over a 24-h period in Btu is")

+//Heat loss for a 100ft2 surface over a 24-h period in Btu 

+Q = (q*100)*24

+

+//Q in SI units i.e. kWh

+Q = (Q*0.2931)/1000;

+

+//At price of 10c/kWh, the total price shall be

+disp("So, the total price in c are")

+//Total price in c

+Price = 10*Q

diff --git a/1244/CH1/EX1.10/Example110.sce b/1244/CH1/EX1.10/Example110.sce
new file mode 100755
index 000000000..6e676efbb
--- /dev/null
+++ b/1244/CH1/EX1.10/Example110.sce
@@ -0,0 +1,33 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.10 ")

+

+//diameter of pipe in m

+d = 0.5;

+//Epsilon is given as

+epsilon = 0.9;

+//sigma(constant) in SI units is

+sigma = 0.0000000567;

+//Temperatures in K are given as

+T1 = 500;

+T2 = 300;

+

+//Radiation heat transfer coefficient in W/m2K

+hr = ((sigma*epsilon)*(T1*T1+T2*T2))*(T1+T2);

+

+//Convection heat transfer coefficient in W/m2K

+hc = 20;

+

+//total heat transfer coefficient in W/m2K

+h = hc+hr;

+

+disp("Rate of heat loss per meter in W/m is")

+//Rate of heat loss per meter in W/m

+q = ((%pi*d)*h)*(T1-T2)

diff --git a/1244/CH1/EX1.11/Example111.sce b/1244/CH1/EX1.11/Example111.sce
new file mode 100755
index 000000000..21b5d1c81
--- /dev/null
+++ b/1244/CH1/EX1.11/Example111.sce
@@ -0,0 +1,26 @@
+

+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.11 ")

+

+//Hot-gas temperature in K

+Tgh = 1300;

+//Heat transfer coefficient on hot side in W/m2K

+h1 = 200;

+//Heat transfer coefficient on cold side in W/m2K

+h3 = 400;

+//Coolant temperature in K

+Tgc = 300;

+//Max temp. in C

+Tsg = 800;

+//Maximum permissible unit thermal resistance per square meter of the metal wall in K/W

+R2 = (Tgh-Tgc)/((Tgh-Tsg)/(1/h1))-1/h1-1/h3;

+disp("Maximum permissible unit thermal resistance per square meter of the metal wall in m2.K/W is")

+R2

diff --git a/1244/CH1/EX1.12/Example112.sce b/1244/CH1/EX1.12/Example112.sce
new file mode 100755
index 000000000..faa910d8a
--- /dev/null
+++ b/1244/CH1/EX1.12/Example112.sce
@@ -0,0 +1,44 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.12 ")

+

+// total length of metal sheet in m

+L = 0.625/39.4;

+//  we estimate the thermal conductivity of the metal sheets to be approximately 43 W/m K

+k = 43;

+// therefore the resistance in K/W offered by metal sheey

+R = L/k;

+

+//heat loss in W/m2 is given as

+q = 1200;

+// overall heat transfer coefficient between the gas and the door is given

+// in W/m2K

+U = 20;

+//The temperature drop between the gas and the interior surface of the door at the specified heat flux is

+deltaT1 = q/U;

+//Hence, the temperature of the Inconel will be in degree C

+T = 1200-deltaT1;

+

+//The heat transfer coefficient between the outer surface of the door and

+//the surroundings at 20°C in W/m2K

+h = 5;

+//The temperature drop at the outer surface in degree C is

+deltaT2 = q/h;

+//Selecting milled alumina-silica chips as insulator (Fig 1.31 on page 48)

+

+// Hence, temperature difference across the insulation is

+deltaT3 = T-deltaT1-deltaT2;

+

+//thermal conductivity for milled alumina-silica chips in W/mK is

+k = 0.27;

+

+disp("The insulation thickness in m is")

+//The insulation thickness in m

+L = (k*deltaT3)/q

diff --git a/1244/CH1/EX1.13/Example113.sce b/1244/CH1/EX1.13/Example113.sce
new file mode 100755
index 000000000..87d4f8f44
--- /dev/null
+++ b/1244/CH1/EX1.13/Example113.sce
@@ -0,0 +1,32 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.13 ")

+

+//Temperature of air in degree K

+Tair = 300;

+//Heat transfer coefficient in W/m2K

+h = 10;

+

+disp("Part a")

+//Radiation solar flux in W/m2

+q = 500;

+//Ambient temperature in K

+Tsurr = 50;

+

+disp("Solving energy balance equaiton by trial and error for the roof temperature, we get temp. in degree K")

+//Room temperature in degree K

+Troof = 303

+

+disp("Part b")

+

+//No heat flux, energy balance equaiton is modified

+disp("Room temperature in degree K")

+//Room temperature in degree K

+Troof = 270

diff --git a/1244/CH1/EX1.14/Example114.sce b/1244/CH1/EX1.14/Example114.sce
new file mode 100755
index 000000000..742a5b17c
--- /dev/null
+++ b/1244/CH1/EX1.14/Example114.sce
@@ -0,0 +1,12 @@
+
+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.14 ")

+

+disp("The given example is theoretical and does not involve any numerical computation")

diff --git a/1244/CH1/EX1.2/Example12.sce b/1244/CH1/EX1.2/Example12.sce
new file mode 100755
index 000000000..56cbb136f
--- /dev/null
+++ b/1244/CH1/EX1.2/Example12.sce
@@ -0,0 +1,39 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.2 ")

+

+//Thermal conductivity of window glass in W/m-K

+k = 0.81;

+//Height of the glass in m

+h = 1;

+//Width of the glass in m

+w = 0.5;

+//Thickness of the glass in m

+t = 0.005;

+//Outside temperature in C

+T2 = 24;

+//Inside temperature in C

+T1 = 24.5;

+

+//Assume that steady state exists and that the temperature is uniform over the inner and outer surfaces

+

+//Cross sectional area in m2

+A = h*w;

+

+disp("Thermal resistance to conduction in K/W is")

+//Thermal resistance to conduction in K/W

+R = t/(k*A)

+

+//The rate of heat loss from the interior to the exterior surface is

+//obtained by dividing temperature difference from the thermal resistence

+

+disp("Heat loss in W from the window glass is")

+//Heat loss in W

+q = (T1-T2)/R

diff --git a/1244/CH1/EX1.3/Example13.sce b/1244/CH1/EX1.3/Example13.sce
new file mode 100755
index 000000000..80652fd93
--- /dev/null
+++ b/1244/CH1/EX1.3/Example13.sce
@@ -0,0 +1,26 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.3 ")

+

+//Area of room in m2 is given as

+A = 20*20;

+//Air temperature in C

+Tair = -3;

+//Roof temperature in C

+Troof = 27;

+//Heat transfer coefficient in W/m2-K

+h = 10;

+

+//Assume that steady state exists and the direction of heat flow is from the

+//roof to the air i.e higher to lower temperature (as it should be).

+

+disp(" The rate of heat transfer by convection from the roof to the air in W")

+//The rate of heat transfer by convection from the roof to the air in W

+q = (h*A)*(Troof-Tair)

diff --git a/1244/CH1/EX1.4/Example14.sce b/1244/CH1/EX1.4/Example14.sce
new file mode 100755
index 000000000..cc26d2a3d
--- /dev/null
+++ b/1244/CH1/EX1.4/Example14.sce
@@ -0,0 +1,36 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.4 ")

+

+//Diameter of rod in m

+d = 0.02;

+// Emissivity and temperautre of rod in K

+epsilon = 0.9;

+T1 = 1000;

+//Temperature of walls of furnace

+T2 = 800;

+

+//Assuming steady state has been reached.

+//Since the walls of the furnace completely enclose the heating rod, all the radiant energy emitted by the surface of the rod is intercepted by the furnace walls

+

+//From eq. 1.17, net heat loss can be given

+

+disp("Net heat loss per unit length considering 1m length in W")

+//Area in m2

+A = (%pi*d)*1;

+//Constant sigma in W/m2-K4

+sigma = 0.0000000567;

+//Net heat loss per unit length considering 1m length in W

+q = ((A*sigma)*epsilon)*(T1^4-T2^4)

+

+//From eq. 1.21 radiation heat transfer coefficient in W/m2-K is

+disp("Radiation heat transfer coefficient in W/m2-K is")

+//Radiation heat transfer coefficient in W/m2-K 

+hr = ((epsilon*sigma)*(T1^4-T2^4))/(T1-T2)

diff --git a/1244/CH1/EX1.5/Example15.sce b/1244/CH1/EX1.5/Example15.sce
new file mode 100755
index 000000000..e9a4a56da
--- /dev/null
+++ b/1244/CH1/EX1.5/Example15.sce
@@ -0,0 +1,31 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.5 ")

+

+//Thickness of inside steel in m and thermal conductivity in W/m-k

+t1 = 0.005;

+k1 = 40;

+//Thickness of outside brick in m and thermal conductivity in W/m-k

+t2 = 0.1;

+k2 = 2.5;

+

+//Inside temperature in C

+T1 = 900;

+//Outside temperature in C

+To = 460;

+

+//Assuming the condition of steady state and using Eq. 1.24

+disp("The rate of heat loss per unit area  in W/m2 is")

+//The rate of heat loss per unit area  in W/m2 

+qk = (T1-To)/(t1/k1+t2/k2)

+

+disp("Temperature at the interface in K is given as")

+//Temperature at the interface in K

+T2 = T1-(qk*t1)/k1

diff --git a/1244/CH1/EX1.6/Example16.sce b/1244/CH1/EX1.6/Example16.sce
new file mode 100755
index 000000000..3dd0aae7c
--- /dev/null
+++ b/1244/CH1/EX1.6/Example16.sce
@@ -0,0 +1,33 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.6 ")

+

+//Thermal conductivity of aluminium in W/m-K

+k = 240;

+//Thickness of each plate in m

+t = 0.01;

+//Temperature at the surfaces of plates in C is given as

+Ts1 = 395;

+Ts3 = 405;

+//From Table 1.6 the contact resistance at the interface in K/W is

+R2 = 0.000275;

+//Thermal resistance of the plates in K/W is

+R1 = t/k;

+R3 = t/k;

+

+disp("Heat flux in W/m2-K is")

+//Heat flux in W/m2-K

+q = (Ts3-Ts1)/(R1+R2+R3)

+

+//Since the temperature drop in each section of this one-dimensional system is propor-tional to the resistance.

+

+disp("Temperature drop due to contact resistance in degree C is")

+//Temperature drop due to contact resistance in degree C

+deltaT = (R2/(R1+R2+R3))*(Ts3-Ts1)

diff --git a/1244/CH1/EX1.7/Example17.sce b/1244/CH1/EX1.7/Example17.sce
new file mode 100755
index 000000000..2f9dfd0f6
--- /dev/null
+++ b/1244/CH1/EX1.7/Example17.sce
@@ -0,0 +1,51 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.7 ")

+

+//Because of symmetry, we need to calculate for only one half of the system

+

+//Thickness of firebrick in inches

+L1 = 1;

+//Thermal conductivity of firebrick in Btu/h-ft-F

+kb = 1;

+//Thickness of steel plate in inches

+L3 = 1/4;

+//Thermal conductivity of steel in Btu/h-ft-F

+ks = 30;

+//Average height of asperities in inches is given as

+L2 = 1/32;

+//Temperature difference between the steel plates in F is

+deltaT = 600;

+

+

+//The thermal resistance of the steel plate is, on the basis of a unit wall area, equal to

+R3 = L3/(12*ks);//12 is added to convert ft to in

+

+//The thermal resistance of the brick asperities is, on the basis of a unit wall area, equal to

+R4 = L2/((0.3*12)*kb);//Considering the 30 percent area

+

+//At temperature of 300F, thermal conductivity of air in Btu/h-ft-F is

+ka = 0.02;

+

+// Thermal resistance of the air trapped between the asperities, is, on the basis of a unit area, equal to

+R5 = L2/((0.7*12)*ka);//Considering the other 70 percent area

+

+//Since R4 and R5 are in parallel, so there combined resistance is

+R2 = (R4*R5)/(R4+R5);

+

+//The thermal resistance of half of the solid brick is

+R1 = L1/(12*kb);

+

+//The overall unit conductance for half the composite wall in Btu/h-ft2-F is then

+kk = 0.5/(R1+R2+R3);

+

+disp("The rate of heat flow per unit area in Btu/h-ft2 is")

+//The rate of heat flow per unit area in Btu/h-ft2

+q = kk*deltaT

diff --git a/1244/CH1/EX1.8/Example18.sce b/1244/CH1/EX1.8/Example18.sce
new file mode 100755
index 000000000..573e393cd
--- /dev/null
+++ b/1244/CH1/EX1.8/Example18.sce
@@ -0,0 +1,47 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.8 ")

+

+//Length for heat transfer for stainless steel in m

+Lss = 0.1;

+

+//Area for heat transfer for stainless steel in m2

+A = 0.01;

+

+//Thermal conductivity for stainless steel in W/m-K

+kss = 144;

+

+//Length for heat transfer for Duralumin in m

+La1 = 0.02;

+

+//Area for heat transfer for Duralumin in m2

+A = 0.01;

+

+//Thermal conductivity for Duralumin in W/m-K

+ka1 = 164;

+

+//Resistance in case of steel  in K/W

+Rk1 = Lss/(A*kss);

+

+//Resistance in case of Duralumin  in K/W

+Rk2 = La1/(A*ka1);

+

+//From Fig. 1.20, contact resistance in K/W

+Ri = 0.05;

+

+//Total resistance to heat transfer in K/W

+Rtotal = Rk1+Rk2+Ri;

+

+//Temperature diff. is given in K

+deltaT = 40;

+

+disp("Maximum allowable rate of heat dissipation in W is")

+//Maximum allowable rate of heat dissipation in W

+q = deltaT/Rtotal

diff --git a/1244/CH1/EX1.9/Example19.sce b/1244/CH1/EX1.9/Example19.sce
new file mode 100755
index 000000000..be6066e01
--- /dev/null
+++ b/1244/CH1/EX1.9/Example19.sce
@@ -0,0 +1,39 @@
+

+

+// Display mode

+mode(0);

+

+// Display warning for floating point exception

+ieee(1);

+

+clc;

+disp("Principles of Heat Transfer, 7th Ed. Frank Kreith et. al Chapter - 1 Example # 1.9 ")

+

+//Cross sectional area in m2

+A = 1;

+//Heat transfer coefficient on hot side in W/m2-K

+hchot = 10;

+//Heat transfer coefficient on cold side in W/m2-K

+hccold = 40;

+

+//Length for heat transfer in m

+L = 0.1;

+//Thermal conductivity in W/m-K

+k = 0.7;

+

+//Resistances in K/w

+R1 = 1/(hchot*A);

+R2 = L/(k*A);

+R3 = 1/(hccold*A);

+

+//Total resistance

+Rtotal = R1+R2+R3;

+

+//Temperature on hot side in K

+T1 = 330;

+//Temperature on cold side in K

+T2 = 270;

+

+disp("Rate of heat transfer per unit area in W is")

+//Rate of heat transfer per unit area in W

+q = (T1-T2)/(R1+R2+R3)