diff options
Diffstat (limited to '3825/CH8')
-rw-r--r-- | 3825/CH8/EX8.1/Ex8_1.sce | 26 | ||||
-rw-r--r-- | 3825/CH8/EX8.10/Ex8_10.sce | 9 | ||||
-rw-r--r-- | 3825/CH8/EX8.11/Ex8_11.sce | 11 | ||||
-rw-r--r-- | 3825/CH8/EX8.12/Ex8_12.sce | 11 | ||||
-rw-r--r-- | 3825/CH8/EX8.2/Ex8_2.sce | 9 | ||||
-rw-r--r-- | 3825/CH8/EX8.3/Ex8_3.sce | 14 | ||||
-rw-r--r-- | 3825/CH8/EX8.4/Ex8_4.sce | 9 | ||||
-rw-r--r-- | 3825/CH8/EX8.9/Ex8_9.sce | 13 |
8 files changed, 102 insertions, 0 deletions
diff --git a/3825/CH8/EX8.1/Ex8_1.sce b/3825/CH8/EX8.1/Ex8_1.sce new file mode 100644 index 000000000..dc426cf3b --- /dev/null +++ b/3825/CH8/EX8.1/Ex8_1.sce @@ -0,0 +1,26 @@ +clc
+P1=10 //pressure in bar
+P2=1 //pressure in bar
+h1=3052.1*10^3
+v1=0.258
+u1=h1-(P1*10^5*v1)
+mprintf("u1=%fkJ/kg\n",u1/1000)//ans vary due to roundoff error
+s1=7.1251 //in kJ/kg K
+s2=s1
+sg=7.3598
+sf=1.3027
+X2=(s2-sf)/(sg-sf)
+mprintf("X2=%f\n",X2)//ans vary due to roundoff error
+hg=2675.4
+hf=417.54
+h2=(X2*hg)+(1-X2)*hf
+mprintf("h2=%fkJ/kg\n",h2)//ans vary due to roundoff error
+vg=1.694
+vf=0.0010434
+v2=(X2*vg)+(1-X2)*vf
+mprintf("v2=%fmetre-cube/kg\n",v2)//ans vary due to roundoff error
+u2=(h2*10^3)-(P2*10^5*v2)
+mprintf("u2=%fkJ/kg\n",u2/1000)//ans vary due to roundoff error
+W=u1-u2
+mprintf("W=%fkJ",W/1000)//ans vary due to roundoff error
+
diff --git a/3825/CH8/EX8.10/Ex8_10.sce b/3825/CH8/EX8.10/Ex8_10.sce new file mode 100644 index 000000000..3ff556e64 --- /dev/null +++ b/3825/CH8/EX8.10/Ex8_10.sce @@ -0,0 +1,9 @@ +clc
+P2=200
+P1=101.325
+hfg=2256.94*18
+R=8.314
+T1=373
+T2=1/((1/T1)-((log(P2/P1))/(hfg/R)))//from the equation formed in book
+mprintf("T2=%fcelsius",T2)//ans vary due to roundoff error
+
diff --git a/3825/CH8/EX8.11/Ex8_11.sce b/3825/CH8/EX8.11/Ex8_11.sce new file mode 100644 index 000000000..99a45614a --- /dev/null +++ b/3825/CH8/EX8.11/Ex8_11.sce @@ -0,0 +1,11 @@ +clc
+T1=100
+T2=120
+CPf=4.23
+CPg=1.55
+h1=2256.94
+deltahfg=(CPg-CPf)*(T2-T1)
+mprintf("hfg2-hfg1=%f\n",deltahfg)
+hfg=h1+deltahfg
+mprintf("hfg at 120 celsius=%fkJ/kg",hfg)
+
diff --git a/3825/CH8/EX8.12/Ex8_12.sce b/3825/CH8/EX8.12/Ex8_12.sce new file mode 100644 index 000000000..8624f1930 --- /dev/null +++ b/3825/CH8/EX8.12/Ex8_12.sce @@ -0,0 +1,11 @@ +clc
+hsf=334.92
+T1=273.15
+T2=263.15
+CPf=4.186
+CPs=2.093
+a=(hsf/T1)+((CPf-CPs)*log(T2/T1))
+mprintf("(hsf/T)2=%f\n",a)//ans vary due to roundoff error
+deltas=-a
+mprintf("deltas=%fkJ/kg K",deltas)//ans vary due to roundoff error
+
diff --git a/3825/CH8/EX8.2/Ex8_2.sce b/3825/CH8/EX8.2/Ex8_2.sce new file mode 100644 index 000000000..c7ffa7d9d --- /dev/null +++ b/3825/CH8/EX8.2/Ex8_2.sce @@ -0,0 +1,9 @@ +clc
+R=8.314
+P2=10^6//pressure in pascals
+P1=100*10^3 //pressure in pascal
+deltas=-R*log(P2/P1)
+mprintf("s2-s1=%fkJ/kmol K\n",deltas)//ans vary due to roundoff error
+TO=300
+Wmin=-TO*deltas
+mprintf("Wmin=%fkJ/kmol K",Wmin)//ans vary due to roundoff error
diff --git a/3825/CH8/EX8.3/Ex8_3.sce b/3825/CH8/EX8.3/Ex8_3.sce new file mode 100644 index 000000000..100957904 --- /dev/null +++ b/3825/CH8/EX8.3/Ex8_3.sce @@ -0,0 +1,14 @@ +clc
+hi=2947.95
+si=7.0248
+se=si
+sg=7.3598
+sf=1.3027
+Xe=(se-sf)/(sg-sf)
+mprintf("Xe=%g\n",Xe)//ans vary due to roundoff error
+hg=2675.4
+hf=417.54
+he=(Xe*hg)+(1-Xe)*hf
+mprintf("he=%fkJ/kg\n",he)//ans vary due to roundoff error
+W=hi-he
+mprintf("W=%fkJ/s",W)//ans vary due to roundoff error
diff --git a/3825/CH8/EX8.4/Ex8_4.sce b/3825/CH8/EX8.4/Ex8_4.sce new file mode 100644 index 000000000..452a3201f --- /dev/null +++ b/3825/CH8/EX8.4/Ex8_4.sce @@ -0,0 +1,9 @@ +clc
+R=8.314
+TO=300
+X1=0.79
+X2=0.21
+deltag=R*TO*((X1*log(X1))+(X2*log(X2)))
+mprintf("ge-gi=%fJ/mol\n",deltag)//ans vary due to roundoff error
+W=(-deltag*10^3)/3600
+mprintf("W=%fW\n",W)//ans vary due to roundoff error
diff --git a/3825/CH8/EX8.9/Ex8_9.sce b/3825/CH8/EX8.9/Ex8_9.sce new file mode 100644 index 000000000..d47fde36e --- /dev/null +++ b/3825/CH8/EX8.9/Ex8_9.sce @@ -0,0 +1,13 @@ +clc
+P1=0.1 //pressure in MPa
+P2=10 //pressure in MPa
+T=298
+Cp=4.2
+beeta=2.07*10^-4
+vf=0.0010029
+a=(T*vf*beeta)/(Cp*10^3)
+mprintf("delT/delP=%E\n",a)//ans vary due to roundoff error
+deltaP=P2-P1
+deltaT=a*deltaP*10^6
+mprintf("deltaT=%fcelsius",deltaT)//ans vary due to roundoff error
+
|