blob: 2ad6d2e058446e9feee7fe23134df2f56f5b2bdb (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
// Example 9_4
clc;funcprot(0);
// Given data
r=10;// The compression ratio
T_1=200+273;// K
P_1=200;// kPa
w_net=1000;// kJ/kg
c_v=0.717;// kJ/kg.K
R=0.287;// kJ/kg.K
k=1.4;// The specific heat ratio
// Calculation
T_2=T_1*(r)^(k-1);// K
function[X]=temperature(y)
X(1)=w_net-((c_v*(T_1-T_2))+(c_v*(y(1)-y(2))))
X(2)=y(1)-(y(2)*(r)^(k-1));
endfunction
y=[1000 1000];
z=fsolve(y,temperature);
T_3=z(1);// K
T_4=z(2);// K
n_carnot=(1-(T_1/T_3));// %
v_1=(R*T_1)/P_1;// m^3/kg
// v_2=v_1/r;
MEP=w_net/(0.9*v_1);// kPa
printf("\nThe maximum possible efficiency,n_carnot=%0.3f or %2.1f percentage.\nMEP=%4.0f kPa",n_carnot,n_carnot*100,MEP);
|
