initial commit / add all books

1 files changed, 24 insertions, 0 deletions

diff --git a/1226/CH20/EX20.18/EX20_18.sce b/1226/CH20/EX20.18/EX20_18.sce
new file mode 100755
index 000000000..61a6fd81a
--- /dev/null
+++ b/1226/CH20/EX20.18/EX20_18.sce
@@ -0,0 +1,24 @@
+clc;funcprot(0);//EXAMPLE 20.18

+// Initialisation of Variables

+ns=2;.............//No of stages

+v1=0.2;.........//Intake volume in m^3/s

+p1=1;..........//Intake pressure in bar

+t1=289;...............//Intake temperature in K

+p3=7;.........//Final pressure in bar

+n=1.25;.......//Compression index

+N=600;........//Compressor rpm

+cp=1.005;.....//Specific heat at constant pressure in kJ/kgK

+R=0.287;......//Gas constant in kJ/kgK

+//Calculations

+p2=sqrt(p1*p3);......//Intermediate pressure in bar

+disp(p2/10,"Intermediate pressure in Mpa:")

+vslp=60*v1/N;..........//Volume of low pressure cylinder in m^3

+vshp=p1*vslp/p2;..........//Volume of high pressure cylinder in m^3

+disp(vslp,"Volume of low pressure cylinder in m^3:")

+disp(vshp,"Volume of high pressure cylinder in m^3:")

+W=(ns*(n/(n-1)))*p1*10^5*(v1/1000)*(((p3/p1)^((n-1)/(ns*n)))-1);...........//Power required to drive the compressor in kW

+disp(W,"Power required to drive the compressor in kW:")

+m=p1*10^5*v1/(R*t1*1000);.........//Mass of air handled in kg/s

+t2=t1*((p2/p1)^((n-1)/n));.........//Temperature at the end of first stage compression in K

+Qr=m*cp*(t2-t1);.......//Heat rejected in the intercooler in kW

+disp(Qr,"Heat rejected in the intercooler in kW:")