4 files changed, 102 insertions, 0 deletions
diff --git a/443/CH20/EX20.2/20_2.sce b/443/CH20/EX20.2/20_2.sce
new file mode 100755
index 000000000..64758211d
--- /dev/null
+++ b/443/CH20/EX20.2/20_2.sce
@@ -0,0 +1,21 @@
+pathname=get_absolute_file_path('20_2.sce')
+filename=pathname+filesep()+'20_2_data.sci'
+exec(filename)
+//Problem no 20.2
+//Mass of air, Ma
+ma=ffr/far
+//Swept Volume Vs
+Vs=(%pi/4)*(d*d)*(120/1000)
+//Volume total Vtot
+Vtot=Vs*(cr/15)
+//Mass of air supplied
+mas=(ma/(60*s))
+//M(ref)
+mref=(Vtot*(ep/(r*t)))
+//Scavanging efficiency nsc
+nsc=mas/mref
+//Scavanging ratio
+scr=(a/(60*s*mref))
+printf("\n\nRESULTS\n\n")
+printf("\nScavanging efficiency:%f\n",nsc)
+printf("\nScavanging ratio:%f\n",scr)
+\ No newline at end of file
diff --git a/443/CH20/EX20.3/20_3.sce b/443/CH20/EX20.3/20_3.sce
new file mode 100755
index 000000000..3f3bca9c8
--- /dev/null
+++ b/443/CH20/EX20.3/20_3.sce
@@ -0,0 +1,23 @@
+pathname=get_absolute_file_path('20_3.sce')
+filename=pathname+filesep()+'20_3_data.sci'
+exec(filename)
+//Question 20.3 
+//hydrogen ratio
+h=(26*1)/(12*12)
+//Fuel to air ratio 
+far=0.33*(1+h)*((CO2+CO)/N) 
+//Mass flow of air
+ma=ff/(far*60)
+//Volume total (in m3/cycle)
+Vtot=(%pi/4*(d*d)*(2*l)*(cr/(cr-1)))
+//Scavanger density
+psc=(ep/(r*t))
+//Theoretical mass flow
+tmf=psc*Vtot*s
+//Scavanger efficiency
+nsc=(ma/tmf)
+//Indicated mean effective pressure (in kN/m2)
+imep = (nsc*psc)*(cr/(cr-1))*(nith*far*CV)
+printf("\n\nRESULTS\n\n")
+printf("\nScavanger efficiency: %f\n",nsc)
+printf("\nIndicated mean effective pressure %f\n",imep)
+\ No newline at end of file
diff --git a/443/CH20/EX20.4/20_4.sce b/443/CH20/EX20.4/20_4.sce
new file mode 100755
index 000000000..68d9ccf78
--- /dev/null
+++ b/443/CH20/EX20.4/20_4.sce
@@ -0,0 +1,22 @@
+pathname=get_absolute_file_path('20_4.sce')
+filename=pathname+filesep()+'20_4_data.sci'
+exec(filename)
+//Question 20.4 
+//psc (in kg/m3)
+psc=ep/(r*t)
+//Mass of air (in kg/cycle)
+ma=nsc*V*psc
+//Scavanging ratio Rsc
+Rsc=(cf/(s*V*psc))
+//Trapping efficiency ntr
+ntr=nsc/Rsc
+//Brake power (in kW)
+bp=((ma*CV*nbth*s*far)/60)
+//the bsfc (in kg/kWh)
+bsfc=(cf*60*far/(bp))
+//Short circuting (in kg/h)
+sc=(cf*(1-ntr)*(60*far))
+printf("\n\nRESULTS\n\n")
+printf("\nbrake power: %f\n",bp)
+printf("\nthe bsfc %f\n",bsfc)
+printf("\nShort Circuting %f\n",sc)
+\ No newline at end of file
diff --git a/443/CH20/EX20.5/20_5.sce b/443/CH20/EX20.5/20_5.sce
new file mode 100755
index 000000000..c197293e3
--- /dev/null
+++ b/443/CH20/EX20.5/20_5.sce
@@ -0,0 +1,36 @@
+pathname=get_absolute_file_path('20_5.sce')
+filename=pathname+filesep()+'20_5_data.sci'
+exec(filename)
+//Question 20.5 page 715
+//Volume total (in m3/cycle)
+Vtot=k*(%pi/4)*(d*d)*(2*l)*(cr/(cr-1))
+//Scavanger Density (psc)
+psc=ep/(r*t)
+//Theortical mass flow rate (in kg/cycle)
+tmf=Vtot*psc
+//Actual mass flow rate (amr) (in kg/cycle)
+amr=nsc*tmf
+//Mass flow rate of air supplied by the blower (mfrb) (in kg/cycle)
+mfrb=Rsc*tmf
+//Actual temperature (dta) (in kelvin)
+dta=((300*((1.15/1)^0.286))-t)/0.75
+//New temperature (T2 in kelvin)
+T2=t+dta
+//Compressor work (in kW)
+Wc=((amr*dta*cp*s)/60)
+//V displacement (in m3/cycle)
+Vdisp=(k*(%pi/4)*(d*d*2*l))
+//Brake mean effective pressure (bmep in bar)
+bmep=bp*1000/(Vdisp*(s/60)*10^(5))
+//Total mechanical loss
+fp = 1.5*20
+ip =bp+fp
+//Fuel consumed per hour (mf)
+mf=bsfc*bp
+//Idicated thermal efficiency (nith)
+nith=(ip/(mf*CV))*100
+//Fuel-air ratio
+far=(mf/(60*720*Vtot))
+printf("\n\nRESULTS\n\n")
+printf("\nIdicated thermal efficiency %f\n",nith)
+printf("\nFuel-air Ratio %f\n",far)
+\ No newline at end of file