4 files changed, 135 insertions, 0 deletions
diff --git a/534/CH4/EX4.1/4_1_Eccentric_Wire.sce b/534/CH4/EX4.1/4_1_Eccentric_Wire.sce
new file mode 100644
index 000000000..97ff4e5de
--- /dev/null
+++ b/534/CH4/EX4.1/4_1_Eccentric_Wire.sce
@@ -0,0 +1,21 @@
+clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 4.1   Page 211 \n'); //Example 4.1
+// Thermal resistance of wire coating associated with peripheral variations in coating thickness
+
+d = .005;        //[m] Diameter of wire
+k = .35;         //[W/m.K] Thermal Conductivity
+h = 15;          //[W/m^2.K] Total coeff with Convection n Radiation
+
+rcr = k/h;        // [m] critical insulation radius
+tcr = rcr - d/2;  // [m] critical insulation Thickness
+
+Rtcond = 2.302*log10(rcr/(d/2))/(2*%pi*k);      //[K/W] Thermal resistance 
+
+//Using Table 4.1 Case 7
+z = .5*tcr;
+D=2*rcr;
+Rtcond2D = (acosh((D^2 + d^2 - 4*z^2)/(2*D*d)))/(2*%pi*k);
+
+printf("\n\n The reduction in thermal resistance of the insulation is %.2f K/W ", Rtcond-Rtcond2D);
+//END
+\ No newline at end of file
diff --git a/534/CH4/EX4.2/4_2_Theoretical_Problem.sce b/534/CH4/EX4.2/4_2_Theoretical_Problem.sce
new file mode 100644
index 000000000..21ba9886b
--- /dev/null
+++ b/534/CH4/EX4.2/4_2_Theoretical_Problem.sce
@@ -0,0 +1,8 @@
+clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 4.2   Page 218 \n')// Example 4.2
+//Theoretical Problem
+
+printf('\n The given example is theoretical and does not involve any numerical computation')
+
+//End
diff --git a/534/CH4/EX4.3/4_3_Column_Matrix.sce b/534/CH4/EX4.3/4_3_Column_Matrix.sce
new file mode 100644
index 000000000..064caa7f9
--- /dev/null
+++ b/534/CH4/EX4.3/4_3_Column_Matrix.sce
@@ -0,0 +1,33 @@
+clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 4.3   Page 224 \n'); //Example 4.2
+// Temperature Distribution and Heat rate per unit length
+
+Ts = 500;        //[K] Temp of surface
+Tsurr = 300;     //[K] Temp of surrounding Air
+h = 10;          //[W/m^2.K] Heat Convection soefficient
+//Support Column
+delx = .25;      //[m]
+dely = .25;      //[m]
+k = 1;         //[W/m.K] From Table A.3, Fireclay Brick at T = 478K
+
+//Applying Eqn 4.42 and 4.48
+A = [-4 1 1 0 0 0 0 0; 
+     2 -4 0 1 0 0 0 0;
+     1 0 -4 1 1 0 0 0;
+     0 1 2 -4 0 1 0 0;
+     0 0 1 0 -4 1 1 0;
+     0 0 0 1 2 -4 0 1;
+     0 0 0 0 2 0 -9 1;
+     0 0 0 0 0 2 2 -9 ];
+     
+C = [-1000; -500; -500; 0; -500; 0; -2000; -1500 ];
+
+T = inv(A)*C;
+
+printf("\n Temp Distribution = ");
+printf("\n    %.2f K ", T);
+
+q = 2*h*[(delx/2)*(Ts-Tsurr)+delx*(T(7)-Tsurr)+delx*(T(8)-Tsurr)/2];
+printf("\n\n Heat rate from column to the airstream %.1f W/m ", q);
+//END
+\ No newline at end of file
diff --git a/534/CH4/EX4.4/4_4_Turbine_Matrix.sce b/534/CH4/EX4.4/4_4_Turbine_Matrix.sce
new file mode 100644
index 000000000..490bed5e7
--- /dev/null
+++ b/534/CH4/EX4.4/4_4_Turbine_Matrix.sce
@@ -0,0 +1,73 @@
+clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 4.4   Page 230 \n'); //Example 4.4
+// Temperature Field and Rate of Heat Transfer
+
+//Operating Conditions
+
+ho = 1000;          //[W/m^2.K] Heat Convection coefficient
+hi = 200;          //[W/m^2.K] Heat Convection coefficient
+Ti = 400;        //[K] Temp of Air
+Tg = 1700;     //[K] Temp of Gas
+h = 10 ;         //[W/m^2.K] Heat Convection coefficient
+
+A = 2*6*10^-6 ;   //[m^2] Cross section of each Channel
+x = .004 ;     //[m] Spacing between joints
+t = .006;      //[m] Thickness
+k = 25;         //[W/m.K] Thermal Conductivity of Blade
+delx = .001 ;     //[m]
+dely = .001 ;     //[m]
+
+//Applying Eqn 4.42 and 4.48
+A = [-(2+ho*delx/k) 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0; 
+     1 -2*(2+ho*delx/k) 1 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0;
+     0 1 -2*(2+ho*delx/k) 1 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 0;
+     0 0 1 -2*(2+ho*delx/k) 1 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0;
+     0 0 0 1 -2*(2+ho*delx/k) 1 0 0 0 0 2 0 0 0 0 0 0 0 0 0 0;
+     0 0 0 0 1 -(2+ho*delx/k) 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0;
+     1 0 0 0 0 0 -4 2 0 0 0 0 1 0 0 0 0 0 0 0 0;
+     0 1 0 0 0 0 1 -4 1 0 0 0 0 1 0 0 0 0 0 0 0;
+     0 0 1 0 0 0 0 1 -4 1 0 0 0 0 1 0 0 0 0 0 0;
+     0 0 0 1 0 0 0 0 1 -4 1 0 0 0 0 1 0 0 0 0 0;
+     0 0 0 0 1 0 0 0 0 1 -4 1 0 0 0 0 1 0 0 0 0;
+     0 0 0 0 0 1 0 0 0 0 2 -4 0 0 0 0 0 1 0 0 0;
+     0 0 0 0 0 0 1 0 0 0 0 0 -4 2 0 0 0 0 1 0 0;
+     0 0 0 0 0 0 0 1 0 0 0 0 1 -4 1 0 0 0 0 1 0;
+     0 0 0 0 0 0 0 0 2 0 0 0 0 2 -2*(3+hi*delx/k) 1 0 0 0 0 1;
+     0 0 0 0 0 0 0 0 0 2 0 0 0 0 1 -2*(2+hi*delx/k) 1 0 0 0 0;
+     0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 1 -2*(2+hi*delx/k) 1 0 0 0;
+     0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 -(2+hi*delx/k) 0 0 0;
+     0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 -2 1 0;
+     0 0 0 0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 1 -4 1;
+     0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 -(2+hi*delx/k)];
+     
+C = [-ho*delx*Tg/k; 
+     -2*ho*delx*Tg/k;
+     -2*ho*delx*Tg/k;
+     -2*ho*delx*Tg/k;
+     -2*ho*delx*Tg/k;
+     -ho*delx*Tg/k;
+     0;
+     0;
+     0;
+     0;
+     0;
+     0;
+     0;
+     0;
+     -2*hi*delx*Ti/k;
+     -2*hi*delx*Ti/k;
+     -2*hi*delx*Ti/k;
+     -hi*delx*Ti/k;
+     0;
+     0;
+     -hi*delx*Ti/k];
+
+T = inv(A)*C;
+
+printf("\n Temp Distribution = ");
+printf("\n    %.1f K ", T);
+
+q = 4*ho*[(delx/2)*(Tg-T(1))+delx*(Tg-T(2))+delx*(Tg-T(3))+ delx*(Tg-T(4))+delx*(Tg-T(5))+delx*(Tg-T(6))/2];
+printf("\n\n Heat rate Transfer %.1f W/m ", q);
+//END
+\ No newline at end of file