summaryrefslogtreecommitdiff
path: root/551/CH9
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /551/CH9
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '551/CH9')
-rwxr-xr-x551/CH9/EX9.1/1.sce35
-rwxr-xr-x551/CH9/EX9.10/10.sce120
-rwxr-xr-x551/CH9/EX9.11/11.sce59
-rwxr-xr-x551/CH9/EX9.12/12.sce63
-rwxr-xr-x551/CH9/EX9.13/13.sce20
-rwxr-xr-x551/CH9/EX9.14/14.sce17
-rwxr-xr-x551/CH9/EX9.15/15.sce34
-rwxr-xr-x551/CH9/EX9.16/16.sce48
-rwxr-xr-x551/CH9/EX9.17/17.sce36
-rwxr-xr-x551/CH9/EX9.18/18.sce69
-rwxr-xr-x551/CH9/EX9.19/19.sce39
-rwxr-xr-x551/CH9/EX9.2/2.sce28
-rwxr-xr-x551/CH9/EX9.20/20.sce26
-rwxr-xr-x551/CH9/EX9.21/21.sce56
-rwxr-xr-x551/CH9/EX9.22/22.sce40
-rwxr-xr-x551/CH9/EX9.3/3.sce69
-rwxr-xr-x551/CH9/EX9.4/4.sce67
-rwxr-xr-x551/CH9/EX9.5/5.sce29
-rwxr-xr-x551/CH9/EX9.6/6.sce27
-rwxr-xr-x551/CH9/EX9.7/7.sce61
-rwxr-xr-x551/CH9/EX9.8/8.sce44
-rwxr-xr-x551/CH9/EX9.9/9.sce58
22 files changed, 1045 insertions, 0 deletions
diff --git a/551/CH9/EX9.1/1.sce b/551/CH9/EX9.1/1.sce
new file mode 100755
index 000000000..3845860c4
--- /dev/null
+++ b/551/CH9/EX9.1/1.sce
@@ -0,0 +1,35 @@
+clc
+V=0.35; //m^3
+m_CO=0.4; //kg
+m_air=1; //kg
+m_O2=0.233; //kg
+m_N2=0.767; //kg
+T=293; //K
+R0=8.314; //kJ/kg K
+M_O2=32; //Molecular mass of O2
+M_N2=28; //Molecular mass of N2
+M_CO=28; //Molecular mass of CO
+
+disp("Partial Pressures=")
+
+p_O2=m_O2*R0*10^3*T/M_O2/V/10^5; //bar
+disp("partial pressure for p_O2")
+disp(p_O2)
+disp("bar")
+
+p_N2=m_N2*R0*10^3*T/M_N2/V/10^5; //bar
+disp("partial pressure for p_N2")
+disp(p_N2)
+disp("bar")
+
+p_CO=m_CO*R0*10^3*T/M_CO/V/10^5; //bar
+disp("partial pressure for p_CO")
+disp(p_CO)
+disp("bar")
+
+
+disp("(ii) Total pressure in the vessel")
+p=p_O2+p_N2+p_CO;
+disp("p=")
+disp(p)
+disp("bar") \ No newline at end of file
diff --git a/551/CH9/EX9.10/10.sce b/551/CH9/EX9.10/10.sce
new file mode 100755
index 000000000..b1d914abf
--- /dev/null
+++ b/551/CH9/EX9.10/10.sce
@@ -0,0 +1,120 @@
+clc
+p=4*10^5; //Pa
+T=293; //K
+R0=8.314;
+
+m_N2=4; //kg
+m_CO2=6; //kg
+
+M_N2=28; //Molecular mass
+M_CO2=44; //Molecular mass
+
+n_N2=m_N2/M_N2; //moles of N2
+n_CO2=m_CO2/M_CO2; //moles of CO2
+
+x_N2=n_N2/(n_N2+n_CO2);
+disp("x_N2=")
+disp(x_N2)
+
+x_CO2=n_CO2/(n_CO2+n_N2);
+disp("x_CO2=")
+disp(x_CO2)
+
+
+disp("(ii) The equivalent molecular weight of the mixture")
+M=x_N2*M_N2 + x_CO2*M_CO2;
+disp("M=")
+disp(M)
+disp("kg/kg-mole")
+
+
+disp("(iii) The equivalent gas constant of the mixture")
+m=m_N2+m_CO2;
+Rmix=(m_N2*(R0/M_N2) + m_CO2*(R0/M_CO2))/m;
+disp("Rmix=")
+disp(Rmix)
+disp("kJ/kg K")
+
+
+disp("(iv) The partial pressures and partial volumes")
+P_N2=x_N2*p/10^5;
+disp("P_N2=")
+disp(P_N2)
+disp("bar")
+
+P_CO2=x_CO2*p/10^5;
+disp("P_CO2=")
+disp(P_CO2)
+disp("bar")
+
+V_N2=m_N2*R0/M_N2*T/p*10^3;
+disp("V_N2")
+disp(V_N2)
+disp("m^3")
+
+V_CO2=m_CO2*R0/M_CO2*T/p*10^3;
+disp("V_CO2")
+disp(V_CO2)
+disp("m^3")
+
+disp("(v) The volume and density of the mixture")
+
+V=m*Rmix*10^3*T/p;
+disp("V=")
+disp(V)
+disp("m^3")
+
+rho_mix=m/V;
+disp("Density of mixture=")
+disp(rho_mix)
+disp("kg/m^3")
+
+
+disp("(vi) cp and cv of the mixture")
+
+y_N2=1.4;
+cv_N2=(R0/M_N2)/(y_N2 - 1);
+cp_N2=cv_N2*y_N2;
+
+y_CO2=1.286;
+cv_CO2=(R0/M_CO2)/(y_CO2 - 1);
+cp_CO2=cv_CO2*y_CO2;
+
+cp=(m_N2*cp_N2 + m_CO2*cp_CO2)/(m_N2+m_CO2);
+disp("cp=")
+disp(cp)
+disp("kJ/kg K")
+
+cv=(m_N2*cv_N2 + m_CO2*cv_CO2)/(m_N2+m_CO2);
+disp("cv=")
+disp(cv)
+disp("kJ/kg K")
+
+
+
+T1=293; //K
+T2=323; //K
+dU=m*cv*(T2-T1);
+disp("Change in internal energy =")
+disp(dU)
+disp("kJ")
+
+dH=m*cp*(T2-T1);
+disp("Change in enthalpy =")
+disp(dH)
+disp("kJ")
+
+dS=m*cv*log(T2/T1); //Constant volume process
+disp("Change in entropy=")
+disp(dS)
+disp("kJ/kg K")
+
+
+disp("When the mixture is heated at constant pressure")
+
+disp("If the mixture is heated at constant pressure ΔU and ΔH will remain the same")
+
+dS=m*cp*log(T2/T1);
+disp("Change in entropy =")
+disp(dS)
+disp("kJ/kg K") \ No newline at end of file
diff --git a/551/CH9/EX9.11/11.sce b/551/CH9/EX9.11/11.sce
new file mode 100755
index 000000000..4e7fac62d
--- /dev/null
+++ b/551/CH9/EX9.11/11.sce
@@ -0,0 +1,59 @@
+clc
+
+Cv_O2=21.07; //kJ/mole K
+Cv_CO=20.86; //kJ/mole K
+
+p_O2=8*10^5; //Pa
+p_CO=1*10^5; //Pa
+
+V_O2=1.8; //m^3
+V_CO=3.6; //m^3
+
+T_O2=323; //K
+T_CO=293; //K
+
+R0=8314;
+
+n_O2=p_O2*V_O2/R0/T_O2;
+n_CO=p_CO*V_CO/R0/T_CO;
+
+n=(n_O2+n_CO);
+V=(V_O2+V_CO);
+
+disp("(i) Final temperature (T) and pressure (p) of the mixture")
+
+//Before mixing
+U1=n_O2*Cv_O2*T_O2 + n_CO*Cv_CO*T_CO;
+
+//After mixing
+//U2=T*(n_O2*Cv_O2 + n_CO*Cv_CO);
+//U1=U2
+
+T=U1/(n_O2*Cv_O2 + n_CO*Cv_CO);
+t=T-273;
+disp("Final temperature =")
+disp(t)
+disp("°C")
+
+p=n*R0*T/V/10^5;
+disp("Final pressure =")
+disp(p)
+disp("bar")
+
+
+disp("(ii) Change of entropy")
+//For oxygen
+dS_O1A=n_O2*R0*log(V/V_O2); //isothermal process
+dS_O2A=n_O2*Cv_O2*log(T_O2/T); //constant volume process
+
+dS_O12=dS_O1A - dS_O2A; // Change of entropy of O2
+
+//For CO
+
+dS_CO12=n_CO*R0*log(V/V_CO) + n_CO*Cv_CO*log(T/T_CO); //Change of entropy of CO
+
+
+dS=(dS_O12 + dS_CO12)/10^3;
+disp("Change of entropy of system =")
+disp(dS)
+disp("kJ/K") \ No newline at end of file
diff --git a/551/CH9/EX9.12/12.sce b/551/CH9/EX9.12/12.sce
new file mode 100755
index 000000000..02d2034fc
--- /dev/null
+++ b/551/CH9/EX9.12/12.sce
@@ -0,0 +1,63 @@
+clc
+
+p_A=16*10^5; //Pa
+p_B=6.4*10^5; //Pa
+
+T_A=328; //K
+T_B=298; //K
+
+n_A=0.6; //kg-mole
+m_B=3; //kg
+
+R0=8314;
+M_A=28;
+y=1.4;
+
+V_A=n_A*R0*T_A/p_A;
+
+m_A=n_A*M_A;
+R=R0/M_A;
+
+V_B=m_B*R*T_B/p_B;
+
+V=V_A+V_B;
+m=m_A+m_B;
+T=303; //K
+
+
+disp("(a) (i) Final equilibrium pressure, p")
+p=m*R*T/V/10^5;
+disp("p=")
+disp(p)
+disp("bar")
+
+cv=R/10^3/(y-1);
+
+disp("(ii) Amount of heat transferred, Q :")
+
+U1=cv*(m_A*T_A + m_B*T_B);
+U2=m*cv*T;
+
+Q=U2-U1;
+disp("Q=")
+disp(Q)
+disp("kJ")
+
+
+disp("(b) If the vessel were insulated :")
+
+disp("(i) Final temperature,")
+
+T=cv*(m_A*T_A + m_B*T_B)/(m*cv);
+t=T-273;
+disp("T=")
+disp(t)
+disp("°C")
+
+
+disp("(ii) Final pressure")
+
+p=m*R*T/V/10^5;
+disp("p=")
+disp(p)
+disp("bar") \ No newline at end of file
diff --git a/551/CH9/EX9.13/13.sce b/551/CH9/EX9.13/13.sce
new file mode 100755
index 000000000..18f783a4a
--- /dev/null
+++ b/551/CH9/EX9.13/13.sce
@@ -0,0 +1,20 @@
+clc
+
+m_O2=3; //kg
+M_O2=32;
+
+m_N2=9; //kg
+M_N2=28;
+
+R0=8.314;
+
+R_O2=R0/M_O2;
+R_N2=R0/M_N2;
+
+x_O2=(m_O2/M_O2)/((m_O2/M_O2) + (m_N2/M_N2));
+x_N2=(m_N2/M_N2)/((m_O2/M_O2) + (m_N2/M_N2));
+
+dS=-m_O2*R_O2*log(x_O2) -m_N2*R_N2*log(x_N2);
+disp("Change in entropy =")
+disp(dS)
+disp("kJ/kg K") \ No newline at end of file
diff --git a/551/CH9/EX9.14/14.sce b/551/CH9/EX9.14/14.sce
new file mode 100755
index 000000000..904df2b43
--- /dev/null
+++ b/551/CH9/EX9.14/14.sce
@@ -0,0 +1,17 @@
+clc
+m_N2=2.5; //kg
+M_N2=28;
+
+p_N2=15; //bar
+p_total=20; //bar
+
+n_N2=m_N2/M_N2;
+p_O2=p_total-p_N2;
+
+n_O2=p_O2/p_N2*n_N2;
+M_O2=32;
+
+m_O2=n_O2*M_O2;
+disp("Mass of O2 added =")
+disp(m_O2)
+disp("kg") \ No newline at end of file
diff --git a/551/CH9/EX9.15/15.sce b/551/CH9/EX9.15/15.sce
new file mode 100755
index 000000000..20845773d
--- /dev/null
+++ b/551/CH9/EX9.15/15.sce
@@ -0,0 +1,34 @@
+clc
+n_O2=1;
+
+//V_O2=0.21*V;
+//V_N2=0.79*V;
+M_N2=28;
+M_O2=32;
+
+disp("(i) Moles of N2 per mole of O2 :")
+n_N2=n_O2*0.79/0.21;
+disp("n_N2=")
+disp(n_N2)
+disp("moles")
+
+n=n_O2+n_N2;
+disp("(ii) p_O2 and p_N2 :")
+p=1; //atm
+
+p_O2=n_O2/n*p;
+disp("p_O2=")
+disp(p_O2)
+disp("atm")
+
+p_N2=n_N2/n*p;
+disp("p_N2=")
+disp(p_N2);
+disp("atm")
+
+
+disp("(iii) The kg of nitrogen per kg of mixture :")
+x=n_N2*M_N2/(n_N2*M_N2+n_O2*M_O2);
+disp(" The kg of nitrogen per kg of mixture =")
+disp(x)
+disp("kg N2/kg mix") \ No newline at end of file
diff --git a/551/CH9/EX9.16/16.sce b/551/CH9/EX9.16/16.sce
new file mode 100755
index 000000000..54d5d5b06
--- /dev/null
+++ b/551/CH9/EX9.16/16.sce
@@ -0,0 +1,48 @@
+clc
+V=0.6; //m^3
+p1=12*10^5; //Pa
+p2=18*10^5; //Pa
+T=298; //K
+R0=8.314;
+x_O2=0.23;
+x_N2=0.77;
+
+n=p1*V/R0/10^3/T;
+
+//Considering 100 kg of air
+m_O2=23; //kg
+m_N2=77; //kg
+M_O2=32;
+M_N2=28;
+m=100; //kg
+
+R=(m_O2/M_O2 + m_N2/M_N2)*R0/m; //for air
+M=R0/R; //for air
+
+m=p1*V/R/T/10^3;
+
+m_O2=x_O2*m;
+disp("Mass of O2=")
+disp(m_O2)
+disp("kg")
+
+m_N2=x_N2*m;
+disp("Mass of N2=")
+disp(m_N2)
+disp("kg")
+
+
+//After adding CO2 in the vessel
+p2=18*10^5; //Pa;
+
+// p_CO2+p_N2+p_O2=18*10^5
+// p_N2 + p_O2=12*10^5
+
+p_CO2=6*10^5; //Pa
+M_CO2=44;
+R_CO2=R0/M_CO2;
+
+m_CO2=p_CO2*V/(R_CO2*10^3*T);
+disp("Mass of CO2 = ")
+disp(m_CO2)
+disp("kg")
diff --git a/551/CH9/EX9.17/17.sce b/551/CH9/EX9.17/17.sce
new file mode 100755
index 000000000..5e53f2cb2
--- /dev/null
+++ b/551/CH9/EX9.17/17.sce
@@ -0,0 +1,36 @@
+clc
+V=6; //m^3
+%A=0.45;
+%B=0.55;
+R_A=0.288; //kJ/kg K
+R_B=0.295; //kJ/kg K
+m=2; //kg
+T=303; //K
+
+disp("(i) The partial pressures")
+m_A=%A*m;
+m_B=%B*m;
+
+p_A=m_A*R_A*10^3*T/V/10^5; //bar
+disp("p_A=")
+disp(p_A)
+disp("bar")
+
+p_B=m_B*R_B*10^3*T/V/10^5; //bar
+disp("p_B=")
+disp(p_B)
+disp("bar")
+
+
+disp("(ii) The total pressure")
+p=p_A+p_B;
+disp("p=")
+disp(p)
+disp("bar")
+
+
+disp("(iii) The mean value of R for the mixture")
+Rm=(m_A*R_A + m_B*R_B)/(m_A + m_B);
+disp("Rm=")
+disp(Rm)
+disp("kJ/kg K") \ No newline at end of file
diff --git a/551/CH9/EX9.18/18.sce b/551/CH9/EX9.18/18.sce
new file mode 100755
index 000000000..a6bff6d7f
--- /dev/null
+++ b/551/CH9/EX9.18/18.sce
@@ -0,0 +1,69 @@
+clc
+m_O2=4; //kg
+m_N2=6; //kg
+p=4*10^5; //Pa
+T=300; //K
+M_O2=32;
+M_N2=28;
+m=10; //kg
+
+disp("(i) The mole fraction of each component")
+n_O2=m_O2/M_O2;
+n_N2=m_N2/M_N2;
+
+x_O2=n_O2/(n_O2+n_N2);
+disp("x_O2=")
+disp(x_O2)
+
+x_N2=n_N2/(n_N2+n_O2);
+disp("x_N2=")
+disp(x_N2)
+
+
+disp("(ii) The average molecular weight")
+M=(n_O2*M_O2 + n_N2*M_N2)/(n_O2 + n_N2);
+disp("M=")
+disp(M)
+
+
+disp("(iii) The specific gas constant")
+R0=8.314;
+R=R0/M;
+disp("R=")
+disp(R)
+disp("kJ/kg K")
+
+disp("(iv) The volume and density")
+
+V=m*R*T*10^3/p;
+disp("V=")
+disp(V)
+disp("m^3")
+
+rho=(m_O2/V) + (m_N2/V);
+disp("density=")
+disp(rho)
+disp("kg/m^3")
+
+
+disp("(v) The partial pressures and partial volumes")
+
+p_O2=n_O2*R0*10^3*T/V/10^5; //bar
+disp("p_O2=")
+disp(p_O2)
+disp("bar")
+
+p_N2=n_N2*R0*10^3*T/V/10^5; //bar
+disp("p_N2=")
+disp(p_N2)
+disp("bar")
+
+V_O2=x_O2*V;
+disp("V_O2=")
+disp(V_O2)
+disp("m^3")
+
+V_N2=x_N2*V;
+disp("V_N2=")
+disp(V_N2)
+disp("m^3")
diff --git a/551/CH9/EX9.19/19.sce b/551/CH9/EX9.19/19.sce
new file mode 100755
index 000000000..8497305c2
--- /dev/null
+++ b/551/CH9/EX9.19/19.sce
@@ -0,0 +1,39 @@
+clc
+cp_CO2=0.85; //kJ/kg K
+cp_N2=1.04; //kJ/kg K
+m_CO2=4; //kg
+T1_CO2=313; //K
+m_N2=8; //kg
+T1_N2=433; //K
+p2=0.7; //bar
+p1_CO2=1.4; //bar
+p1_N2=1;
+R=8.314;
+M_CO2=44;
+M_N2=28;
+R_CO2=R/M_CO2;
+R_N2=R/M_N2;
+
+disp("(i) Final temperature, T2")
+T2=(m_CO2*cp_CO2*T1_CO2 + m_N2*cp_N2*T1_N2)/(m_CO2*cp_CO2 + m_N2*cp_N2);
+disp("T2=")
+disp(T2)
+disp("K")
+
+
+disp("(ii) Change in entropy")
+n_CO2=0.0909;
+n_N2=0.2857;
+n=n_CO2 + n_N2;
+
+x_CO2=n_CO2/n;
+x_N2=n_N2/n;
+
+p2_CO2=x_CO2*p2;
+p2_N2=x_N2*p2;
+
+
+dS=m_CO2*cp_CO2*log(T2/T1_CO2) - m_CO2*R_CO2*log(p2_CO2/p1_CO2) + m_N2*cp_N2*log(T2/T1_N2) - m_N2*R_N2*log(p2_N2/p1_N2);
+disp("dS=")
+disp(dS)
+disp("kJ/K") \ No newline at end of file
diff --git a/551/CH9/EX9.2/2.sce b/551/CH9/EX9.2/2.sce
new file mode 100755
index 000000000..c39d55eaf
--- /dev/null
+++ b/551/CH9/EX9.2/2.sce
@@ -0,0 +1,28 @@
+clc
+R0=8.314;
+
+M_O2=32;
+M_N2=28;
+M_Ar=40;
+M_CO2=44;
+
+R_O2=R0/M_O2; //kJ/kg K
+R_N2=R0/M_N2; //kJ/kg K
+R_Ar=R0/M_Ar; //kJ/kg K
+R_CO2=R0/M_CO2; //kJ/kg K
+
+%O2=0.2314;
+%N2=0.7553;
+%Ar=0.0128;
+%CO2=0.0005;
+
+disp("(i) Gas constant for air")
+R=%O2*R_O2 + %N2*R_N2 + %Ar*R_Ar + %CO2*R_CO2;
+disp("R=")
+disp(R)
+disp("kJ/kg K")
+
+disp("(ii) Apparent molecular weight.")
+M=R0/R;
+disp("M=")
+disp(M) \ No newline at end of file
diff --git a/551/CH9/EX9.20/20.sce b/551/CH9/EX9.20/20.sce
new file mode 100755
index 000000000..b5c900bd5
--- /dev/null
+++ b/551/CH9/EX9.20/20.sce
@@ -0,0 +1,26 @@
+clc
+
+cv_O2=0.39; //kJ/kg K
+cv_N2=0.446; //kJ/kg K
+n_O2=1;
+n_N2=2;
+M_O2=32;
+M_N2=28;
+m_O2=32; //kg
+m_N2=2*28; //kg
+T_O2=293; //K
+T_N2=301; //K
+R0=8.314;
+p_O2=2.5*10^5; //Pa
+p_N2=1.5*10^5; //Pa
+
+T2=(m_O2*cv_O2*T_O2 + m_N2*cv_N2*T_N2)/(m_O2*cv_O2 + m_N2*cv_N2);
+
+V_O2=n_O2*R0*10^5*T_O2/p_O2;
+V_N2=n_N2*R0*10^5*T_N2/p_N2;
+V=V_O2+V_N2;
+
+dS=m_O2*[cv_O2*log(T2/T_O2) + R0/M_O2*log(V/V_O2)] + m_N2*[cv_N2*log(T2/T_N2) + R0/M_N2*log(V/V_N2)];
+disp("dS=")
+disp(dS)
+disp("kJ") \ No newline at end of file
diff --git a/551/CH9/EX9.21/21.sce b/551/CH9/EX9.21/21.sce
new file mode 100755
index 000000000..8550d57c6
--- /dev/null
+++ b/551/CH9/EX9.21/21.sce
@@ -0,0 +1,56 @@
+clc
+cv_N2=0.744; //kJ/kg K
+cv_H2=10.352; //kJ/kg K
+cp_N2=1.041; //kJ/kg K
+cp_H2=14.476; //kJ/kg K
+
+V=0.45; //m^3
+V_H2=0.3; //m^3
+V_N2=0.15; //m^3
+
+p_H2=3*10^5; //Pa
+p_N2=6*10^5; //Pa
+
+T_H2=403; //K
+T_N2=303; //K
+
+R_H2=8.314/2;
+R_N2=8.314/28;
+
+disp("(i) Temperature of equilibrium mixture")
+
+m_H2=p_H2*V_H2/(R_H2*10^3)/T_H2;
+m_N2=p_N2*V_N2/(R_N2*10^3)/T_N2;
+
+T2=(m_H2*cv_H2*T_H2 + m_N2*cv_N2*T_N2)/(m_H2*cv_H2 + m_N2*cv_N2);
+disp("T2=")
+disp(T2)
+disp("K")
+
+
+disp("(ii) Pressure of the mixture")
+p2_H2=m_H2*R_H2*10^3*T2/V;
+p2_N2=m_N2*R_N2*10^3*T2/V;
+
+p2=p2_H2+p2_N2;
+disp("p2=")
+disp(p2/10^5)
+disp("bar")
+
+disp("(iii) Change in entropy :")
+
+dS_H2=m_H2*[cp_H2*log(T2/T_H2) - R_H2*log(p2_H2/p_H2)];
+disp("Change in entropy of H2 =")
+disp(dS_H2)
+disp("kJ/K")
+
+dS_N2=m_N2*[cp_N2*log(T2/T_N2) - R_N2*log(p2_N2/p_N2)];
+disp("Change in entropy of N2 =")
+disp(dS_N2)
+disp("kJ/K")
+
+dS=dS_H2+dS_N2;
+
+disp("Total change in entropy =")
+disp(dS)
+disp("kJ/K") \ No newline at end of file
diff --git a/551/CH9/EX9.22/22.sce b/551/CH9/EX9.22/22.sce
new file mode 100755
index 000000000..8887a3044
--- /dev/null
+++ b/551/CH9/EX9.22/22.sce
@@ -0,0 +1,40 @@
+clc
+
+cv_N2=0.745; //kJ/kg K
+cv_CO2=0.653; //kJ/kg K
+
+cp_N2=1.041; //kJ/kg K
+cp_CO2=0.842; //kJ/kg K
+
+m_N2=4; //kg
+m_CO2=6; //kg
+pmix=4; //bar
+m=m_N2+m_CO2;
+
+T1=298; //K
+T2=323; //K
+
+cv_mix=(m_N2*cv_N2 + m_CO2*cv_CO2)/(m_N2+m_CO2);
+disp("cv_mix=")
+disp(cv_mix)
+disp("kJ/kg K")
+
+cp_mix=(m_N2*cp_N2 + m_CO2*cp_CO2)/(m_N2+m_CO2);
+disp("cp_mix=")
+disp(cp_mix)
+disp("kJ/kg K")
+
+dU=m*cv_mix*(T2-T1);
+disp("Change in internal energy=")
+disp(dU)
+disp("kJ")
+
+dH=m*cp_mix*(T2-T1);
+disp("Change in enthalpy=")
+disp(dH)
+disp("kJ")
+
+dS=m_N2*cv_N2*log(T2/T1) + m_CO2*cv_CO2*log(T2/T1);
+disp("Change in entropy=")
+disp(dS)
+disp("kJ/K") \ No newline at end of file
diff --git a/551/CH9/EX9.3/3.sce b/551/CH9/EX9.3/3.sce
new file mode 100755
index 000000000..314e74cd9
--- /dev/null
+++ b/551/CH9/EX9.3/3.sce
@@ -0,0 +1,69 @@
+clc
+p=1; //bar
+
+//For oxygen
+m_O2=0.2314;
+M_O2=32;
+n_O2=m_O2/M_O2;
+
+//For Nitrogen
+m_N2=0.7553;
+M_N2=28;
+n_N2=m_N2/M_N2;
+
+//For Argon
+m_Ar=0.0128;
+M_Ar=40;
+n_Ar=m_Ar/M_Ar;
+
+//For CO2
+m_CO2=0.0005;
+M_CO2=44;
+n_CO2=m_CO2/M_CO2;
+
+
+n=n_O2 + n_N2 + n_Ar + n_CO2;
+
+//Let Vi/V be A
+
+A_O2=n_O2/n * 100;
+disp("Vi/V of O2=")
+disp(A_O2)
+disp("%")
+
+A_N2=n_N2/n * 100;
+disp("Vi/V of N2=")
+disp(A_N2)
+disp("%")
+
+A_Ar=n_Ar/n *100;
+disp("Vi/V of Ar")
+disp(A_Ar)
+disp("%")
+
+A_CO2=n_CO2/n * 100;
+disp("Vi/V of CO2=")
+disp(A_CO2)
+disp("%")
+
+
+P_O2=n_O2/n*p;
+disp("Partial pressure of O2=")
+disp(P_O2)
+disp("bar")
+
+P_N2=n_N2/n*p;
+disp("Partial pressure of N2=")
+disp(P_N2)
+disp("bar")
+
+P_Ar=n_Ar/n*p;
+disp("Partial pressure of Ar=")
+disp(P_Ar)
+disp("bar")
+
+P_CO2=n_CO2/n*p;
+disp("Partial pressure of CO2=")
+disp(P_CO2)
+disp("bar")
+
diff --git a/551/CH9/EX9.4/4.sce b/551/CH9/EX9.4/4.sce
new file mode 100755
index 000000000..81928213a
--- /dev/null
+++ b/551/CH9/EX9.4/4.sce
@@ -0,0 +1,67 @@
+clc
+p=1*10^5; //Pa
+T=293; //K
+n_CO2=1; //moles of CO2
+n=4; //moles of air
+M_CO2=44;
+M_N2=28;
+M_O2=32;
+
+//Let A be the volumeetric analysis
+A_O2=0.21;
+A_N2=0.79;
+
+n_O2=A_O2*n;
+n_N2=A_N2*n;
+
+disp("(i) The masses of CO2, O2 and N2, and the total mass")
+
+m_CO2=n_CO2*M_CO2;
+disp("Mass of CO2=")
+disp(m_CO2)
+disp("kg")
+
+m_O2=n_O2*M_O2;
+disp("Mass of O2=")
+disp(m_O2)
+disp("kg")
+
+m_N2=n_N2*M_N2;
+disp("Mass of N2=")
+disp(m_N2)
+disp("kg")
+
+m=m_CO2 + m_O2 + m_N2;
+disp("Total mass =")
+disp(m)
+disp("kg")
+
+
+disp("(ii) The percentage carbon content by mass")
+//Since the molecular weight of carbon is 12, therefore, there are 12 kg of carbon present for every mole of CO2
+m_C=12; //kg
+
+%C=m_C/m*100;
+disp("Percentage carbon in mixture")
+disp(%C)
+disp("%")
+
+
+disp("(iii) The apparent molecular weight and the gas constant for the mixture")
+n=n_CO2 + n_O2 + n_N2;
+M=n_CO2/n*M_CO2 + n_O2/n*M_O2 + n_N2/n*M_N2;
+disp("Apparent Molecular weight")
+disp(M)
+
+R0=8.314;
+R=R0/M;
+disp("Gas constant for the mixture=")
+disp(R)
+disp("kJ/kg K")
+
+
+disp("(iv) The specific volume of the mixture")
+v=R*10^3*T/p;
+disp("specific volume=")
+disp(v)
+disp("m^3/kg")
diff --git a/551/CH9/EX9.5/5.sce b/551/CH9/EX9.5/5.sce
new file mode 100755
index 000000000..350304629
--- /dev/null
+++ b/551/CH9/EX9.5/5.sce
@@ -0,0 +1,29 @@
+clc
+p=1*10^5; //Pa
+T=298; //K
+M_H2=2;
+M_O2=32;
+R0=8314;
+
+// ratio = V_H2/V_O2=2;
+ratio=2;
+
+disp("(i) The mass of O2 required")
+//Let the mass of O2 per kg of H2 = x kg
+m_H2=1; //kg
+n_H2=m_H2/M_H2;
+
+// n_O2=x/M_O2
+x=M_O2*n_H2/ratio;
+disp("Mass of O2 per kg of H2=")
+disp(x)
+disp("kg")
+
+
+disp("(ii) The volume of the container")
+n_O2=x/M_O2;
+n=n_H2 + n_O2;
+V=n*R0*T/p;
+disp("V=")
+disp(V)
+disp("m^3") \ No newline at end of file
diff --git a/551/CH9/EX9.6/6.sce b/551/CH9/EX9.6/6.sce
new file mode 100755
index 000000000..cdb8a3ac8
--- /dev/null
+++ b/551/CH9/EX9.6/6.sce
@@ -0,0 +1,27 @@
+clc
+
+//Let composition of mixture by volume be denoted by c1
+//Let Final composition desired be denoted by c2
+
+c1_H2=0.78;
+c1_CO=0.22;
+
+c2_H2=0.52;
+c2_CO=0.48;
+
+M_H2=2;
+M_CO=28;
+
+M=c1_H2*M_H2 + c1_CO*M_CO;
+
+// Let x kg of mixture be removed and y kg of CO be added.
+
+x=(c1_H2 - c2_H2)/c1_H2*M;
+disp("Mass of mixture removed =")
+disp(x)
+disp("kg")
+
+y=M_CO/M*x;
+disp("Mass of CO added=")
+disp(y)
+disp("kg") \ No newline at end of file
diff --git a/551/CH9/EX9.7/7.sce b/551/CH9/EX9.7/7.sce
new file mode 100755
index 000000000..bf9d035b8
--- /dev/null
+++ b/551/CH9/EX9.7/7.sce
@@ -0,0 +1,61 @@
+clc
+
+ratio=1/8; //volume ratio; v1/v2
+T1=1223; //K
+
+cp_CO2=1.235; //kJ/kg K
+cp_O2=1.088; //kJ/kg K
+cp_N2=1.172; //kJ/kg K
+
+n_CO2=0.13;
+n_O2=0.125;
+n_N2=0.745;
+
+M_CO2=44;
+M_O2=32;
+M_N2=28;
+
+m_CO2=M_CO2*n_CO2;
+m_O2=M_O2*n_O2;
+m_N2=M_N2*n_N2;
+
+m=m_CO2 + m_O2 + m_N2;
+
+// Let Fraction by mass be denoted by F
+F_CO2=m_CO2/m;
+F_O2=m_O2/m;
+F_N2=m_N2/m;
+
+
+cp=F_CO2*cp_CO2 + F_O2*cp_O2 + F_N2*cp_N2;
+
+R0=8.314;
+R=F_CO2*R0/M_CO2 + F_O2*R0/M_O2 + F_N2*R0/M_N2;
+
+cv=cp - R;
+n=1.2;
+
+disp("(i) The workdone")
+T2=T1*(ratio)^(n-1);
+W=R*(T1-T2)/(n-1);
+disp("W=")
+disp(W)
+disp("kJ/kg")
+
+
+disp("(ii) The heat flow")
+du=cv*(T2-T1);
+Q=du + W;
+disp("Q=")
+disp(Q)
+disp("kJ/kg")
+
+
+disp("(iii) Change of entropy per kg of mixture")
+ds_1A=R*log(1/ratio); //isothermal process
+ds_2A=cv*log(T1/T2);
+
+ds_12=ds_1A - ds_2A;
+disp("change of entropy=")
+disp(ds_12)
+disp("kJ/kg K") \ No newline at end of file
diff --git a/551/CH9/EX9.8/8.sce b/551/CH9/EX9.8/8.sce
new file mode 100755
index 000000000..97140a003
--- /dev/null
+++ b/551/CH9/EX9.8/8.sce
@@ -0,0 +1,44 @@
+clc
+
+M_CO2=44;
+M_H2=2;
+M_N2=28;
+M_CH4=16;
+M_CO=28;
+
+// Let volumetric analysis be denoted by V
+V_CO=0.28;
+V_H2=0.13;
+V_CH4=0.04;
+V_CO2=0.04;
+V_N2=0.51;
+
+Cp_CO=29.27; //kJ/mole K
+Cp_H2=28.89; //kJ/mole K
+Cp_CH4=35.8; //kJ/mole K
+Cp_CO2=37.22; //kJ/mole K
+Cp_N2=29.14; //kJ/mole K
+
+R0=8.314;
+
+Cp=V_CO*Cp_CO + V_H2*Cp_H2 + V_CO2*Cp_CO2 + V_CH4*Cp_CH4 + V_N2*Cp_N2;
+disp("Cp=")
+disp(Cp)
+disp("kJ/mole K")
+
+Cv=Cp-R0;
+disp("Cv=")
+disp(Cv)
+disp("kJ/mole K")
+
+M=V_CO*M_CO + V_H2*M_H2 + V_CO2*M_CO2 + V_CH4*M_CH4 + V_N2*M_N2;
+
+cp=Cp/M;
+disp("cp=")
+disp(cp)
+disp("kJ/kg K")
+
+cv=Cv/M;
+disp("cv")
+disp(cv)
+disp("kJ/kg K") \ No newline at end of file
diff --git a/551/CH9/EX9.9/9.sce b/551/CH9/EX9.9/9.sce
new file mode 100755
index 000000000..25948018c
--- /dev/null
+++ b/551/CH9/EX9.9/9.sce
@@ -0,0 +1,58 @@
+clc
+
+p=1.3; //bar
+R0=8.314;
+
+M_CO2=44;
+M_O2=32;
+M_N2=28;
+M_CO=28;
+
+m_O2=0.1;
+m_N2=0.7;
+m_CO2=0.15;
+m_CO=0.05;
+
+//Considering 1 kg of mixture
+m=1; //kg
+
+//let moles be denoted by n
+n_O2=m_O2/M_O2;
+n_N2=m_N2/M_N2;
+n_CO2=m_CO2/M_CO2;
+n_CO=m_CO/M_CO;
+
+M=1/(m_O2/M_O2 + m_N2/M_N2 + m_CO2/M_CO2 + m_CO/M_CO);
+
+n=m/M;
+
+x_O2=n_O2/n;
+x_N2=n_N2/n;
+x_CO2=n_CO2/n;
+x_CO=n_CO/n;
+
+disp("(i) Partial pressures of the constituents")
+P_O2=x_O2*p;
+disp("Partial pressure of O2=")
+disp(P_O2)
+disp("bar")
+
+P_N2=x_N2*p;
+disp("Partial pressure of N2=")
+disp(P_N2)
+disp("bar")
+
+P_CO2=x_CO2*p;
+disp("Partial pressure of CO2=")
+disp(P_CO2)
+disp("bar")
+
+P_CO=x_CO*p;
+disp("Partial pressure of CO=")
+disp(P_CO)
+disp("bar")
+
+disp("Gas constant of mixture =")
+R_mix=R0/M;
+disp(R_mix)
+disp("kJ/kg K") \ No newline at end of file
generated by cgit (git 2.34.1) at 2024-12-20 08:08:19 +0000