initial commit / add all books

5 files changed, 150 insertions, 0 deletions

diff --git a/1379/CH12/EX12.1.1/example12_1.sce b/1379/CH12/EX12.1.1/example12_1.sce
new file mode 100755
index 000000000..dcb03c38a
--- /dev/null
+++ b/1379/CH12/EX12.1.1/example12_1.sce
@@ -0,0 +1,26 @@
+

+

+//example 12.1 

+clc; funcprot(0);

+// Initialization of Variable

+rho=1.22;

+pi=3.1428;

+rhos=518;

+rhoav=321;

+mu=1.73/10^5;

+g=9.81;

+d=0.65/1000;

+d2=25.5/100;//dia of duct

+ms=22.7/60;//mass flow rate

+//calculation

+e=(rhos-rhoav)/(rhos-rho);

+//coeff of quadratic eqn in U

+//a*x^2+b*x+c=0

+c=-(1-e)*(rhos-rho)*g;

+b=150*(1-e)^2*mu/d^2/e^3;

+a=1.75*(1-e)*rho/d/e^3;

+y=poly([c b a],'U','coeff');

+U=roots(y);

+Us=ms*4/pi/d2^2/rhos;//superficial speed

+Ua=e/e*(U(2)/e+Us/(1-e));

+disp(Ua,"the actual linear flow rate through duct in (m/s):")

diff --git a/1379/CH12/EX12.1.2/example12_2.sce b/1379/CH12/EX12.1.2/example12_2.sce
new file mode 100755
index 000000000..aa0097a1c
--- /dev/null
+++ b/1379/CH12/EX12.1.2/example12_2.sce
@@ -0,0 +1,30 @@
+

+

+//example 12.2 

+clc; funcprot(0);

+// Initialization of Variable

+rho=1.22;//density of air

+pi=3.1428;

+rhos=910;//density of polyethene

+d=3.4/1000;//dia of particles

+mu=1.73/10^5;

+g=9.81;

+dt=3.54/100;//dia of duct

+//calculation

+a=2*d^3*rho*g*(rhos-rho)/3/mu^2;

+disp(a,"R/rho/U^2*(Re^2)=");

+//using Chart

+Re=2*10^3;

+U=mu*Re/d/rho;

+b=U/(g*dt)^.5;

+if b>0.35 then

+    disp("choking can occur of this pipe system");

+else

+    disp("choking can not occur of this pipe system");

+end

+//part 2

+Uc=15;//actual gas velocity

+e=((Uc-U)^2/2/g/dt/100+1)^(1/-4.7);

+Usc=(Uc-U)*(1-e);//superficial speed of solid

+Cmax=Usc*rhos*pi*dt^2/4;

+disp(Cmax,"the maximum carrying capacity of polythene particles in (kg/s)");

diff --git a/1379/CH12/EX12.1.3/example12_3.sce b/1379/CH12/EX12.1.3/example12_3.sce
new file mode 100755
index 000000000..10aa51a2e
--- /dev/null
+++ b/1379/CH12/EX12.1.3/example12_3.sce
@@ -0,0 +1,40 @@
+

+

+//example 12.3 

+clc; funcprot(0);

+// Initialization of Variable

+rho=1.22;//density of air

+pi=3.1428;

+rhos=1400;//density of coal

+mu=1.73/10^5;

+g=9.81;

+U=25;

+Ut=2.80;

+l=50;

+ms=1.2;//mass flow rate

+mg=ms/10;//mass flow of gas

+//calculation

+Qs=ms/rhos;//flow of solid

+Qg=mg/rho;//flow of gas

+us=U-Ut;//actual linear velocity

+A=Qg/U;

+Us=Qs/A;//solid velocity

+e=(us-Us)/us;

+d=sqrt(4*A/pi);

+function [y ]= fround(x,n)

+// fround(x,n)

+// Round the floating point numbers x to n decimal places

+// x may be a vector or matrix// n is the integer number of places to round to

+y=round(x*10^n)/10^n;

+endfunction

+[d]=fround(d,4);

+Re=d*rho*U/mu;

+//using moody's chart

+phi=2.1/1000;//friction factor

+P1=2*phi*U^2*l*rho/d*2;

+f=0.05/us;

+P2=2*l*f*(0.0098)*rhos*us^2/d;

+P2=fround(P2/1000,1)*1000

+delP=rho*e*U^2+rhos*(0.0098)*us^2+P1+P2;

+//disp(delP,"the pressure difference in kN/m^2 ");

+printf('The Pressure value in (kN/m^2) is %.1f',delP/1000);

diff --git a/1379/CH12/EX12.1.4/example12_4.sce b/1379/CH12/EX12.1.4/example12_4.sce
new file mode 100755
index 000000000..78c4537f3
--- /dev/null
+++ b/1379/CH12/EX12.1.4/example12_4.sce
@@ -0,0 +1,39 @@
+

+

+//example 12.4 

+clc; funcprot(0);

+// Initialization of Variable

+rho=1.22;//density of air

+pi=3.1428;

+rhos=1090;//density of steel

+mu=1.73/10^5;

+g=9.81;

+d=14.5/100;

+Qg=0.4;

+Qs=5000/3600/1090;

+Ut=6.5;

+ar=0.046/1000;//absolute roughness

+l=18.5;//length

+//calculation

+function [y ]= fround(x,n)

+// fround(x,n)

+// Round the floating point numbers x to n decimal places

+// x may be a vector or matrix// n is the integer number of places to round to

+y=round(x*10^n)/10^n;

+endfunction

+Us=Qs/pi/d^2*4;//solid velocity

+U=Qg/pi/d^2*4;

+us=U-Ut;//actual linear velocity

+e=1-Us/us;

+e=fround(e,4);

+Re=rho*U*d/mu;

+rr=ar/d;//relative roughness

+//using moody's diagram

+phi=2.08/1000;

+P1=2*phi*U^2*l*rho/d*2;

+f=0.05/us;

+P2=2*l*f*(1-e)*rhos*us^2/d;

+P2=fround(P2/1000,2)*1000;

+delP=rhos*(1-e)*us^2+rhos*(1-e)*g*l+P1+P2;

+//disp(delP,"the pressure difference in kN/m^2 ");

+printf('The Pressure value in (kN/m^2) is %.2f',delP/1000)

diff --git a/1379/CH12/EX12.1.5/example12_5.sce b/1379/CH12/EX12.1.5/example12_5.sce
new file mode 100755
index 000000000..f99f24ae6
--- /dev/null
+++ b/1379/CH12/EX12.1.5/example12_5.sce
@@ -0,0 +1,15 @@
+

+

+//example 12.5 

+clc; funcprot(0);

+// Initialization of Variable

+l=25;

+pi=3.1428;

+rhos=2690;//density of ore

+emin=0.6;

+emax=0.8;

+//calculation

+Pmax=rhos*(1-emin)*g*l;

+disp(Pmax,"The maximum pressure drop in (N/m^2):");

+Pmin=rhos*(1-emax)*g*l;

+disp(Pmin,"The minimum pressure drop in (N/m^2):");