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

---
 2762/CH3/EX3.1.1/3_1_1.sce | 16 ++++++++++++++++
 2762/CH3/EX3.1.2/3_1_2.sce | 15 +++++++++++++++
 2762/CH3/EX3.1.3/3_1_3.sce | 19 +++++++++++++++++++
 2762/CH3/EX3.1.4/3_1_4.sce | 27 +++++++++++++++++++++++++++
 2762/CH3/EX3.1.5/3_1_5.sce | 15 +++++++++++++++
 2762/CH3/EX3.2.1/3_2_1.sce | 26 ++++++++++++++++++++++++++
 2762/CH3/EX3.2.2/3_2_2.sce | 16 ++++++++++++++++
 2762/CH3/EX3.3.2/3_3_2.sce | 25 +++++++++++++++++++++++++
 2762/CH3/EX3.3.3/3_3_3.sce | 20 ++++++++++++++++++++
 2762/CH3/EX3.4.3/3_4_3.sce | 45 +++++++++++++++++++++++++++++++++++++++++++++
 2762/CH3/EX3.5.1/3_5_1.sce | 18 ++++++++++++++++++
 11 files changed, 242 insertions(+)
 create mode 100755 2762/CH3/EX3.1.1/3_1_1.sce
 create mode 100755 2762/CH3/EX3.1.2/3_1_2.sce
 create mode 100755 2762/CH3/EX3.1.3/3_1_3.sce
 create mode 100755 2762/CH3/EX3.1.4/3_1_4.sce
 create mode 100755 2762/CH3/EX3.1.5/3_1_5.sce
 create mode 100755 2762/CH3/EX3.2.1/3_2_1.sce
 create mode 100755 2762/CH3/EX3.2.2/3_2_2.sce
 create mode 100755 2762/CH3/EX3.3.2/3_3_2.sce
 create mode 100755 2762/CH3/EX3.3.3/3_3_3.sce
 create mode 100755 2762/CH3/EX3.4.3/3_4_3.sce
 create mode 100755 2762/CH3/EX3.5.1/3_5_1.sce

(limited to '2762/CH3')

diff --git a/2762/CH3/EX3.1.1/3_1_1.sce b/2762/CH3/EX3.1.1/3_1_1.sce
new file mode 100755
index 000000000..80d14caec
--- /dev/null
+++ b/2762/CH3/EX3.1.1/3_1_1.sce
@@ -0,0 +1,16 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3
+//Example 3.1-1
+//Principles of Momentum Transfer and Applications
+//given data
+rho=1.137;
+mu=1.9e-5;
+Dp=0.042;
+v0=23;
+Re=(Dp*v0*rho)/mu
+//from the mentioned graph,
+Cd=0.47;//drag coefficient as seen from the graph
+Ap=(%pi*Dp*Dp)/4;//surface area of sphere
+Fd=Cd*(v0*v0/2)*rho*Ap;//drag force
+mprintf("drag coefficient= %f",Cd)
+mprintf("drag force= %f N",Fd)
diff --git a/2762/CH3/EX3.1.2/3_1_2.sce b/2762/CH3/EX3.1.2/3_1_2.sce
new file mode 100755
index 000000000..4ef052c2b
--- /dev/null
+++ b/2762/CH3/EX3.1.2/3_1_2.sce
@@ -0,0 +1,15 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3
+//Example 3.1-2
+//Principles of Momentum Transfer and Applications
+//given data
+rho=997.2;//density of water at 300K in kg/m3
+v=1;//velocity of air in m/s
+Dp=0.09;//diameter of a cylinder
+mu=0.9142/1000;//viscosity of water Pa.s
+Re=Dp*v*rho/mu;//Reynolds Number
+Cd=1.4;//drag coefficient; found from a graph by the american chemical society
+L=1;//length of the tube in m
+Ap=L*Dp;
+Fd=Cd*(v*v/2)*rho*Ap;//drag force in newtons
+mprintf("the drag force is %f N",Fd)
diff --git a/2762/CH3/EX3.1.3/3_1_3.sce b/2762/CH3/EX3.1.3/3_1_3.sce
new file mode 100755
index 000000000..3829dd2c6
--- /dev/null
+++ b/2762/CH3/EX3.1.3/3_1_3.sce
@@ -0,0 +1,19 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3
+//Example 3.1-3
+//Principles of Momentum Transfer and Applications
+//given data
+basis=1;//taking basis as 1 m3 of packed bed
+rho=962;//bulk density of packed bed
+m=rho*basis;//total mass
+rho2=1600;//density of solid cylinders 
+V=m/rho2;//volume of the cylinder
+E=(basis-V)/(basis);//void fraction
+mprintf("void fraction = %f",E);
+D=0.02;//diameter of cylinder
+Av=6/D;// Av= Sp/Vp where Sp is the surface area of the particle and D is the diameter of the particle
+Dp=6/Av;//effectice diameter
+mprintf(" ii) effectice diameter = %f m",Dp);
+a=(6/Dp)*(1-E);//value of a
+mprintf(" iii) value of a= %f m-1",a)
+//end
diff --git a/2762/CH3/EX3.1.4/3_1_4.sce b/2762/CH3/EX3.1.4/3_1_4.sce
new file mode 100755
index 000000000..87325f335
--- /dev/null
+++ b/2762/CH3/EX3.1.4/3_1_4.sce
@@ -0,0 +1,27 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3
+//Example 3.1-4
+//Principles of Momentum Transfer and Applications
+//given data
+D=0.61;//diameter of bed in m
+h=2.44;//height of bed
+A=(3.14*D*D)/4;//cross section area in m2
+mdot=0.358;//mass flow rate of air in kg/s
+E=0.38//void fraction
+G=mdot/A;
+Dp=0.0127;//diameter of spheres in m
+mu=1.9e-5;//viscosity of air Pa.s
+delL=2.44;
+Re= (Dp*G/((1-E)*mu));//Reynolds Number
+delP=0.05e+5;//assumed pressure difference in pascal
+p1=1.115*101325;//air entering at this pressure in Pa
+p2=p1-delP;
+avgP=(p1+p2)/2;//average pressure
+M=28.97;//molecular weight of air in SI units 
+R=8314.34;//gas constant
+T=311;//temp of air in K
+avgrho=(M/(R*T))*avgP;//Avg density 
+//(delP*rho/G*G)*(Dp/delL)*(E^3/(1-E))= (150/Re)+1.75: erguns equation in dimensionless groups
+delPn=((150/Re)+1.75)*((G*G*delL*(1-E))/(avgrho*E*E*E*Dp));//calculated pressure drop
+mprintf("calculated pressure drop= %f Pa",delPn)
+//end
diff --git a/2762/CH3/EX3.1.5/3_1_5.sce b/2762/CH3/EX3.1.5/3_1_5.sce
new file mode 100755
index 000000000..eba51209a
--- /dev/null
+++ b/2762/CH3/EX3.1.5/3_1_5.sce
@@ -0,0 +1,15 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3
+//Example 3.1-5
+//Principles of Momentum Transfer and Applications
+//given data
+p=[0.25 0.40 0.35];//percentages of diff particle sizes
+s=[25 50 75];//sizes of the particles in mm
+phi=0.68;//sphericity
+sump=0
+for i=1:3
+    term=p(1,i)/(phi*s(1,i))
+    sump=sump+term
+end
+Dpm=1/sump;
+mprintf("mean diameter= %f mm",Dpm)
diff --git a/2762/CH3/EX3.2.1/3_2_1.sce b/2762/CH3/EX3.2.1/3_2_1.sce
new file mode 100755
index 000000000..89c3c9e1e
--- /dev/null
+++ b/2762/CH3/EX3.2.1/3_2_1.sce
@@ -0,0 +1,26 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3
+//Example 3.2-1
+//Principles of Momentum Transfer and Applications
+//given data
+rho=1.043;//density of air at 328.5K in kg/m3
+v=23;//velocity of air in m/s
+D=0.6;//diameter of a cylinder
+mu=2.03/100000;//viscosity of air Pa.s
+delh=0.205;// 0.205m of water pitot tube reading
+rhow=1000;//density of water 
+delP=delh*(rhow-rho)*9.80665;//pressure diff and g=9.80655 m/s2
+patm=101325;//atm pressure in pascals
+p1=patm+0.02008*100000;//absolute pressure+ pressure diff
+rhoc=(p1/patm)*1.043;//corrected air density
+delH=10.7/1000;//manometer reading, m of water
+Cp=0.98;
+delP=delH*(rhow-rhoc)*9.80655;//pressure diff in Pa
+v=Cp*((2*delP)/rhoc)^0.5;//max vel at center
+Re=D*v*rhoc/mu;//Reynolds Number
+vr=0.85;//from the given graph the ratio of avg vel/max vel is 0.85
+vavg=vr*v;//the average velcity in m/s
+mprintf(" average velcity = %f m/s",vavg)
+A=(3.14/4)*(D*D);//cross sec area in m2
+V=A*vavg;//volumetric flow rate in m3/s
+mprintf("volumetric flow rate = %f m3/s",V)
diff --git a/2762/CH3/EX3.2.2/3_2_2.sce b/2762/CH3/EX3.2.2/3_2_2.sce
new file mode 100755
index 000000000..8f18271b8
--- /dev/null
+++ b/2762/CH3/EX3.2.2/3_2_2.sce
@@ -0,0 +1,16 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3
+//Example 3.2-2
+//Principles of Momentum Transfer and Applications
+//given data
+delP=9.32e+4;//pressure diff in N/m2
+D1=0.1541;//external diameter in m
+D0=0.0566;//internal diameter in m
+Dr=D0/D1;
+Co=0.61;
+rho=878; //oil density in kg/m3
+v0=(Co/(sqrt(1-(Dr^4))))*sqrt((2*delP)/rho);//velocity calculation in m/s
+A=(%pi/4)*D0*D0;//cross section area
+V=A*v0;//volumetric flow rate
+mprintf("the volumetric flow rate is %f m3/s",V);
+//end
diff --git a/2762/CH3/EX3.3.2/3_3_2.sce b/2762/CH3/EX3.3.2/3_3_2.sce
new file mode 100755
index 000000000..723c71c89
--- /dev/null
+++ b/2762/CH3/EX3.3.2/3_3_2.sce
@@ -0,0 +1,25 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3
+//Example 3.3-2
+//Principles of Momentum Transfer and Applications
+//given data
+Ps=741.7;//suction pressure in mm hg
+Pd=769.6;//discharge pressure in mm hg
+Patm=760;//atmospheric pressure in mm Hg
+rho1=28.97*(1/22.414)*(273.2/366.3)*(Ps/Patm);//air density at suction: mol wt= 28.97 kg air/kg mol for 22.414 m3/kg mol at 101.3 kPa and 273.2 K
+rho2=rho1*(Pd/Ps);
+rhoavg=(rho1+rho2)/2
+V=28.32;//volumetric flow rate in m3/s
+Ts=294.1;//temp at suction
+mdot=V*(1/60)*(1/22.414)**(273.2/Ts)*28.97;//mass flow rate of gas
+Patm=760;//atm pressure in mm Hg
+Hp=((Pd-Ps)/Patm)*(101325/rhoavg);//pressure head in J/kg
+v1=0;//air is stationary
+v2=45.7;//discharge velocity in m/s
+vd=(((v2^2)-(v1^2))/2);//developed velocity
+z1=0;
+sumF=0;
+Ws=Hp+vd;//substituting and solving for Ws by mechanical energy balance equation for a closed system in J/kg
+n=60/100;//efficiency given is 60%
+bkW= (Ws*mdot)/(n*1000);// brake kW
+mprintf(" brake kW= %f hP",bkW)
diff --git a/2762/CH3/EX3.3.3/3_3_3.sce b/2762/CH3/EX3.3.3/3_3_3.sce
new file mode 100755
index 000000000..6106624c6
--- /dev/null
+++ b/2762/CH3/EX3.3.3/3_3_3.sce
@@ -0,0 +1,20 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3
+//Example 3.3-3
+//Principles of Momentum Transfer and Applications
+//given data
+p1=137.9*1000;
+p2=551.6*1000;
+T1=26.7+273.2;
+mmol=7.56/1000;
+M=16;
+mdot=mmol*M;
+gam=1.31;
+R=8314.3;
+nWs1=(gam/(gam-1))*(R*T1/M)*((p2/p1)^((gam-1)/gam)-1);
+n=80/100;
+bkW1=(nWs1*mdot)/(n*1000);
+mprintf("i) brake power= %f kW",bkW1)
+nWs2=(R*T1/M)*log(p2/p1);
+bkW2=(nWs2*mdot)/(n*1000);
+mprintf("ii) brake power= %f kW",bkW2)
diff --git a/2762/CH3/EX3.4.3/3_4_3.sce b/2762/CH3/EX3.4.3/3_4_3.sce
new file mode 100755
index 000000000..77f1c3057
--- /dev/null
+++ b/2762/CH3/EX3.4.3/3_4_3.sce
@@ -0,0 +1,45 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3
+//Example 3.4-3
+//Principles of Momentum Transfer and Applications
+//given data
+H1=1.83;
+DT1=H1;
+V1=(%pi*DT1*DT1*H1)/4;
+V2=3*V1;//given
+R=(V2/V1)^(1/3);
+DT2=R*DT1;
+Da1=0.61;
+Da2=R*Da1;
+W1=0.122;
+W2=R*W1;
+J1=0.15;
+J2=R*J1;
+N1=1.5;//no. of revs
+N2=N1*((1/R)^(2/3))
+rho=929;
+mu=0.01;
+Re=(Da2*Da2*N2*rho)/(mu)
+Np=5;
+P2=Np*rho*(N2^3)*(Da2^5);
+P1=Np*rho*(N1^3)*(Da1^5);
+//a)
+N2=N1*((1/R)^(2/3));
+sP1=P1/V1;
+sP2=P2/V2;
+
+mprintf("scaled up no. of revs %f rev/s",N2);
+mprintf("scaled up Power %f W",P2);
+mprintf(" power per unit volume= %f kW/m3",sP1/1000)
+if (sP1/1000)<0.8  then
+    disp(" Value of power is less than permissible condition(0.8 kW/m3 for mass transfer)")
+end
+mprintf(" scaled up Power %f m3",P2);
+mprintf(" power per unit volume %f W/m3",(P2/(V2*1000)));
+//b)
+N2b=N1*(1/R);
+mprintf(" scaled up revolutions %f rev/s",N2b);
+P2b=Np*rho*(N2b^3)*(Da2^5);
+mprintf(" scaled up Power %f kW",P2b);
+mprintf(" power per unit volume %f W/m3",(P2b/V2));
+
diff --git a/2762/CH3/EX3.5.1/3_5_1.sce b/2762/CH3/EX3.5.1/3_5_1.sce
new file mode 100755
index 000000000..fb7720c55
--- /dev/null
+++ b/2762/CH3/EX3.5.1/3_5_1.sce
@@ -0,0 +1,18 @@
+//Transport Processes and Seperation Process Principles
+//Chapter 3//Example 3.5-1
+//Principles of Momentum Transfer and Applications
+//given data
+Kd=15.23;
+nd=0.4;
+D=0.0524;
+V=0.0728;
+L=14.9;
+rho=1041;
+delP=(Kd*4*L/D)*((8*V/D)^nd);//pressure drop
+Ff=delP/rho;//friction loss
+nd=0.4;
+g=8;
+Re=((D^nd)*(V^(2-nd))*rho)/(Kd*(g^(nd-1)));
+f=16/Re;//friction factor
+delP=4*f*rho*(L/D)*(V*V/2);
+mprintf("pressure drop= %f kN/m2",delP/1000)
-- 
cgit