diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /1226/CH20/EX20.4/EX20_4.sce | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '1226/CH20/EX20.4/EX20_4.sce')
| -rwxr-xr-x | 1226/CH20/EX20.4/EX20_4.sce | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/1226/CH20/EX20.4/EX20_4.sce b/1226/CH20/EX20.4/EX20_4.sce new file mode 100755 index 000000000..3e1b8e0ca --- /dev/null +++ b/1226/CH20/EX20.4/EX20_4.sce @@ -0,0 +1,28 @@ +clc;funcprot(0);//EXAMPLE 20.4
+// Initialisation of Variables
+p1=1;........//Suction pressure in bar
+t1=293;.......//Suction temperature in K
+p2=6;..........//Discharge pressure in bar
+t2=453;.......//Discharge temperature in K
+N=1200;.........//Compressor rpm
+Ps=6.25;........//Shaft power in kW
+ma=1.7;........//Mass of air delivered in kg/min
+D=0.14;......//Engine bore in m
+L=0.10;.......//Engine stroke in m
+R=287;..........//Gas constant in kJ/kgK
+//Calculations
+Vd=(%pi/4)*D*D*L*N;.........//Displlacement volume in m^3/min
+FAD=ma*R*t1/(p1*10^5);........//Free air delivered
+etav=FAD/Vd;.....//Volumetric efficiency
+disp(etav*100,"Volumetric Efficiency in %:")
+n=1/(1-((log(t2/t1))/(log(p2/p1))));........//Index of compression
+IP=(n/(n-1))*(ma/60)*(R/1000)*t1*(((p2/p1)^((n-1)/n))-1);..........//Indicated power in kW
+disp(IP,"Indicated power in kW:")
+Piso=((ma/60)*(R/1000)*t1*(log(p2/p1)));..........//Isothermal power
+etaiso=Piso/IP;..............//Isothermal efficiency
+disp(etaiso*100,"Isothermal efficiency in %:")
+etamech=IP/Ps;...........//Mechanical efficiency
+disp(etamech*100,"Mechanical efficiency in %:")
+etao=Piso/Ps;........//Overall efficiency
+disp(etao*100,"Overall efficiency in %:")
+
|