initial commit / add all books

5 files changed, 133 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
diff --git a/914/CH13/EX13.10/ex13_10.sce b/914/CH13/EX13.10/ex13_10.sce
new file mode 100755
index 000000000..6c71a88ff
--- /dev/null
+++ b/914/CH13/EX13.10/ex13_10.sce
@@ -0,0 +1,17 @@
+clc;

+warning("off");

+printf("\n\n example13.10 - pg701");

+// given

+d=0.01;  //[m] - diameter of cyclindrical porous plug

+D=2*10^-9;  //[m^2/sec] - diffusion coefficient

+t=60*60;  //[sec]

+r=d/2;

+m=0;

+Ca_inf=0;

+Ca_0=10;

+X=(D*t)/((r)^2);

+// from fig 13.14 the ordinate is

+Y=0.7;

+Ca_c=Ca_inf-Y*(Ca_inf-Ca_0);

+printf("\n\n the concentration of KCL at the centre after 60 min is \n Ca = %f kg/m^3",Ca_c);

+

diff --git a/914/CH13/EX13.6/ex13_6.sce b/914/CH13/EX13.6/ex13_6.sce
new file mode 100755
index 000000000..ab238cbe1
--- /dev/null
+++ b/914/CH13/EX13.6/ex13_6.sce
@@ -0,0 +1,29 @@
+clc;

+warning("off");

+printf("\n\n example13_6 - pg684");

+// given

+d=1*0.0254;  //[m]

+Lr=d/2;  //[m];

+Lz=(1.2/2)*(0.0254);

+x=Lz;

+r=Lr;

+k=0.481;

+h=20;

+mr=k/(h*Lr);

+mz=k/(h*Lz);

+nr=r/Lr;

+nz=x/Lz;

+t=1.2;  //[sec]

+alpha=1.454*10^-4;

+Xr=(alpha*t)/(Lr^2);

+Xz=(alpha*t)/(Lz^2);

+// using the above value of m,n,X the value for Ycz and Ycr from fig 13.14 is

+Ycr=0.42;

+Ycz=0.75;

+Yc=Ycr*Ycz;

+T_infinity=400;  //[K]

+To=295;

+Tc=T_infinity-(Yc*(T_infinity-To));

+printf("\n\n The temperature t the centre is \n Tc = %f K",Tc);

+

+

diff --git a/914/CH13/EX13.7/ex13_7.sce b/914/CH13/EX13.7/ex13_7.sce
new file mode 100755
index 000000000..52ea2973e
--- /dev/null
+++ b/914/CH13/EX13.7/ex13_7.sce
@@ -0,0 +1,43 @@
+clc;

+warning("off");

+printf("\n\n example13_7 - pg684");

+// given

+T_x0=300;  //[K]

+Tw=400;  //[K]

+L=0.013;  //[m]

+alpha=2.476*(10^-5);  //[m^/sec]

+h=600;  //[W/m^2*K]

+pcp=3.393*(10^6);  //[J/m^3*K]

+L=0.013;  //[m]

+deltax=L/10;

+betaa=0.5;

+deltat=0.03;

+deltat=betaa*((deltax)^2)*(1/alpha);

+T_infinity=400;  //[K]

+// to be sure that the solution is stable, it is customary to truncate this number

+deltat=0.03;  //[sec]

+// betaa=alpha*deltat*((1/deltax)^2);

+    for i=1:11

+    Told(i)=300;

+end

+a=((2*h*deltat)/(pcp*deltax));

+b=((2*alpha*deltat)/(pcp*((deltax)^2)));

+for j=1:11

+Tnew(1)=(T_infinity*0.08162)+(Told(1)*(1-0.08162-0.8791))+(Told(2)*0.8791)

+for k=1:9

+    Tnew(k+1)=(betaa*Told(k+2))+((1-2*betaa)*(Told(k+1)))+(betaa*Told(k));

+end

+Tnew(11)=((2*betaa)*(Told(10)))

+Told=Tnew;

+end

+disp(Told);

+

+

+

+

+

+

+

+

+

+

diff --git a/914/CH13/EX13.9/ex13_9.sce b/914/CH13/EX13.9/ex13_9.sce
new file mode 100755
index 000000000..40627205a
--- /dev/null
+++ b/914/CH13/EX13.9/ex13_9.sce
@@ -0,0 +1,20 @@
+clc;

+warning("off");

+printf("\n\n example 13_9 - pg700");

+// given

+p=2050;  //[kg/m^3] - density of soil

+cp=1840;  //[J/kg*K] - heat cpapacity of soil

+k=0.52;  //[W/m*K] - thermal conductivity of soil

+alpha=0.138*10^-6;  //[m^2/sec]

+t=4*30*24*3600;  //[sec] - no. of seconds in 4 months

+Tx=-5;  //[degC]

+Tinf=-20;  //[degC]

+T0=20;  //[degC]

+// from the fig 13.24 the dimensionless distance Z is 

+Z=0.46;

+// then the depth is

+x=2*((alpha*t)^(1/2))*Z

+printf("\n\n the depth is \n x = %f m = %f ft",x,x*(3.6/1.10));

+

+

+