From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 1379/CH12/EX12.1.1/example12_1.sce | 26 +++++++++++++++++++++++++
 1379/CH12/EX12.1.2/example12_2.sce | 30 ++++++++++++++++++++++++++++
 1379/CH12/EX12.1.3/example12_3.sce | 40 ++++++++++++++++++++++++++++++++++++++
 1379/CH12/EX12.1.4/example12_4.sce | 39 +++++++++++++++++++++++++++++++++++++
 1379/CH12/EX12.1.5/example12_5.sce | 15 ++++++++++++++
 5 files changed, 150 insertions(+)
 create mode 100755 1379/CH12/EX12.1.1/example12_1.sce
 create mode 100755 1379/CH12/EX12.1.2/example12_2.sce
 create mode 100755 1379/CH12/EX12.1.3/example12_3.sce
 create mode 100755 1379/CH12/EX12.1.4/example12_4.sce
 create mode 100755 1379/CH12/EX12.1.5/example12_5.sce

(limited to '1379/CH12')

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):");
-- 
cgit