diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3825/CH10/EX10.6 | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '3825/CH10/EX10.6')
| -rw-r--r-- | 3825/CH10/EX10.6/Ex10_6.sce | 37 |
1 files changed, 37 insertions, 0 deletions
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
+
+
+
+
|