diff options
Diffstat (limited to '534/CH1')
-rw-r--r-- | 534/CH1/EX1.1/1_1_Wall_Heat_Loss.sce | 20 | ||||
-rw-r--r-- | 534/CH1/EX1.2/1_2_Emissive_Power_Irradiation.sce | 26 | ||||
-rw-r--r-- | 534/CH1/EX1.3/1_3_Theoretical_Problem.sce | 8 | ||||
-rw-r--r-- | 534/CH1/EX1.4/1_4_Coolant_Fuid_Velocity.sce | 30 | ||||
-rw-r--r-- | 534/CH1/EX1.5/1_5_Theoretical_Problem.sce | 8 | ||||
-rw-r--r-- | 534/CH1/EX1.6/1_6_Human_Body_Heat_Loss.sce | 48 | ||||
-rw-r--r-- | 534/CH1/EX1.7/1_7_Cure_Temperature.sce | 58 | ||||
-rw-r--r-- | 534/CH1/EX1.8/1_8_Theoretical_Problem.sce | 8 |
8 files changed, 206 insertions, 0 deletions
diff --git a/534/CH1/EX1.1/1_1_Wall_Heat_Loss.sce b/534/CH1/EX1.1/1_1_Wall_Heat_Loss.sce new file mode 100644 index 000000000..e92a71d6a --- /dev/null +++ b/534/CH1/EX1.1/1_1_Wall_Heat_Loss.sce @@ -0,0 +1,20 @@ +clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 1.1 Page 5 ')//Example 1.1
+// Find Wall Heat Loss - Problem of Pure Conduction Unidimensional Heat
+
+L=.15; //[m] - Thickness of conducting wall
+delT = 1400 - 1150; //[K] - Temperature Difference across the Wall
+A=.5*1.2; //[m^2] - Cross sectional Area of wall = H*W
+k=1.7; //[W/m.k] - Thermal Conductivity of Wall Material
+
+//Using Fourier's Law eq 1.2
+Q = k*delT/L; //[W/m^2] - Heat Flux
+
+q = A*Q; //[W] - Rate of Heat Transfer
+
+printf("\n \n Heat Loss through the Wall = %.2f W",q);
+//END
+
+
+
diff --git a/534/CH1/EX1.2/1_2_Emissive_Power_Irradiation.sce b/534/CH1/EX1.2/1_2_Emissive_Power_Irradiation.sce new file mode 100644 index 000000000..20610d7e4 --- /dev/null +++ b/534/CH1/EX1.2/1_2_Emissive_Power_Irradiation.sce @@ -0,0 +1,26 @@ +clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 1.2 Page 11 \n')// Example 1.2
+// Find a) Emissive Power & Irradiation b)Total Heat Loss per unit length
+
+d=.07; //[m] - Outside Diameter of Pipe
+Ts = 200+273.15; //[K] - Surface Temperature of Steam
+Tsurr = 25+273.15; //[K] - Temperature outside the pipe
+e=.8; // Emissivity of Surface
+h=15; //[W/m^2.k] - Thermal Convectivity from surface to air
+stfncnstt=5.67*10^(-8); // [W/m^2.K^4] - Stefan Boltzmann Constant
+//Using Eq 1.5
+E = e*stfncnstt*Ts^4; //[W/m^2] - Emissive Power
+G = stfncnstt*Tsurr^4; //[W/m^2] - Irradiation falling on surface
+
+printf("\n (a) Surface Emissive Power = %.2f W/m^2",E);
+printf("\n Irradiation Falling on Surface = %.2f W/m^2",G);
+
+//Using Eq 1.10 Total Rate of Heat Transfer Q = Q by convection + Q by radiation
+q = h*(%pi*d)*(Ts-Tsurr)+e*(%pi*d)*stfncnstt*(Ts^4-Tsurr^4); //[W]
+
+printf("\n\n (b) Total Heat Loss per unit Length of Pipe= %.2f W",q);
+//END
+
+
+
diff --git a/534/CH1/EX1.3/1_3_Theoretical_Problem.sce b/534/CH1/EX1.3/1_3_Theoretical_Problem.sce new file mode 100644 index 000000000..64bc95f58 --- /dev/null +++ b/534/CH1/EX1.3/1_3_Theoretical_Problem.sce @@ -0,0 +1,8 @@ +clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 1.3 Page 18 \n')// Example 1.3
+//Theoretical Problem
+
+printf('\n The given example is theoretical and does not involve any numerical computation')
+
+//End
diff --git a/534/CH1/EX1.4/1_4_Coolant_Fuid_Velocity.sce b/534/CH1/EX1.4/1_4_Coolant_Fuid_Velocity.sce new file mode 100644 index 000000000..814df67ef --- /dev/null +++ b/534/CH1/EX1.4/1_4_Coolant_Fuid_Velocity.sce @@ -0,0 +1,30 @@ +clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 1.4 Page 20 \n')// Example 1.4
+// Find Velocity of Coolant Fluid
+
+Ts = 56.4+273.15; //[K] - Surface Temperature of Steam
+Tsurr = 25+273.15; //[K] - Temperature of Surroundings
+e=.88; // Emissivity of Surface
+
+//As h=(10.9*V^.8)[W/m^2.k] - Thermal Convectivity from surface to air
+stfncnstt=5.67*10^(-8); // [W/m^2.K^4] - Stefan Boltzmann Constant
+
+A=2*.05*.05; // [m^2] Area for Heat transfer i.e. both surfaces
+
+E = 11.25; //[W] Net heat to be removed by cooling air
+
+Qrad = e*stfncnstt*A*(Ts^4-Tsurr^4);
+
+//Using Eq 1.10 Total Rate of Heat Transfer Q = Q by convection + Q by radiation
+Qconv = E - Qrad;//[W]
+
+//As Qconv = h*A*(Ts-Tsurr) & h=10.9 Ws^(.8)/m^(-.8)K.V^(.8)
+
+V = [Qconv/(10.9*A*(Ts-Tsurr))]^(1/0.8);
+
+printf("\n\n Velocity of Cooling Air flowing= %.2f m/s",V);
+//END
+
+
+
diff --git a/534/CH1/EX1.5/1_5_Theoretical_Problem.sce b/534/CH1/EX1.5/1_5_Theoretical_Problem.sce new file mode 100644 index 000000000..e7156c2eb --- /dev/null +++ b/534/CH1/EX1.5/1_5_Theoretical_Problem.sce @@ -0,0 +1,8 @@ +clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 1.5 Page 23 \n')// Example 1.5
+//Theoretical Problem
+
+printf('\n The given example is theoretical and does not involve any numerical computation')
+
+//End
diff --git a/534/CH1/EX1.6/1_6_Human_Body_Heat_Loss.sce b/534/CH1/EX1.6/1_6_Human_Body_Heat_Loss.sce new file mode 100644 index 000000000..9c8c40257 --- /dev/null +++ b/534/CH1/EX1.6/1_6_Human_Body_Heat_Loss.sce @@ -0,0 +1,48 @@ +clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 1.6 Page 26 ')// Example 1.6
+// Find Skin Temperature & Heat loss rate
+
+A=1.8; // [m^2] Area for Heat transfer i.e. both surfaces
+Ti = 35+273; //[K] - Inside Surface Temperature of Body
+Tsurr = 297; //[K] - Temperature of surrounding
+Tf = 297; //[K] - Temperature of Fluid Flow
+e=.95; // Emissivity of Surface
+L=.003; //[m] - Thickness of Skin
+k=.3; // Effective Thermal Conductivity
+h=2; //[W/m^2.k] - Natural Thermal Convectivity from body to air
+stfncnstt=5.67*10^(-8); // [W/m^2.K^4] - Stefan Boltzmann Constant
+//Using Eq 1.5
+
+Tsa=305; //[K] Body Temperature Assumed
+
+i=-1;
+while(i==-1)
+ hr = e*stfncnstt*(Tsa+Tsurr)*(Tsa^2+Tsurr^2); //[W/m^2.K] - Radiative Heat transfer Coeff on assumption
+
+ //Using Eq 1.8 & Eq 1.9 k(Ti-Ts)/L = h(Ts - Tf) + hr(Ts - Tsurr)
+Ts = (k*Ti/L + (h+hr)*Tf)/(k/L +(h+hr));
+ c=abs(Ts-Tsa);
+ if(c<=0.0001)
+ i=1;
+ break;
+ end
+ Tsa=Ts;
+end
+
+q = k*A*(Ti-Ts)/L; //[W]
+
+printf("\n\n (I) In presence of Air")
+printf("\n (a) Temperature of Skin = %.2f K",Ts);
+printf("\n (b) Total Heat Loss = %.2f W",q);
+
+//When person is in Water
+h = 200; //[W/m^2.k] - Thermal Convectivity from body to water
+hr = 0; // As Water is Opaque for Thermal Radiation
+Ts = (k*Ti/L + (h+hr)*Tf)/(k/L +(h+hr)); //[K] Body Temperature
+q = k*A*(Ti-Ts)/L; //[W]
+printf("\n\n (II) In presence of Water")
+printf("\n (a) Temperature of Skin = %.2f K",Ts);
+printf("\n (b) Total Heat Loss = %.2f W",q);
+
+//END
\ No newline at end of file diff --git a/534/CH1/EX1.7/1_7_Cure_Temperature.sce b/534/CH1/EX1.7/1_7_Cure_Temperature.sce new file mode 100644 index 000000000..83187249c --- /dev/null +++ b/534/CH1/EX1.7/1_7_Cure_Temperature.sce @@ -0,0 +1,58 @@ +clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 1.7 Page 30 \n')//Example 1.7
+// (a) Cure Temperature for h = 15 W/m^2
+// (b) Value of h for cure temp = 50 deg C
+
+Tsurr = 30+273; //[K] - Temperature of surrounding
+Tf = 20+273; //[K] - Temperature of Fluid Flow
+e=.5; // Emissivity of Surface
+a = .8; // Absorptivity of Surface
+G = 2000; //[W/m^2] - Irradiation falling on surface
+h=15; //[W/m^2.k] - Thermal Convectivity from plate to air
+stfncnstt=5.67*10^(-8); // [W/m^2.K^4] - Stefan Boltzmann Constant
+T=375; //[K] Value initially assumed for trial-error approach
+//Using Eq 1.3a & 1.7 and trial-and error approach of Newton Raphson
+while(1>0)
+f=((a*G)-(h*(T-Tf)+e*stfncnstt*(T^4 - Tsurr^4)));
+fd=(-h*T-4*e*stfncnstt*T^3);
+Tn=T-f/fd;
+if(((a*G)-(h*(Tn-Tf)+e*stfncnstt*(Tn^4 - Tsurr^4)))<=.01)
+ break;
+end;
+T=Tn;
+end
+
+printf("\n (a) Cure Temperature of Plate = %i degC\n",T-273);
+//solution (b)
+Treq=50+273;
+function[T]=Tvalue(h)
+ T=240;
+ while(1>0)
+ f=((a*G)-(h*(T-Tf)+e*stfncnstt*(T^4 - Tsurr^4)));
+ fd=(-h*T-4*e*stfncnstt*T^3);
+ Tn=T-f/fd;
+ if(((a*G)-(h*(Tn-Tf)+e*stfncnstt*(Tn^4 - Tsurr^4)))<=.01)
+ break;
+ end;
+ T=Tn;
+ end
+ funcprot(0)
+endfunction
+
+h = [2:.5:100];
+Tm = [1:1:197];
+for i=1:1:197;
+ Tm(i)=Tvalue(h(i));
+end
+
+T=Treq;
+hnew=((a*G)-(e*stfncnstt*(T^4 - Tsurr^4)))/(T-Tf);
+clf()
+xtitle("Graph Temp vs Convection Coeff", "h (W/m^2/K)", "T (degC)");
+x=[0 hnew hnew];
+y=[Treq-273 Treq-273 0];
+plot(h,Tm-273,x,y);
+legend("Plot","h at T = 50 degC");
+printf("\n (b) Air flow must provide a convection of = %i W/m^2.K", hnew);
+//END
\ No newline at end of file diff --git a/534/CH1/EX1.8/1_8_Theoretical_Problem.sce b/534/CH1/EX1.8/1_8_Theoretical_Problem.sce new file mode 100644 index 000000000..2d458cad6 --- /dev/null +++ b/534/CH1/EX1.8/1_8_Theoretical_Problem.sce @@ -0,0 +1,8 @@ +clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 1.8 Page 40 \n')// Example 1.8
+//Theoretical Problem
+
+printf('\n The given example is theoretical and does not involve any numerical computation')
+
+//End
|