summaryrefslogtreecommitdiff
path: root/534/CH10
diff options
context:
space:
mode:
Diffstat (limited to '534/CH10')
-rw-r--r--534/CH10/EX10.1/10_1_Boiling_Water_pan.sce35
-rw-r--r--534/CH10/EX10.2/10_2_Horizontal_cylinder.sce43
-rw-r--r--534/CH10/EX10.3/10_3_Condensation_Chimney.sce36
-rw-r--r--534/CH10/EX10.4/10_4_Steam_Condenser.sce32
4 files changed, 146 insertions, 0 deletions
diff --git a/534/CH10/EX10.1/10_1_Boiling_Water_pan.sce b/534/CH10/EX10.1/10_1_Boiling_Water_pan.sce
new file mode 100644
index 000000000..93db94d87
--- /dev/null
+++ b/534/CH10/EX10.1/10_1_Boiling_Water_pan.sce
@@ -0,0 +1,35 @@
+clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 10.1 Page 632 \n'); //Example 10.1
+// Power Required by electruc heater to cause boiling
+// Rate of water evaporation due to boiling
+// Critical Heat flux corresponding to the burnout point
+
+//Operating Conditions
+Ts = 118+273 ;//[K] Surface Temperature
+Tsat = 100+273 ;//[K] Saturated Temperature
+D = .3 ;//[m] Diameter of pan
+g = 9.81 ;//[m^2/s] gravitaional constant
+//Table A.6 Saturated water Liquid Properties T = 373 K
+rhol = 957.9 ;//[kg/m^3] Density
+cp = 4.217*10^3 ;//[J/kg] Specific Heat
+u = 279*10^-6 ;//[N.s/m^2] Viscosity
+Pr = 1.76 ;// Prandtl Number
+hfg = 2257*10^3 ;//[J/kg] Specific Heat
+si = 58.9*10^-3 ;//[N/m]
+//Table A.6 Saturated water Vapor Properties T = 373 K
+rhov = .5956 ;//[kg/m^3] Density
+
+Te = Ts-Tsat;
+//From Table 10.1
+C = .0128;
+n = 1;
+q = u*hfg*[g*(rhol-rhov)/si]^.5*(cp*Te/(C*hfg*Pr^n))^3;
+qs = q*%pi*D^2/4;
+
+m = qs/hfg;
+
+qmax = .149*hfg*rhov*[si*g*(rhol-rhov)/rhov^2]^.25;
+
+printf("\n Boiling Heat transfer rate = %.1f kW \n Rate of water evaporation due to boiling = %i kg/h \n Critical Heat flux corresponding to the burnout point = %.2f MW/m^2",qs/1000,m*3600,qmax/10^6);
+//END \ No newline at end of file
diff --git a/534/CH10/EX10.2/10_2_Horizontal_cylinder.sce b/534/CH10/EX10.2/10_2_Horizontal_cylinder.sce
new file mode 100644
index 000000000..97bf09ef6
--- /dev/null
+++ b/534/CH10/EX10.2/10_2_Horizontal_cylinder.sce
@@ -0,0 +1,43 @@
+clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 10.2 Page 635 \n'); //Example 10.2
+// Power Dissipation per unith length for the cylinder, qs
+
+//Operating Conditions
+Ts = 255+273 ;//[K] Surface Temperature
+Tsat = 100+273 ;//[K] Saturated Temperature
+D = 6*10^-3 ;//[m] Diameter of pan
+e = 1 ;// eimssivity
+stfncnstt=5.67*10^(-8) ;// [W/m^2.K^4] - Stefan Boltzmann Constant
+g = 9.81 ;//[m^2/s] gravitaional constant
+//Table A.6 Saturated water Liquid Properties T = 373 K
+rhol = 957.9 ;//[kg/m^3] Density
+hfg = 2257*10^3 ;//[J/kg] Specific Heat
+//Table A.4 Water Vapor Properties T = 450 K
+rhov = .4902 ;//[kg/m^3] Density
+cpv = 1.98*10^3 ;//[J/kg.K] Specific Heat
+kv = 0.0299 ;//[W/m.K] Conductivity
+uv = 15.25*10^-6 ;//[N.s/m^2] Viscosity
+
+Te = Ts-Tsat;
+
+hconv = .62*[kv^3*rhov*(rhol-rhov)*g*(hfg+.8*cpv*Te)/(uv*D*Te)]^.25;
+hrad = e*stfncnstt*(Ts^4-Tsat^4)/(Ts-Tsat);
+
+//From eqn 10.9 h^(4/3) = hconv^(4/3) + hrad*h^(1/3)
+//Newton Raphson
+h=250; //Initial Assumption
+while(1>0)
+f = h^(4/3) - [hconv^(4/3) + hrad*h^(1/3)];
+fd = (4/3)*h^(1/3) - [(1/3)*hrad*h^(-2/3)];
+hn=h-f/fd;
+if((hn^(4/3) - [hconv^(4/3) + hrad*hn^(1/3)])<=.01)
+ break;
+end;
+h=hn;
+end
+
+q = h*%pi*D*Te;
+
+printf("\n Power Dissipation per unith length for the cylinder, qs= %i W/m",q);
+//END \ No newline at end of file
diff --git a/534/CH10/EX10.3/10_3_Condensation_Chimney.sce b/534/CH10/EX10.3/10_3_Condensation_Chimney.sce
new file mode 100644
index 000000000..1153e9ae8
--- /dev/null
+++ b/534/CH10/EX10.3/10_3_Condensation_Chimney.sce
@@ -0,0 +1,36 @@
+clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 10.3 Page 648 \n'); //Example 10.3
+// Heat Transfer and Condensation Rates
+
+//Operating Conditions
+Ts = 50+273 ;//[K] Surface Temperature
+Tsat = 100+273 ;//[K] Saturated Temperature
+D = .08 ;//[m] Diameter of pan
+g = 9.81 ;//[m^2/s] gravitaional constant
+L = 1 //[m] Length
+//Table A.6 Saturated Vapor Properties p = 1.0133 bars
+rhov = .596 ;//[kg/m^3] Density
+hfg = 2257*10^3 ;//[J/kg] Specific Heat
+//Table A.6 Saturated water Liquid Properties T = 348 K
+rhol = 975 ;//[kg/m^3] Density
+cpl = 4193 ; //[J/kg.K] Specific Heat
+kl = 0.668 ;//[W/m.K] Conductivity
+ul = 375*10^-6 ;//[N.s/m^2] Viscosity
+uvl = ul/rhol; ;//[N.s.m/Kg] Kinematic viscosity
+Ja = cpl*(Tsat-Ts)/hfg;
+hfg2 = hfg*(1+.68*Ja);
+//Equation 10.43
+Re = [3.70*kl*L*(Tsat-Ts)/(ul*hfg2*(uvl^2/g)^.33334)+4.8]^.82;
+
+//From equation 10.41
+hL = Re*ul*hfg2/(4*L*(Tsat-Ts));
+q = hL*(%pi*D*L)*(Tsat-Ts);
+
+m = q/hfg;
+//Using Equation 10.26
+del = [4*kl*ul*(Tsat-Ts)*L/(g*rhol*(rhol-rhov)*hfg2)]^.25;
+
+
+printf("\n Heat Transfer Rate = %.1f kW and Condensation Rates= %.4f kg/s \n And as del(L) %.3f mm << (D/2) %.2f m use of vertical cylinder correlation is justified",q/1000,m,del*1000,D/2);
+//END \ No newline at end of file
diff --git a/534/CH10/EX10.4/10_4_Steam_Condenser.sce b/534/CH10/EX10.4/10_4_Steam_Condenser.sce
new file mode 100644
index 000000000..cb6a4ca11
--- /dev/null
+++ b/534/CH10/EX10.4/10_4_Steam_Condenser.sce
@@ -0,0 +1,32 @@
+clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 10.4 Page 652 \n'); //Example 10.4
+// Condensation rate per unit length of tubes
+
+//Operating Conditions
+Ts = 25+273 ;//[K] Surface Temperature
+Tsat = 54+273 ;//[K] Saturated Temperature
+D = .006 ; //[m] Diameter of pan
+g = 9.81 ;//[m^2/s] gravitaional constant
+N = 20 // No of tubes
+
+//Table A.6 Saturated Vapor Properties p = 1.015 bar
+rhov = .098 ;//[kg/m^3] Density
+hfg = 2373*10^3 ;//[J/kg] Specific Heat
+//Table A.6 Saturated water Liquid Properties Tf = 312.5 K
+rhol = 992 ;//[kg/m^3] Density
+cpl = 4178 ;//[J/kg.K] Specific Heat
+kl = 0.631 ; //[W/m.K] Conductivity
+ul = 663*10^-6 ; //[N.s/m^2] Viscosity
+
+Ja = cpl*(Tsat-Ts)/hfg;
+hfg2 = hfg*(1+.68*Ja);
+//Equation 10.46
+h = .729*[g*rhol*(rhol-rhov)*kl^3*hfg2/(N*ul*(Tsat-Ts)*D)]^.25;
+//Equation 10.34
+m1 = h*(%pi*D)*(Tsat-Ts)/hfg2;
+
+m = N^2*m1;
+
+printf("\n For the complete array of tubes, the condensation per unit length is %.3f kg/s.m",m);
+//END \ No newline at end of file
generated by cgit (git 2.34.1) at 2024-12-19 12:45:01 +0000