7 files changed, 252 insertions, 0 deletions
diff --git a/542/CH1/EX1.1/Example_1_1.sce b/542/CH1/EX1.1/Example_1_1.sce
new file mode 100755
index 000000000..51fa689c1
--- /dev/null
+++ b/542/CH1/EX1.1/Example_1_1.sce
@@ -0,0 +1,18 @@
+clear; 
+clc;
+printf("\n Example 1.1");
+//Given size analysis of a powdered material
+d=[1,101];//diameter of the powdered particles
+x=[0,1];//mass fractions of the particles
+plot2d(d,x,style=2,rect=[0,0,120,1])
+xtitle("size analysis of powder","particle size(um)","mass fraction(x)")
+d=100*x+1; // from the given plot
+//calculation of surface mean diameter
+function[ds]=surface_mean_diameter(x0,x1)
+    ds=1/(integrate('1/(100*x+1)','x',x0,x1))
+    funcprot(0)
+endfunction
+ds=surface_mean_diameter(0,1);//deduced surface mean diameter according to def.
+printf("\n The surface mean diameter is %fum",ds);
+
+ 
+\ No newline at end of file
diff --git a/542/CH1/EX1.2/Example_1_2.sce b/542/CH1/EX1.2/Example_1_2.sce
new file mode 100755
index 000000000..efda889b0
--- /dev/null
+++ b/542/CH1/EX1.2/Example_1_2.sce
@@ -0,0 +1,87 @@
+clear;
+clc;
+printf("\n Example 1.2");
+//from given differential eq we get these functions
+//particle number distribution for the size range 0-10um
+
+
+//n=0.5*d^2;
+//const of integration is0 since at n=0,d=0
+
+//particle number distribution for the size range 10-100um
+//n=83-(0.33*(10^(5))*d^(-3))
+//c2=83,since at d=10um,n=50
+
+//number distribution plot for the powdered material of size range 0-100um
+function[n]= number_distribution(d)
+    if(d<=10) then
+        n=0.5*d^2;
+    else
+        n=83-(0.33*(10^(5))*d^(-3));
+        end
+        funcprot(0)
+endfunction
+d=0;
+while(d<=100)
+  n=number_distribution(d);
+    plot(d,n,"+-");
+    d=d+1; 
+end 
+xtitle("number_distribution_plot","diameter(um)","number distribution");
+ps=[0 6.2 9.0 10.0 11.4 12.1 13.6 14.7 16.0 17.5 19.7 22.7 25.5 31.5 100];
+function[n1]=difference(i)
+//ps=[0 6.2 9.0 10.0 11.4 12.1 13.6 14.7 16.0 17.5 19.7 22.7 25.5 31.5 10];
+//according to the given particle sizes particle sizes are in um
+    n1=number_distribution(ps(i+1))-number_distribution(ps(i));
+    funcprot(0);
+endfunction
+function[da]=average(i)
+    da= (ps(i+1)+ps(i))/2;
+    funcprot(0);
+endfunction 
+tot_n1d12=0;
+tot_n1d13=0;
+i=1;
+for i=1:14 
+        tot_n1d12=tot_n1d12+difference(i)*(average(i))^2;
+        tot_n1d13=tot_n1d13+difference(i)*(average(i))^3;
+end
+printf("\n tot_n1d12 =%d \n tot_n1d13=%d",tot_n1d12,tot_n1d13);
+function[s]=surface_area(j)
+    s=(difference(j)*(average(j))^2)/tot_n1d12;
+    funcprot(0);
+endfunction
+su=0;
+j=0;
+xset('window',1);
+
+plot(0,0,"o-");
+for j=1:14
+    su=su+surface_area(j);
+    plot(ps(j+1),su,"o-");
+end
+xtitle("surface area and mass distribution plot","diameter(um)","surface area or mass distribution");
+//mass distribution plot
+function[x]=mass_distribution(k)
+    x=(difference(k)*(average(k))^3)/tot_n1d13;
+    funcprot(0);
+endfunction
+ma=0;
+k=0;
+plot(0,0,"+-");
+for k=1:14
+    ma=ma+mass_distribution(k);
+    plot(ps(k+1),ma,"+-");
+end
+//evaluating surface mean diameter
+function[d]=surface_mean_diameter(l)
+    e=0;
+    for l=1:14
+       n=(mass_distribution(l)/average(l));
+       e=e+n;
+    end
+d=1/e;
+    funcprot(0);
+endfunction
+printf("\nthe surface mean diameter is: %fum",surface_mean_diameter());
+    
diff --git a/542/CH1/EX1.3/Example_1_3.sci b/542/CH1/EX1.3/Example_1_3.sci
new file mode 100755
index 000000000..06c11b05f
--- /dev/null
+++ b/542/CH1/EX1.3/Example_1_3.sci
@@ -0,0 +1,30 @@
+clear all;
+clc;
+printf("\n Example 1.3");
+p=0.20;//components analysed represents 20 percent of the mixture by mass
+//for a completely unmixed system
+so=p*(1-p);
+//for a completely random mixture :
+n=100;//Each of the sample removed contains 100 particles
+sr=p*(1-p)/n;
+s=[0.025 0.006 0.015 0.018 0.019];
+time_secs=[30 60 90 120 150];
+printf("\n degree of mixing is :\n")
+function[b]=degree_of_mixing()
+for i=1:5
+    b(i)=(so-s(i))/(so-sr);
+    disp(b(i));//b is the degree of mixing
+end
+    return b;
+funcprot(0)
+endfunction
+plot2d(time_secs,degree_of_mixing(),style=3)
+xtitle("degree of mixing curve","time_secs","degree_of_mixing")
+//plot of sample variance vs time(secs)
+xset('window',1)
+plot2d(time_secs,s,style=2)
+xtitle("sample variance curve","time_secs","sample variance")
+//from the graph the maxima is at 60 secs
+
+
+
diff --git a/542/CH1/EX1.4/Example_1_4.sci b/542/CH1/EX1.4/Example_1_4.sci
new file mode 100755
index 000000000..6ad45e484
--- /dev/null
+++ b/542/CH1/EX1.4/Example_1_4.sci
@@ -0,0 +1,30 @@
+//minimum size of the particle in the mixture of quartz and galena(mm)
+clear all;
+clc;
+printf("\n Example 1.4");
+
+//maximum size of the particle(mm)
+d_max=0.065;
+//minimum size of the particle(mm)
+d_min=0.015;
+//density of quartz(kg/m^3)
+p_quartz=2650;
+//density of galena (kg/m^3)
+p_galena=7500;
+//minimum density of the particle which will give this seperation
+//When stoke's law is applied the required density is as given below
+function[d]=stoke_required_density()
+    p=poly([0],'p');
+    d=roots((p-7500)-(p-2650)*(d_max/d_min)^2);
+    funcprot(0);
+endfunction
+d=stoke_required_density();
+printf("\n required density is = %d kg/m^3",d);
+//When Newton's law is applied then the required density is as given below
+function[e]=newton_required_density()
+    r=poly([0],'r');
+    e=roots((r-7500)-(r-2650)*(d_max/d_min));
+    funcprot(0);
+endfunction
+e=newton_required_density();
+printf("\nrequired density is by newton law =%d kg/m^3",e);
+\ No newline at end of file
diff --git a/542/CH1/EX1.5/Example_1_5.sce b/542/CH1/EX1.5/Example_1_5.sce
new file mode 100755
index 000000000..36b935f56
--- /dev/null
+++ b/542/CH1/EX1.5/Example_1_5.sce
@@ -0,0 +1,18 @@
+clear;
+clc;
+printf("\n Example 1.5");
+//efficiency of the collector for different size ranges
+efficiency_1=45;//in percentage for the size range of 0-5um
+efficiency_2=80;//in percentage for the size range of 5-10um
+efficiency_3=96;//in percentage for the size range greater than 10um
+
+//mass percent of the ndust for various size range
+mass_1=50; //in percentage for the size range of 0-5um
+mass_2=30; //in percetage for the size range of 5-10um
+mass_3=20; //in percentage for the size range greater than 10um
+// on the basis of 100kg dust
+mass_retained_1=0.45*50;//mass retained(kg) in the size range of 0-5um
+mass_retained_2=0.80*30;//mass retained(kg) in the size range of 5-10um
+mass_retained_3=0.96*20;//mass retained(kg) in the size range greater than10um
+overall_efficiency=0.45*50+0.80*30+0.96*20;
+printf("\n the overall efficiency is =%f",overall_efficiency);
diff --git a/542/CH1/EX1.6/Example_1_6.sce b/542/CH1/EX1.6/Example_1_6.sce
new file mode 100755
index 000000000..79d69bf1f
--- /dev/null
+++ b/542/CH1/EX1.6/Example_1_6.sce
@@ -0,0 +1,35 @@
+clear;
+clc;
+printf("\n Example 1.6");
+//To calculate mass flow of the dust emitted
+mass_1=10;//in percentage in the size range of 0-5um
+mass_2=15;//in percentage in the size range of 5-10um
+mass_3=35;//in percentage in the size range of 10-20um
+mass_4=20;//in percentage in the size range of 20-40um
+mass_5=10;//in percentage in the size range of 40-80um
+mass_6=10;//in percentage in the size range of 80-160um
+efficiecny_1=20;//in percentage in the size range of 0-5um
+efficiency_2=40;//in percentage in the size range of 5-10um
+efficiency_3=80;//in percentage in the size range of 10-20um
+efficiency_4=90;//in percentage in the size range of 20-40um
+efficiency_5=95;//in percentage in the size range of 40-80um
+efficiency_6=100;//in percentage in the size range of 80-160um
+dust_burden=18;//in g/m^3 at the entrance
+//taking 1m^3 as the basis of calculation
+total_mass_retained=18*(0.1*0.20+0.15*0.40+0.35*0.80+0.2*0.9+0.1*0.95+0.1*1);
+printf("\ntotal mass retained =%fg",total_mass_retained);
+total_efficiency=(total_mass_retained/18)*100;
+printf("\ntotal efficiency is =%f",total_efficiency);
+total_mass_emitted=18-total_mass_retained;
+printf("total mass emitted is:%fg",total_mass_emitted);
+t=18*(0.1*0.80+0.15*0.60+0.35*0.20);
+printf("\ntotal mass emitted less than 20um is %fg",t);
+e=t*100/total_mass_emitted;
+printf("\nThe efficiency of particles emitted is %f",e);
+//gas flow is 0.3m^3/sec
+f=0.3*total_mass_emitted;
+printf("\nmass flow rate is:%fkg/sec",f);
+ 
+
+
+
diff --git a/542/CH1/EX1.7/Example_1_7.sce b/542/CH1/EX1.7/Example_1_7.sce
new file mode 100755
index 000000000..233d03d0d
--- /dev/null
+++ b/542/CH1/EX1.7/Example_1_7.sce
@@ -0,0 +1,34 @@
+clear;
+clc;
+printf("\n Example 1.7");
+Ai=(%pi/4)*(0.075)^2;//cross sectional area at the gas inlet in m^2
+do=0.075;//gas outlet diameter in m
+p=1.3;//gas density in kg/m^3
+Z=1.2;//height of the seperator in m
+dt=0.3;//seperator diameter in m
+v=1.5;//gas entry velocity in m/sec
+G=(Ai*v*p);//mass flow rate of the gas in kg/sec
+printf("\n cross sectional area at the gas inlet is %fm^2",Ai);
+printf("\ngas outlet diameter is %fm",do);
+printf("\ngas density is %fkg/m^3",p);
+printf("\nheight of the seperator is %fm",Z);
+printf("\nseperator diameter is%fm",dt);
+printf("\nmass flow rate of the gas is %fkg/sec",G);
+function[u]=terminal_vel()
+    u=0.2*(Ai)^2*(do)*p*9.8/(%pi*Z*(dt)*G);//velocity is in m/sec
+    funcprot(0);
+endfunction
+u=terminal_vel();
+printf("\nthe terminal velocity of the smallest particle retained by the seperator =%fm/sec",u);
+function[d]=particle_diameter(u)
+    u=terminal_vel();
+    n=0.018*10^(-3);//viscosity in mNs/m^2
+    ps=2700;//density of the particle in kg/m^3
+    d=((u*18*n)/(9.8*(ps-p)))^(0.5);//particle size in um
+    funcprot(0);
+endfunction
+u=terminal_vel();
+d=particle_diameter(u);
+do=d*10^6;
+printf("\n particle diameter by the stoke law is %fum",do);
+