initial commit / add all books

23 files changed, 511 insertions, 0 deletions

diff --git a/32/CH3/EX3.01/3_01.sce b/32/CH3/EX3.01/3_01.sce
new file mode 100755
index 000000000..c32e266fd
--- /dev/null
+++ b/32/CH3/EX3.01/3_01.sce
@@ -0,0 +1,18 @@
+//pathname=get_absolute_file_path('3.01.sce')

+//filename=pathname+filesep()+'3.01-data.sci'

+//exec(filename)

+//Pressure in the gas cylinder(in kPa):

+p=689

+//Final volume(in m^3):

+v2=0.045

+//Initial volume(in m^3):

+v1=0.04

+//Work done by the paddle(in kJ):

+Pw=-4.88

+//Work done by the system on the piston(in kJ):

+w=p*(v2-v1)

+//Net Work of the system(in kJ):

+wn=w+Pw

+printf("\nRESULTS\n")

+printf("\nWork done on the piston=%f kJ",w)

+printf("\nWork done on the system=%f kJ",-wn)
\ No newline at end of file
diff --git a/32/CH3/EX3.02/3_02.sce b/32/CH3/EX3.02/3_02.sce
new file mode 100755
index 000000000..224cc36f7
--- /dev/null
+++ b/32/CH3/EX3.02/3_02.sce
@@ -0,0 +1,13 @@
+//pathname=get_absolute_file_path('3.02.sce')

+//filename=pathname+filesep()+'3.02-data.sci'

+//exec(filename)

+//Mass of the gas(in kg):

+m=0.5

+//Initial internal energy(in kJ/kg):

+u1=26.6

+//Final internal energy(in kJ/kg):

+u2=37.8

+//Heat required(in kJ):

+Q=(u2-u1)*m

+printf("\nRESULT\n")

+printf("Heat required= %f kJ",Q)
\ No newline at end of file
diff --git a/32/CH3/EX3.03/3_03.sce b/32/CH3/EX3.03/3_03.sce
new file mode 100755
index 000000000..fac03498f
--- /dev/null
+++ b/32/CH3/EX3.03/3_03.sce
@@ -0,0 +1,15 @@
+//pathname=get_absolute_file_path('3.03.sce')

+//filename=pathname+filesep()+'3.03-data.sci'

+//exec(filename)

+//Mass flow rate(in kg/hr):

+m=50

+//Initial temp(in C):

+t1=800

+//Final temp(in C):

+t2=50

+//Heat capacity at const pressure(in kJ/kg.K):

+Cp=1.08

+//Heat to be removed(in kJ/hr):

+Q=m*Cp*(t2-t1)

+printf("\nRESULT\n")

+printf("Heat should be removed at %d kJ/hr",-Q)
\ No newline at end of file
diff --git a/32/CH3/EX3.04/3_04.sce b/32/CH3/EX3.04/3_04.sce
new file mode 100755
index 000000000..506944b21
--- /dev/null
+++ b/32/CH3/EX3.04/3_04.sce
@@ -0,0 +1,12 @@
+//pathname=get_absolute_file_path('3.04.sce')

+//filename=pathname+filesep()+'3.04-data.sci'

+//exec(filename)

+//Volume of the cylinnder(in m^3):

+v=0.78

+//Atmospheric pressure(in kPa):

+p=101.325

+//Work done(in kJ):

+w=p*v

+printf("\nRESULT\n")

+printf("\nWork done by air= %f",-w)

+printf("\nWork done by surroundings= %f",w)
\ No newline at end of file
diff --git a/32/CH3/EX3.05/3_05.sce b/32/CH3/EX3.05/3_05.sce
new file mode 100755
index 000000000..c27e622b1
--- /dev/null
+++ b/32/CH3/EX3.05/3_05.sce
@@ -0,0 +1,25 @@
+//pathname=get_absolute_file_path('3.05.sce')

+//filename=pathname+filesep()+'3.05-data.sci'

+//exec(filename)

+//Mass of the gas(in kg):

+m=5

+//Value of n in P*(V^n)=const:

+n=1.3

+//Initial pressure(in MPa):

+p1=1

+//Initial volume(in m^3):

+v1=0.5

+//Final pressure(in MPa):

+p2=0.5

+//Final volume(in m^3):

+v2=v1*((p1/p2)^(1/n))

+//Work done(in kJ):

+w=(p2*v2-p1*v1)*10^3/(1-n)

+//Change in internal energy(in kJ/kg):

+du=1.8*(p2*v2-p1*v1)*10^3

+//Heat interaction(in kJ):

+Q=du+w

+printf("\nRESULT\n")

+printf("\nHeat interaction = %f kJ",Q)

+printf("\nWork interaction = %f kJ",w)

+printf("\nChange in internal energy = %f kJ",du)
\ No newline at end of file
diff --git a/32/CH3/EX3.06/3_06.sce b/32/CH3/EX3.06/3_06.sce
new file mode 100755
index 000000000..f6a724ee8
--- /dev/null
+++ b/32/CH3/EX3.06/3_06.sce
@@ -0,0 +1,27 @@
+//pathname=get_absolute_file_path('3.06.sce')

+//filename=pathname+filesep()+'3.06-data.sci'

+//exec(filename)

+//Initial pressure(in MPa):

+p1=1

+//Final pressure(in MPa):

+p2=2

+//Initial volume(in m^3):

+v1=0.05

+//Value of n:

+n=1.4

+//Final volume(in m^3):

+v2=v1*((p1/p2)^(1/n))

+//Change in internal energy(in kJ/kg):

+du=7.5*(p2*v2-p1*v1)*10^3

+//Work done(in kJ):

+w=(p2*v2-p1*v1)*10^3/(1-n)

+//Heat interaction(in kJ):

+Q=du+w

+printf("\nRESULT\n")

+printf("\nHeat interaction = %f kJ",Q)

+printf("\nWork interaction = %f kJ",w)

+printf("\nChange in internal energy = %f kJ",du)

+//If 180 kJ heat transfer takes place:

+//Work done(in kJ):

+w2=180-du

+printf("\nNew work = %f kJ",w2)
\ No newline at end of file
diff --git a/32/CH3/EX3.07/3_08.sce b/32/CH3/EX3.07/3_08.sce
new file mode 100755
index 000000000..88763f520
--- /dev/null
+++ b/32/CH3/EX3.07/3_08.sce
@@ -0,0 +1,13 @@
+//pathname=get_absolute_file_path('3.08.sce')

+//filename=pathname+filesep()+'3.08-data.sci'

+//exec(filename)

+//Initial temperature(in K):

+t1=627+273

+//Final temperature(in K):

+t2=27+273

+//Specific heat at const pressure(in kJ/kg.K):

+Cp=1.005

+//Exit velocity(in m/s):

+c2=sqrt(2*Cp*10^3*(t1-t2))

+printf("\nRESULT\n")

+printf("\nExit Velocity = %f m/s",c2)
\ No newline at end of file
diff --git a/32/CH3/EX3.08/3_08.sce b/32/CH3/EX3.08/3_08.sce
new file mode 100755
index 000000000..88763f520
--- /dev/null
+++ b/32/CH3/EX3.08/3_08.sce
@@ -0,0 +1,13 @@
+//pathname=get_absolute_file_path('3.08.sce')

+//filename=pathname+filesep()+'3.08-data.sci'

+//exec(filename)

+//Initial temperature(in K):

+t1=627+273

+//Final temperature(in K):

+t2=27+273

+//Specific heat at const pressure(in kJ/kg.K):

+Cp=1.005

+//Exit velocity(in m/s):

+c2=sqrt(2*Cp*10^3*(t1-t2))

+printf("\nRESULT\n")

+printf("\nExit Velocity = %f m/s",c2)
\ No newline at end of file
diff --git a/32/CH3/EX3.09/3_09.sce b/32/CH3/EX3.09/3_09.sce
new file mode 100755
index 000000000..f814a1168
--- /dev/null
+++ b/32/CH3/EX3.09/3_09.sce
@@ -0,0 +1,15 @@
+//pathname=get_absolute_file_path('3.09.sce')

+//filename=pathname+filesep()+'3.09-data.sci'

+//exec(filename)

+//Work interaction(in kJ):

+w=-200

+//Increase in enthalpy(in kJ/kg):

+dh=100

+//Heat picked up by the cooling water(in kJ/kg):

+qc=-90

+//Heat flow(in kJ/kg):

+Q=dh+w

+//Heat transferred to atmosphere(in kJ/kg):

+Qa=Q-qc

+printf("\nRESULT\n")

+printf("\nHeat transferred to atmosphere = %d kJ/kg",Qa)
\ No newline at end of file
diff --git a/32/CH3/EX3.10/3_10.sce b/32/CH3/EX3.10/3_10.sce
new file mode 100755
index 000000000..0e9c52b76
--- /dev/null
+++ b/32/CH3/EX3.10/3_10.sce
@@ -0,0 +1,24 @@
+//pathname=get_absolute_file_path('3.10.sce')

+//filename=pathname+filesep()+'3.10-data.sci'

+//exec(filename)

+//Seating capacity:

+c=500

+//Heat requirement per person(in kcal/hr):

+q=50

+//Enthalpy of water entering the pipe(in kcal/kg):

+h1=80

+//Enthalpy of water leaving the pipe(in kcal/kg):

+h2=45

+//Difference in elevation of inlet and exit pipe(in m):

+z=10

+//Acceleration due to gravity(in m/s^2):

+g=9.81

+//Heat to be supplied(in kcal/hr):

+Q=c*q

+//Heat lost by water(in kcal/kg):

+Ql=-Q

+//By SFEE:

+//Quantity of water circulated(in kg/hr):

+m=(Ql*10^3*4.18)/(g*z+(h2-h1)*10^3*4.18)

+printf("\nRESULT\n")

+printf("\nWater circulation rate = %f kg/min",m/60)
\ No newline at end of file
diff --git a/32/CH3/EX3.11/3_11.sce b/32/CH3/EX3.11/3_11.sce
new file mode 100755
index 000000000..c9afa91a8
--- /dev/null
+++ b/32/CH3/EX3.11/3_11.sce
@@ -0,0 +1,22 @@
+//pathname=get_absolute_file_path('3.11.sce')

+//filename=pathname+filesep()+'3.11-data.sci'

+//exec(filename)

+//Enthalpy of steam entering the injector(in kcal/kg):

+h1=720

+//Enthalpy of water entering(in kcal/kg):

+h2=24.6

+//Enthalpy of water and steam mixture leaving the injector(in kcal/kg):

+h3=100

+//Depth of water injector from steam injector(in m):

+z=2

+//Velocity of steam entering the injector(in m/s):

+v1=50

+//Velocity of mixture leaving the injector(in m/s):

+v3=25

+//Heat loss from injector to surroundings(in kcal/kg):

+q=12

+//By applying SFEE:

+//Steam supply rate(in kg/s):

+m=(((v3^2)/2+h3*10^3*4.18)-(h2*10^3*4.18+g*z))/(((v1^2)/2+h1*10^3*4.18)-((v3^2)/2+h3*10^3*4.18)-(q*10^3*4.18))

+printf("\nRESULT\n")

+printf("\nSteam supply rate = %f kg/s",m)
\ No newline at end of file
diff --git a/32/CH3/EX3.12/3_12.sce b/32/CH3/EX3.12/3_12.sce
new file mode 100755
index 000000000..f64e975d4
--- /dev/null
+++ b/32/CH3/EX3.12/3_12.sce
@@ -0,0 +1,16 @@
+//pathname=get_absolute_file_path('3.12.sce')

+//filename=pathname+filesep()+'3.12-data.sci'

+//exec(filename)

+//Atmospheric pressure(in bar):

+p=1.013

+//Volume to which the baloon is inflated(in m^3):

+v=0.4

+//Work done by cylinder(in kJ):

+w1=0

+//Work done by the balloon(in kJ):

+w2=p*10^5*v

+//Total work(in kJ):

+w=w1+w2

+printf("\nRESULT\n")

+printf("\nWork done by the system upon atmoshere = %f kJ",w/(10^3))

+printf("\nWork done by the atmoshere = %f kJ",-w/(10^3))
\ No newline at end of file
diff --git a/32/CH3/EX3.13/3_13.sce b/32/CH3/EX3.13/3_13.sce
new file mode 100755
index 000000000..d63d78fda
--- /dev/null
+++ b/32/CH3/EX3.13/3_13.sce
@@ -0,0 +1,15 @@
+//pathname=get_absolute_file_path('3.13.sce')

+//filename=pathname+filesep()+'3.13-data.sci'

+//exec(filename)

+//Heat added(in J/s):

+Qa=5000

+//Turbine work(in J/s):

+Wt=0.25*Qa

+//Heat rejected(in J/s):

+Qr=0.75*Qa

+//Work by feed pump(in J/s):

+Wp=0.002*Qa

+//Capacity of generator(in W):

+C=Wt-Wp

+printf("\nRESULT\n")

+printf("\nCapacity of generator = %f kW ",C/(10^3))
\ No newline at end of file
diff --git a/32/CH3/EX3.14/3_14.sce b/32/CH3/EX3.14/3_14.sce
new file mode 100755
index 000000000..8018ac062
--- /dev/null
+++ b/32/CH3/EX3.14/3_14.sce
@@ -0,0 +1,30 @@
+//pathname=get_absolute_file_path('3.14.sce')

+//filename=pathname+filesep()+'3.14-data.sci'

+//exec(filename)

+//Ambient temperature(in K):

+T1=27+273

+//Temperature of air inside heat exchanger(in K):

+T2=750+273

+//Temperature of air leaving turbine(in K):

+T3=600+273

+//Temperature of air leaving the nozzle(in K):

+T4=500+273

+//Velocity of air entering turbine(in m/s):

+c2=50

+//Velocity of air entering the nozzle(in m/s):

+c3=60

+//Specific heat at constant pressure(in kj?kg.K):

+Cp=1.005

+//By applying SFEE between points 1 & 2:

+//Heat transfer to air in heat exchanger(in kJ):

+Q12=Cp*(T2-T1)

+printf("\nRESULT\n")

+printf("\nHeat transfer to air in heat exchanger =%f kJ",Q12)

+//By applying SFEE between points 2 & 3:

+//Power output from turbine(in kJ/s):

+Wt=Cp*(T2-T3)+(c2^2-c3^2)*10^(-3)/2

+printf("\nPower output from turbine = %f kJ/s",Wt)

+//By applying SFEE between points 3 & 4:

+//Velocity at exit of the nozzle(in m/s):

+c4=sqrt(2*(Cp*(T3-T4)+(c3^2)*10^(-3)/2))

+printf("\nVelocity at exit of the nozzle = %f m/s",c4)
\ No newline at end of file
diff --git a/32/CH3/EX3.15/3_15.sce b/32/CH3/EX3.15/3_15.sce
new file mode 100755
index 000000000..08cb1b911
--- /dev/null
+++ b/32/CH3/EX3.15/3_15.sce
@@ -0,0 +1,23 @@
+//pathname=get_absolute_file_path('3.15.sce')

+//filename=pathname+filesep()+'3.15-data.sci'

+//exec(filename)

+//Initial pressure(in MPa):

+p1=0.5

+//Initial temperature(in K):

+T1=400

+//Ratio of v2 to v1:

+r1=2

+//Ratio of v3 to v1:

+r2=6

+//Universal gas constant(in kJ/kg):

+R=8.314

+//Work from state 1 to 2(in kJ):

+Wa=R*T1

+//Temperature at point 2(in K):

+T2=2*T1

+//Work done from state 2 to 3(in kJ):

+Wb=R*T2*log(r2/r1)

+//Total work done by air(in kJ):

+W=Wa+Wb

+printf("\nRESULT\n")

+printf("\nWork done = %f kJ",W)
\ No newline at end of file
diff --git a/32/CH3/EX3.16/3_16.sce b/32/CH3/EX3.16/3_16.sce
new file mode 100755
index 000000000..cb31d0299
--- /dev/null
+++ b/32/CH3/EX3.16/3_16.sce
@@ -0,0 +1,17 @@
+//pathname=get_absolute_file_path('3.16.sce')

+//filename=pathname+filesep()+'3.16-data.sci'

+//exec(filename)

+//Initial pressure(in MPa):

+pi=0.5

+//Initial volume(in m^3):

+vi=0.5

+//Final pressure(in MPa):

+pf=1

+//Atmospheric pressure(in Pa):

+patm=1.013*10^5

+//Final volume(in m^3):

+vf=3*vi

+//Work done(in J):

+W=(vf-vi)*(pi+pf)*10^5/2

+printf("\nRESULT\n")

+printf("\nWork done = %d J",W)
\ No newline at end of file
diff --git a/32/CH3/EX3.17/3_17.sce b/32/CH3/EX3.17/3_17.sce
new file mode 100755
index 000000000..7473d3036
--- /dev/null
+++ b/32/CH3/EX3.17/3_17.sce
@@ -0,0 +1,37 @@
+//pathname=get_absolute_file_path('3.17.sce')

+//filename=pathname+filesep()+'3.17-data.sci'

+//exec(filename)

+//Initial pressure(in MPa):

+pi=0.5

+//Initial volume(in m^3):

+vi=0.5

+//Final pressure(in MPa):

+pf=1

+//Atmospheric pressure(in Pa):

+patm=1.013*10^5

+//Adiabatic index of compression for H2:

+rH2=CpH2/(CpH2-RH2)

+//Adiabatic index of compression for N2:

+rN2=CpN2/(CpN2-RN2)

+//Final pressure of hydrogen(in Pa):

+p2=p1*(v1/v2)^rH2

+printf("\nRESULT\n")

+printf("\nFinal pressure of hydrogen = %f MPa",p2/(10^6))

+//Partition work:

+Pw=0

+printf("\nPartition work = %d",Pw)

+//Work done upon H2(in J):

+WH2=(p1*v1-p2*v2)/(rH2-1)

+//Work done by nitrogen(in J):

+WN2=-WH2

+printf("\nWork done by hyrogen = %d J",WH2)

+printf("\nWork done by nitrogen = %d J",WN2)

+//Mass of N2(in kg):

+mN2=p1*v1/(RN2*10^3*T1)

+//Final temperature of N2(in K):

+T2=p2*vN2*T1/(p1*v1)

+//Cv of N2(in kJ/kg):

+CvN2=CpN2-RN2

+//Heat added to N2(in kJ):

+QN2=mN2*CvN2*10^3*(T2-T1)+WN2

+printf("\nHeat added to nitrogen = %f kJ",QN2/(10^3))
\ No newline at end of file
diff --git a/32/CH3/EX3.18/3_18.sce b/32/CH3/EX3.18/3_18.sce
new file mode 100755
index 000000000..edd2e1be3
--- /dev/null
+++ b/32/CH3/EX3.18/3_18.sce
@@ -0,0 +1,29 @@
+//pathname=get_absolute_file_path('3.18.sce')

+//filename=pathname+filesep()+'3.18-data.sci'

+//exec(filename)

+//Volume of the cylinder(in m^3):

+v1=2

+//Pressure in the cylinder(in Pa):

+p1=0.5*10^6

+//Temperature of the cylinder(in K):

+T1=375

+//Specific heat at const pressure(in kJ/kg.K):

+Cp=1.003

+//Specific heat at const volume(in kJ/kg.K):

+Cv=0.716

+//Gas constant for air(in kJ/kg.K):

+Ra=0.287

+//Atmospheric pressure(in Pa):

+patm=1.013*10^5

+//Compression ratio:

+r=1.4

+//Initial mass of air(in kg):

+m1=p1*v1/(Ra*T1)

+//Final temperature(in K):

+T2=T1*(patm/p1)^((r-1)/r)

+//Final mass of air left in tank(in kg):

+m2=patm*v1/(Ra*T2)

+//Kinetic energy available(in kJ):

+KE=m1*Cv*T1-m2*Cv*T2-(m1-m2)*Cp*T2

+printf("\nRESULT\n")

+printf("\nAmount of work available = %f J",KE)
\ No newline at end of file
diff --git a/32/CH3/EX3.19/3_19.sce b/32/CH3/EX3.19/3_19.sce
new file mode 100755
index 000000000..552410c9f
--- /dev/null
+++ b/32/CH3/EX3.19/3_19.sce
@@ -0,0 +1,28 @@
+//pathname=get_absolute_file_path('3.19.sce')

+//filename=pathname+filesep()+'3.19-data.sci'

+//exec(filename)

+//Pressure in the vessel(in Pa):

+p1=0.5*10^6

+//Volume of 1st chamber(in m^3):

+v1=0.5

+//Temperature in the vessel(in K):

+T1=300

+//Final pressure(in Pa):

+p2=10^6

+//Volume of 2nd chamber(in m^3):

+v2=0.5

+//Final temperature(in K):

+T2=500

+//Universal gas constant(in J/kg.K):

+R=8314

+//Moles in chamber 1:

+n1=p1*v1/(R*T1)

+//Moles in chamber 2:

+n2=p2*v2/(R*T2)

+//Temperature of the mixture(in K):

+T3=(n1*T1+n2*T2)/(n1+n2)

+//Final pressure(in MPa):

+p3=(n1+n2)*R*T3/(v1+v2)

+printf("\nRESULT\n")

+printf("\nFinal pressure = %f MPa",p3/(10^6))

+printf("\nFinal temperature = %f K",T3)
\ No newline at end of file
diff --git a/32/CH3/EX3.20/3_20.sce b/32/CH3/EX3.20/3_20.sce
new file mode 100755
index 000000000..cf910f4be
--- /dev/null
+++ b/32/CH3/EX3.20/3_20.sce
@@ -0,0 +1,13 @@
+//pathname=get_absolute_file_path('3.20.sce')

+//filename=pathname+filesep()+'3.20-data.sci'

+//exec(filename)

+//Volume of the bottle(in m^3):

+v=0.5

+//Pressure in the bottle(in Bar):

+p=1.0135

+//Displacement work(in N-m):

+W=p*10^5*(0-v)

+//Heat transfer(in N-m):

+Q=-W

+printf("\nRESULT\n")

+printf("\nHeat transferred = %d N-m",Q)
\ No newline at end of file
diff --git a/32/CH3/EX3.21/3_21.sce b/32/CH3/EX3.21/3_21.sce
new file mode 100755
index 000000000..29d8b510c
--- /dev/null
+++ b/32/CH3/EX3.21/3_21.sce
@@ -0,0 +1,38 @@
+//pathname=get_absolute_file_path('3.21.sce')

+//filename=pathname+filesep()+'3.21-data.sci'

+//exec(filename)

+//Volume of bottle(in m^3):

+v1=0.3

+//Pressure in the bottle(in bar):

+p1=35

+//Temperature in the bottle(in K):

+T1=40+273

+//Turbo generator's actual output(in kJ/s):

+w1=5

+//Final prssure(in bar):

+p2=1

+//Final volume(in m^3):

+v2=v1

+//Gas constant for air(in kJ/kg.K):

+Ra=0.287

+//Compression ratio:

+r=1.4

+//% of output which is consumed= 60%

+//Specific heat at const volume(in kJ/kg):

+Cv=0.718

+//Specific heat at const pressure(in kJ/kg):

+Cp=1.005

+//Final temperature(in K):

+T2=T1*(p2/p1)^((r-1)/r)

+//Initial mass in the bottle(in kg):

+m1=p1*10^2*v1/(Ra*T1)

+//Final mass in the bottle(in kg):

+m2=p2*10^2*v2/(Ra*T2)

+//Maximum work that can be obtained(in kJ):

+W=(m1*Cv*T1-m2*Cv*T2)-(m1-m2)*Cp*T2

+//Input to the turbo generator(in kJ/s):

+i=w1/0.6

+//Time duration(in s):

+t=W/i

+printf("\nRESULT\n")

+printf("\nDuration = %f seconds",t)
\ No newline at end of file
diff --git a/32/CH3/EX3.22/3_22.sce b/32/CH3/EX3.22/3_22.sce
new file mode 100755
index 000000000..92c1229b6
--- /dev/null
+++ b/32/CH3/EX3.22/3_22.sce
@@ -0,0 +1,37 @@
+//pathname=get_absolute_file_path('3.22.sce')

+//filename=pathname+filesep()+'3.22-data.sci'

+//exec(filename)

+//Pressure at state 1(in bar):

+p1=1.5

+//Temperature at state 1(in K):

+T1=77+273

+//Pressure at state 2(in bar):

+p2=7.5

+//Mass of the air(in kg):

+m=3

+//Value of n:

+n=1.2

+//Gas constant for air(in kJ/kg.K):

+Ra=0.287

+//Temperature at state 2(in K):

+T2=T1*(p2/p1)^((n-1)/n)

+//Initial volume(in m^3):

+v1=m*Ra*T1/(p1*10^2)

+//Volume at state 2(in m^3):

+v2=(p1*(v1^n)/p2)^(1/n)

+//Temperature at state 3(in K):

+T3=T1

+//Volume at state 3(in m^3):

+v3=v2*T3/T2

+//Pressure at state 3(in bar):

+p3=7.5

+//Compression work during process 1-2(in kJ):

+W12=m*Ra*(T2-T1)/(1-n)

+//Work during process 2-3(in kJ):

+W23=p2*(10^2)*(v3-v2)

+//Work during process 3-1(in kJ):

+W31=p3*10^2*v3*log(v1/v3)

+//Net work(in kJ):

+Wn=W12+W23+W31

+printf("\nRESULT\n")

+printf("\nHeat transferred from the system = %f kJ",-Wn)
\ No newline at end of file
diff --git a/32/CH3/EX3.23/3_23.sce b/32/CH3/EX3.23/3_23.sce
new file mode 100755
index 000000000..51ca87370
--- /dev/null
+++ b/32/CH3/EX3.23/3_23.sce
@@ -0,0 +1,31 @@
+//pathname=get_absolute_file_path('3.23.sce')

+//filename=pathname+filesep()+'3.23-data.sci'

+//exec(filename)

+//Volume of air bottle(in m^3):

+v1=0.15

+//Initial pressure(in bar):

+p1=40

+//Initial temperature(in K):

+T1=27+273

+//Final presure(in bar):

+p2=2

+//Final volume(in m^3):

+v2=v1

+//Gas constant for air(in kJ/kg):

+Ra=0.287

+//Specific heat at const pressure(in kJ/kg):

+Cp=1.005

+//Specific heat at const volume(in kJ/kg):

+Cv=0.718

+//Compression ratio:

+r=1.4

+//Initial mass of air in bottle(in kg):

+m1=p1*10^2*v1/(Ra*T1)

+//Final temperature(in K):

+T2=T1*(p2/p1)^((r-1)/r)

+//Final mass of air in bottle(in kg):

+m2=p2*10^2*v2/(Ra*T2)

+//Energy available for running of turbine due to emptying of bottle(in kJ):

+E=m1*Cv*T1-m2*Cv*T2-(m1-m2)*Cp*T2

+printf("\nRESULT\n")

+printf("\nWorj available from turbine = %f",E)
\ No newline at end of file