initial commit / add all books

7 files changed, 141 insertions, 0 deletions

diff --git a/965/CH3/EX3.1/1.sci b/965/CH3/EX3.1/1.sci
new file mode 100644
index 000000000..5518f928a
--- /dev/null
+++ b/965/CH3/EX3.1/1.sci
@@ -0,0 +1,40 @@
+clc;

+clear all;

+disp("steady state temperature distribution")

+disp("Let th = t-ta")

+disp("the controllign differential equation for the given problem is given by")

+disp("d2th/dx2+d2th/dy2 =0------(1)")

+disp("the boundary conditions are :")

+disp("i) at x = infinity, th =0")

+disp("ii) at x = 0, th =th0")

+disp("iii) at y =L, th =0")

+disp("iv) at y = 0, th =0")

+disp("The solution of eq. 1 is  th = X(x)Y(y) ------ (2)")

+disp("substituting the solution in controlling equation, we get ")

+disp("1/X*d2X/dx2 =-1/Y*d2Y/dy2 =  + or - lambda^2")

+disp("The required equations are :")

+disp("d2X^2/dx2-lambda^2*X =0 ------(iii)")

+disp("d2Y^2/dy2+lambda^2*Y =0 ------(iv)")

+disp("the solutions of eqns are :")

+disp("X = A*exp(lambda*x)+B*exp(-lambda*x)")

+disp("Y = C*cos(lambda*y)+D*sin(lambda*y)")

+disp("th = (A*exp(lambda*x)+B*exp(-lambda*x))*(C*cos(lambda*y)+D*sin(lambda*y))")

+disp(" from boundary condition i), we have ")

+disp("0 = (A*exp(lambda*x)+B*exp(-lambda*x))*(C*cos(lambda*y)+D*sin(lambda*y)")

+disp("A = 0 and  th =B*(C*cos(lambda*y)+D*sin(lambda*y)")

+disp("from boundary condition iv), we have")

+disp("0 = C*B*exp(-lambda*x)")

+disp("hence C = 0 and equation reduces to  th = B*D*sin(lambda*y)*exp(-lambda*x)")

+disp("from boundary condition iii) we get, 0 = E*exp(-lambda*x)*sin(lambda*L), where E = B*D")

+disp("since E is not 0, sin (lambda*L)=0")

+disp("lambda = 0, %pi/L,2*%pi/L.....")

+disp("lambdan = n*%pi/L, where n = 0,1,2....")

+disp("hence , th = E*exp(-lambdan*x)*sin(lambdan*y)")

+disp("from boundary eqn ii) we have ")

+disp("th = sum(En*sin(lambdan*y), 1, infinity)")

+disp("This is an expression of th0 in terms of Fourier series, where En are Fourier coefficients.")

+disp("by integrating we get")

+disp("th = 2*th0/L*(sum(((1-(-1)^n)/lambdan *exp(-lambdan*x)*sin(lambdan*y))")

+

+,

+

diff --git a/965/CH3/EX3.3/3.sci b/965/CH3/EX3.3/3.sci
new file mode 100644
index 000000000..f5ec0484a
--- /dev/null
+++ b/965/CH3/EX3.3/3.sci
@@ -0,0 +1,12 @@
+clc;

+clear all;

+disp("Heat loss calculation")

+r=0.6/2;//m

+L=1;//m

+H=1.8;//m

+k=0.51;// W/(m*C)

+tp=95;// degree C

+te=25;// degree C

+Sfc=2*3.1416*L/(log(2*H/r));

+Q=k*Sfc*(tp-te);

+disp("W",Q,"Heat loss from the pipe meter length, Q =")

diff --git a/965/CH3/EX3.4/4.sci b/965/CH3/EX3.4/4.sci
new file mode 100644
index 000000000..545ca9a36
--- /dev/null
+++ b/965/CH3/EX3.4/4.sci
@@ -0,0 +1,14 @@
+clc;

+clear all;

+disp("Sphere Surface temperature")

+r=1.6/2;//m

+

+H=5.5;//m

+k=0.51;// W/(m*C)

+Qg=580;// W

+

+te=6;// degree C

+Sfc=4*3.1416*r/(1-r/(2*H));

+//Qg=k*Sfc*(t-te);

+t=Qg/(k*Sfc)+te;

+disp("degree C",t,"surface temperature of sphere, t =")

diff --git a/965/CH3/EX3.5/5.sci b/965/CH3/EX3.5/5.sci
new file mode 100644
index 000000000..52758e6cf
--- /dev/null
+++ b/965/CH3/EX3.5/5.sci
@@ -0,0 +1,18 @@
+clc;

+clear all;

+disp("Heat loss/length")

+a=0.6;// m

+b=0.45;//m

+H=1.5;//m

+k=0.51;// W/(m*C)

+tp=105;// degree C

+ts=5;// degree C

+

+x=log(1+H/a);

+y=H/b;

+

+Sfc=2.756*((x^(-.59))*(y^(-.078)));

+

+Q=k*Sfc*(tp-ts);

+

+disp("W",Q,"Heat loss per meter length, Q =")

diff --git a/965/CH3/EX3.6/6.sci b/965/CH3/EX3.6/6.sci
new file mode 100644
index 000000000..b95cbc9e2
--- /dev/null
+++ b/965/CH3/EX3.6/6.sci
@@ -0,0 +1,16 @@
+clc;

+clear all;

+disp("heat transfer rate, temperature")

+r1=0.600/2;// m

+r2=0.18/2;//m

+D=1.8;//m

+L=90;//m

+tp1=180;// degree C

+tp2=12;// degree C

+k=0.45;// W/(m*C)

+

+Sfc=2*3.1416*L/(acosh((D^2-r1^2-r2^2)/(2*r1*r2)));

+Sfc

+Q=k*Sfc*(tp1-tp2);

+

+disp("W",Q,"Heat loss per meter length, Q =")

diff --git a/965/CH3/EX3.7/7.sci b/965/CH3/EX3.7/7.sci
new file mode 100644
index 000000000..2332ec997
--- /dev/null
+++ b/965/CH3/EX3.7/7.sci
@@ -0,0 +1,19 @@
+clc;

+clear all;

+disp("Heat loss calculations")

+k=1.05;// W/(m*C)

+ti=600;//degree C

+to=70;// degree C

+L=0.12;//m

+A=0.6*0.6;//m^2

+Sfcwall=A/L;

+

+D=0.6;//m

+Sfcedge=0.54*D;

+

+Sfccorner=0.15*L;

+Sfct=6*Sfcwall+12*Sfcedge+8*Sfccorner;

+

+Q=k*Sfct*(ti-to);

+

+disp("W",Q,"Heat loss through walls, Q =")

diff --git a/965/CH3/EX3.8/8.sci b/965/CH3/EX3.8/8.sci
new file mode 100644
index 000000000..839c75bee
--- /dev/null
+++ b/965/CH3/EX3.8/8.sci
@@ -0,0 +1,22 @@
+clc;

+clear all;

+disp("Heat loss rate")

+ti=350;// degree C

+to=50;// degree C

+k=1.5;// W/(m*C)

+x=0.3;// m

+l=4.5;//m

+b=3.75;//m

+h=3;//m

+Ai=2*(l*b+l*h+b*h);//m^2

+y=4*(l+b+h);// length of edges

+Sfce=y;

+Sfcc=0.15*x;

+

+Am=Ai+Sfce*.54*x+8*Sfcc*x;

+Am

+disp("m^2")

+Q=k*Am*(ti-to)/x;

+disp("W",Q,"rate of heat loss Q =")

+

+