1 files changed, 33 insertions, 0 deletions
diff --git a/2144/CH3/EX3.27/ex3_27.sce b/2144/CH3/EX3.27/ex3_27.sce
new file mode 100755
index 000000000..4111273d6
--- /dev/null
+++ b/2144/CH3/EX3.27/ex3_27.sce
@@ -0,0 +1,33 @@
+// Exa 3.27
+clc;
+clear;
+close;
+// Given data
+P1 = 1 * 10^5;// in N/m^2
+V1 = 0.1;// in m^3
+V2 = 0.01;// in m^3
+T1 = 90;// in degree C
+T1 = T1  +273;// in K
+R = 0.287;// in kJ/kg-K
+R = R  *10^3;
+C_v = 0.717;// in kJ/kg-K
+C_P = 1.005;// in kJ/kg-K
+m = (P1 * V1)/(R*T1);// in kg
+Gamma = 1.4;
+T2 = T1 * ((V1/V2)^(Gamma - 1));// in K
+del_U = m*C_v*(T1-T2);// in kJ
+disp(del_U,"The change in internal energy in kJ is : ")
+del_E = m * C_P*(T2-T1);// in kJ
+disp(del_E,"The change in enthalpy in kJ is : ")
+U2 = m*C_v*T2;//Internal energy at 2 in kJ
+T= 473;// temp. of entering air
+E = V1*C_P*T;//Enthalpy of entering air in kJ
+// U3= (m+V1)*C_v*T3 ; (internal energy at 3)
+// U3= U2+E
+T3 = (E+U2)/( (m+V1)*C_v );// in K
+disp(T3,"Temperature in K is");
+m=m+.1;
+P3 =m* R*T3/V2;// in N/m^2
+disp(P3*10^-6,"The pressure in MN/m^2 is");
+
+// Note: There is a calculation error to evaluating the value of P3. So the answer in the book of P3 is wrong.