initial commit / add all books

8 files changed, 212 insertions, 0 deletions

diff --git a/1052/CH23/EX23.1/231.sce b/1052/CH23/EX23.1/231.sce
new file mode 100755
index 000000000..d2170661e
--- /dev/null
+++ b/1052/CH23/EX23.1/231.sce
@@ -0,0 +1,11 @@
+clc;

+//Example 23.1 page no 323

+printf("Example 23.1 page no 323\n\n");

+//calculation of aerodynamic diameter for the following particles

+d_es=1.4//equivalent dia of solid sphere,micrometer

+sg_s=2//specific gravity of solid sphere

+d_eh=2.8//equivalent diameter of hollow sphere, mirometer

+sg_h=0.51//specific gravity of hollow sphere

+d_pa1=d_es*sqrt(sg_s)//aerodynamic dia for solid sphere

+d_pa2=round(d_eh*sqrt(sg_h))//aerodynamic dia for hollow sphere

+printf("\n d_pa1=%f micron\nd_pa2=%f micron",d_pa1,d_pa2);

diff --git a/1052/CH23/EX23.2/232.sce b/1052/CH23/EX23.2/232.sce
new file mode 100755
index 000000000..1782871f5
--- /dev/null
+++ b/1052/CH23/EX23.2/232.sce
@@ -0,0 +1,8 @@
+clc;

+//Example 23.2 page no 323

+printf("Example 23.2 page no 323\n\n");

+//calculation of aerodynamic diameter of irregular saped sphere

+d_e=1.3//eq. diameter,micron

+sg=2.35

+d_pa=d_e*sqrt(sg)//aerodynamic diameter

+printf("\n aerodynamic diameter d_pa=%f micron",d_pa);

diff --git a/1052/CH23/EX23.3/233.sce b/1052/CH23/EX23.3/233.sce
new file mode 100755
index 000000000..b393cf7c7
--- /dev/null
+++ b/1052/CH23/EX23.3/233.sce
@@ -0,0 +1,9 @@
+clc;

+//Example 23.3 page no 335

+printf("Example 23.3 page no 335\n\n");

+//calculation of cunningham correction factor 

+dp=0.4//particle diameter

+lemda=6.53e-2

+A=1.257 + 0.40*exp(-1.10*dp/(2*lemda))

+C= 1 + 2*A*lemda/dp//cunningham correction factor(CCF)

+printf("CCF C=%f ",C);

diff --git a/1052/CH23/EX23.4/234.sce b/1052/CH23/EX23.4/234.sce
new file mode 100755
index 000000000..367df8b74
--- /dev/null
+++ b/1052/CH23/EX23.4/234.sce
@@ -0,0 +1,50 @@
+clc;

+//Example 23.4

+//page no 336

+printf("Example 23.4 page no 336\n\n");

+//three different diameter sized fly ash particls settle through air 

+//we have to calculate the particle terminal velocity and determine how far each will fall in 30 seconds 

+//assume the particles are speherical 

+SG=2.31//specific gravity of fly ash

+rho_w=62.4//density of water 

+rho_p=SG*rho_w//density of particles

+//properties of air

+R=0.7302//gas constant

+T=698//temperature,R

+P=1//pressure ,atm

+Mw=29//mol. wt of air

+rho_a=P*Mw/(R*T)//density of air,lb/ft^3

+meu=1.41e-5//viscosity of air,lb/ft.s

+g=32.174//grav. acc

+D1=0.4//diameter of particle 1,microns

+D2=40//diameter of particle 2,microns

+D3=400//diameter of particle 3,microns

+K1=(D1/(25400*12))*(g*rho_p*rho_a/(meu^2))^(1/3)//dimensionless constant for particle 1

+K2=(D2/(25400*12))*(g*rho_p*rho_a/(meu^2))^(1/3)//dimensionless constant for particle 2

+K3=(D3/(25400*12))*(g*rho_p*rho_a/(meu^2))^(1/3)//dimensionless constant for particle 3

+printf("\n dimensionless constant K1=%f \n K2=%f \n K3=%f ",K1,K2,K3);

+//first we determine which fluid particle dynamic law applies for the above values of K

+//for particle 1,strokes law applies

+//for particle 2,strokes law applies

+//for particle 3,intermediate law applies

+//terminal settling velocity for each particle

+v1=(D1/(25400*12))^2*g*rho_p/(18*meu)

+printf("\n terminal settling velocity for particle 1 v1=%f ft/s",v1);

+v2=(D2/(25400*12))^2*g*rho_p/(18*meu)

+printf("\n terminal settling velocity v2=%f ft/s",v2);

+v3=(D3/(25400*12))^1.14*0.153*g^0.71*rho_p^0.71/(rho_a^0.29*meu^0.43)

+printf("\n terminal settling velocity v3=%f ft/s ",v3);

+//calculation of how far x,the fly ash particles will fall in 30 seconds

+t=30//time,sec

+x2=v2*t//distance travel by 2 particle

+x3=v3*t//distance travel by 3 particle

+printf("\n distance by 2 particle x2=%f ft\n distance by 3 particle x3=%f ft",x2,x3);

+//for 1 particle K1 and v1 value are without the CCF.With the correction factor lemda=6.53e-8,gives

+lemda=6.53e-8//correction factor

+y=-1.10*(D1/(25400*12))/(2*lemda)

+A =1.257 + 0.40*exp(y)

+C=1 + 2*A*lemda/(D1/(25400*12))//cunningham correction factor(ccf)

+//now equation 23.36 can be employed 

+v1_corrected=v1*C//corrected velocity of 1 particle

+x1=v1_corrected*t//distance travel by 1 particle

+printf("\n distance travel by 1 particle x1=%f ft",x1); 

diff --git a/1052/CH23/EX23.5/235.sce b/1052/CH23/EX23.5/235.sce
new file mode 100755
index 000000000..df2750757
--- /dev/null
+++ b/1052/CH23/EX23.5/235.sce
@@ -0,0 +1,24 @@
+clc;

+//Example 23.5

+//page no 338

+printf("\n Example 23.5 page no 338\n\n");

+//refer to example 23.5

+//we have to calculate size of a flyash particle that will settle with a velocity of 1.384 ft/s

+SG=2.31//specific gravity of fly ash

+rho_w=62.4//density of water 

+rho_p=SG*rho_w//density of particles

+//properties of air

+R=0.7302//gas constant

+T=698//temperature,R

+P=1//pressure ,atm

+Mw=29//mol. wt of air

+rho_a=P*Mw/(R*T)//density of air,lb/ft^3

+meu=1.41e-5//viscosity of air,lb/ft.s

+g=32.174//grav. acc

+v=1.384//velocity at which particle settle down,ft/s

+W= v^3*rho_a^2/(g*rho_p*meu)//dimensionless constant

+printf("\n dimensionless constant W=%f ",W);

+//since W < 0.2222 stokes' law applies

+D_p=sqrt(18*meu*v/(g*rho_p))//diameter of particle

+printf("\n diameter of particle D_p=%f ft",D_p);

+  

diff --git a/1052/CH23/EX23.7/237.sce b/1052/CH23/EX23.7/237.sce
new file mode 100755
index 000000000..b7d224c84
--- /dev/null
+++ b/1052/CH23/EX23.7/237.sce
@@ -0,0 +1,50 @@
+clc;

+//Example 23.7

+//page no 340

+printf("\n Example 23.7 page no 340\n\n");

+// In a plant manufacturing ivory soap detergent explodes one windy day

+//we have to calculate the distance from the plant where the soap particles will start to deposit and where they will cease to deposit

+//the smallest particle wll travel the greatest distance while the largest will travel the least distance 

+//for the minimumdistance ,we use largest particle 

+D_l=3.28e-3//largest diameter,ft

+g=32.174//grav. acc.

+SG=0.8//specific gravity of soap particle

+rho_w=62.4

+rho_p=SG*rho_w//density of particle

+rho_a=0.0752//density of given atmosphere,lb/ft^3

+meu=1.18e-5//viscosity 

+K_l = D_l*(g*(rho_p-rho_a)*rho_p/(meu^2))^(1/3)//dimensionless constant

+printf("\n dimensionless constant K_l=%f ",K_l);

+//value of K indicates the intermediate range applies 

+//the settling velocity is given by 

+v_l=0.153*g^0.71*D_l^1.14*rho_p^0.71/(meu^0.43*rho_a^0.29)

+printf("\n settling velocity v_l=%f ft/s",v_l);

+H=400//vertical height blowen by particle,ft

+t_l=H/v_l//descent time

+v_w=20//wind velocity in miles/h

+L=t_l*v_w*(5280/3600)//horizontal distance travelled by particles

+printf("\n descent time t_l=%f second\n horizontal distance L=%f ft",t_l,L);

+//for the minimum distance we use smallest particle

+D_s=6.89e-6//diameter of smallest particle,ft

+K_s=D_s*(g*(rho_p-rho_a)*rho_a/(meu^2))^(1/3)

+printf("\n dimensionless constant K_s=%f ",K_s);

+//velocity is in the stokes regime and is given by

+v_s=g*D_s^2*rho_p/(18*meu)

+printf("\n settling velocity v_s=%f ft/s",v_s);

+t_s=H/v_s//descent time 

+L_s=t_s*v_w*(5280/3600)//horizontal distance travelled 

+printf("\n descent time t_s=%f s\nhorizontal distance travelled by smallest particle L_s=%f ft",t_s,L_s);

+m=100*2000//mass of particles

+V_act=m/rho_p//actual volume of particles

+e=0.5//void fraction

+V_b=V_act/e//bulk volume

+printf("\ actual volume V_act=%f ft^3\nbulk volume V_b=%f ",V_act,V_b);

+L_d=L_s-L//length of drop area

+printf("\n L_d=%f ",L_d);

+W=100//width ,ft

+A_d=L_d*W//deposition area

+H_d=V_b/A_d//deposition height

+printf("\n deposition height H_d=%f ft",H_d);

+//deposition height can be ,at bestt, described asa sprinkling 

+

+

diff --git a/1052/CH23/EX23.8/238.sce b/1052/CH23/EX23.8/238.sce
new file mode 100755
index 000000000..6e9f4fec9
--- /dev/null
+++ b/1052/CH23/EX23.8/238.sce
@@ -0,0 +1,42 @@
+clc;

+//Example 23.8

+//page no 342

+printf("Example 23.8 page no 342\n\n");

+//a small sphere is observed to fall through caster oil 

+v_t=0.042//terminal velocity of particle 

+meu_f=0.9//viscosity of oil

+rho_f=970//density of oil

+g=9.807//grav. acc.

+D_p=0.006//diameter of particle

+rho_p=(18*meu_f*v_t)/(g*D_p^2)  + rho_f

+printf("\n density of particle rho_p=%f kg/m^3",rho_p);

+neu_f=9.28e-4//dynamic viscosity of fluid

+R_e=D_p*v_t/neu_f//reynolds no

+printf("\n reynolds no R_e=%f ",R_e);

+//since R_e < 0.3

+//calculation of the settling criterion factor ,K

+K=D_p*(g*rho_f*(rho_p-rho_f)/(meu_f^2))^(1/3)//the settling criterion factor

+printf("\n K=%f ",K);

+//since K <3.3, stokes law applies 

+//the drag coeff. C_d 

+C_d=24/R_e

+printf("\n drag coeff C_d=%f ",C_d);

+F_d=3*%pi*meu_f*D_p*v_t//drag force

+printf("\n drag force F_d=%f N",F_d);

+F_b=(%pi/6)*D_p^3*rho_f*g//buoyancy force 

+printf("\n buoyancy force F_b=%f N",F_b);

+//Consider the case when same sphere is dropped in water 

+rho_w=1000//density of water,kg/m^3

+meu_w=0.001//viscosity of water,kg/m.s

+//the particle will move faster because of the lower viscosity of water ,stokes law will almost definietly not apply

+K_w=D_p*(g*rho_w*(rho_p-rho_w)/(meu_w^2))^(1/3)//the settling criterion factor

+printf("\n k_w settling factor =%f ",K_w);

+//since K_w = 158 > 43.6,the flow is in the Newton's law regime 

+//employ eq. 23.31 but include the (buoyant) density ratio factor

+v_t_w=1.75*sqrt((rho_p-rho_w)/(rho_w)*g*D_p)//terminal velocity 

+printf("\n terminal velocity in water v_t_w=%f m/s",v_t_w);  

+

+

+

+

+ 

diff --git a/1052/CH23/EX23.9/239.sce b/1052/CH23/EX23.9/239.sce
new file mode 100755
index 000000000..daf1e3f19
--- /dev/null
+++ b/1052/CH23/EX23.9/239.sce
@@ -0,0 +1,18 @@
+clc;

+//Example 23.9

+//page no 344

+printf("Example 23.9 page no 344\n\n");

+//the bottom of a ship,moving in water

+rho=1000//density of water

+v=12//velocity of boat,m/s

+L=20//length,m

+W=5//width ,m

+meu=1e-3//viscosity

+R_e=rho*v*L/meu//reynolds no

+printf("Reynolds no R_e=%f ",R_e);

+//from reynolds no  flow is turbulent

+C_d=0.031/(R_e^(1/7))//coeff. discharge\

+printf("\ncoeff. discharge C_d=%f ",C_d);

+//calculation of the drag on area LW

+F_d=(1/2)*C_d*rho*v^2*L*W//drag force

+printf("\n drag force F_d=%f N",F_d);