initial commit / add all books

11 files changed, 242 insertions, 0 deletions

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)