13 files changed, 332 insertions, 0 deletions
diff --git a/3825/CH10/EX10.1/Ex10_1.sce b/3825/CH10/EX10.1/Ex10_1.sce
new file mode 100644
index 000000000..d135bf559
--- /dev/null
+++ b/3825/CH10/EX10.1/Ex10_1.sce
@@ -0,0 +1,32 @@
+clc
+y1=0.4 //CH4
+y2=0.2 //C2H6
+y3=0.25 //H2
+y4=0.15 //N2
+M1=16
+M2=30
+M3=2
+M4=28
+m1=0.64
+m2=0.6
+m3=0.05
+m4=0.42
+M=(y1*M1)+(y2*M2)+(y3*M3)+(y4*M4)
+mprintf("M=%fkg/kmol\n",M)
+m=m1+m2+m3+m4
+mprintf("m=%fkg\n",m)
+yCH4=m1/m
+mprintf("phi1=%f\n",yCH4)//ans vary due to roundoff error
+yC2H6=m2/m
+mprintf("phi2=%f\n",yC2H6)//ans vary due to roundoff error
+yH2=m3/m
+mprintf("phi3=%f\n",yH2)//ans vary due to roundoff error
+yN2=m4/m
+mprintf("phi4=%f\n",yN2)//ans vary due to roundoff error
+N=100
+N1=N*y1
+N2=N*y2
+N3=N*y3
+N4=N*y4
+mprintf("Ni=.\n%f\n%f\n%f\n%f",N1,N2,N3,N4)
+
diff --git a/3825/CH10/EX10.10/Ex10_10.sce b/3825/CH10/EX10.10/Ex10_10.sce
new file mode 100644
index 000000000..cf18cc4bf
--- /dev/null
+++ b/3825/CH10/EX10.10/Ex10_10.sce
@@ -0,0 +1,36 @@
+clc
+phi1=0.9
+Ps=2.337
+pw=phi1*Ps
+mprintf("pw=%fkPa\n",pw)//ans vary due to roundoff error
+P=101.325
+Pa=P-pw
+w1=0.622*(pw/Pa)
+mprintf("w1=%f kg H2O/kg dry air\n",w1)//ans vary due to roundoff error
+Mw=pw/P
+mprintf("Mole fraction of water=%f\n",Mw)//ans vary due to roundoff error
+R=8.314*10^3
+T=293
+v1=(R*T)/(P*10^3)
+mprintf("v1=%fmetre-cube/kmol\n",v1)//ans vary due to roundoff error
+ma1=((1-Mw)*(10/v1))*28.97
+mprintf("Flow rate of dry air=%f kg/min\n",ma1)//ans vary due to roundoff error
+TDB=40
+phi2=0.2
+w2=0.009
+pw=1.3 //in kPa
+Mw=pw/P
+mprintf("Mole fraction of water vapor=%f\n",Mw)//ans vary due to roundoff error
+v2=(R*(TDB+273))/(P*10^3)
+mprintf("v2=%fmetre-cube/kmol\n",v2)//ans vary due to roundoff error
+ma2=((1-Mw)*(30/v2))*28.97
+mprintf("ma2=%fkg/min\n",ma2)//ans in textbook is wrong
+w3=(w2+w1*(ma1/ma2))/(1+(ma1/ma2))
+mprintf("w3=%f kg H20/kg dry air\n",w3)//ans vary due to roundoff error
+h1star=54
+h2star=64
+h3star=(h2star+h1star*(ma1/ma2))/(1+(ma1/ma2))
+mprintf("h3star=%fkJ/kg dry air",h3star)//ans vary due to roundoff error
+
+
+
diff --git a/3825/CH10/EX10.11/Ex10_11.sce b/3825/CH10/EX10.11/Ex10_11.sce
new file mode 100644
index 000000000..38f4ea15c
--- /dev/null
+++ b/3825/CH10/EX10.11/Ex10_11.sce
@@ -0,0 +1,16 @@
+clc
+w3=0.011
+w2=0.009
+w1=0.0118
+P=101.325 //in bar
+pw1=1.8
+pw2=1.3
+a=(w2-w3)/(w3-w1)
+mprintf("ma1/ma2=%f\n",a)//ans vary due to roundoff error
+pa1=P-pw1
+pa2=P-pw2
+V1=10
+T1=293.15
+T2=313.15
+V2=(pa1*V1*T2)/(a*pa2*T1)
+mprintf("V2=%fmetre-cube/min",V2)//ans vary due to roundoff error
diff --git a/3825/CH10/EX10.12/Ex10_12.sce b/3825/CH10/EX10.12/Ex10_12.sce
new file mode 100644
index 000000000..5458a09d0
--- /dev/null
+++ b/3825/CH10/EX10.12/Ex10_12.sce
@@ -0,0 +1,26 @@
+clc
+phi1=0.8
+TDB=40 //in celsius
+w1=0.0284
+pw1=4.4 //in kPa
+h1star=114
+P=101.325
+R=8.314*10^3
+V=1
+M=28.97
+T1=313.15
+pa1=P-pw1
+ma1=(pa1*10^3*V*M)/(R*T1)
+mprintf("ma1=%fkg/s dry air\n",ma1)//ans vary due to roundoff error
+w4=0.0056
+h4star=34.5
+w2=w4
+h2star=20.5
+mw=ma1*(w1-w2)
+mprintf("mw=%fkg/s\n",mw)//ans vary due to roundoff error
+hw3=25.21
+Qc=ma1*(h1star-h2star)-(mw*hw3)
+mprintf("Qc=%fkJ/s\n",Qc)//ans vary due to roundoff error
+ma2=ma1
+Q=ma2*(h4star-h2star)
+mprintf("Q=%fkJ/s",Q)//ans vary due to roundoff error
diff --git a/3825/CH10/EX10.14/Ex10_14.sce b/3825/CH10/EX10.14/Ex10_14.sce
new file mode 100644
index 000000000..6f07d9d04
--- /dev/null
+++ b/3825/CH10/EX10.14/Ex10_14.sce
@@ -0,0 +1,10 @@
+clc
+TDB=49
+phi1=0.2
+w1=0.0147
+w2=0.023
+ma1=1
+mw=ma1*(w2-w1)
+mprintf("mw=%fkg/s\n",mw)
+mprintf("Water required to operate coler for 5 hours=%fkg",mw*5*3600)
+
diff --git a/3825/CH10/EX10.15/Ex10_15.sce b/3825/CH10/EX10.15/Ex10_15.sce
new file mode 100644
index 000000000..d8a86710f
--- /dev/null
+++ b/3825/CH10/EX10.15/Ex10_15.sce
@@ -0,0 +1,11 @@
+clc
+w1=0.009
+h1star=64.5
+w2=0.029
+h2star=109.5
+h3=167.45
+h4=125.66
+mw3=1000
+ma=(mw3*(h3-h4))/((h2star-h1star)-((w2-w1)*h4))
+mprintf("ma=%fkg/s\n",ma)//ans vary due to roundoff error
+mprintf("Make up water required=%fkg/s",ma*(w2-w1))//ans vary due to roundoff error
diff --git a/3825/CH10/EX10.2/Ex10_2.sce b/3825/CH10/EX10.2/Ex10_2.sce
new file mode 100644
index 000000000..e86812a96
--- /dev/null
+++ b/3825/CH10/EX10.2/Ex10_2.sce
@@ -0,0 +1,44 @@
+clc
+y1=0.2 //O2
+y2=0.5 //N2
+y3=0.1 //H2O
+y4=0.2 //CO2
+P=100 //in kPa
+p1=P*y1
+p2=P*y2
+p3=P*y3
+p4=P*y4
+mprintf("pi=.(in kPa)\n%f\n%f\n%f\n%f\n",p1,p2,p3,p4)
+M1=32
+M2=28
+M3=18
+M4=44
+M=(y1*M1)+(y2*M2)+(y3*M3)+(y4*M4)
+mprintf("Molar mass of mixture=%ikg/kmol\n",M)
+m1=M1*y1
+m2=M2*y2
+m3=M3*y3
+m4=M4*y4
+mprintf("mi=.(kg*10^3)\n%f\n%f\n%f\n%f\n",m1,m2,m3,m4)
+phi1=m1/M
+phi2=m2/M
+phi3=m3/M
+phi4=m4/M
+mprintf("phi i=.\n%f\n%f\n%f\n%f\n",phi1,phi2,phi3,phi4)
+P=p1+p2+p3
+mprintf("Final pressure=%ikPa\n",P)
+y1=p1/P
+y2=p2/P
+y3=p3/P
+mprintf("yi=.\n%f\n%f\n%f\n",y1,y2,y3)
+m1=M1*y1
+m2=M2*y2
+m3=M3*y3
+mprintf("mi=.\n%f\n%f\n%f\n",m1,m2,m3)
+phi1=m1/M
+phi2=m2/M
+phi3=m3/M
+mprintf("phi i=.\n%f\n%f\n%f\n",phi1,phi2,phi3)
+M=(y1*M1)+(y2*M2)+(y3*M3)
+mprintf("Molar mass=%f\n",M)//ans vary due to roundoff error
+
diff --git a/3825/CH10/EX10.4/Ex10_4.sce b/3825/CH10/EX10.4/Ex10_4.sce
new file mode 100644
index 000000000..079368bcd
--- /dev/null
+++ b/3825/CH10/EX10.4/Ex10_4.sce
@@ -0,0 +1,34 @@
+clc
+y1=0.75 //H2
+y2=0.25 //N2
+CP1=28.6455
+CP2=29.1783
+CP=(y1*CP1)+(y2*CP2)
+mprintf("CP=%fkJ/kmol K\n",CP)//ans vary due to roundoff error
+Cv1=20.3311
+Cv2=20.8641
+Cv=(y1*Cv1)+(y2*Cv2)
+mprintf("Cv=%fkJ/kmol K\n",Cv)//ans vary due to roundoff error
+gama=CP/Cv
+mprintf("gamma=%f\n",gama)//ans vary due to roundoff error
+P1=100 //pressure in kPa
+P2=500 //pressure in kPa
+T1=300
+T2=T1*((P2/P1)^((gama-1)/gama))
+mprintf("T2=%fK\n",T2)//ans vary due to roundoff error
+ws=-CP*(T2-T1)
+mprintf("-ws=%fkJ/kmol\n",-ws)//ans vary due to roundoff error
+M1=2.016
+M2=28.013
+M=(y1*M1)+(y2*M2)
+mprintf("Molar mass=%fkg/kmol\n",M)//ans vary due to roundoff error
+Ws=-(-ws/M)
+mprintf("-Ws=%fkJ/kg of mixture\n",-Ws)//ans vary due to roundoff error
+R=8.314
+deltas1=(CP1*log(T2/T1))-(R*log(P2/P1))
+mprintf("s2-s1=%fkJ/kmol K\n",deltas1)//ans vary due to roundoff error
+deltas2=(CP2*log(T2/T1))-(R*log(P2/P1))
+mprintf("s2-s1=%fkJ/kmol K\n",deltas2)//ans vary due to roundoff error
+deltas=(y1*deltas1)+(y2*deltas2)
+mprintf("s2-s1=%fkJ/kmol K",deltas)//ans vary due to roundoff error
+
diff --git a/3825/CH10/EX10.5/Ex10_5.sce b/3825/CH10/EX10.5/Ex10_5.sce
new file mode 100644
index 000000000..5eba311f5
--- /dev/null
+++ b/3825/CH10/EX10.5/Ex10_5.sce
@@ -0,0 +1,23 @@
+clc
+y1=0.6 //H2
+y2=0.4 //CH4
+p1=60 //in kPa
+p2=40 //in kPa
+T=300 //in kelvin
+V=5 //in metre-cube
+R=8.314*10^3
+N1=(p1*10^3*V)/(R*T)
+mprintf("NH2=%fkmol\n",N1)//ans vary due to roundoff error
+N2=(p2*10^3*V)/(R*T)
+mprintf("NCH4=%fkmol\n",N2)//ans vary due to roundoff error
+y1=0.5
+y2=0.5
+mprintf("Total methane to be added=%fkmol\n",(N1-N2))//ans vary due to roundoff error
+N2=N1
+N=N1+N2
+Pf=(N*R*T)/V
+mprintf("Pf=%fkPa",Pf/1000)//ans vary due to roundoff error
+
+
+
+
diff --git a/3825/CH10/EX10.6/Ex10_6.sce b/3825/CH10/EX10.6/Ex10_6.sce
new file mode 100644
index 000000000..197b2eac0
--- /dev/null
+++ b/3825/CH10/EX10.6/Ex10_6.sce
@@ -0,0 +1,37 @@
+clc
+P1=100 //pressure in kPa
+V=2 //in metre-cube
+R=8.314*10^3
+T0=300 //in kelvin
+N1=(P1*10^3*V)/(R*T0)
+mprintf("N1=%fkmol\n",N1)//ans vary due to roundoff error
+NHe0=0.0401
+Nair0=0.0401
+CvHe=12.4717
+Cvair=20.7889
+U0=((NHe0*CvHe)+(Nair0*Cvair))*T0
+mprintf("Uo=%fkJ\n",U0)//ans vary due to roundoff error
+Pf=4 //pressure in MPa
+a=(Pf*10^6*V)/R
+CPHe=20.786
+THe=600
+Nairf=0.0401
+b=a/(CPHe*THe)
+c=-NHe0+(U0/(CPHe*THe))
+d=c+Nairf-(CvHe*b)
+e=(Nairf*c)-(Nairf*Cvair*b)
+f=1
+//the above are the coefficients of quadratic equation formed for NHEf obtained from equation formed in book
+NHef=(-d+sqrt(d^2-4*f*e))/(2*f)
+mprintf("NHe=%fkmol\n",NHef)//ans vary due to roundoff error
+Tf=a/(NHef+Nairf)
+mprintf("Tf=%fK\n",Tf)//ans vary due to roundoff error
+yHe=NHef/(NHef+Nairf)
+mprintf("yHe=%f\n",yHe)//ans vary due to roundoff error
+yair=Nairf/(NHef+Nairf)
+mprintf("yair=%f\n",yair)//ans vary due to roundoff error
+mprintf("Helium that enters tank=%fkmol",NHef-Nairf)//ans vary due to roundoff error
+
+
+
+
diff --git a/3825/CH10/EX10.7/Ex10_7.sce b/3825/CH10/EX10.7/Ex10_7.sce
new file mode 100644
index 000000000..7d0fc73f4
--- /dev/null
+++ b/3825/CH10/EX10.7/Ex10_7.sce
@@ -0,0 +1,22 @@
+clc
+Ps=4.241 //pressure in kPa
+phi=0.6
+pw=Ps*phi
+mprintf("Partial pressure of water vapor=%fkPa\n",pw)//ans vary due to roundoff error
+P=101.325 //in kPa
+pa=P-pw
+mprintf("pa=%fkPa\n",pa)//ans vary due to roundoff error
+w=0.622*(pw/pa)
+mprintf("w=%f kg water/kg dry air\n",w)//ans vary due to roundoff error
+Mw=18
+R=8.314*10^3
+T=303 //temperature in kelvin
+V=100 //in metre-cube
+mprintf("Mass of water vapor=%fkg\n",(pw*10^3*V*Mw)/(R*T))//ans vary due to roundoff error
+Ma=28.97
+mprintf("Mass of dry air=%fkg\n",(pa*10^3*V*Ma)/(R*T))//ans vary due to roundoff error
+Ps=8.2578
+phi=0.9
+pw=Ps*phi
+mprintf("pw=%fkPa",pw)//ans vary due to roundoff error
+
diff --git a/3825/CH10/EX10.8/Ex10_8.sce b/3825/CH10/EX10.8/Ex10_8.sce
new file mode 100644
index 000000000..916e6bd9d
--- /dev/null
+++ b/3825/CH10/EX10.8/Ex10_8.sce
@@ -0,0 +1,21 @@
+clc
+Ps=4.241 //in kPa
+pw2=Ps
+P=101.325 //in kPa
+pa=P-pw2
+w2=0.622*(pw2/pa)
+mprintf("w2=%f kg H2O/kg dry air\n",w2)//ans vary due to roundoff error
+hfg2=2430.74
+hg1=2574.4
+hf2=125.66
+Cp=1.005
+T1=40 //in celsius
+T2=30 //in celsius
+w1=(Cp*(T2-T1)+w2*hfg2)/(hg1-hf2)
+mprintf("w1=%f kg H2O/kg dry air\n",w1)//ans vary due to roundoff error
+pw=P/(1+(0.622/w1))
+mprintf("pw=%fkPa\n",pw)//ans vary due to roundoff error
+Ps=7.375
+phi1=pw/Ps
+mprintf("phi1=%f",phi1)//ans vary due to roundoff error
+
diff --git a/3825/CH10/EX10.9/Ex10_9.sce b/3825/CH10/EX10.9/Ex10_9.sce
new file mode 100644
index 000000000..df3ea3d97
--- /dev/null
+++ b/3825/CH10/EX10.9/Ex10_9.sce
@@ -0,0 +1,20 @@
+clc
+//Ex uses data from Ex10_8.sce
+//we first write Ex10_7 for it
+Ps=4.241 //in kPa
+pw2=Ps
+P=101.325 //in kPa
+pa=P-pw2
+w2=0.622*(pw2/pa)
+hfg2=2430.74
+hg1=2574.4
+hf2=125.66
+Cp=1.005
+T1=40 //in celsius
+T2=30 //in celsius
+w1=(Cp*(T2-T1)+w2*hfg2)/(hg1-hf2)
+pw=P/(1+(0.622/w1))
+ha1=Cp*T1
+h1star=ha1+(w1*hg1)
+mprintf("h1star=%fkJ/kg dry air",h1star)//ans vary due to roundoff error
+