summaryrefslogtreecommitdiff
path: root/3825/CH3
diff options
context:
space:
mode:
Diffstat (limited to '3825/CH3')
-rw-r--r--3825/CH3/EX3.1/Ex3_1.sce41
-rw-r--r--3825/CH3/EX3.10/Ex3_10.sce6
-rw-r--r--3825/CH3/EX3.11/Ex3_11.sce13
-rw-r--r--3825/CH3/EX3.12/Ex3_12.sce16
-rw-r--r--3825/CH3/EX3.13/Ex3_13.sce12
-rw-r--r--3825/CH3/EX3.14/Ex3_14.sce6
-rw-r--r--3825/CH3/EX3.15/Ex3_15.sce25
-rw-r--r--3825/CH3/EX3.17/Ex3_17.sce11
-rw-r--r--3825/CH3/EX3.2/Ex3_2.sce30
-rw-r--r--3825/CH3/EX3.3/Ex3_3.sce13
-rw-r--r--3825/CH3/EX3.4/Ex3_4.sce10
-rw-r--r--3825/CH3/EX3.5/Ex3_5.sce16
-rw-r--r--3825/CH3/EX3.6/Ex3_6.sce21
13 files changed, 220 insertions, 0 deletions
diff --git a/3825/CH3/EX3.1/Ex3_1.sce b/3825/CH3/EX3.1/Ex3_1.sce
new file mode 100644
index 000000000..f570506a2
--- /dev/null
+++ b/3825/CH3/EX3.1/Ex3_1.sce
@@ -0,0 +1,41 @@
+clc
+//following is the function to solve cubic equation
+function[x1]=cub(a,b,c,d)
+ funcprot(0)
+ r=b/a;
+ s=c/a;
+ t=d/a;
+ p=(3*s-r^2)/3;
+ q=2*r^3/27-r*s/3+t;
+ D=(p/3)^3+(q/2)^2;
+ u=(-q/2+sqrt(D))^(1/3);
+ v=(-q/2-sqrt(D))^(1/3);
+ y1=u+v;
+ x1=y1-(r/3);
+ funcprot(0)
+endfunction
+R=8.314*(10^3) //universal gas constant
+T=300 //temperature in kelvin
+P=60*(10^6) //pressure in Pascals
+v=(R*T)/P //volume of cylinder in metre cube
+mprintf("v=%fmetre cube\n",v)//ans may vary due to roundoff error
+a=228.296*(10^(-3))//vander waals constant for methane in Pa(metre cube/mol)square
+b=0.043*(10^(-3))//vander waaals constant for methane in metre cube/mol
+R=8.314 //universal gas constant
+v=cub(P,(-(P*b)-(R*T)),a,(-a*b))//(P+a/v*v)*(v-b)=RT(van der Waals equation)
+mprintf("v=%fmetre-cube/mol\n",v)//ans may vary due to roundoff error
+Tc=190.7 //Tc for methane in kelvin
+Pc=46.41 //Pc for methane in bar
+a=0.42748*R*R*(Tc^2.5)/(Pc*(10^5)*(T^0.5)) //Redlich-Kwong equation
+b=0.0867*R*Tc/(Pc*10^5) //Redlich-Kwong equation
+mprintf("a=%fmetre^6Pa/mol-square\n",a)//ans may vary due to roundoff error
+mprintf("b=%fmetre-cube/mol\n",b)//ans may vary due to roundoff error
+v=cub(P,(-R*T),(-P*b*b-R*T*b+a),(-a*b))//P=(RT/v-b)-(a/v(v+b))
+mprintf("v=%fmetre cube/kmol\n",v*1000) //ans may vary due to roundoff error
+
+
+
+
+
+
+
diff --git a/3825/CH3/EX3.10/Ex3_10.sce b/3825/CH3/EX3.10/Ex3_10.sce
new file mode 100644
index 000000000..5887d29f0
--- /dev/null
+++ b/3825/CH3/EX3.10/Ex3_10.sce
@@ -0,0 +1,6 @@
+clc
+//Ex 3_7,3_8,3_9 and 3_10 use Molier Diagram
+h1=3275
+h2=2725
+deltah=h2-h1
+mprintf("deltah=%fkJ/kg",deltah)
diff --git a/3825/CH3/EX3.11/Ex3_11.sce b/3825/CH3/EX3.11/Ex3_11.sce
new file mode 100644
index 000000000..88f23ed95
--- /dev/null
+++ b/3825/CH3/EX3.11/Ex3_11.sce
@@ -0,0 +1,13 @@
+clc
+T=200//temperature in degree celsius
+P=15.549 //pressure in bars
+hf=852.37
+hg=2790.9
+hfg=hg-hf
+mprintf("Saturation pressure=%f bar\n",P)
+mprintf("Latent heat of vaporisation=%f kJ/kg\n",hfg)
+vg=0.1272 //in metre-cube/kg
+ug=(hg*10^3)-(P*10^5*vg)
+mprintf("ug=%fkJ/kg",ug/1000)//ans may vary due to roundoff error
+
+
diff --git a/3825/CH3/EX3.12/Ex3_12.sce b/3825/CH3/EX3.12/Ex3_12.sce
new file mode 100644
index 000000000..2311307f9
--- /dev/null
+++ b/3825/CH3/EX3.12/Ex3_12.sce
@@ -0,0 +1,16 @@
+clc
+X=0.85
+vg=0.8854 //in metre-cube/kg
+vf=0.001060 //in metre-cube/kg
+hf=504.7 //in kJ/kg
+hg=2706.3 //in kJ/kg
+v=(X*vg)+(1-X)*vf
+mprintf("v=%f metre-cube/kg\n",v)//ans may vary due to roundoff error
+h=(X*hg)+(1-X)*hf
+mprintf("h=%f kJ/kg\n",h)
+P=2 //pressure in bar
+u=(h*10^3)-(P*10^5*v)
+mprintf("u=%fkJ/kg",u/1000)//ans varies due to roundoff error
+
+
+
diff --git a/3825/CH3/EX3.13/Ex3_13.sce b/3825/CH3/EX3.13/Ex3_13.sce
new file mode 100644
index 000000000..b5ccaee59
--- /dev/null
+++ b/3825/CH3/EX3.13/Ex3_13.sce
@@ -0,0 +1,12 @@
+clc
+P1=30 //pressure in bar
+P2=35 //pressure in bar
+P3=32 //pressure in bar
+vg1=0.06663 //at P1
+vg2=0.05703 //at P2
+hg1=2802.3 //at P1
+hg2=2802 //at P2
+vg=vg1+((vg2-vg1)*(P3-P1)/(P2-P1))
+hg=hg1+((hg2-hg1)*(P3-P1)/(P2-P1))
+mprintf("vg=%f metre-cube/kg\n",vg)
+mprintf("hg=%f kJ/kg\n",hg)
diff --git a/3825/CH3/EX3.14/Ex3_14.sce b/3825/CH3/EX3.14/Ex3_14.sce
new file mode 100644
index 000000000..6577ba256
--- /dev/null
+++ b/3825/CH3/EX3.14/Ex3_14.sce
@@ -0,0 +1,6 @@
+clc
+hf=504.7
+hg=2706.3
+h=2600
+X=(h-hf)/(hg-hf)//from relation h=Xhg+(1-X)hf
+mprintf("X=%f",X)//ans may vary due to roundoff error
diff --git a/3825/CH3/EX3.15/Ex3_15.sce b/3825/CH3/EX3.15/Ex3_15.sce
new file mode 100644
index 000000000..1dc459050
--- /dev/null
+++ b/3825/CH3/EX3.15/Ex3_15.sce
@@ -0,0 +1,25 @@
+clc
+P1=24
+P2=26
+T1=300
+T2=400
+v1=0.10336 //at P1 and T1
+v2=0.09483 //at P2 and T1
+v3=0.12522 //at P1 and T2
+v4=0.11526 //at P2 and T2
+h1=3013.4
+h2=3007.4
+h3=3242.3
+h4=3239
+P3=25
+h5=h1+(((h2-h1)*(P3-P1))/(P2-P1))//interpolation at T=300
+v5=v1+(((v2-v1)*(P3-P1))/(P2-P1))//interpolation at T=300
+h6=h3+(((h4-h3)*(P3-P1))/(P2-P1))//interpolation at T=400
+v6=v3+(((v4-v3)*(P3-P1))/(P2-P1))//interpolation at T=400
+T3=350
+h7=h5+(((h6-h5)*(T3-T1))/(T2-T1))//interpolation at T=350
+v7=v5+(((v6-v5)*(T3-T1))/(T2-T1))//interpolation at T=350
+mprintf("v=%fmetre-cube/kg\n",v7)//ans may vary due to roundoff error
+mprintf("h=%fkJ/kg\n",h7)//ans may vary due to roundoff error
+
+
diff --git a/3825/CH3/EX3.17/Ex3_17.sce b/3825/CH3/EX3.17/Ex3_17.sce
new file mode 100644
index 000000000..33d58b906
--- /dev/null
+++ b/3825/CH3/EX3.17/Ex3_17.sce
@@ -0,0 +1,11 @@
+clc
+vc=0.00317
+vf=0.0010434
+vg=1.694
+X=(vc-vf)/(vg-vf)
+mprintf("X=%f\n",X)
+mprintf("volume of liquid=%fmetre-cube\n",(1-X)*vf)
+mprintf("Volume of vapor=%fmetre-cube\n",X*vg)
+mprintf("Ratio of liquid to vapor by volume=%f\n",((1-X)*vf)/(X*vg))//ans may vary due to roundoff error
+
+
diff --git a/3825/CH3/EX3.2/Ex3_2.sce b/3825/CH3/EX3.2/Ex3_2.sce
new file mode 100644
index 000000000..28e043a2d
--- /dev/null
+++ b/3825/CH3/EX3.2/Ex3_2.sce
@@ -0,0 +1,30 @@
+clc
+R=8.314*(10^3)//universal gas constant
+T=473 //temperature of bath in kelvin
+v1=0.6 //volume of steel vessel in metre-cube
+P=R*T/v1 //pressure developed by ideal gas law
+mprintf("P=%fMPa\n",P*10^-6)//ans may vary due to roundoff error
+a=453.046*(10^-3) //vander waals constant in Pa(metre-cube/mol)^2
+b=0.057*(10^-3)//vander waals constant in metre cube/mol
+P=((R*T)/(v1-b))-(a/(v1*v1))//pressure by vander waals equation
+mprintf("P=%fMPa\n",P*10^-6)//ans may vary due to roundofff error
+Pc=51.17 //pressure in bars
+Tc=283.1 //temperature in kelvin
+a=(0.42748*R*R*(Tc^2.5))/((Pc*10^5)*(T^0.5))//Redlich-Kwong equation
+b=0.0867*R*Tc/(Pc*(10^5))//Redlich-Kwong equation
+mprintf("a=%fPam^6/mol square\n",a*10^-6)//ans may vary due to roundoff error
+mprintf("b=%fm^3/mol\n",b*10^-3)//ans may vary due to roundoff error
+P=(R*T/(v1-b))-(a/(v1*(v1+b))) //pressure by Redlich-Kwong euation
+mprintf("P=%fMPa",P*10^-6)//ans may vary due to roundoff error
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/3825/CH3/EX3.3/Ex3_3.sce b/3825/CH3/EX3.3/Ex3_3.sce
new file mode 100644
index 000000000..27f7c2b44
--- /dev/null
+++ b/3825/CH3/EX3.3/Ex3_3.sce
@@ -0,0 +1,13 @@
+clc
+Pc=46.41*10^5 //pressure for methane in Pascals
+Tc=190.7 //temperature for methane in kelvin
+P=60*10^6 //methane pressure in pascals
+T=300 //methane temperature in kelvins
+Pr=P/Pc
+mprintf("Pr=%f\n",Pr)//ans may vary due to roundoff error
+Tr=T/Tc
+mprintf("Tr=%f\n",Tr)//ans may vary due to roundoff error
+Z=1.34
+R=8.314*10^3
+v=Z*R*T/P //volume of cylinder for storage
+mprintf("v=%fmetre-cube/kmol",v)//ans varies due to roundoff error
diff --git a/3825/CH3/EX3.4/Ex3_4.sce b/3825/CH3/EX3.4/Ex3_4.sce
new file mode 100644
index 000000000..1262f800b
--- /dev/null
+++ b/3825/CH3/EX3.4/Ex3_4.sce
@@ -0,0 +1,10 @@
+clc
+Tc=405.5 //temperature in kelvin
+Pc=112.77 //pressure in bar
+T=473 //temperature in kelvin
+Tr=T/Tc
+mprintf("Tr=%f\n",Tr)//ans vary due to roundoff error
+Pr=1.9 //obtained from compressibility chart in the book with given conditions
+P=Pc*Pr
+mprintf("P=%fMPa",P/10)//ans vary due to roundoff error
+
diff --git a/3825/CH3/EX3.5/Ex3_5.sce b/3825/CH3/EX3.5/Ex3_5.sce
new file mode 100644
index 000000000..5d403b83a
--- /dev/null
+++ b/3825/CH3/EX3.5/Ex3_5.sce
@@ -0,0 +1,16 @@
+clc
+Tc=405.5 //temperature in kelvin
+Pc=112.77 //pressure in bar
+T=473 //temperature in kelvin
+Tr=T/Tc //reduced temperature
+mprintf("Tr=%f\n",Tr)//ans vary due to roundoff error
+V=0.1
+n=10^3
+R=8.314
+vr=Pc*10^5*(V/n)/(R*Tc)//pseudo-reduced volume using compressibility chart
+mprintf("vr=%f\n",vr)//ans may vary due to roundoff error
+Pr=1.9
+P=Pr*Pc//equation for reduced pressure
+mprintf("P=%fMPa",P/10)//ans vary due to roundoff error
+
+
diff --git a/3825/CH3/EX3.6/Ex3_6.sce b/3825/CH3/EX3.6/Ex3_6.sce
new file mode 100644
index 000000000..8c87e5be7
--- /dev/null
+++ b/3825/CH3/EX3.6/Ex3_6.sce
@@ -0,0 +1,21 @@
+clc
+M=30 //molar mass for ethane in kg/kmol
+Tc=305.43 //temperature in kelvin
+Pc=4.884 //pressure in MPa
+P=70 //pressure in bar
+Pr=(P*10^5)/(Pc*10^6)
+mprintf("Pr=%f\n",Pr)//ans vary due to roundoff error
+M1=7 //mass of ethane present in kg
+Moles=(M1/M)*1000 //moles of ethane
+mprintf("Moles of ethane=%fmol\n",Moles)
+V=0.1 //volume of ethane in metre-cube
+v=V/Moles
+R=8.314
+vr=Pc*10^6*v/(R*Tc)
+mprintf("vr=%f\n",vr)//ans may vary due to roundoff error
+Tr=1.4
+T=Tr*Tc
+mprintf("T=%f K",T)//ans may vary due to roundoff error
+
+
+
generated by cgit (git 2.34.1) at 2024-12-22 09:46:24 +0000