diff options
Diffstat (limited to '914/CH13/EX13.1/ex13_1.sce')
-rwxr-xr-x | 914/CH13/EX13.1/ex13_1.sce | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/914/CH13/EX13.1/ex13_1.sce b/914/CH13/EX13.1/ex13_1.sce new file mode 100755 index 000000000..9319d61b5 --- /dev/null +++ b/914/CH13/EX13.1/ex13_1.sce @@ -0,0 +1,24 @@ +clc;
+warning("off");
+printf("\n\n example13.1 - pg651");
+// given
+h=12; //[W/m^2*K] - heat transfer coefficeint
+k=400; //[W/m*K] - thermal conductivity
+// (a) for sphere
+r=5*10^-2; //[m] - radius of copper sphere
+Lc=((4*%pi*((r)^3))/3)/(4*%pi*((r)^2));
+Nbi=h*Lc*(1/k);
+printf("\n\n (a) The biot no. is \n Nbi = %e",Nbi);
+// (b) for cyclinder
+r=0.05; //[m] - radius of cyclinder
+L=0.3; //[m] - height of cyclinder
+Lc=(%pi*((r)^2)*L)/(2*%pi*r*L);
+Nbi=h*Lc*(1/k);
+printf("\n\n (b) The biot no. is \n Nbi = %e",Nbi);
+// (c) for a long square rod
+L=.4; //[m] - length of copper rod
+r=0.05; //[m] - radius of a cyclinder havimg same cross sectional area as that of square
+x=((%pi*r^2)^(1/2));
+Lc=((x^2)*L)/(4*x*L);
+Nbi=h*Lc*(1/k);
+printf("\n\n (c) The biot no. is \n Nbi = %e",Nbi);
\ No newline at end of file |