initial commit / add all books

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