diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /542/CH9 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '542/CH9')
| -rwxr-xr-x | 542/CH9/EX9.1/Example_9_1.sci | 19 | ||||
| -rwxr-xr-x | 542/CH9/EX9.3/Example_9_3.sci | 37 |
2 files changed, 56 insertions, 0 deletions
diff --git a/542/CH9/EX9.1/Example_9_1.sci b/542/CH9/EX9.1/Example_9_1.sci new file mode 100755 index 000000000..0e380f822 --- /dev/null +++ b/542/CH9/EX9.1/Example_9_1.sci @@ -0,0 +1,19 @@ +clear;
+clc;
+printf("\n Example 9.1");
+d_particle = 5; //particle size is in um
+p = 1000; //density of water in kg/m^3
+ps = 2800; //density of solids in kg/m^3
+viscosity = 10^(-3); //viscosity is in Ns/m^2
+uo = ((d_particle*10^(-6))^2)*(ps-p)*9.81/(18*viscosity);
+printf("\n Terminal falling velocity of particles of diameter = %.2f m/sec",uo*10^5);
+Q = 0.25; //volumetric flow rate is in m^3/sec
+printf("\n E = %.2f*10^(4) m^2",(Q/uo)*10^(-4));
+
+printf("\n For coal-in-oil mixture");
+uo1 = 0.04/(Q/uo);
+printf("\n uo = %.2f*10^-6 m/sec",uo1*10^6);
+
+d = sqrt((18*0.01*uo1)/((1300-850)*9.81));
+printf("\n d = %d um",(d/3)*10^6);
+
diff --git a/542/CH9/EX9.3/Example_9_3.sci b/542/CH9/EX9.3/Example_9_3.sci new file mode 100755 index 000000000..83b3e3310 --- /dev/null +++ b/542/CH9/EX9.3/Example_9_3.sci @@ -0,0 +1,37 @@ +clear;
+clc;
+printf("\n Example 9.3");
+//In the filter press
+// V^2 + 2(AL/v)V = 2(-deltaP)A^2*t/(ruv)
+
+l = 0.025; //l is in meters
+L = 0.003; //L is in meters
+deltaP = 350; //it is in N/m^2
+t = 3600; //t is in secs
+
+// x = v/ru
+x = poly([0],'x');
+x1 = roots(0.025^2 + 2*0.003*0.025 - 2*3.5*10^(5)*3600*x);
+printf("\n the value of ru/v = %.2f*10^12",(1/x1)*10^(-12));
+
+//In the centrifuge
+R = 0.15; //R is in meters
+H = 0.20; //H is in meters
+V = 0.00225; //V is in m^3
+r = poly([0],'r');
+r1 = roots(%pi*(R^2 - r^2)*H-V);
+printf("\n Value of ro = %f mm",r1(1)/2);
+printf("\n angular frequency = %.1frad/s",(r1(1)/2*10^3)*2*(%pi));
+
+//(R^2 -r^2)(1+2L/R)+2r^2ln(r/R) = 2vtpw^2/ru(R^2-ro^2)
+t = poly([0],'t');
+t1 = roots((R^2 - r1(1)^2)*(1+2*(L/R))+2*(r1(1)^2)*log(r1(1)/R)-2*t*1000*408.4^(2)/(3.25*10^12)*(R^2-(r1(1)/2)^2));
+printf(" \n time required = %f secs",t1);
+
+
+
+
+
+
+
+
|