summaryrefslogtreecommitdiff
path: root/914/CH13/EX13.1/ex13_1.sce
diff options
context:
space:
mode:
Diffstat (limited to '914/CH13/EX13.1/ex13_1.sce')
-rwxr-xr-x914/CH13/EX13.1/ex13_1.sce24
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