"Homework 20"


gamma = 1.35
AF = 15
x_r = 0.06
r_c = 9.2
V_d = 0.0022
eta_m = 0.84
Q_HV = 44000
mfm = (1-x_r)/(1+AF)

"Volume Calculations"

V_1 = V_d / (1-1/r_c)
V_2 = V_1 -V_d

"Point 1"

T_1 = 300
P_1 = 100
u_1 = intenergy(Air,T=T_1)
s_1 = entropy(Air,T=T_1,P=P_1)
sv_1 = volume(Air,T=T_1,P=P_1)
mass = V_1 / sv_1

"Point 2"

T_2a = T_1 * r_c^(gamma-1)
P_2a = P_1 * r_c^gamma
u_2 = intenergy(Air,T=T_2b)
s_2 = entropy(Air,T=T_2b,P=P_2b)
sv_2 = volume(Air,T=T_2b,P=P_2b)
s_2=s_1
sv_2=sv_1/r_c
c_v = CV(Air,T=T_2a)
W_12b = -mass*(u_2-u_1)
W_12a = (P_2a*V_2 - P_1*V_1)/(1-gamma)


"Point 3"

mf = mass * mfm
mfm * Q_HV = c_v * (T_3a - T_2a)
P_3a / T_3a = P_2a / T_2a
u_3 = u_2 + mfm* Q_HV
u_3 = intenergy(Air,T=T_3b)
s_3 = entropy(Air,T=T_3b,P=P_3b)
sv_3 = volume(Air,T=T_3b,P=P_3b)
sv_3=sv_2


"Point 4"

P_4a = P_3a * (1/r_c^gamma)
T_4a = T_3a * ( 1/r_c^(gamma-1))
u_4 = intenergy(Air,T=T_4b)
s_4 = entropy(Air,T=T_4b,P=P_4b)
sv_4 = volume(Air,T=T_4b,P=P_4b)
sv_4=sv_1
s_4 = s_3
w_34a = (P_4a * V_1 - P_3a*V_2)/(1-gamma)
w_34b = mass*(u_3-u_4)

"Net work per cycle"

W_a = w_34a + w_12a
W_b = w_34b + w_12b

"Fuel Conversion Efficiency"

eta_a = W_a / (mf*Q_HV)
eta_b = W_b /(mf*Q_HV)