7 files changed, 171 insertions, 0 deletions
diff --git a/1379/CH9/EX9.1.1/example9_1.sce b/1379/CH9/EX9.1.1/example9_1.sce
new file mode 100755
index 000000000..011f27a93
--- /dev/null
+++ b/1379/CH9/EX9.1.1/example9_1.sce
@@ -0,0 +1,28 @@
+
+
+//exapple 9.1 
+clc; funcprot(0);
+// Initialization of Variable
+rho=1.2;
+mu=1.85/100000;
+pi=3.1428;
+d=3;
+v=50*1000/3600;
+//calculation part 1
+Re=d*rho*v/mu;
+//from chart of drag coeff. vs Re
+Cd=0.2;//coeff. of drag
+Ad=pi*d^2/4;//projected area
+Fd=Ad*Cd*rho*v^2/2;
+disp(Fd , "The drag force on sphere in N");
+//part 2
+v=2;
+l=0.25;
+Re=l*v*rho/mu;
+zi=4*pi*(l^3*3/4/pi)^(2/3)/6/l^2;//sphericity
+//using graph
+Cd=2;
+Ad=l^2;
+Fd=Ad*Cd*rho*v^2/2;
+disp(Fd , "The drag force on cube in N");
+
diff --git a/1379/CH9/EX9.1.2/example9_2.sce b/1379/CH9/EX9.1.2/example9_2.sce
new file mode 100755
index 000000000..3dec8c4c6
--- /dev/null
+++ b/1379/CH9/EX9.1.2/example9_2.sce
@@ -0,0 +1,29 @@
+
+
+//exapple 9.2 
+clc; funcprot(0);
+// Initialization of Variable
+rho=1.2;
+mu=1.85/100000;
+pi=3.1428;
+g=9.81;
+d=1.38;
+t=0.1;//thickness
+v=30*1000/3600;
+T=26.2;//Tension
+m=0.51//mass
+theta=60*pi/180;
+//calculation
+Fd=T*cos(theta);
+disp(Fd,"Drag force in N:");
+A=pi*d^2/4;
+Ad=A*cos(theta);//area component to drag
+Cd=2*Fd/Ad/rho/v^2;//coeff of drag
+disp(Cd , "The drag coefficient:")
+Fg=m*g;//force of gravity
+Fb=rho*pi*d^2/4*t*g;//buoyant force
+Fl=Fg-Fb+T*sin(theta);
+disp(Fl , "The lift force in N :");
+Al=A*sin(theta);
+Cl=2*Fl/Al/rho/v^2;
+disp(Cl ,"The coefficient of lift:")
diff --git a/1379/CH9/EX9.1.3/example9_3.sce b/1379/CH9/EX9.1.3/example9_3.sce
new file mode 100755
index 000000000..b66ec9b82
--- /dev/null
+++ b/1379/CH9/EX9.1.3/example9_3.sce
@@ -0,0 +1,16 @@
+
+
+//exapple 9.3 
+clc; funcprot(0);
+// Initialization of Variable
+rhog=1200;//density of glycerol
+mu=1.45;
+pi=3.1428;
+g=9.81;
+rhos=2280;//density of sphere
+v=0.04;//terminal velocity;
+a=2*mu*g*(rhos-rhog)/v^3/3/rhog^2;//a=Cd/2/Re
+//using graph of Cd/2/Re vs Re
+Re=0.32;
+d=Re*mu/v/rhog;
+disp(d , "Diameter of sphere in (m):");
diff --git a/1379/CH9/EX9.1.4/example9_4.sce b/1379/CH9/EX9.1.4/example9_4.sce
new file mode 100755
index 000000000..1485c18e6
--- /dev/null
+++ b/1379/CH9/EX9.1.4/example9_4.sce
@@ -0,0 +1,16 @@
+
+
+//exapple 9.4 
+clc; funcprot(0);
+// Initialization of Variable
+rhoa=1.218;//density of air
+mu=1.73/100000;
+pi=3.1428;
+g=9.81;
+rhop=2280;//density of polythene
+d=0.0034;//diameter
+a=4*d^3*(rhop-rhoa)*rhoa*g/3/mu^2;//a=Cd*Re^2
+//using graph of Cd*Re^2 vs Re
+Re=2200;
+v=Re*mu/d/rhog;
+disp(v , "The terminal vrlocity in (m/s)");
diff --git a/1379/CH9/EX9.1.5/example9_5.sce b/1379/CH9/EX9.1.5/example9_5.sce
new file mode 100755
index 000000000..7b3009a8a
--- /dev/null
+++ b/1379/CH9/EX9.1.5/example9_5.sce
@@ -0,0 +1,30 @@
+
+
+//exapple 9.2 
+clc; funcprot(0);
+// Initialization of Variable
+pi=3.1428;
+rho=825;
+mu=1.21;
+g=9.81;
+l=0.02;
+de=0.02;//dia exterior
+di=0.012;//dia interior
+//calculation
+//part 1
+zi=pi*(6*(pi*de^2/4-pi*di^2/4)*l/pi)^(2/3)/(pi*l*(di+de)+2*pi*(de^2/4-di^2/4));
+disp(zi, "sphericity of Raschig ring is:");
+//part 2
+u=0.04;
+ds=0.003//diameter of each sphere
+zi=pi*(6*pi*ds^3/pi)^(2/3)/6/pi/ds^2;//sphericity
+disp(zi, "sphericity of given object is:");
+Ap=4*ds^2-4*3/4*(ds^2-pi*ds^2/4);//projected area
+dp=sqrt(4*Ap/pi);//projected dia
+Re=dp*u*rho/mu;
+disp(Re, "Reynolds no. for the object:");
+//using graph b/w Re and zi and Cd
+Cd=105;//coeff. of drag
+Fd=Ap*Cd*u^2*rho/2;
+disp(Fd,"The drag force on object in (N):")
+
diff --git a/1379/CH9/EX9.1.6/example9_6.sce b/1379/CH9/EX9.1.6/example9_6.sce
new file mode 100755
index 000000000..5c80aeddf
--- /dev/null
+++ b/1379/CH9/EX9.1.6/example9_6.sce
@@ -0,0 +1,28 @@
+
+
+//exapple 9.6 
+clc; funcprot(0);
+// Initialization of Variable
+rho=998;//density of water
+mu=1.25/1000;//viscosity of water
+w=100;//mass of water
+pi=3.1428;
+g=9.81;
+rhog=2280;//density of glass
+wg=60;//mass of glass
+d=45*10^-6;//diameter of glass sphere
+//claculation
+rhom=(w+wg)/(w/rho+wg/rhog);//density of mixure
+e=w/rho/(w/rho+wg/rhog);//volume fraction of watter
+//using charts
+zi=exp(-4.19*(1-e));
+
+K=d*(g*rho*(rhog-rho)*zi^2/mu^2)^(1/3);//stoke's law coeff.
+disp(K);
+if K<3.3 then
+    disp("settling occurs in stoke-s law range");
+    U=g*d^2*e*zi*(rhog-rhom)/18/mu;
+    disp(U,"settling velocity in m/s:")
+else
+    disp("settling does not occurs in stoke-s law range");
+end
diff --git a/1379/CH9/EX9.1.7/example9_7.sce b/1379/CH9/EX9.1.7/example9_7.sce
new file mode 100755
index 000000000..50f11f30f
--- /dev/null
+++ b/1379/CH9/EX9.1.7/example9_7.sce
@@ -0,0 +1,24 @@
+
+
+//exapple 9.7 
+clc; funcprot(0);
+// Initialization of Variable
+rhog=1200;//density of glycerol
+mu=1.45;//viscosity of glycerol
+pi=3.1428;
+g=9.81;
+rhos=2280;//density of sphere
+d=8/1000;
+s=0;
+uf=0.8*0.026;
+//calculation
+function[a]=intre()
+    u=linspace(0,uf,1000);
+    for i=1:1000
+        y=((pi/6*d^3*rhos*g-pi*d^3/6*rhog*g-0.5*pi*d^2/4*24*mu/d/rhog*rhog*u(i))/pi*6/d^3/rhos)^(-1)*uf/1000;
+        s=s+y;
+    end
+    a=s;
+endfunction
+[t]=intre();
+disp(t,"Time taken by particle to reach 80% of its velocity in (s):");