7 files changed, 247 insertions, 0 deletions
diff --git a/1379/CH10/EX10.1.1/example10_1.sce b/1379/CH10/EX10.1.1/example10_1.sce
new file mode 100755
index 000000000..3a3ff7d1e
--- /dev/null
+++ b/1379/CH10/EX10.1.1/example10_1.sce
@@ -0,0 +1,28 @@
+
+
+//example 10.1
+clc; funcprot(0);
+clf()
+//exapple 10.1 
+// Initialization of Variable
+t=[0 0.5 1 2 3 4 5 6 7 8 9 10];//time
+h=[1.10 1.03 .96 .82 .68 .54 .42 .35 .31 .28 .27 .27];
+Cl=[0 0 0 0 0 0 0 0 0 0 0];
+m=0.05;
+V=1/1000;//volume
+//calculations
+Co=m/V;//concentration at t=0
+v(1)=(h(1)-h(2))/(t(2)-t(1));
+Cl(1)=Co;
+for i=2:11
+  
+        v(i)=(h(i-1)-h(i+1))/(t(i+1)-t(i-1));//slope or settling velocity
+        Cl(i)=Co*h(1)/(h(i)+v(i)*t(i));
+
+    
+end
+plot(t,h,'r--d');
+clf();
+plot(Cl,v,'r->');
+xtitle("Concentration vs Settling veocity" , "Concentration(kg/m^3)" , "Settling velocity (m/h)");
+
diff --git a/1379/CH10/EX10.1.2/example10_2.sce b/1379/CH10/EX10.1.2/example10_2.sce
new file mode 100755
index 000000000..582813ef8
--- /dev/null
+++ b/1379/CH10/EX10.1.2/example10_2.sce
@@ -0,0 +1,31 @@
+
+
+//example 10.2
+clc; funcprot(0);
+clf()
+//exapple 10.2 
+// Initialization of Variable
+t=[0 0.5 1 2 3 4 5 6 7 8 9 10];//time
+h=[1.10 1.03 .96 .82 .68 .54 .42 .35 .31 .28 .27 .27];
+Cl=50:5:100;
+U=[19.53 17.71 16.20 14.92 13.82 12.87 12.04 11.31 10.65 9.55];//mass ratio of liquid to solid
+v=[0.139 0.115 0.098 0.083 0.071 0.062 0.055 0.049 0.043 0.034];//terminal velocity
+//above value taken from graph given with ques.
+C=130;//conc. of solids
+Q=0.06;//slurry rate
+Cmax=130//maximum solid conc.
+rhos=2300;//density of solid
+rho=998;//density of water
+V=rho*(1/C-1/rhos);
+F=Q*Cl(1)*3600;
+for i=1:10
+A(i)=F*(U(i)-V)/rho/v(i);
+end
+plot(v,A,'r-');
+xtitle("","Settling Velocity(m/h)", "Area(m^2)")
+//maxima finding using datatraveller in the graph
+disp(A,"the area for each settling velocity");
+disp("1005 m^2 is the maximum area found out from the plot");
+Qu=Q-F/3600/Cmax;
+disp(Qu, "Volumetric flow rate of clarified water in (m^3/s):")
+ 
diff --git a/1379/CH10/EX10.1.3/example10_3.sce b/1379/CH10/EX10.1.3/example10_3.sce
new file mode 100755
index 000000000..b0ebc539f
--- /dev/null
+++ b/1379/CH10/EX10.1.3/example10_3.sce
@@ -0,0 +1,58 @@
+
+
+//example 10.3
+clc; funcprot(0);
+//exapple 10.3 
+// Initialization of Variable
+rho1=2600;//density lighter
+rho2=5100;//density heavier
+pd1=0.000015:0.000010:0.000095;//particle diameter lighter
+pd2=0.000025:0.00001:0.000095;//particle diameter heavier
+wp1=[0 22 35 47 59 68 75 81 100];//weight distribution lighter
+wp2=[0 21 33.5 48 57.5 67 75 100];//weight distribution heavier
+rho=998.6;//density water
+mu=1.03/1000;//viscosity water
+g=9.81;
+u=0.004;//velocity of water
+d=95/1000000;//paeticle diameter maximum
+//calculation
+//part 1
+Re=d*u*rho/mu;
+d1=sqrt(18*mu*u/g/(rho1-rho));
+d2=sqrt(18*mu*u/g/(rho2-rho));
+function[a]=inter(d,f,g,b);//interpolation linear
+    for i=1:b
+        if d<=f(i+1)& d>f(i)  then
+            break
+        else 
+            continue
+        end
+        break
+    end
+    a=(d-f(i))/(f(i+1)-f(i))*(g(i+1)-g(i))+g(i);
+endfunction
+[a]=inter(d1,pd1,wp1,9);
+[b]=inter(d2,pd2,wp2,8);
+v2=1/(1+5)*100-b/100*1/(1+5)*100;
+v1=5/(1+5)*100-a/100*5/(1+5)*100;
+pl2=(v2)/(v2+v1);
+disp(pl2, "The fraction of heavy ore remained in bottom");
+ //part 2
+ rho=1500;
+ mu=6.25/10000;
+ a=log10(2*d^3*rho*g*(rho1-rho)*3*mu^2);//log10(Re^2(R/rho/mu^2))
+  //using value from chart(graph)
+Re=10^0.2136;
+u=Re*mu/rho/d;
+d2=sqrt(18*mu*u/g/(rho1-rho));
+[b]=inter(d2,pd2,wp2,8);
+disp(100-b+3.5,"The percentage of heavy ore left in this case");
+//part 3
+a=0.75//% of heavy ore in overhead product
+s=100*5/6/(100*5/6+0.75*100/6);
+disp(s,"the fraction of light ore in overhead product:");
+//part 4
+da=pd2(1);
+db=pd1(9);
+rho=(da^2*rho2-db^2*rho1)/(-db^2+da^2);
+  disp(rho,"The minimum density required to seperate 2 ores in kg/m^3:")
diff --git a/1379/CH10/EX10.1.4/example10_4.sce b/1379/CH10/EX10.1.4/example10_4.sce
new file mode 100755
index 000000000..5cf847dbb
--- /dev/null
+++ b/1379/CH10/EX10.1.4/example10_4.sce
@@ -0,0 +1,33 @@
+
+
+//example 10.4
+clc; funcprot(0);
+//exapple 10.4
+// Initialization of Variable
+rho=998;
+w0=40;//density of slurry
+mu=1.01/1000;
+g=9.81;
+rho1=2660;//density quartz
+h=0.25;
+t=18.5*60;
+mp=[5 11.8 20.2 24.2 28.5 37.6 61.8];
+d=[30.2 21.4 17.4 16.2 15.2 12.3 8.8]/1000000;
+u=h/t;
+d1=sqrt(18*mu*u/g/(rho1-rho));
+function[a]=inter(d,f,g,b);//interpolation linear
+    for i=1:b
+        if d>f(i+1)& d<=f(i)  then
+            break
+        else 
+            continue
+        end
+        break
+    end
+    
+    a=-(d-f(i+1))/(f(i)-f(i+1))*(g(i+1)-g(i))+g(i+1);
+endfunction
+[a]=inter(d1,d,mp,6);
+phi=1-a/100;
+rhot=phi*(rho1-rho)/rho1*w0+rho;
+disp(rhot,"the density of suspension at depth 25cm in kg/m^3 is")
diff --git a/1379/CH10/EX10.1.5/example10_5.sce b/1379/CH10/EX10.1.5/example10_5.sce
new file mode 100755
index 000000000..9f1d2a896
--- /dev/null
+++ b/1379/CH10/EX10.1.5/example10_5.sce
@@ -0,0 +1,34 @@
+
+
+//example 10.5
+clc; funcprot(0);
+clf()
+//exapple 10.5 
+// Initialization of Variable
+t=[0 45 135 495 1875 6900 66600 86400];//time
+m=[0.1911 0.1586 0.1388 0.1109 0.0805 0.0568 0.0372 0.0359];//mass total
+rho1=3100;//density of cement
+mu=1.2/1000;//viscosity of desperant liquid
+rho=790;//density of desperant liquid
+h=0.2;
+V=10;
+s=0;
+d(1)=100/1000000;//assumed value
+for i=1:7
+    d(i+1)=sqrt(18*mu*h/g/t(i+1)/(rho1-rho));//dia of particles
+    mc(i+1)=m(i+1)-0.2/100*V;//mass of cement
+    s=s+mc(i+1);  
+end
+mc(1)=m(1)-0.2*V/100;
+s=s+mc(1);
+mp(1)=100;
+for i=1:7
+    mp(i+1)=mc(i+1)/mc(1)*100;//mass percent below size
+end
+plot(mp,d);
+xtitle("", "%undersize", "Particle Size(m)");
+u=h/t(2);
+Re=d(2)*u*rho/mu;
+if Re<2 then
+    disp("since Re<2 for 81% of particles so settlement occurs mainly by stoke-s law")
+end
diff --git a/1379/CH10/EX10.1.6/example10_6.sce b/1379/CH10/EX10.1.6/example10_6.sce
new file mode 100755
index 000000000..b9c228dca
--- /dev/null
+++ b/1379/CH10/EX10.1.6/example10_6.sce
@@ -0,0 +1,38 @@
+
+
+//example 10.6
+clc; funcprot(0);
+//exapple 10.6
+clf() 
+// Initialization of Variable
+rho=998;
+rho1=2398;//density of ore
+mu=1.01/1000;
+g=9.81;
+h=25/100;
+t=[114 150 185 276 338 396 456 582 714 960];
+m=[0.1429 0.2010 0.2500 0.3564 0.4208 0.4781 0.5354 0.6139 0.6563 0.7277];
+for i=1:10
+ms=0.0573+m(10);//total mass setteled
+d(i)=sqrt(18*mu*h/g/(rho1-rho)/t(i));
+P(i)=m(i)/ms*100;//mass percent of sample
+end
+plot(t,P);
+xtitle("","Settling time (s)","mass percent in (%)");
+disp(P,d,"& its percentage mass distribution respectively" ,"the particle size distribution in (m)" );
+for i=2:9 
+        del(i)=(P(i+1)-P(i-1))/(t(i+1)-t(i-1));//slope 
+        W(i)=P(i)-t(i)*del(i);
+        W(1)=P(1)-P(1);
+           
+end
+W(10)=P(10)-t(10)*0.025;
+disp("mass% and diameter(m)respectively with serial no:")
+for i=4:10
+    disp(i-4);
+    disp("mass% is")
+    disp( "for diameter in(m) of",W(i));
+    disp(d(i));
+   
+end
+
diff --git a/1379/CH10/EX10.1.7/example10_7.sce b/1379/CH10/EX10.1.7/example10_7.sce
new file mode 100755
index 000000000..1a0709af9
--- /dev/null
+++ b/1379/CH10/EX10.1.7/example10_7.sce
@@ -0,0 +1,25 @@
+
+
+//example 10.7
+clc; funcprot(0);
+//exapple 10.7 
+// Initialization of Variable
+rho=1002;//density of disperant
+rho1=2240;//density of kaolin
+mu=1.01/1000;//viscosity
+g=9.81;
+t=600;
+h2=0.2;
+h1=0.4;
+dg=15*10^-6;//particle size to be removed
+//calculations
+//part 1
+d=sqrt(18*mu*h2/g/(rho1-rho)/t);
+x=dg/d;
+f=h2/h1*(1-x^2);//fraction separated after first decanting
+g=f*(1-f);
+disp(g,"fraction of particles separated after second decanting");
+disp(f+g,"total fraction of particles separated after decanting")
+//part 2
+h=(1-20/40*(1-x^2))^6;
+disp(h,"fraction of particles separated after sixth decanting");