1 files changed, 36 insertions, 0 deletions
diff --git a/1226/CH17/EX17.32/EX17_32.sce b/1226/CH17/EX17.32/EX17_32.sce
new file mode 100755
index 000000000..21640613a
--- /dev/null
+++ b/1226/CH17/EX17.32/EX17_32.sce
@@ -0,0 +1,36 @@
+clc;funcprot(0);//EXAMPLE 17.32
+// Initialisation of Variables
+n=1;...................//No of cylinders
+d=0.032;................//Diameter of circular orifice in m
+Cd=0.62;.............//Co efficient of discharge
+hw=150;.................//Pressure across orfice in mm of water
+t=20+273;..............//Temperature of air in the room in Kelvin
+p=1.0132;.................//Ambient pressure in bar
+pd=0.00178;............//Piston displacement in m^3
+R=287;....................//Gas constant in J/kg
+r=6.5;..................//Compression ratio
+fc=0.135;................//Fuel consumption in kg/min
+C=43900;.................//Calorific value of fuel in kJ/kg
+BP=28;................//Brake power in kW
+N=2500;...................//Engine rpm
+k=0.5;....................//Four stroke engine
+g=9.81;.......................//Acceleration due to gravity in m/s^2
+rhow=1000;....................//Density of water in kg/m^3
+ga=1.4;........................//Degree of freedom
+//calculations
+mbyv=(p*10^5)/(R*t);
+pw=(hw/rhow)*rhow;....................//Pressure across orifice in kg/m^2
+H=pw/mbyv;........................//Head of air column causing the flow in m
+ma=Cd*(%pi/4)*d*d*sqrt(2*g*H);................//Air flow through orifice in m^3/s
+maps=(ma*60)/(N/2);........................//Air consumption per stroke
+etav=maps/pd;.................//Volumetric efficiency
+disp(etav*100,"Volumetric efficiency (in %):")
+ac=ma*60*mbyv;...............//Mass of air drawn into cylinder per min in kg
+afr=ac/fc;...................//Air fuel ratio
+disp(afr,"Air fuel ratio:")
+pmb=(6*BP)/(n*pd*N*k*10);...................//Mean effective pressure in bar
+disp(pmb,"Mean effective pressure (in bar):")
+etast=1-(1/(r^(ga-1)));...............//Air standard efficiency
+etabth=BP/((fc/60)*C);...............//Brake thermal efficiency
+etarel=etabth/etast;.................//Relative efficiency
+disp(etarel*100,"Relative efficiency (in %):")