16 files changed, 138 insertions, 0 deletions
diff --git a/1309/CH1/EX1.1/Result1_1.pdf b/1309/CH1/EX1.1/Result1_1.pdf
new file mode 100755
index 000000000..4a2d7d865
--- /dev/null
+++ b/1309/CH1/EX1.1/Result1_1.pdf
diff --git a/1309/CH1/EX1.1/ch1_1.sce b/1309/CH1/EX1.1/ch1_1.sce
new file mode 100755
index 000000000..9c567fa6a
--- /dev/null
+++ b/1309/CH1/EX1.1/ch1_1.sce
@@ -0,0 +1,17 @@
+clear;
+clc;
+printf("\t\t\tChapter1_Example1\n\n\n");
+// determination of surface temperature on one side of a firewall
+k=9.4; // thermal conductivity in [BTU/hr.ft. ˚Rankine]
+q=6.3; // heat flux in [BTU/s. sq.ft]
+T1=350; // the outside surface temperature of one aide of the wall [˚F]
+// converting heat flux into BTU/hr  sq.ft
+Q=6.3*3600 // [BTU/hr.sq.ft]
+printf("\nThe heat flux is %.2f BTU/hr. sq.ft",Q);
+dx=0.5; // thickness in [inch]
+//converting distance into ft
+Dx=0.5/12 // thickness in [ft]
+printf("\nThe thickness of firewall is %.2f ft",Dx);
+// solving for temeprature T2
+T2=T1-(Q*Dx/k); // [˚F]
+printf("\nThe required temperature on the other side of the firewall is %.1f degree Fahrenheit",T2);
diff --git a/1309/CH1/EX1.2/Result1_2.pdf b/1309/CH1/EX1.2/Result1_2.pdf
new file mode 100755
index 000000000..f3ad05de8
--- /dev/null
+++ b/1309/CH1/EX1.2/Result1_2.pdf
diff --git a/1309/CH1/EX1.2/ch1_2.sce b/1309/CH1/EX1.2/ch1_2.sce
new file mode 100755
index 000000000..83cb01171
--- /dev/null
+++ b/1309/CH1/EX1.2/ch1_2.sce
@@ -0,0 +1,12 @@
+clear;
+clc;
+printf("\t\t\tchapter1_example2\n\n\n");
+// determination of thermal conductivity of aluminium
+k_ss=14.4; // thermal conductivity of stainless steel in [W/m.K]
+printf("\nThe thermal conductivity of stainless steel is %.1f W/m.K",k_ss);
+dt_ss=40; // [K]
+dt_al=8.65; // [K]
+dz_ss=1; // [cm]
+dz_al=3; // [cm]
+k_al=k_ss*dt_ss*dz_al/(dt_al*dz_ss);// thermal conductivity of Al in [W/m.K]
+printf("\nThe thermal conductivity of aluminium is %d W/m.K",k_al);
diff --git a/1309/CH1/EX1.3/Result1_3.pdf b/1309/CH1/EX1.3/Result1_3.pdf
new file mode 100755
index 000000000..720c31a67
--- /dev/null
+++ b/1309/CH1/EX1.3/Result1_3.pdf
diff --git a/1309/CH1/EX1.3/ch1_3.sce b/1309/CH1/EX1.3/ch1_3.sce
new file mode 100755
index 000000000..91be46380
--- /dev/null
+++ b/1309/CH1/EX1.3/ch1_3.sce
@@ -0,0 +1,11 @@
+clear;
+clc;
+printf("\t\t\tchapter1_example3\n\n\n");
+// determination of heat transferred by convection
+h_c=3; // convective coefficient in [BTU/hr.ft^2
+A=30*18; // Cross sectional area in ft^2
+T_w=140; // Roof surface temperature in degree Fahrenheit
+T_inf=85; // Ambient temperature in degree Fahrenheit
+dT= (T_w-T_inf);
+Q_c=h_c*A*dT; // Convective heat transfer in BTU/hr
+printf("\nThe heat transferred by convection is %d BTU/hr",Q_c);
diff --git a/1309/CH1/EX1.4/Result1_4.pdf b/1309/CH1/EX1.4/Result1_4.pdf
new file mode 100755
index 000000000..9e3585de3
--- /dev/null
+++ b/1309/CH1/EX1.4/Result1_4.pdf
diff --git a/1309/CH1/EX1.4/ch1_4.sce b/1309/CH1/EX1.4/ch1_4.sce
new file mode 100755
index 000000000..ad120cb5b
--- /dev/null
+++ b/1309/CH1/EX1.4/ch1_4.sce
@@ -0,0 +1,18 @@
+clear;
+clc;
+printf("\t\t\tchapter1_example4\n\n\n");
+// determining average film conductance
+D=2.43/100; // diameter in meter
+L=20/100; // length in meter
+A=3.14*D*L; // cross-sectional area in sq.m
+cp=4200; // specific heat of water in J/kg.K
+T_b2=21.4; // temperature of bulk fluid in degree celsius
+T_in=20; // temperature of inlet water in degree celsius
+T_w=75; // temperature of wall in degree celsius
+Q=500; // volumetric flow rate in cc/s
+density=1000; // density of water in kg/cu.m
+m=Q*density/10^6; // mass flowa rate in kg/s
+printf("\nThe mass flow rate is %.1f kg/s",m);
+// using definition of specific heat and Newton's law of cooling
+hc=m*cp*(T_b2-T_in)/(A*(T_w-T_in));
+printf("\nThe average film conductance is %d W/sq.m. K",hc);
diff --git a/1309/CH1/EX1.5/Result1_5.pdf b/1309/CH1/EX1.5/Result1_5.pdf
new file mode 100755
index 000000000..dde281729
--- /dev/null
+++ b/1309/CH1/EX1.5/Result1_5.pdf
diff --git a/1309/CH1/EX1.5/ch1_5.sce b/1309/CH1/EX1.5/ch1_5.sce
new file mode 100755
index 000000000..c4d0e950f
--- /dev/null
+++ b/1309/CH1/EX1.5/ch1_5.sce
@@ -0,0 +1,15 @@
+clear;
+clc;
+printf("\t\t\tchapter1_example5\n\n\n");
+// determination of heat loss rate by radiation
+W=14; // width in ft
+L=30; // length in ft
+A=W*L; // area in ft^2
+F_12=1; // view factor assumed to be 1
+T1=120+460; // driveway surface temperature  in degree Rankine
+printf("\nThe driveway surface temperature is %d degree Rankine",T1);
+T2=0; // space temperature assumed to be 0 degree Rankine
+sigma=0.1714*10^(-8); // value of Stefan-Boltzmann's constant in BTU/(hr.ft^2.(degree Rankine)^4)
+e=0.9; // surface emissivity
+q=sigma*A*e*F_12*((T1)^4-(T2)^4);
+printf("\nThe heat loss rate by radiation is %d BTU/hr",q);
diff --git a/1309/CH1/EX1.6/Result1_6.pdf b/1309/CH1/EX1.6/Result1_6.pdf
new file mode 100755
index 000000000..32a4b1bae
--- /dev/null
+++ b/1309/CH1/EX1.6/Result1_6.pdf
diff --git a/1309/CH1/EX1.6/ch1_6.sce b/1309/CH1/EX1.6/ch1_6.sce
new file mode 100755
index 000000000..f25808fc1
--- /dev/null
+++ b/1309/CH1/EX1.6/ch1_6.sce
@@ -0,0 +1,10 @@
+clear;
+clc;
+printf("\t\t\tchapter1_example6\n\n\n");
+// determination of radiation thermal conductance
+A=14*30; // area in sq.ft
+T1=120+460; // driveway surface temperature in degree Rankine
+T2=0; // surface temperature assumed to be 0 degree Rankine
+Qr=73320; // heat loss rate in BTU/hr
+hr=Qr/(A*(T1-T2)); // radiation thermal conductance in BTU/(hr.ft^2.(degree Rankine)
+printf("\nthe radiation thermal conductance is %.2f BTU/(hr. sq.ft.(degree Rankine))",hr);
diff --git a/1309/CH1/EX1.7/Result1_7.pdf b/1309/CH1/EX1.7/Result1_7.pdf
new file mode 100755
index 000000000..13c21f68a
--- /dev/null
+++ b/1309/CH1/EX1.7/Result1_7.pdf
diff --git a/1309/CH1/EX1.7/ch1_7.sce b/1309/CH1/EX1.7/ch1_7.sce
new file mode 100755
index 000000000..07e10a8c1
--- /dev/null
+++ b/1309/CH1/EX1.7/ch1_7.sce
@@ -0,0 +1,27 @@
+clear;
+clc;
+printf("\t\t\tchapter1_example7\n\n\n");
+// Identification of all resistances and their values
+// Estimation of heat transfer per unit area
+// Determination of the inside and outside wall temperatures
+printf("\n\t\t\tSolution to part (b)\n");
+A=1; // assuming A=1 m^2 for convenience
+hc1_avg=(5+25)/2; // taking average of extreme values for hc [W/m^2.K]
+Rc1=1/(hc1_avg*A); // resistance on left side of wall [K/W]
+printf("\nThe resistance on left side of wall is %.3f K/W",Rc1);
+k=(0.38+0.52)/2; // thermal conductivity of common brick in W/M.k
+L=0.1; //10 cm converted into m
+Rk=(L/(k*A));// resistance of construction material, assume common brick
+printf("\nThe resistance of construction material of wall is %.3f K/W",Rk);
+Rc2=Rc1;
+printf("\nThe resistance on right side of wall is %.3f K/W",Rc2);
+printf("\n\n\t\t\tSolution to part (c)\n");
+T_inf1=1000; // temperature of exhaust gases in K
+T_inf2=283; // temperature of ambient air in K
+q=(T_inf1-T_inf2)/(Rc1+Rk+Rc2); // heat transferred per unit area
+printf("\nThe Heat transferred per unit area is %d W = %.3f kW",q,q/1000);
+printf("\n\n\t\t\tSolution to part (d)\n");
+T_in=T_inf1-Rc1*q; //
+T_out=T_inf2+Rc2*q;
+printf("\nThe inside wall temperature is %d K",T_in);
+printf("\nThe outside wall temperature is %d K",T_out);
diff --git a/1309/CH1/EX1.8/Result1_8.pdf b/1309/CH1/EX1.8/Result1_8.pdf
new file mode 100755
index 000000000..5d538e1e0
--- /dev/null
+++ b/1309/CH1/EX1.8/Result1_8.pdf
diff --git a/1309/CH1/EX1.8/ch1_8.sce b/1309/CH1/EX1.8/ch1_8.sce
new file mode 100755
index 000000000..a7f679141
--- /dev/null
+++ b/1309/CH1/EX1.8/ch1_8.sce
@@ -0,0 +1,28 @@
+clear;
+clc;
+printf("\t\t\tchapter1_example8\n\n\n");
+// determination of surface temperature
+k=0.604; // [BTU/(hr.ft.degree Rankine)]
+hc=3; // average value for natural convection in BTU/(hr.ft^2.degree Rankine)
+ew=0.93;
+f_wr=1; // shape factor
+sigma= 0.1714*10^(-8) // BTU/(hr.ft^2.degree Rankine).
+L=4/12; // length in ft
+T1=80+460; // temperature of side-walk in degree Rankine
+T_inf=20+460; // temperature of ambient air in degree Rankine
+T_r=0; // assuming space temperature to be 0 degree Rankine
+// LHS of the form a*Tw+b*Tw^4=c
+a=((k/L)+hc);
+b=(sigma*ew*f_wr);
+c=(k*T1/L)+(hc*T_inf)+(sigma*f_wr*ew*T_r^4);
+printf("\nRHS=%d",c);
+    Tw=[470;480;490;485;484.5];
+for i=1:5
+    LHS(i)=a*Tw(i)+b*Tw(i)^4;
+end
+printf("\nSolving by trial and error yields the following, where LHS is the left-hand side of the equation");
+printf("\n\tTw\tLHS");
+for i=1:5
+    printf("\n\t%.1f\t%d",Tw(i),LHS(i));
+end
+printf("\nThe Surface temperature is %.1f degree R = %.1f degree F",Tw(5),Tw(i)-460);