initial commit / add all books

17 files changed, 363 insertions, 0 deletions

diff --git a/2015/CH9/EX9.1/9_1.sce b/2015/CH9/EX9.1/9_1.sce
new file mode 100755
index 000000000..130be9bbb
--- /dev/null
+++ b/2015/CH9/EX9.1/9_1.sce
@@ -0,0 +1,25 @@
+clc

+//initialisation of variables

+t1=305 //temp in k

+r=0.287 //kj/kg

+p2=6 //pressure in bar

+p1=1.013 //pressure in bar

+g=1.4 //const value

+n=1.28

+v1=100 //volume

+//CALCULATIONS

+rp=(p2/p1)

+wiso=r*t1*log(p2/p1)

+wadia=(g/(g-1))*r*t1*0.6623

+wpoly=(n/(n-1))*r*t1*0.4756

+m=(p1*v1*100)/(r*t1)

+ipr=(wiso*m)/60

+apr=(wadia*m)/60

+//RESULTS

+printf('work for isthermal compression is %2fknm/kg',wiso)

+printf('\nwork for adiabatic compression is %2fknm/kg',wadia)

+printf('\nwork for polytropic compression is %2fknm/kg',wpoly)

+printf('\nmass of air compressed is %2fkg/min',m)

+printf('\nisothermal power required is %2fkw',ipr)

+printf('\nadiabatic power required is %2fkw',apr)

+    

diff --git a/2015/CH9/EX9.10/9_10.sce b/2015/CH9/EX9.10/9_10.sce
new file mode 100755
index 000000000..a31bb3c41
--- /dev/null
+++ b/2015/CH9/EX9.10/9_10.sce
@@ -0,0 +1,12 @@
+clc

+//initialisation of variables

+sp=200 //mean speed m/s

+//CALCULATIONS

+d=(21/(0.7773*1.18*200))^0.5

+l=1.5*d

+s=200/(3*d)

+//RESULTS

+disp('volumetric efficiency is 0.7773')

+printf('\ndiameter is %2fm',d)

+printf('\nstroke is %2fm',l)

+printf('\nspeed of compressor is %2frev/min',s)

diff --git a/2015/CH9/EX9.11/9_11.sce b/2015/CH9/EX9.11/9_11.sce
new file mode 100755
index 000000000..4be13c8d8
--- /dev/null
+++ b/2015/CH9/EX9.11/9_11.sce
@@ -0,0 +1,30 @@
+

+clc

+//initialisation of variables

+r=0.287

+p=1.01325 //pressure in bar

+v=5 //volume in m*m*m

+t=288 //temparature in k

+t1=303 //temparature in k

+t2=403 //temparature in k

+p2=4.08 //pressure in bar

+p1=0.98 //pressure in bar

+p3=17 //pressure in bar

+n=1.25

+c=0.06 //clearance volume by swept volume

+//CALCULATIONS

+m=(p*v)/(r*t)

+rp=p2/p1

+t2s=(t1*(p2/p1)^((n-1)/n))

+wr=(n/n-1)*r*(t2-t1)

+wc=2*wr

+veff=(1+c-c*(rp)^(1/n))

+x=(p*100*v*t1)/(p1*100*t) //x=(v1-v4)

+vs=x/veff

+vsc=vs/125

+d1=((4*vsc)/%pi)^(1/3)

+//RESULTS

+printf('volumetric efficiency is %2f',veff)

+printf('\nstoke volume is %2fm*m*m/min',vs)

+printf('\nstroke volume per cycle is %2fm*m*m',vsc)

+printf('\nstoke of piston is %2f',d1)

diff --git a/2015/CH9/EX9.12/9_12.sce b/2015/CH9/EX9.12/9_12.sce
new file mode 100755
index 000000000..62aed33a9
--- /dev/null
+++ b/2015/CH9/EX9.12/9_12.sce
@@ -0,0 +1,30 @@
+clc

+//initialisation of variables

+t1=303 //temparature in k

+p2=4.08 //pressure in bar

+p1=1 //pressure in bar

+t5=303 //temparature in k

+x=0.3247 //x=v2/v1 where the relation is v2=v1*(1/rp)^1/n

+y=0.0385 //y=v3/v1

+vo=0.2862 //vo=volume of air delivered/v1

+vf=0.8299 //vf=vome of free air /v1

+n=1.25

+p3=17.5 //pressure in bar

+r=0.287

+tatm=2911 //temp in k

+patm=1.02 //pressure in bar

+w=291 

+//CALCULATIONS

+t2=(t1*(p2/p1)^((n-1)/n))

+veff=vf/(1-y)

+a=(r*(t2-t1)*5)

+t3=(t1*(p3/p2)^((n-1)/n))

+hp=(5*r*(t3-t1))

+iso=(r*tatm*log(p3/patm))/10 //its ln

+ieff=iso/w

+//RESULTS

+printf('volumetric efficiency is %2f',veff)

+printf('\nwork required for lp cyclinder is %2f',a)

+printf('\nwork required for hp cyclinder is %2f',hp)

+printf('\nwork required for isothermal is %2f',iso)

+printf('\nisothermal efficiency is %2f',ieff)

diff --git a/2015/CH9/EX9.13/9_13.sce b/2015/CH9/EX9.13/9_13.sce
new file mode 100755
index 000000000..60e626f9a
--- /dev/null
+++ b/2015/CH9/EX9.13/9_13.sce
@@ -0,0 +1,18 @@
+clc

+//initialisation of variables

+p2=1.5 //pressure in bar

+p1=1 //pressure in bar

+v=0.05 //volume in m*m*m

+g=1.4

+r=1.4

+n=120 //number of cycles

+//CALCULATIONS

+wa=v*(p2-p1)*100

+wi=3.5*100*p1*v*(((p2/p1)^((r-1)/r))-1)

+reff=wi/wa

+vo=v/4

+pr=wa*n/60

+//RESULTS

+printf('roots efficiency is %2f',reff)

+printf('\nvolume of air is %2fm*m*m/cycle',vo)

+printf('\npower required is %2fkw',pr)

diff --git a/2015/CH9/EX9.14/9_14.sce b/2015/CH9/EX9.14/9_14.sce
new file mode 100755
index 000000000..2cdab6379
--- /dev/null
+++ b/2015/CH9/EX9.14/9_14.sce
@@ -0,0 +1,29 @@
+clc

+//initialisation of variables

+p2=1.5 //pressure in bar

+p1=1 //pressure in bar

+v=0.05 //volume in m*m*m

+x=0.35 //increse in pressure

+g=1.4

+r=1.4

+n=120 //number of cycles

+//CALCULATIONS

+wa=v*(p2-p1)*100

+wi1=3.5*100*p1*v*(((p2/p1)^((r-1)/r))-1)

+ceff=wi1/wa

+vo=v/4

+pr=wa*n/60

+prs=x*(p2-p1)

+p3=p1+prs

+wi2=3.5*100*p1*v*(((p3/p1)^((r-1)/r))-1)

+vi=v*(p1/p3)^(1/g)

+w2=vi*(p2-p3)*100

+tw=w2+wi2

+comeff=wi1/tw

+po=tw*2

+//RESULTS

+printf('compressor efficiency is %2f',ceff)

+printf('\nwork required for internal compression is %2fknm/rev',wi2)

+printf('\npower required is %2fkw',pr)

+printf('\ncompressor efficiency 2 is %2f',comeff)

+printf('\npower required 2 is %2fkw',po)

diff --git a/2015/CH9/EX9.15/9_15.sce b/2015/CH9/EX9.15/9_15.sce
new file mode 100755
index 000000000..668a5a42e
--- /dev/null
+++ b/2015/CH9/EX9.15/9_15.sce
@@ -0,0 +1,18 @@
+clc

+//initialisation of variables

+t1=295 //temp in k

+p1=1.02 //pressure in bar

+p2=7.14 //pressure in bar

+cp=1.005 //kj/kg

+g=1.4

+wr=250 //kj/kg

+//CALCULATIONS

+t2s=t1*(p2/p1)^((g-1)/g)

+wi=cp*(t2s-t1)

+ieff=wi/wr

+t2=(wr/cp)+t1

+//RESULTS

+printf('isentropic work is %2fkj/kg',wi)

+printf('\nisentropic efficiency is %2f',ieff)

+printf('\ntemparature 2 is %2fk',t2)

+disp('index of compression is 1.46')

diff --git a/2015/CH9/EX9.16/9_16.sce b/2015/CH9/EX9.16/9_16.sce
new file mode 100755
index 000000000..e9846f560
--- /dev/null
+++ b/2015/CH9/EX9.16/9_16.sce
@@ -0,0 +1,21 @@
+

+clc

+//initialisation of variables

+t1=310 //temp in k

+p1=1 //pressure in bar

+p2=4 //pressure in bar

+cp=1.005 //kj/kg

+v1=28 //m*m*m volume

+r=0.287

+ce=0.7 //copression efficiency

+g= 32.2 //ft/sec^2

+//CALCULATIONS

+t2s=t1*(p2/p1)^((g-1)/g)

+wi=cp*(t2s-t1)

+m=(p1*v1*100)/(r*t1)

+apr=(m*wi)/60

+iei=wi/ce

+//RESULTS

+printf('isentropic work is %2f',apr)

+printf('\nadiabatic power required is %2f',m)

+printf('\nindicated enthalpy increase is %2f',iei)

diff --git a/2015/CH9/EX9.17/9_17.sce b/2015/CH9/EX9.17/9_17.sce
new file mode 100755
index 000000000..097441e0e
--- /dev/null
+++ b/2015/CH9/EX9.17/9_17.sce
@@ -0,0 +1,24 @@
+

+clc

+//initialisation of variables

+p2=6 //prressure in bar

+p1=1 //pressure in bar

+t1=313 //temp in k

+a1=45 //angle in degrees

+a2=10 //angle in degrees

+a3=55 //angle in degrees

+r=1.4

+cp=1.005 //kj/kg

+ieff=0.85 //isentropic efficiency

+c=200 //m/s

+//CALCULATIONS

+t2s=(t1*(p2/p1)^((r-1)/r))

+t2=(((t2s-t1)/ieff)+t1)

+w=cp*(t2-t1)

+cro=(c*(sin(45*(%pi/180))/sin(55*(%pi/180))))

+cv=c-cro

+n=w/cv

+//RESULTS

+printf('actual work is %2fkj/kg',w)

+printf('\nchange in whirl velocities is %2fkj/kg/stage',cv)

+printf('\nnumber of stages is %2fstages',n)

diff --git a/2015/CH9/EX9.2/9_2.sce b/2015/CH9/EX9.2/9_2.sce
new file mode 100755
index 000000000..807023bf5
--- /dev/null
+++ b/2015/CH9/EX9.2/9_2.sce
@@ -0,0 +1,15 @@
+clc

+//initialisation of variables

+p2=135 //bar pressure

+p1=1 //bar pressure

+x=5 //x=p2/p1

+//CALCULATIONS

+s=log(p2)/log(x)

+rp=(p2/p1)^0.25

+//RESULTS

+printf('s is %2f',s)

+printf('\nrp is %2f',rp)

+disp('number of stages are 4')

+disp('1st intermediate pressure is 3.4087 bar')

+disp('2nd intermediate pressure is 11.619 bar')

+disp('3rd intermediate pressure is 39.605 bar')

diff --git a/2015/CH9/EX9.3/9_3.sce b/2015/CH9/EX9.3/9_3.sce
new file mode 100755
index 000000000..9cf6e6116
--- /dev/null
+++ b/2015/CH9/EX9.3/9_3.sce
@@ -0,0 +1,28 @@
+clc

+//initialisation of variables

+p2=3.24 //pressure in bar

+p1=1 //pressure in bar

+v1=16 //volume in m*m*m

+n=1.35

+rp=3.24 //pressure

+r=10.5

+t1=294 //temparature in k

+t2=294 //temparature in k

+cp=1.005 //kj/kg 

+rx=0.287

+//CALCULATIONS

+w1=(2*n/(n-1))*p1*v1*100*0.35630 //(3.24)^0.2592-1

+w2=(n/(n-1))*p1*v1*100*0.8396 //(10.5)^0.2592-1

+pr1=w1/60

+pr2=w2/60

+tb=t1*(r)^(n-1/n)

+t3=t2*(rp)^((n-1)/n)

+m=(p1*v1*100)/(rx*t1)

+hr=m*cp*(t3-t2)

+ma=hr/(4.18*25)

+//RESULTS

+printf('minimum power required are %2fkw and %2fkw',pr1,pr2)

+printf('\nmass of air compressed is %2fkg/min',m)

+printf('\nheat rejected by air compressor is %2fkj/min',hr)

+printf('\nmass of water is %2fkg/min',ma)

+

diff --git a/2015/CH9/EX9.4/9_4.sce b/2015/CH9/EX9.4/9_4.sce
new file mode 100755
index 000000000..2e073139c
--- /dev/null
+++ b/2015/CH9/EX9.4/9_4.sce
@@ -0,0 +1,25 @@
+clc

+//initialisation of variables

+p2=4.08 //pressure in bar

+p1=1 //pressure in bar

+n=1.22

+r=0.287

+p=1.01325 //pressure in bar

+v=145 //volume

+t=288 //temparature in k

+p3=17.5 //pressure in bar

+t1=307 //temp in k

+t2=313 //temp in k

+//CALCULATIONS

+wlp=5.54*r*t1*(((p2/p1)^((n-1)/n))-1)

+whp=5.54*r*t2*(((p2/p1)^((n-1)/n))-1)

+w=wlp+whp

+m=(p*v)/(r*t)

+pr=(w*m)/60

+p2=(p1*p3)^0.5

+x=(p2)^0.5 //x=d1/d2

+//RESULTS

+printf('total work required is %2fknm/kg',w)

+printf('\nmass of free air is %2fkg/min',m)

+printf('\npower required to drive the compressor is %2fkw',pr)

+printf('\nratio of cylinder diameters is %2f',x)

diff --git a/2015/CH9/EX9.5/9_5.sce b/2015/CH9/EX9.5/9_5.sce
new file mode 100755
index 000000000..64d32e976
--- /dev/null
+++ b/2015/CH9/EX9.5/9_5.sce
@@ -0,0 +1,14 @@
+clc

+//initialisation of variables

+c1=0.05 //percentage

+c2=0.10 //percentage

+c3=0.20 //percentage

+rp=10

+//CALCULATIONS

+eff1=(1+c1-c1*(rp)^(0.78125))

+eff2=(1+c2-c2*(rp)^(0.78125))

+eff3=(1+c3-c3*(rp)^(0.78125))

+//RESULTS

+printf('volumetric effiency 1 is %2f',eff1)

+printf('\nvolumetric effiency 2 is %2f',eff2)

+printf('\nvolumetric effiency 3 is %2f',eff3)

diff --git a/2015/CH9/EX9.6/9_6.sce b/2015/CH9/EX9.6/9_6.sce
new file mode 100755
index 000000000..040b582be
--- /dev/null
+++ b/2015/CH9/EX9.6/9_6.sce
@@ -0,0 +1,18 @@
+clc

+//initialisation of variables

+d=0.2 //diameter in m

+lc=0.01 //linear clearance

+l=0.3 //lenght

+rp=7

+n=1.25

+pi=(22/7)

+//CALCULATIONS

+cv=((pi/4)*((d)^2)*lc)

+sv=((pi/4)*(d)^2*l)

+cr=cv/sv

+veff=(1+cr-cr*(rp)^(1/n))

+x=veff*sv

+//RESULTS

+printf('clearance ratio is %2f',cr)

+printf('\nvolumetric efficiency is %2f',veff)

+printf('\nvolume of air taken in is %2fm*m*/stroke',x)

diff --git a/2015/CH9/EX9.7/9_7.sce b/2015/CH9/EX9.7/9_7.sce
new file mode 100755
index 000000000..42862ed22
--- /dev/null
+++ b/2015/CH9/EX9.7/9_7.sce
@@ -0,0 +1,14 @@
+clc

+//initialisation of variables

+n=1.2

+r=0.287

+t1=310 //temparature in degrees

+p2=7 //pressure in bar

+p1=1 //pressure in bar

+//CALCULATIONS

+rp=(p2/p1)

+wr=((n/(n-1))*r*t1*((rp)^((n-1)/n)-1))

+//RESULTS

+disp('volumetric efficiency is 0.797')

+disp('volumetric efficiency referred to atmospheric conditions is 0.73')

+printf('work required is %2fknm/kg',wr)

diff --git a/2015/CH9/EX9.8/9_8.sce b/2015/CH9/EX9.8/9_8.sce
new file mode 100755
index 000000000..b849a9c99
--- /dev/null
+++ b/2015/CH9/EX9.8/9_8.sce
@@ -0,0 +1,13 @@
+clc

+//initialisation of variables

+veff=0.8 //efficiency

+rp=7 

+n=1.2 //constant value

+pi=(22/7)

+//CALCULATIONS

+c=(veff-1)/(1-(rp)^(1/n))

+vs=2/c

+d=((4*vs)/pi)^(1/3)

+//RESULTS

+printf('stroke volume is %2fm*m*m',vs)

+printf('\nlenght of stroke is %2fm',d)

diff --git a/2015/CH9/EX9.9/9_9.sce b/2015/CH9/EX9.9/9_9.sce
new file mode 100755
index 000000000..e01e0fdbd
--- /dev/null
+++ b/2015/CH9/EX9.9/9_9.sce
@@ -0,0 +1,29 @@
+clc

+//initialisation of variables

+sp=1400 //speed in revolutions per min

+ma=15 //mass in kgs

+r=0.287

+p1=1 //pressure in bar

+t1=303 //temparature in k

+p2=7 //pressure in bar

+c=0.05 //clearance volume/stoke volume

+pi=(22/7)

+n=1.2

+m1=15

+meff=0.85 //mechanical efficinecy

+//CALCULATIONS

+rp=(p2/p1)

+m=ma/sp

+va=(m1*r*t1)/(p1*100)

+eff1=(1+c-c*(rp)^(1/n))

+vs=va/eff1

+d1=((4*vs)/pi)^(1/3)

+pr=((n/(n-1))*m1*r*t1*((rp)^((n-1)/n)-1))/60

+prs=pr/meff

+d2=((prs*4)/(7*100*pi*700))^0.333

+//RESULTS

+printf('volumetric efficiency is %2f',eff1)

+printf('\nlengh of the stroke is %2fm',d1)

+printf('\nindicated power is %2fkw',pr)

+printf('\npower required at the shaft of the compressor is %2fkw',prs)

+printf('\ndiameter of the piston is %2fm',d2)