1 files changed, 50 insertions, 0 deletions
diff --git a/1309/CH4/EX4.7/ch4_7.sce b/1309/CH4/EX4.7/ch4_7.sce
new file mode 100755
index 000000000..99d2b30c3
--- /dev/null
+++ b/1309/CH4/EX4.7/ch4_7.sce
@@ -0,0 +1,50 @@
+clc;
+clear;
+printf("\t\t\tChapter4_example7\n\n\n");
+// properties of water at 68 degree fahrenheit from appendix table C11
+rou=62.46;
+cp=0.9988;
+k=0.345;
+alpha=k/(rou*cp);
+printf("\nThe diffusivity at 68 degree fahrenheit is %.2e sq.ft/hr",alpha);
+D=2.5/12;
+L=4.75/12;
+Vs=%pi*D^2*L/4;
+As=(%pi*D*L)+(%pi*D^2)/2;
+Lc=Vs/As;
+printf("\nThe volume of the can is %.4f cu.ft",Vs);
+printf("\nThe surface area of the can is %.3f sq.ft",As);
+printf("\nThe characteristic length is %.3f ft",Lc);
+hc=1.7;
+Bi=hc*Lc/k;
+printf("\nThe Biot number is %.3f",Bi);
+t=4;
+// for the cylinder solution
+Fo_cylinder=alpha*t/(D/2)^2;
+Bi_cylinder=hc*(D/2)/k;
+printf("\nFor the cylinder, The Fourier number is %.2f and Biot Number is %.3f",Fo_cylinder,Bi_cylinder);
+reciprocal_Bi_cylinder=1/Bi_cylinder;
+printf("\nThe reciprocal for Biot number for cylinder is %.2f",reciprocal_Bi_cylinder);
+dim_T_cylinder=0.175; //The value of dimensionless temperature of cylinder from figure 4.7a at corresponding values of Fo and 1/Bi
+// for the infinite plate solution
+Fo_plate=alpha*t/(L/2)^2;
+Bi_plate=hc*L/(2*k);
+printf("\nFor the infinite plate, The Fourier number is %.3f and Biot Number is %.2f",Fo_plate,Bi_plate);
+reciprocal_Bi_plate=1/Bi_plate;
+printf("\nThe reciprocal for Biot number for infinite plate is %.2f",reciprocal_Bi_plate);
+dim_T_plate=0.55; //The value of dimensionless temperature of infinite plate from figure 4.7a at corresponding values of Fo and 1/Bi
+// Table 4. I, for the short-cylinder problem, indicates that the solution is the product of the infinite-cylinder problem (Figure 4.7) and the infinite-plate problem (Figure 4.6).
+// For short cylinder problem
+dim_T_shortcylinder=dim_T_cylinder*dim_T_plate;
+printf("\nThe value of dimensionless temperature for short cylinder is %.3f ",dim_T_shortcylinder);
+T_inf=30;
+T_i=72;
+Tc=dim_T_shortcylinder*(T_i-T_inf)+T_inf;
+printf("\nThe temperature at centre of can is %.1f degree celsius",Tc);
+dim_Tw_cylinder=0.77; //The dimensionless temperature from figure 4.7b corresponding to the value of 1/Bi and r/R=1
+dim_Tw_plate=0.65; //The dimensionless temperature from figure 4.6b corresponding to the value of 1/Bi and x/L=1
+dim_Tw_shortcylinder=dim_Tw_cylinder*dim_Tw_plate;
+printf("\nThe value of dimensionless temperature  at the wall for short cylinder is %.2f ",dim_Tw_shortcylinder);
+Tw=dim_Tw_shortcylinder*(Tc-T_inf)+T_inf;
+printf("\nThe wall temperature is %.1f degree F",Tw);
+