Turbojet Engine With T-S Diagram
Turbojet Engine With T-S Diagram
The diffuser placed before the compressor decelerates the incoming air relative to
the engine. A pressure rise known as the ram effect is associated with this
deceleration.
The gas generator section consists of a compressor, combustor, and turbine, with
the same functions as the corresponding components of a stationary gas turbine
power plant. In a turbojet engine, the turbine power output need only be sufficient
to drive the compressor and auxiliary equipment, combustion gas leaving turbine
expand through the nozzle to a high velocity before being discharged to the
surroundings.
Goal of report to get familiar with matlab code and implementation on the
thermodynamic equation on turbojet engine.
Analysis
Assumption
Steady state
Each component is analyzed as a control volume at steady state.
The diffuser, compressor, turbine, and nozzle processes are isentropic.
There is no pressure drop for flow through the combustor.
The turbine work output equals the work required to drive the compressor.
Except at the inlet and exit of the engine, kinetic energy effects can be
ignored. Potential energy effects are negligible throughout.
The working fluid is air modeled as an ideal gas.
Thermodynamic analysis
T2=T1+(Va^2)/(2000)
T3=(rp)^(.4/1.4)*T2
T5=T4-(T3-T2)
T6=T5/((P5/P6)^(k-1)/k)
P2=P1*(T2/T1)^(k/(k-1))
P5=P4/((T4/T5)^(k/(k-1)))
rho1=(Patm/(R*Tatm))
m_dot=rho1*A1*Va
Q_dotin = m_dot*Cp*(T4-T3)
rho_6=(P6)/(R*T6)
A6=m_dot/(rho_6*V6)
F=m_dot*(V6-Va)
W_dotp=(F*Va)/1000
eta_th=(W_dotp/Q_dotin)*100
Code Development
Result of matlab and hand calculation are same
Matlab result
Va=285; %m/s
Patm=17; %kpa
Cp=1.005;
k=1.4;
R=(8.314/29);
Tatm=214 ;%K
d=0.5 ;%m
rp=18;
T1=Tatm
T2=Tatm+(Va^2)/(2000*Cp) %K
T3=(rp)^(.4/1.4)*T2
T4=1275; %K
T5=T4-(T3-T2)
T6=T5/((P5/P6)^(k-1)/k)
P1=Patm
STATE=[1;2;3;4;5;6];
TEMPERATURE=[T1;T2;T3;T4;T5;T6];
PRESSURE=[P1;P2;P3;P4;P5;P6];
t=table(STATE,PRESSURE,TEMPERATURE)
P2=P1*(T2/T1)^(k/(k-1))
P3=18*P2
P4=P3
P5=P4/((T4/T5)^(k/(k-1)))
P6=P1
rho1=(Patm/(R*Tatm));
A1=(pi*d^2)/4;
m_dot=rho1*A1*Va
Q_dotin = m_dot*Cp*(T4-T3) %Kj
V6=sqrt(2000*Cp*(T5-T6)) %m/s
rho_6=(P6)/(R*T6);
A6=m_dot/(rho_6*V6) %m^2
d6=sqrt((4*A6)/pi)
F=m_dot*(V6-Va)
W_dotp=(F*Va)/1000
eta_th=(W_dotp/Q_dotin)*100
Result
[2] https://www.mathworks.com/
[4] https://ecourses.ou.edu/cgi-
bin/ebook.cgi?doc&topic=th&chap_sec=09.1&page=theory