diff options
Diffstat (limited to '3843/CH9/EX9.6')
| -rw-r--r-- | 3843/CH9/EX9.6/Ex9_6.sce | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/3843/CH9/EX9.6/Ex9_6.sce b/3843/CH9/EX9.6/Ex9_6.sce new file mode 100644 index 000000000..100fe7f7b --- /dev/null +++ b/3843/CH9/EX9.6/Ex9_6.sce @@ -0,0 +1,23 @@ +// Example 9_6
+clc;funcprot(0);
+// Given data
+r=16;// The compression ratio
+T_1=200+273;// K
+P_1=200;// kPa
+r_c=2;// The cut off ratio
+r_p=1.3;// The pressure ratio
+c_p=1.00;// kJ/kg.K
+c_v=0.717;// kJ/kg.K
+R=0.287;// kJ/kg.K
+k=1.4;// The specific heat ratio
+
+// Calculation
+n=1-((1/(r^(k-1)))*(((r_p*r_c^k)-1)/((k*r_p*(r_c-1))+(r_p-1))));// The thermal efficiency
+T_2=T_1*(r)^(k-1);// K
+T_3=T_2*r_p;// K
+T_4=T_3*r_c;// K
+q_in=(c_v*(T_3-T_2))+(c_p*(T_4-T_3));// kJ/kg
+w_out=n*q_in;// kJ/kg
+v_1=(R*T_1)/P_1;// m^3/kg
+MEP=w_out/(v_1*(1-(1/r)));// kPa
+printf("\nThe thermal efficiency,n=%0.3f or %2.1f percentage. \nThe heat input,q_in=%4.0f kJ/kg. \nThe work output,w_out=%4.0f kJ/kg. \nThe MEP=%4.0f kPa",n,n*100,q_in,w_out,MEP);
|