4 files changed, 68 insertions, 0 deletions
diff --git a/551/CH7/EX7.17/17.sce b/551/CH7/EX7.17/17.sce
new file mode 100755
index 000000000..f5e80e1ab
--- /dev/null
+++ b/551/CH7/EX7.17/17.sce
@@ -0,0 +1,31 @@
+clc
+B=5*10^(-5); // /K
+K=8.6*10^(-12); // m^2/N
+v=0.114*10^(-3); //m^3/kg
+p2=800*10^5; //Pa
+p1=20*10^5; //Pa
+T=288; //K
+disp("(i) Work done on the copper = ")
+W=-v*K/2*(p2^2-p1^2);
+disp(W)
+disp("J/kg")
+
+disp("(ii) Change in entropy =")
+ds=-v*B*(p2-p1);
+disp(ds)
+disp("J/kg K")
+
+disp("(iii) The heat transfer =")
+Q=T*ds;
+disp(Q)
+disp("J/kg")
+
+disp("(iv) Change in internal energy =")
+du=Q-W;
+disp(du)
+disp("J/kg")
+
+disp("(v) cp – cv =")
+R=B^2*T*v/K;
+disp(R)
+disp("J/kg K")
+\ No newline at end of file
diff --git a/551/CH7/EX7.18/18.sce b/551/CH7/EX7.18/18.sce
new file mode 100755
index 000000000..b0a666801
--- /dev/null
+++ b/551/CH7/EX7.18/18.sce
@@ -0,0 +1,10 @@
+clc
+vg=0.1274; //m^3/kg
+vf=0.001157; //m^3/kg
+// dp/dT=32; //kPa/K
+T3=473; //K
+
+h_fg=32*10^3*T3*(vg-vf)/10^3;
+disp("h_fg=")
+disp(h_fg)
+disp("kJ/kg")
+\ No newline at end of file
diff --git a/551/CH7/EX7.19/19.sce b/551/CH7/EX7.19/19.sce
new file mode 100755
index 000000000..341aee2f9
--- /dev/null
+++ b/551/CH7/EX7.19/19.sce
@@ -0,0 +1,12 @@
+clc
+h_fg=334; //kJ/kg
+v_liq=1; //m^3/kg
+v_ice=1.01; //m^3/kg
+T1=273; //K
+T2=263; //K
+p1=1.013*10^5; //Pa
+
+p2=(p1+h_fg*10^3/(v_ice-v_liq)*log(T1/T2))/10^5;
+disp("p2=")
+disp(p2)
+disp("bar")
+\ No newline at end of file
diff --git a/551/CH7/EX7.20/20.sce b/551/CH7/EX7.20/20.sce
new file mode 100755
index 000000000..d997b33c0
--- /dev/null
+++ b/551/CH7/EX7.20/20.sce
@@ -0,0 +1,15 @@
+clc
+h_fg=294.54; //kJ/kg
+
+// log(p)=7.0323 - 3276.6/T - 0.652*log(T)
+// Differentiating both sides, we get
+// 1/2.302/p*dp/dT=3276.6/T^2-0.652/2.302/T
+
+//Putting p=0.1 bar , we get
+p=0.1; //bar
+T=523; //K
+
+vg=h_fg*10^3/T/(2.302*3276.6*p*10^5/T^2 - 0.652*p*10^5/T);
+disp("vg=")
+disp(vg)
+disp("m^3/kg")
+\ No newline at end of file