diff options
Diffstat (limited to '3825/CH5')
-rw-r--r-- | 3825/CH5/EX5.1/Ex5_1.sce | 15 | ||||
-rw-r--r-- | 3825/CH5/EX5.10/Ex5_10.sce | 19 | ||||
-rw-r--r-- | 3825/CH5/EX5.11/Ex5_11.sce | 23 | ||||
-rw-r--r-- | 3825/CH5/EX5.12/Ex5_12.sce | 26 | ||||
-rw-r--r-- | 3825/CH5/EX5.13/Ex5_13.sce | 36 | ||||
-rw-r--r-- | 3825/CH5/EX5.14/Ex5_14.sce | 38 | ||||
-rw-r--r-- | 3825/CH5/EX5.2/Ex5_2.sce | 18 | ||||
-rw-r--r-- | 3825/CH5/EX5.3/Ex5_3.sce | 11 | ||||
-rw-r--r-- | 3825/CH5/EX5.4/Ex5_4.sce | 16 | ||||
-rw-r--r-- | 3825/CH5/EX5.5/Ex5_5.sce | 18 | ||||
-rw-r--r-- | 3825/CH5/EX5.6/Ex5_6.sce | 22 | ||||
-rw-r--r-- | 3825/CH5/EX5.7/Ex5_7.sce | 16 | ||||
-rw-r--r-- | 3825/CH5/EX5.8/Ex5_8.sce | 31 | ||||
-rw-r--r-- | 3825/CH5/EX5.9/Ex5_9.sce | 11 |
14 files changed, 300 insertions, 0 deletions
diff --git a/3825/CH5/EX5.1/Ex5_1.sce b/3825/CH5/EX5.1/Ex5_1.sce new file mode 100644 index 000000000..e515b2005 --- /dev/null +++ b/3825/CH5/EX5.1/Ex5_1.sce @@ -0,0 +1,15 @@ +clc
+v1=0.8854 //in metre-cube/kg
+V=0.1 //in metre-cube
+m=V/v1 //total mass
+mprintf("m=%fkg\n",m)//ans vary due to roundoff error
+v2=v1
+vf=0.001053 //in metre-cube/kg
+vg=1.159 //in metre-cube/kg
+T=111.37 //in degree celsius
+X2=(v2-vf)/(vg-vf)//quality of steam
+mprintf("X2=%f\n",X2)//ans vary due to roundoff error
+mvapor=X2*m
+mprintf("Mass of vapor=%fkg\n",mvapor)//ans may vary due to roundoff error
+mliquid=m-mvapor
+mprintf("Mass of liquid=%fkg",mliquid)//ans may vary due to roundoff error
diff --git a/3825/CH5/EX5.10/Ex5_10.sce b/3825/CH5/EX5.10/Ex5_10.sce new file mode 100644 index 000000000..37ea59978 --- /dev/null +++ b/3825/CH5/EX5.10/Ex5_10.sce @@ -0,0 +1,19 @@ +clc
+P1=100*10^3
+v1=1
+v2=1/16
+n=1.25
+T1=300
+P2=P1*((v1/v2)^n)
+mprintf("P2=%fMPa\n",P2/(10^6))
+T2=(T1*P2*v2)/(P1*v1)
+mprintf("T2=%fK\n",T2)
+R=8.314
+W=(R*(T1-T2))/(n-1)
+mprintf("W=%fkJ/mol\n",W/1000)
+gama=1.4
+q=((R*(T2-T1))/(gama-1))+W
+mprintf("q=%fkJ/mol",q/1000)
+
+
+
diff --git a/3825/CH5/EX5.11/Ex5_11.sce b/3825/CH5/EX5.11/Ex5_11.sce new file mode 100644 index 000000000..f051cc492 --- /dev/null +++ b/3825/CH5/EX5.11/Ex5_11.sce @@ -0,0 +1,23 @@ +clc
+P1=200//pressure in kPa
+V1=1//volume in metre-cube
+R=8.314//universal gas constant
+T1=127 //temperature in degree celsius
+NHe=(P1*10^3*V1)/(R*(T1+273))//ideal gas law
+mprintf("NHe=%fmol\n",NHe)//ans vary due to roundoff error
+P2=400 //pressure in kPa
+V2=1//volume in metre-cube
+T2=227 //temperature in degree celsius
+NN2=(P2*10^3*V2)/(R*(T2+273))//ideal gas law
+mprintf("NN2=%fmol\n",NN2)//ans vary due to roundoff error
+CvHe=1.5*R
+CvN2=2.5*R
+Tf=((NHe*CvHe*(T1+273))+(NN2*CvN2*(T2+273)))/((NN2*CvN2)+(NHe*CvHe))//temperature in kelvin from linear equation
+mprintf("Tf=%fK\n",Tf)//ans vary due to roundoff error
+Vf=2//volume in metre-cube
+Pf=((NHe+NN2)*R*Tf)/Vf//ideal gas law
+mprintf("Pf=%fkPa",Pf/1000)//ans vary due to roundoff error
+
+
+
+
diff --git a/3825/CH5/EX5.12/Ex5_12.sce b/3825/CH5/EX5.12/Ex5_12.sce new file mode 100644 index 000000000..8cdeddbe6 --- /dev/null +++ b/3825/CH5/EX5.12/Ex5_12.sce @@ -0,0 +1,26 @@ +clc
+v1=0.0011145
+V=0.05
+mliquid=V/v1
+mprintf("Mass of liquid=%fkg\n",mliquid)//ans may vary due to roundoff error
+h1=719.12
+P1=0.792*10^3
+u1=h1-(P1*v1)
+mprintf("u1=%fkJ/kg\n",u1)//ans may vary due to roundoff error
+v2=1/mliquid
+mprintf("v2=%fmetre-cube/kg\n",v2)//ans may vary due to roundoff error
+u2=u1
+T=147.73 //temperature in Kelvin
+vf=1.088*10^-3
+hf=622.4
+vg=0.41845
+hg=2742.55
+X2=(v2-vf)/(vg-vf)
+mprintf("X2=%f\n",X2)//ans vary due to roundoff error
+h2=(X2*hg)+(1-X2)*hf
+mprintf("h2=%fkJ/kg\n",h2)//ans vary due to roundoff error
+P2=4.5*10^5
+u2=(h2*10^3)-(P2*v2)
+mprintf("u2=%fkJ/kg\n",u2/1000)//ans may vary due to roundoff error
+
+
diff --git a/3825/CH5/EX5.13/Ex5_13.sce b/3825/CH5/EX5.13/Ex5_13.sce new file mode 100644 index 000000000..54cf48385 --- /dev/null +++ b/3825/CH5/EX5.13/Ex5_13.sce @@ -0,0 +1,36 @@ +clc
+Pa=500*10^3
+T1=300 //temperature in kelvin
+T2=300
+P1=100*10^3
+Ta=(Pa*T1)/P1
+mprintf("Ta=%fK\n",Ta)
+Cv=20.93
+Cp=29.302
+ua=Ta*Cv
+u1=T1*Cv
+mprintf("ua-u1=%fkJ/mol\n",(ua-u1)/1000)//ans vary due to roundoff error
+q1a=ua-u1
+qa2=Cp*(T2-Ta)
+q=qa2
+mprintf("qa2=%fkJ/mol\n",qa2/1000)//ans vary due to roundoff error
+deltau=Cv*(T2-Ta)
+mprintf("u2-ua=%fkJ/mol\n",deltau/1000)//ans vary due to roundoff error
+W=q-deltau
+mprintf("W=%fkJ/mol\n",W/1000)//aans vary due to roundoff error
+R=8.314
+P1=100
+P2=500
+T=T1
+W=R*T*log(P1/P2)
+mprintf("W=%fkJ/mol\n",W/1000)//ans vary due to roundoff error
+gama=1.4
+Tb=T1*((P2/P1)^(gama-1))
+mprintf("Tb=%fK\n",Tb)//ans vary due to roundoff error
+ub=Cv*Tb
+u1=Cv*T1
+deltau=ub-u1
+mprintf("ub-u1=%fkJ/mol\n",deltau/1000)//ans vary due to roundoff error
+W=-deltau
+mprintf("Total Work done=%fkJ/mol\n",W/1000)//ans vary due to roundoff error
+
diff --git a/3825/CH5/EX5.14/Ex5_14.sce b/3825/CH5/EX5.14/Ex5_14.sce new file mode 100644 index 000000000..d1878fcd1 --- /dev/null +++ b/3825/CH5/EX5.14/Ex5_14.sce @@ -0,0 +1,38 @@ +clc
+P1=2//pressure in bar
+vg=0.8854//in metre-cube/kg
+Pa=1//pressure in bar
+K=250//spring constant in kN/m
+A=0.05//cross-sectional area in metre-square
+h1=2706.3//in kJ/kg
+v1=vg
+V1=vg/10
+P2=4 //pressure in bar
+V0=V1-((((P1-Pa)*10^5)*(A*A))/(K*10^3))//from force balance equation
+mprintf("V0=%fmetre-cube\n",V0)//ans vary due to roundoff error
+V2=((((P2-Pa)*10^5)*(A*A))/(K*10^3))+V0//from force balance equation
+mprintf("V2=%fmeter-cube\n",V2)//ans may vary due to roundoff error
+Tc=500//temperature in degree celsius
+Td=600//temperature in degree celsius
+vc=0.8892//in metre-cube/kg
+vd=1.0054//in metre-cube/kg
+v2=V2*10
+T=Tc+(((Td-Tc)/(vd-vc))*(v2-vc))//by interpolation
+mprintf("T=%fdegree-celsius\n",T)//ans vary due to roundoff error
+h2=3515.2//in kJ/kg
+u2=(h2*10^3)-(P2*10^5*v2)//first law of thermodynamics
+mprintf("u2=%fkJ/kg\n",u2/1000)//ans vary due to roundoff error
+W=(((P1+P2)*10^5)*(V2-V1))/2//from P-V diagram in textbook
+mprintf("W=%fkJ\n",W/1000)//ans vary due to roundoff error
+u1=(h1*10^3)-(P1*10^5*v1)//first law of thermodynamics
+deltaU=u2-u1
+Q=(deltaU/10^4)+(W/1000)//first law of thermodynamics
+mprintf("Q=%fkJ",Q)//ans vary due to roundoff error
+
+
+
+
+
+
+
+
diff --git a/3825/CH5/EX5.2/Ex5_2.sce b/3825/CH5/EX5.2/Ex5_2.sce new file mode 100644 index 000000000..b4c309aa3 --- /dev/null +++ b/3825/CH5/EX5.2/Ex5_2.sce @@ -0,0 +1,18 @@ +clc
+P1=100 //pressure in kPa
+V1=1 //volume in metre-cube
+T1=20 //temperature in celsius
+R=8.314*10^3
+T2=60 //temperature in celsius
+N=(P1*10^3*V1)/(R*(T1+273.15))//conversion of celsius to kelvin
+mprintf("N=%fkmol\n",N)//ans ma vary due to roundoff error
+Cv=12.4717 //in kJ/kmolK
+deltaT=T2-T1
+deltaU=N*Cv*deltaT
+mprintf("DeltaU=%fkJ\n",deltaU)//ans vary due to roundoff error
+W=-deltaU
+mprintf("w=%fkJ\n",W)//ans vary due to roundoff error
+P2=P1*(T2+273.15)/(T1+273.15)//conversion of degree to kelvin
+mprintf("P2=%fkPa",P2)//ans in the textbook is wrong
+
+
diff --git a/3825/CH5/EX5.3/Ex5_3.sce b/3825/CH5/EX5.3/Ex5_3.sce new file mode 100644 index 000000000..ae6b82b61 --- /dev/null +++ b/3825/CH5/EX5.3/Ex5_3.sce @@ -0,0 +1,11 @@ +clc
+hf=417.54 //in kJ/kg
+hg=2675.4 //in kJ/kg
+vf=0.0010434 //in metre-cube/kg
+vg=1.694 //in metre-cube/kg
+P=100 //pressure in kPa
+W=P*10^3*(vg-vf)
+mprintf("W=%fkJ\n",W*10^-3)//ans vary due to roundoff error
+q=hg-hf
+mprintf("q=%fkJ\n",q)
+
diff --git a/3825/CH5/EX5.4/Ex5_4.sce b/3825/CH5/EX5.4/Ex5_4.sce new file mode 100644 index 000000000..f11f7a95f --- /dev/null +++ b/3825/CH5/EX5.4/Ex5_4.sce @@ -0,0 +1,16 @@ +clc
+dw=999.8 //density of water in kg/metre-cube
+dice=916.23 //density of ice in kg/metre-cube
+P=100 //pressure in kPa
+vw=1/dw
+mprintf("vw=%fmetre-cube/kg\n",vw)//ans vary due to roundoff error
+vi=1/dice
+mprintf("vi=%fmetre-cube/kg\n",vi)
+W=P*10^3*(vw-vi)
+mprintf("W=%fJ\n",W)//ans may vary due to roundoff error
+deltaU=334.98 //in kJ
+q=deltaU+(W*10^-3)
+mprintf("q=%fkJ\n",q)//ans may vary due to roundoff error
+
+
+
diff --git a/3825/CH5/EX5.5/Ex5_5.sce b/3825/CH5/EX5.5/Ex5_5.sce new file mode 100644 index 000000000..efe819357 --- /dev/null +++ b/3825/CH5/EX5.5/Ex5_5.sce @@ -0,0 +1,18 @@ +clc
+hf=640.12 //in kJ/kg
+hg=2747.5 //in kJ/kg
+X1=0.75
+h1=(X1*hg)+(1-X1)*hf
+mprintf("h1=%fkJ/kg\n",h1)//ans vary due to roundoff error
+h2=2855.1 //in kJ/kg
+q=h2-h1
+mprintf("q=%fkJ/kg\n",q)//ans vary due to roundoff error
+vf=0.0010928 //in metre-cube/kg
+vg=0.3747 //in metre-cube/kg
+v1=(X1*vg)+(1-X1)*vf
+mprintf("v1=%fmetre-cube/kg\n",v1)//ans vary due to roundoff error
+P=5*10^5
+v2=0.425 //in metre-cub/kg
+W=P*(v2-v1)
+mprintf("W=%fkJ/kg\n",W*10^-3)//ans vary due to roundof error
+
diff --git a/3825/CH5/EX5.6/Ex5_6.sce b/3825/CH5/EX5.6/Ex5_6.sce new file mode 100644 index 000000000..a82734526 --- /dev/null +++ b/3825/CH5/EX5.6/Ex5_6.sce @@ -0,0 +1,22 @@ +clc
+R=8.314
+T=300
+a=422.546
+b=0.0373
+v1=30
+v2=5
+W=integrate('((R*T)/(v-b))-(a/(v.^2))','v',v1,v2)
+mprintf("W=%fkJ/kmol\n",W)//ans may vary due to roundoff error
+deltaU=a*(1/v1-1/v2)
+mprintf("U2-U1=%fkJ/kmol\n",deltaU)//ans vary due to roundoff error
+Q=deltaU+W
+mprintf("Q=%fkJ/kmol\n",Q)//ans vary due to roundoff error
+
+
+
+
+
+
+
+
+
diff --git a/3825/CH5/EX5.7/Ex5_7.sce b/3825/CH5/EX5.7/Ex5_7.sce new file mode 100644 index 000000000..3f01fd7d5 --- /dev/null +++ b/3825/CH5/EX5.7/Ex5_7.sce @@ -0,0 +1,16 @@ +clc
+T1=300 //temperature in kelvin
+P1=100 //pressure in kPa
+P2=2 //pressure in MPa
+gama=1.4 //Cp/Cv ratio
+T2=T1*(((P2*10^6)/(P1*10^3))^((gama-1)/gama))
+mprintf("T2=%fK\n",T2)//ans vary due to roundoff error
+Cr=(T1*P2*10^6)/(P1*10^3*T2)
+mprintf("Compression ratio=%f\n",Cr)//ans vary due to roundoff error
+R=8.314
+W=R*(T1-T2)/(gama-1)
+mprintf("W=%fkJ/mol",W/1000)//ans vary due to roundoff error
+
+
+
+
diff --git a/3825/CH5/EX5.8/Ex5_8.sce b/3825/CH5/EX5.8/Ex5_8.sce new file mode 100644 index 000000000..b457bebb2 --- /dev/null +++ b/3825/CH5/EX5.8/Ex5_8.sce @@ -0,0 +1,31 @@ +clc
+h1=3456.2//in kJ/kg
+v1=0.11608//in metre-cube
+P1=3//in MPa
+u1=(h1*10^3)-(P1*10^6*v1)
+mprintf("u1=%fkJ/kg\n",u1/1000)
+P2=100//in kPa
+sf=1.3027//in kJ/kg K
+sg=7.3598//in kJ/kg K
+s2=7.2345//in kJ/kg K
+s1=s2//isentropic process
+X2=(s1-sf)/(sg-sf)//entropy equation using quality factor
+mprintf("X2=%f\n",X2)//ans vary due to roundoff error
+hg=2675.4//in kJ/kg
+hf=417.54//in kJ/kg
+vg=1.694//in metre-cube/kg
+vf=0.0010434//in metre-cube/kg
+h2=(hg*X2)+(1-X2)*hf//enthalpy equation using quality factor
+mprintf("h2=%fkJ/kg\n",h2)//ans vary due to roundoff error
+v2=(vg*X2)+(1-X2)*vf//specific volume eqaution using quality factor
+mprintf("v2=%fmetre-cube/kg\n",v2)//ans vary due to roundoff error
+u2=(h2*10^3)-(P2*10^3*v2)//first law of thermodynamics
+mprintf("u2=%fkJ/kg\n",u2/1000)//ans vary due to roundoff error
+W=u1-u2//first law of thermodynamics
+mprintf("W=%fkJ/kg",W/1000)//ans vary due to roundoff errror
+
+
+
+
+
+
diff --git a/3825/CH5/EX5.9/Ex5_9.sce b/3825/CH5/EX5.9/Ex5_9.sce new file mode 100644 index 000000000..10b9833be --- /dev/null +++ b/3825/CH5/EX5.9/Ex5_9.sce @@ -0,0 +1,11 @@ +clc
+P1=100*10^3
+P2=10^6
+gama=1.4
+T1=300 //temperature in kelvin
+T2=T1*((P2/P1)^((gama-1)/gama))
+mprintf("T2=%fK\n",T2)//ans may vary due to roundoff error
+R=8.314
+W=(R*(T1-T2))/(gama-1)
+mprintf("W=%fJ/mol",W)//ans vary due to roundoff error
+
|