1 files changed, 46 insertions, 0 deletions
diff --git a/32/CH12/EX12.13/12_13.sce b/32/CH12/EX12.13/12_13.sce
new file mode 100755
index 000000000..0704bda41
--- /dev/null
+++ b/32/CH12/EX12.13/12_13.sce
@@ -0,0 +1,46 @@
+//pathname=get_absolute_file_path('12.13.sce')
+//filename=pathname+filesep()+'12.13-data.sci'
+//exec(filename)
+//Diameter of HP, LP and IP cylinder(in m):
+dhp=0.25
+dip=0.40
+dlp=0.85
+//MEPs of the cylinders(in kPa):
+mephp=0.5*10^3
+mepip=0.3*10^3
+meplp=0.1*10^3
+//Pressure at which steam is supplied(in kPa):
+p1=1.5*10^3
+//Pressure at exhaust(in kPa):
+p4=25
+//Cut-off occurs at:
+r1=0.60
+//Area of HP cylinder(in m^2):
+AHP=%pi*dhp^2/4
+//Area of IP cylinder(in m^2):
+AIP=%pi*dip^2/4
+//Area of LP cylinder(in m^2):
+ALP=%pi*dlp^2/4
+//Mep of HP referred to LP cylinder(in kPa):
+mep1=mephp*AHP/ALP
+//Mep of IP referred to LP cylinder(in kPa):
+mep2=mepip*AIP/ALP
+//Overall mep referred to LP cylinder(in kPa):
+mept=mep1+mep2+meplp
+//Overall expansion ratio:
+r=ALP/(r1*AHP)
+//Hypothetical mep(in kPa):
+mep=p1/r*(1+log(r))-p4
+//Overall diagram factor:    
+d1=mept/mep
+//% of HP cylinder output:
+P1=mep1/mept*100
+//% of HP cylinder output:
+P2=mep2/mept*100
+//% of HP cylinder output:
+P3=meplp/mept*100
+printf("\nRESULT\n")
+printf("\nActual mep referred to LP = %f kPa",mept)
+printf("\nHypothetical mep referred to LP = %f kPa",mep)
+printf("\nOverall diagram factor = %f",d1)
+printf("\nPercentage of HP, IP and LP cylinder outputs = %f, %f and %f percent",P1,P2,P3) 
+\ No newline at end of file