diff options
Diffstat (limited to '2705/CH12/EX12.2/Ex12_2.sce')
-rwxr-xr-x | 2705/CH12/EX12.2/Ex12_2.sce | 33 |
1 files changed, 33 insertions, 0 deletions
diff --git a/2705/CH12/EX12.2/Ex12_2.sce b/2705/CH12/EX12.2/Ex12_2.sce new file mode 100755 index 000000000..fd3cd191b --- /dev/null +++ b/2705/CH12/EX12.2/Ex12_2.sce @@ -0,0 +1,33 @@ +clear;
+clc;
+disp('Example 12.2');
+
+// aim : To determine the increases in pressure, temperature and internal energy per kg of air
+
+// Given values
+T1 = 273;// [K]
+P1 = 140;// [kN/m^2]
+C1 = 900;// [m/s]
+C2 = 300;// [m/s]
+cp = 1.006;// [kJ/kg K]
+cv =.717;// [kJ/kg K]
+
+// solution
+R = cp-cv;// [kJ/kg K]
+Gamma = cp/cv;// heat capacity ratio
+// for frictionless adiabatic flow, (C2^2-C1^2)/2=Gamma/(Gamma-1)*R*(T1-T2)
+
+T2 =T1-((C2^2-C1^2)*(Gamma-1)/(2*Gamma*R))*10^-3; // [K]
+T_inc = T2-T1;// increase in temperature [K]
+
+P2 = P1*(T2/T1)^(Gamma/(Gamma-1));// [MN/m^2]
+P_inc = (P2-P1)*10^-3;// increase in pressure,[MN/m^2]
+
+U_inc = cv*(T2-T1);// Increase in internal energy per kg,[kJ/kg]
+mprintf('\n The increase in pressure is = %f MN/m^2\n',P_inc);
+mprintf('\n Increase in temperature is = %f K\n',T_inc);
+mprintf('\n Increase in internal energy is = %f kJ/kg\n',U_inc);
+
+// there is minor variation in result
+
+// End
|