Analyse thermodynamique du moteur turbojet
Analyse thermodynamique du moteur turbojet
Le moteur à turboréacteur est couramment utilisé pour la propulsion des avions en raison de leur
décélération.
La section du générateur de gaz se compose d'un compresseur, d'un brûleur et d'une turbine, avec
les mêmes fonctions que les composants correspondants d'une turbine à gaz stationnaire
centrale électrique. Dans un moteur turbojet, la puissance de sortie de la turbine doit seulement être suffisante
pour entraîner le compresseur et les équipements auxiliaires, les gaz de combustion sortant de la turbine
s'étendre à travers la buse à une vitesse élevée avant d'être déchargé dans le
alentours.
Objectif du rapport pour se familiariser avec le code Matlab et son implémentation sur le
État stable
Chaque composant est analysé comme un volume de contrôle à l'état stable.
La puissance de sortie de la turbine est égale au travail nécessaire pour entraîner le compresseur.
À part à l'entrée et à la sortie du moteur, les effets de l'énergie cinétique peuvent être
ignored. Potential energy effects are negligible throughout.
Le fluide de travail est l'air modélisé comme un gaz idéal.
Analyse thermodynamique
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
Développement de code
Le résultat de MATLAB et le calcul manuel sont les mêmes
Résultat Matlab
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)^(0.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] [Link]
[4][Link]
bin/[Link]?doc&topic=th&chap_sec=09.1&page=théorie