initial commit / add all books

6 files changed, 250 insertions, 0 deletions

diff --git a/1379/CH6/EX6.1.1/Example6_1.sce b/1379/CH6/EX6.1.1/Example6_1.sce
new file mode 100755
index 000000000..883a9dc99
--- /dev/null
+++ b/1379/CH6/EX6.1.1/Example6_1.sce
@@ -0,0 +1,58 @@
+

+//example 6.1

+clc; funcprot(0);

+//exapple 6.1 

+// Initialization of Variable

+atp=100.2*1000;

+g=9.81;

+rho_w=996;

+rho_toluene=867;

+vap_pre_toluene=4.535*1000;

+viscosity_toluene=5.26/10000;

+//calculation

+m=(atp-vap_pre_toluene)/rho_toluene/g;

+disp(m,"Max. height of toluene supported by atm. pressure (in m):");

+//part(1)

+hopw=0.650;//head of pump in terms of water

+hopt=hopw*rho_w/rho_toluene;//head of pump in terms of toluene

+Q=1.8*10^-3;//flow in m^3/s

+d=2.3*10^-2;//diameter of pipe

+pi=3.14127;

+//u=4*Q/pi/d^2

+//substituting this for reynolds no.

+Re=4*Q*rho_toluene/pi/d/viscosity_toluene;//reynolds no.

+disp(Re ,"reynolds no :");

+phi=0.0396*Re^-0.25;

+//since both LHS and RHS are function of x(max. ht. ab. toluene) 

+//we define a new variable to solve the eqn

+//y=(atp/rho_toluene/g)-(vap_pre_toluene/rho_toluene/g)-(4*phi*16*Q^2*x/pi^2/d^5/g)-hopt;

+//y=x 

+//these are two equations

+b=[0;((atp/rho_toluene/g)-(vap_pre_toluene/rho_toluene/g)-hopt)];

+A=[1 -1;1 4*phi*16*Q^2/pi^2/d^5/g];

+x=A\b;

+disp(x(2,1), "the maximum height above toulene in the tank the pump can be located without risk while flow rate is 1.80dm^3/s (in m):");

+//solution of part(2)

+l=9//length 

+u=sqrt(((atp/rho_toluene/g)-(vap_pre_toluene/rho_toluene/g)-hopt-l)*d*g/4/phi/l);//fluid velocity in pipes

+Q=pi*d^2*u/4;

+disp(Q,"Maximum delivery rate if pump is located 9m above toluene tank(in m^3/s)")

+//solution of part(3)

+//clubing d together we get

+Q=1.8/1000;

+a=(atp/rho_toluene/g)-(vap_pre_toluene/rho_toluene/g)-hopt-l;

+b=a*pi^2*g/4/9/16/Q^2/0.0396/(4*Q*rho_toluene/pi/viscosity_toluene)^-0.25;

+d=(1/b)^(1/4.75);

+disp(d , "minimum smooth diameter of suction pipe which will have flow rate as (1.8 dm^3/s) for pump kept at 9 m high (in m):");

+

+

+

+

+

+

+

+

+

+

+

+

diff --git a/1379/CH6/EX6.1.2/example6_2.sce b/1379/CH6/EX6.1.2/example6_2.sce
new file mode 100755
index 000000000..11f5d991e
--- /dev/null
+++ b/1379/CH6/EX6.1.2/example6_2.sce
@@ -0,0 +1,18 @@
+

+//example 6.2

+clc; funcprot(0);

+//exapple 6.2 

+// Initialization of Variable

+Q1=24.8/1000;//flow in pump 1

+d1=11.8/100;//diameter of impeller 1

+H1=14.7//head of pump 1

+N1=1450//frequency of motor 1

+Q2=48/1000//flow in pump 2

+//calculation

+H2=1.15*H1;//head of pump 2

+specific_speed=N1*Q1^0.5/H1^0.75;

+N2=specific_speed*H2^0.75/Q2^0.5;//frequency of motor 2

+disp(N2 ,"frequency of motor 2 in rpm");

+d2=sqrt(N2^2*H1/H2/N1^2/d1^2);

+disp(1/d2 , "diametr of impeller 2 (in m)");

+

diff --git a/1379/CH6/EX6.1.3/example6_3.sce b/1379/CH6/EX6.1.3/example6_3.sce
new file mode 100755
index 000000000..aa5501287
--- /dev/null
+++ b/1379/CH6/EX6.1.3/example6_3.sce
@@ -0,0 +1,56 @@
+

+//example 6.3

+clc; funcprot(0);

+clf()

+//exapple 6.3 

+// Initialization of Variable

+Q=[0 0.01 0.02 0.03 0.04 0.05];//discharge

+effi_hyd=[65.4 71 71.9 67.7 57.5 39.2];

+effi_over=[0 36.1 56.0 61.0 54.1 37.0];

+H_sys=[0 0 0 0 0 0]

+d=0.114;//diameter of pipe

+d_o=0.096;//diameter of impeller

+h=8.75;//elevation

+g=9.81;//acc. of gravity

+rho=999;//denisity of water

+l=60;//length of pipe

+theta=0.611;//angle in radians

+B=0.0125;//width of blades

+pi=3.1412

+mu=1.109/1000;//viscosity of water

+omega=2*pi*1750/60;

+// calculation

+ for i=1:6

+     if i==1 then

+     H_sys(i)=h;

+ else

+     

+    H_sys(i)=h+8*Q(i)^2/pi^2/d^4/g*(1+8*l*0.0396/d*(4*rho*Q(i)/pi/d/mu)^-0.25);

+end,

+end;

+H_theor=omega^2*d_o^2/g-omega*Q/2/pi/g/B/tan(theta);

+//disp(H_sys"head of system (in m)");

+//disp(H_theor);

+for i=1:6

+    H_eff(i)=effi_hyd(i)*H_theor(i)/100;

+end

+//disp(H_eff);

+plot(Q,effi_hyd, 'r--d');

+plot(Q,effi_over, 'g');

+plot(Q,H_eff,'k');

+plot(Q,H_theor);

+plot(Q,H_sys ,'c-');

+title('system characteritics');

+ylabel('Head(m)or Efficiency(%)');

+xlabel('volumetric flow rate(m^3/s)');

+//calculation of power

+//at intersecting point using datatrip b/w H_sys &H_eff

+Q=0.0336

+effi_over=59.9

+H_eff=13.10

+P=H_eff*rho*g*Q/effi_over/10;

+disp(P ,"Power required to pump fluid at this rate(in KW):")

+

+

+

+

diff --git a/1379/CH6/EX6.1.4/example6_4.sce b/1379/CH6/EX6.1.4/example6_4.sce
new file mode 100755
index 000000000..197cdf003
--- /dev/null
+++ b/1379/CH6/EX6.1.4/example6_4.sce
@@ -0,0 +1,35 @@
+
+
+clc; funcprot(0);
+clf()
+//exapple 6.4 
+// Initialization of Variable
+//each is increased by five units to make each compatible for graph plotting
+Q=[0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1];//flow rate
+HeffA=[20.63 19.99 17.80 14.46 10.33 5.71 0 0 0 0 0 ];//Heff of pump A
+HeffB=[18 17 14.95 11.90 8.10 3.90 0 0 0 0 0];//Heff of pump B
+alpha=1;
+h=10.4;
+d=0.14;
+l=98;
+pi=3.1412;
+g=9.81;
+rho=999;
+for i=1:11
+     if i==1 then
+     H_sys(i)=h;
+ else
+     
+    H_sys(i)=h+8*Q(i)^2/pi^2/d^4/g*(1+8*l*0.0396/d*(4*rho*Q(i)/pi/d/mu)^-0.25);
+end,
+end;
+//H_sys is head of the system
+disp(H_sys, "the head of system in terms of height of water :");
+plot(Q,H_sys,'r--d');
+plot(Q,HeffA ,'-c');
+plot(Q,HeffB);
+//at intersecting point using datatrip b/w H_sys &H_effA
+disp(0.03339,"the flow rate at which H_sys takes over HeffA"); 
+
+
+
diff --git a/1379/CH6/EX6.1.5/example6_5.sce b/1379/CH6/EX6.1.5/example6_5.sce
new file mode 100755
index 000000000..efd9ca2f3
--- /dev/null
+++ b/1379/CH6/EX6.1.5/example6_5.sce
@@ -0,0 +1,57 @@
+

+//example 6.5

+clc; funcprot(0);

+//exapple 6.5 

+// Initialization of Variable

+rho=1000;

+dc=.15;

+l=7.8;

+g=9.81;

+pi=3.1428;

+atp=105.4*1000;

+vap_pre=10.85*1000;

+sl=.22;

+dp=0.045;

+h=4.6;

+//("x(t)=sl/2*cos(2*pi*N*t)"  "the function of displcement");

+//"since we have to maximize the acceleration double derivate the terms");

+//since double derivation have the term cos(kt) 

+//finding it maxima

+t=linspace(0,5,100);

+k=1;

+function[m,v]= maximacheckerforcosine()

+h=0.00001;

+a=0.00;

+for i=1:400 

+    if (cos(a+h)-cos(a-h))/2*h==0 & cos(i-1)>0 then

+break;

+else 

+    a=0.01+a;

+end 

+break;   

+end

+m=i-1;

+v=cos(i-1);

+endfunction;

+[a, b]= maximacheckerforcosine();

+disp(a,"time t when the acceleration will be maximum(s)");

+//double derivative will result in a square of value of N

+//lets consider its coefficient all will be devoid of N^2 

+k=sl/2*(2*pi)^2//accn max of piston

+kp=k*1/4*pi*dc^2/1*4/pi/dp^2;//accn coeff. ofsuction pipe

+f=1/4*pi*dp^2*l*rho*kp;//force exerted by piston

+p=f/1*4/pi/dp^2;//pressure exerted by piston

+//calculation

+o=atp-h*rho*g-vap_pre;

+//constant term of quadratic eqn

+y=poly([o 0 -p],'N', 'coeff')

+a=roots(y);

+disp(abs(a(1,1)),"Maximum frequency of oscillation if cavitation o be avoided(in Hz)");

+

+

+

+

+

+

+

+

diff --git a/1379/CH6/EX6.1.6/example6_6.sce b/1379/CH6/EX6.1.6/example6_6.sce
new file mode 100755
index 000000000..154994bb6
--- /dev/null
+++ b/1379/CH6/EX6.1.6/example6_6.sce
@@ -0,0 +1,26 @@
+

+//example 6.6

+clc; funcprot(0);

+//exapple 6.6 

+// Initialization of Variable

+rhos=1830;//density of acid

+atp=104.2*1000;//atmospheric pressure

+temp=11+273;//temp in kelvin

+M=28.8/1000;//molar mass of air

+R=8.314;//universal gas constant

+g=9.81;//acceleration of gravity

+pi=3.14;

+d=2.45;//diameter of tank

+l=10.5;//length of tank

+h_s=1.65;//height of surface of acid from below

+effi=0.93//efficiency

+//calculation

+mliq=pi*d^2*l*rhos/4;

+h_atm=atp/rhos/g;//height conversion of atp

+h_r=4.3-1.65;//height difference

+mair=g*h_r*mliq*M/(effi*R*temp*log(h_atm/(h_atm+h_s)));//mass of air

+disp(mair,"mass of air required to lift the sulphuric acid tank");

+disp("The negative sign indicates air is expanding & work done is magnitude of value in kg:");

+m=abs(mair/mliq);

+disp(m, "The mass of air required for per kilo of acid transferred:");

+