initial commit / add all books

5 files changed, 153 insertions, 0 deletions

diff --git a/534/CH9/EX9.1/9_1_Vertical_Plate.sce b/534/CH9/EX9.1/9_1_Vertical_Plate.sce
new file mode 100644
index 000000000..5030e4211
--- /dev/null
+++ b/534/CH9/EX9.1/9_1_Vertical_Plate.sce
@@ -0,0 +1,27 @@
+clear;

+clc;

+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 9.1   Page 569 \n'); //Example 9.1

+// Boundary Layer thickness at trailing edge.

+

+//Operating Conditions

+Ts = 70+273;    //[K] Surface Temperature

+Tsurr = 25+273;    //[K] Surrounding Temperature

+v1 = 0;            //[m/s] Velocity of free air

+v2 = 5;            //[m/s] Velocity of free air

+L = .25;            //[m] length

+

+//Table A.4 Air Properties T = 320 K

+uv = 17.95*10^-6;         //[m^2/s] Kinematic Viscosity

+be = 3.12*10^-3;          //[K^-1]  Tf^-1

+Pr = 269;               // Prandtl number 

+g = 9.81;        //[m^2/s]gravitational constt

+

+Gr = g*be*(Ts-Tsurr)*L^3/uv^2;

+del = 6*L/(Gr/4)^.25;

+printf("\n Boundary Layer thickness at trailing edge for no air stream %.3f m",del);

+

+Re = v2*L/uv;

+printf("\n\n For air stream at 5 m/s As the Reynolds Number is %.2e the free convection boundary layer is Laminar",Re);

+del2 = 5*L/(Re)^.5;

+printf("\n Boundary Layer thickness at trailing edge for air stream at 5 m/s is %.4f m",del2);

+//END
\ No newline at end of file
diff --git a/534/CH9/EX9.2/9_2_Glass_door.sce b/534/CH9/EX9.2/9_2_Glass_door.sce
new file mode 100644
index 000000000..b0fb82b80
--- /dev/null
+++ b/534/CH9/EX9.2/9_2_Glass_door.sce
@@ -0,0 +1,28 @@
+clear;

+clc;

+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 9.2   Page 572 \n'); //Example 9.2

+// Heat transfer by convection between screen and room air.

+

+//Operating Conditions

+Ts = 232+273;    //[K] Surface Temperature

+Tsurr = 23+273;    //[K] Surrounding Temperature

+L = .71;            //[m] length

+w = 1.02;           //[m] Width

+

+//Table A.4 Air Properties T = 400 K

+k = 33.8*10^-3            ;//[W/m.K]

+uv = 26.4*10^-6         ;//[m^2/s] Kinematic Viscosity

+al = 38.3*10^-6        ;//[m^2/s]

+be = 2.5*10^-3          ;//[K^-1]  Tf^-1

+Pr = .69               ;// Prandtl number 

+g = 9.81                ;//[m^2/s] gravitational constt

+

+Ra = g*be*(Ts-Tsurr)/al*L^3/uv;

+printf("\n\n As the Rayleigh Number is %.2e the free convection boundary layer is turbulent",Ra);

+//From equatiom 9.23

+Nu = [.825 + .387*Ra^.16667/[1+(.492/Pr)^(9/16)]^(8/27)]^2;

+h = Nu*k/L;

+q = h*L*w*(Ts-Tsurr);

+

+printf("\n Heat transfer by convection between screen and room air is %i W",q);

+//END
\ No newline at end of file
diff --git a/534/CH9/EX9.3/9_3_Rectangular_Duct.sce b/534/CH9/EX9.3/9_3_Rectangular_Duct.sce
new file mode 100644
index 000000000..587ded94a
--- /dev/null
+++ b/534/CH9/EX9.3/9_3_Rectangular_Duct.sce
@@ -0,0 +1,35 @@
+clear;

+clc;

+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 9.3   Page 577 \n'); //Example 9.3

+// Heat Loss from duct per meter of length

+

+//Operating Conditions

+Ts = 45+273;    //[K] Surface Temperature

+Tsurr = 15+273    ;//[K] Surrounding Temperature

+H = .3            ;//[m] Height 

+w = .75           ;//[m] Width

+

+//Table A.4 Air Properties T = 303 K

+k = 26.5*10^-3            ;//[W/m.K]

+uv = 16.2*10^-6         ;//[m^2/s] Kinematic Viscosity

+al = 22.9*10^-6        ;//[m^2/s] alpha

+be = 3.3*10^-3          ;//[K^-1]  Tf^-1

+Pr = .71               ;// Prandtl number 

+g = 9.81                ;//[m^2/s] gravitational constt

+

+Ra = g*be*(Ts-Tsurr)/al*H^3/uv;    //Length = Height

+//From equatiom 9.27

+Nu = [.68 + .67*Ra^.25/[1+(.492/Pr)^(9/16)]^(4/9)];

+//for Sides

+hs = Nu*k/H;

+

+Ra2 = g*be*(Ts-Tsurr)/al*(w/2)^3/uv;       //Length = w/2

+//For top eq 9.31

+ht = [k/(w/2)]*.15*Ra2^.3334;

+//For bottom Eq 9.32

+hb = [k/(w/2)]*.27*Ra2^.25;

+

+q = (2*hs*H+ht*w+hb*w)*(Ts-Tsurr);

+

+printf("\n Rate of heat loss per unit length of duct is %i W/m",q);

+//END
\ No newline at end of file
diff --git a/534/CH9/EX9.4/9_4_Steam_Pipe.sce b/534/CH9/EX9.4/9_4_Steam_Pipe.sce
new file mode 100644
index 000000000..4653eb1cd
--- /dev/null
+++ b/534/CH9/EX9.4/9_4_Steam_Pipe.sce
@@ -0,0 +1,30 @@
+clear;

+clc;

+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 9.4   Page 580 \n'); //Example 9.4

+// Heat Loss from pipe per meter of length

+

+//Operating Conditions

+Ts = 165+273;    //[K] Surface Temperature

+Tsurr = 23+273;    //[K] Surrounding Temperature

+D = .1            ;//[m] Diameter

+e = .85          ;// emissivity

+stfncnstt=5.67*10^(-8)    ;// [W/m^2.K^4] - Stefan Boltzmann Constant 

+

+//Table A.4 Air Properties T = 303 K

+k = 31.3*10^-3            ;//[W/m.K] Conductivity

+uv = 22.8*10^-6         ;//[m^2/s] Kinematic Viscosity

+al = 32.8*10^-6        ;//[m^2/s] alpha

+be = 2.725*10^-3          ;//[K^-1]  Tf^-1

+Pr = .697               ;// Prandtl number 

+g = 9.81                ;//[m^2/s] gravitational constt

+

+Ra = g*be*(Ts-Tsurr)/al*D^3/uv;    

+//From equatiom 9.34

+Nu = [.60 + .387*Ra^(1/6)/[1+(.559/Pr)^(9/16)]^(8/27)]^2;

+h = Nu*k/D;

+

+qconv = h*%pi*D*(Ts-Tsurr);

+qrad = e*%pi*D*stfncnstt*(Ts^4-Tsurr^4);

+

+printf("\n Rate of heat loss per unit length of pipe is %i W/m",qconv+qrad);

+//END
\ No newline at end of file
diff --git a/534/CH9/EX9.5/9_5_Radiation_Shield.sce b/534/CH9/EX9.5/9_5_Radiation_Shield.sce
new file mode 100644
index 000000000..614664f2b
--- /dev/null
+++ b/534/CH9/EX9.5/9_5_Radiation_Shield.sce
@@ -0,0 +1,33 @@
+clear;

+clc;

+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 9.5   Page 592 \n'); //Example 9.5

+// Heat Loss from pipe per unit of length

+// Heat Loss if air is filled with glass-fiber blanket insulation

+

+//Operating Conditions

+To = 35+273    ;//[K] Shield Temperature

+Ti = 120+273    ;//[K] Tube Temperature

+Di = .1            ;//[m] Diameter inner

+Do = .12            ;//[m] Diameter outer

+L = .01            ;//[m] air gap insulation

+

+//Table A.4 Air Properties T = 350 K

+k = 30*10^-3            ;//[W/m.K] Conductivity

+uv = 20.92*10^-6         ;//[m^2/s] Kinematic Viscosity

+al = 29.9*10^-6        ;//[m^2/s] alpha

+be = 2.85*10^-3          ;//[K^-1]  Tf^-1

+Pr = .7               ;// Prandtl number 

+g = 9.81                ;//[m^2/s] gravitational constt

+//Table A.3 Insulation glass fiber T=300K

+kins = .038               ;//[W/m.K] Conductivity

+

+Lc = 2*[2.303*log10(Do/Di)]^(4/3)/((Di/2)^-(3/5)+(Do/2)^-(3/5))^(5/3);

+Ra = g*be*(Ti-To)/al*Lc^3/uv;    

+keff = .386*k*(Pr/(.861+Pr))^.25*Ra^.25;

+q = 2*%pi*keff*(Ti-To)/(2.303*log10(Do/Di));

+

+//From equatiom 9.58 and 3.27

+qin = q*kins/keff;

+

+printf("\n Heat Loss from pipe per unit of length is %i W/m \n Heat Loss if air is filled with glass-fiber blanket insulation %i W/m",q,qin);

+//END
\ No newline at end of file