Comprendre le modèle de cycle Otto de la norme aérienne

Le cycle Otto est l'idéalisation thermodynamique du moteur à quatre temps à allumage par étincelles. Il remplace les processus complexes d'admission, de compression, de combustion et d'échappement par quatre étapes réversibles qui isolent les mécanismes de conversion d'énergie essentiels. Dans le modèle standard de l'air, le fluide de travail est l'air (traité comme un gaz idéal avec des chaleurs spécifiques constantes), la combustion est approximative comme addition de chaleur instantanée à volume constant, et l'échappement est remplacé par un rejet de chaleur en volume constant.

Les quatre processus qui forment le cycle sont les suivants :

  • Compresse isentropique (1→2): Le piston compresse la charge sans transfert de chaleur, avec élévation de la température et de la pression.
  • Ajout de chaleur en volume constant (2→3): Un dégagement instantané d'énergie (représentant la combustion) au centre mort supérieur augmente la température et la pression à leurs valeurs maximales.
  • Dilatation isentropique (3→4): Les gaz à haute pression poussent le piston vers le bas, livrant le travail net au vilebrequin.
  • Rejet de chaleur en volume constant (4→1): Le contenu du cylindre perd de la chaleur à l'état initial, fermant la boucle.

Comme le modèle utilise un rapport thermique spécifique constant et ignore l'échange de gaz, il fournit un repère rapide qui peut être systématiquement affiné. La mise en œuvre de ce cycle dans MATLAB vous offre un outil flexible et interactif pour explorer comment le rapport de compression, l'ajout de chaleur et les propriétés des fluides affectent les performances. L'hypothèse air-standard fournit également une limite supérieure pratique sur l'efficacité; les moteurs réels sont toujours en deçà du niveau de frottement, de pertes de chaleur et de combustion incomplète.

Mise en place d'un environnement MATLAB structuré

Commencez par créer un nouveau fichier de script ou de fonction. Un espace de travail propre rend le débogage et les modifications futures simples. Utilisez un en-tête descriptif qui documente le but, les paramètres d'entrée et les sorties. Pré-alternez les tableaux pour les variables d'état lorsque vous balayez les paramètres pour éviter un redimensionnement dynamique des frais généraux. Adoptez des noms de variables significatifs tels que plutôt que simplement (bien que ] soit standard dans les textes thermodynamiques) et maintenez les unités cohérentes dans l'ensemble (unités SI recommandées). Un bloc d'en-tête typique peut lire :

% otto_cycle_analysis.m
% Air‑standard Otto cycle analysis
% Inputs:
% r – compression ratio (dimensionless)
% T1, P1 – intake state (K, Pa)
% gamma – specific heat ratio (cp/cv)
% Q_in – specific heat addition (J/kg)
% Outputs:
% state_table – array of T, P, v at four corners
% W_net, eta – specific net work (J/kg) and thermal efficiency

En utilisant les tableaux MATLAB=2 intégrés ou pour stocker les quatre états, vous améliorez la lisibilité lorsque vous tracez ou exportez les résultats plus tard. Par exemple, créez un tableau de structure avec des champs , , indexés de 1 à 4. Cela rend intuitives les propriétés de l'état de référence : donne la température après combustion.

Définition des paramètres critiques du moteur

La précision de la simulation dépend des données d'entrée. Au minimum, spécifiez le rapport de compression, la pression et la température d'entrée, le rapport de chaleur spécifique et l'énergie ajoutée par unité de masse. Les valeurs typiques pour un moteur à essence sont:

  • Rapport de compression r: 8-12 pour les moteurs modernes à allumage par étincelles; des rapports plus élevés améliorent l'efficacité mais augmentent le risque de collision.
  • Température d'admission T1: 300 K (ambient) ou 350 K si l'air d'admission est chauffé par le collecteur et les gaz résiduels.
  • Pression d'admission P1: 101325 Pa (1 atm) pour les moteurs à aspiration naturelle; jusqu'à 200 kPa ou plus pour les moteurs turbocompressés.
  • Rapport thermique spécifique γ: 1,4 pour l'air sec; 1,3–1,35 pour les mélanges d'air combustible et les produits de combustion en raison de la présence de molécules triatomiques comme le CO2 et le H2O.
  • Ajout de chaleur spécifique Q in: environ 2,5–3 MJ/kg pour l'essence, en fonction de la valeur de chauffage inférieure et du rapport stœchiométrique air-carburant. Pour une essence typique, la valeur de chauffage inférieure est d'environ 44 MJ/kg, mais mélangée avec l'air (rapport air-carburant d'environ 14,7:1), la densité énergétique du mélange devient d'environ 2,9 MJ/kg de mélange.

Définissez ces variables comme des variables MATLAB en haut du script afin qu'elles puissent être facilement modifiées pour les études de sensibilité.

r = 9.0; % compression ratio
T1 = 300; % K
P1 = 101325; % Pa
gamma = 1.35; % approximate for fuel‑air mixture
Q_in = 2500e3; % J/kg
R_air = 287; % J/(kg·K)

Vous pouvez également inclure la masse molaire d'air (28,97 g/mol) si vous passez ultérieurement aux calculs basés sur la mole.

L'informatique dans les quatre coins

Avec les paramètres définis, les relations air-standard donnent directement les points d'état. L'analyse est par unité de masse, donc nous travaillons avec un volume spécifique. La détermination du volume de dégagement à l'unité (v2 = 1) simplifie l'échelle relative; puis v1 = r parce que le rapport de compression v1/v2 = r. Cette normalisation est pratique parce que tous les volumes sont exprimés par rapport au volume de dégagement.

État 1 (entrée, centre de mortalité du bas)

v1 = r; % relative to clearance volume
% T1 and P1 are inputs

État 2 (Compression isentropique vers le haut du centre mort)

Pour un processus isoentropique avec constante γ:

v2 = 1;
T2 = T1 * r^(gamma-1);
P2 = P1 * r^gamma;

Ces relations proviennent de la loi idéale sur les gaz et de la relation polytropique P v^γ = constante. La hausse de température à travers la compression est significative: pour r=9 et γ=1.35, T2 γ 300 × 9^(0.35) ↓ 690 K.

État 3 (ajout de chaleur continue-volume)

Calculer d'abord la chaleur spécifique à volume constant:

cv = R_air / (gamma-1);
T3 = T2 + Q_in / cv;
P3 = P2 * (T3 / T2);
v3 = v2; % no volume change

La température à l'état 3 peut dépasser 2500 K pour la combustion typique d'essence. La pression maximale P3 atteint souvent 60 à 80 bar, ce qui établit les exigences structurelles pour le bloc moteur et la tête.

État 4 (Extension isentropique vers le volume original)

Le rapport d'expansion entre TDC et BDC est r, mais en termes de volumes spécifiques v3 = v2 = 1 et v4 = v1 = r, donc le rapport de volume est 1/r.

v4 = v1;
T4 = T3 * (1/r)^(gamma-1);
P4 = P3 * (1/r)^gamma;

Ces quatre points définissent le cycle en termes de température, de pression et de volume spécifique. Conservez-les dans des tableaux pour le tracé et le calcul ultérieurs. Une validation utile est de vérifier que P4 est légèrement au-dessus de P1; la différence entraîne la chute des gaz d'échappement.

Calcul du travail net, efficacité thermique et pression efficace moyenne

Le travail spécifique net est la différence entre la chaleur ajoutée et la chaleur rejetée. Pour le modèle de la norme air, la chaleur rejetée dans le processus 4→1 est:

Q_out = cv * (T4 - T1);
W_net = Q_in - Q_out;

L'efficacité thermique est:

eta = W_net / Q_in; % from energy balance
eta_direct = 1 - 1/r^(gamma-1); % direct formula
% Both should match within numerical round-off

La pression effective moyenne (PEM)[ est une mesure de performance qui compare le travail net par cycle au volume de déplacement. Pour le cycle de vol normalisé de l'air, le volume de déplacement spécifique est v1 – v2 = r – 1.

mep = W_net / (v1 - v2); % Pa

Les valeurs typiques pour les moteurs à essence à aspiration naturelle sont de 8-12 bar (0,8–1,2 MPa). Sortie MEP avec efficacité pour fournir une image plus complète. Vous pouvez également convertir en barre pour la lisibilité: .

Y compris les chaleurs spécifiques variables pour une meilleure exactitude

Pour tenir compte de cette hypothèse, remplacer la constante γ par une corrélation dépendante de la température. Une méthode robuste consiste à utiliser les coefficients polynômes de la NASA pour cp[] des principales espèces (O2, N2, CO2, H2O). Le NIST Chemistry WebBook[ fournit ces coefficients pour les substances pures. Pour le modèle de standard air, vous pouvez approximativement le fluide de travail comme air sec ou comme produit de combustion - avec une composition fixe.

La mise en œuvre de la variable γ nécessite une solution itérative car l'exposant isentropique lui-même dépend de la température inconnue. Utilisez MATLABs ou une simple boucle. Par exemple, pendant la compression isentropique, la relation devient :

% Use an anonymous function to find T2 such that
% integral of cp dT from T1 to T2 = R * ln(v1/v2)
% This replaces the constant‑gamma formula.
% For air, cp(T) can be approximated by:
% cp = a0 + a1*T + a2*T^2 + a3*T^3 + a4*T^4
% Using NASA 7-coefficient polynomials.
fun = @(T2) integral(@(T) polyval(cp_coeffs, T), T1, T2) - R_air * log(r);
T2 = fzero(fun, T1 * r^0.4);

La ressource NASA polynôme offre des coefficients téléchargeables pour les espèces d'air et de combustion. Modifiez votre script pour charger ces coefficients et les appeler lors des calculs d'état. Vous observerez que les chaleurs spécifiques variables réduisent l'efficacité prévue de 2 à 5 points de pourcentage aux rapports de compression typiques, ce qui s'harmonise étroitement avec les données réelles du moteur.

Pour une approximation rapide sans intégration complète des polynômes, utilisez un rapport thermique spécifique moyen qui varie linéairement avec la température : γ(T) = 1,4 – 0,00005*(T – 300). Cette corrélation simple capture la tendance et peut être mise en œuvre dans une boucle de temps jusqu'à la convergence.

Construire une analyse modulaire de la fonction

Plutôt qu'un script monolithique, encapsuler la logique de cycle dans une fonction MATLAB. Cela la rend réutilisable pour les balayages de paramètres et réduit la duplication de code. Un squelette de fonction est montré ci-dessous:

function [states, performance] = ottoCycle(r, T1, P1, Q_in, gamma, R_air)
 % states: structure with T, P, v for states 1-4
 % performance: structure with W_net, eta, mep
 % (calculation code from previous section goes here)
end

Pour les chaleurs spécifiques variables, ajouter un drapeau de fluide:

function [states, performance] = ottoCycle(r, T1, P1, Q_in, fluid_type)
 % fluid_type: 'constant_gamma' or 'variable_cp'
 % (logic branches based on flag)
end

Avec la fonction en place, un script pilote peut balayer les rapports de compression:

r_values = 6:1:14;
for i = 1:length(r_values)
 [~, perf] = ottoCycle(r_values(i), 300, 101325, 2500e3, 'constant_gamma');
 eta(i) = perf.eta;
 mep(i) = perf.mep;
end
plot(r_values, eta*100, 'b-', r_values, mep/1e5, 'r--');

Cette approche modulaire montre immédiatement la diminution du rendement de l'augmentation du taux de compression au-delà de 12. La courbe met en évidence un compromis clé : l'efficacité continue d'augmenter mais à un rythme décroissant, tandis que le MEP atteint souvent des sommets autour de r=10 pour les apports thermiques typiques.

Générer des diagrammes de diagnostic P‐V et T‐S

La visualisation aide à vérifier le modèle et à communiquer les résultats. Les fonctions de tracé MATLAB , produisent des diagrammes de qualité de publication.

Diagramme de pression-volume (P‐V)

Les courbes de compression et d'expansion suivent P v^γ = constante. Générer une série de points intermédiaires pour les courbes lisses:

n = 50;
v_comp = linspace(v2, v1, n); % reverse direction for plotting
P_comp = P2 * (v_comp ./ v2).^(-gamma);
v_exp = linspace(v3, v4, n);
P_exp = P3 * (v_exp ./ v3).^(-gamma);
figure;
hold on;
plot(v_comp, P_comp, 'b', 'LineWidth', 1.5);
plot(v_exp, P_exp, 'b', 'LineWidth', 1.5);
plot([v2 v3], [P2 P3], 'r--', 'LineWidth', 1.2);
plot([v4 v1], [P4 P1], 'r--', 'LineWidth', 1.2);
xlabel('Specific volume (relative)'); ylabel('Pressure (Pa)');
title('Air‑Standard Otto Cycle P–V Diagram');
legend('Compression', 'Expansion', 'Heat addition', 'Heat rejection');
grid on;

La zone fermée représente le travail net. Vous pouvez intégrer numériquement la boucle en utilisant pour vérifier le W net calculé. Par exemple, intégrer la courbe supérieure (compression + expansion) et soustraire la courbe inférieure (ligne de rejet de chaleur).

Diagramme température-entropie (T-S)

Pour un gaz idéal avec des chaleurs spécifiques constantes, le changement d'entropie par rapport à l'état 1 est :

% s = cp * log(T / T1) - R * log(P / P1)
cp = gamma * R_air / (gamma - 1);
s1 = 0;
s2 = cp * log(T2/T1) - R_air * log(P2/P1);
s3 = cp * log(T3/T1) - R_air * log(P3/P1);
s4 = cp * log(T4/T1) - R_air * log(P4/P1);

Pour la variable γ, utilisez l'intégration numérique de dQ/T pour calculer l'entropie. Le diagramme T‐S révèle la nature isotrope des processus 1‐2 et 3‐4 (lignes verticales) et du chauffage/rejet à volume constant (chemins courbes). La zone sous la courbe de chauffage représente la chaleur ajoutée, tandis que la zone sous la courbe de refroidissement représente la chaleur rejetée; leur différence est le travail net.

Pour améliorer la courbe T‐S, on place les lignes de pression constante (isobars) à partir de la relation de gaz idéale. Cela permet de contextualiser le cycle dans un cadre thermodynamique plus large.

Effectuer des études de sensibilité

Une fois le modèle de base fonctionnel, varier systématiquement les paramètres d'entrée pour comprendre leur influence.

  • : L'efficacité du lot et l'efficacité du MEP par rapport à r. augmente rapidement à des niveaux r bas puis à des plateaux; le MEP atteint souvent des sommets à un r modéré en raison de l'écart entre l'augmentation de la pression et l'ajout de chaleur.
  • Variation de l'ajout de chaleur : changer les modèles Q in différentes énergies de carburant ou rapports air-carburant. L'augmentation de Q in augmente la pression maximale et le travail, mais peut dépasser les limites de matériaux; inclure une contrainte que P3 < 150 bar (15 MPa) pour refléter la résistance typique du moteur.
  • Effet spécifique du rapport thermique[: γ réduit le rendement. Les résultats de la constante de recouvrement‐γ et de la variable‐γ montrent la correction. Ceci est particulièrement important pour l'analyse de carburants alternatifs comme le gaz naturel ou l'hydrogène.

Pour un balayage à deux paramètres, utilisez des boucles imbriquées et entreposez les résultats dans une matrice. Puis créez un tracé de contour ou une carte couleur pour visualiser l'effet combiné :

r_vec = 6:1:14;
Q_vec = 2e6:0.2e6:3e6;
for i = 1:length(r_vec)
 for j = 1:length(Q_vec)
 [~, perf] = ottoCycle(r_vec(i), 300, 101325, Q_vec(j), 1.35, 287);
 eta_mat(i,j) = perf.eta;
 end
end
contourf(r_vec, Q_vec/1e6, eta_mat', 20);
xlabel('Compression ratio'); ylabel('Q_in (MJ/kg)');
colorbar; title('Thermal Efficiency (%)');

Ces cartes de sensibilité sont précieuses pour l'analyse des compromis dans la conception des moteurs. Elles aident à identifier les régions d'exploitation où l'efficacité est élevée alors que la pression maximale demeure acceptable.

Valider vos résultats

Avant d'utiliser le modèle pour les décisions d'ingénierie, vérifiez par rapport aux repères connus.La formule d'efficacité standard de l'air η = 1 - 1/r^(γ−1) fournit une validation directe pour le cas de constante-γ. En outre, comparez la pression de crête calculée P3 avec les données des cartes de moteurs ou exemples de manuels publiés (p. ex. Stone, Introduction aux moteurs à combustion interne.Les écarts sont souvent attribuables à la supposition de γ ou de Q in. Pour une validation approfondie, téléchargez une trace de pression de cylindre publique à partir de ressources comme la base de données des moteurs des Laboratoires nationaux de Sandia et superposez votre boucle P-V. Cet exercice révèle que le modèle de standard de l'air se démarque de la réalité, motivant une modélisation plus complexe.

Une autre technique de validation consiste à vérifier le bilan énergétique : le travail net plus la chaleur rejetée doit égaler la chaleur d'entrée dans la précision du point flottant. Vérifiez également que les processus de compression et d'expansion satisfont à la relation isotrope T v^(γ‐1) = constante pour détecter les erreurs de codage.

Élargir le modèle pour les effets réalistes

Le cadre de base peut être élargi pour inclure les phénomènes du monde réel:

  • Durée de combustion finale[: Remplacer l'ajout instantané de chaleur par une fonction Wiebe qui distribue la libération d'énergie sur un intervalle d'angle de vilebrequin. Utilisez MATLAB= pour intégrer l'équation d'énergie avec le mouvement du piston. La fonction Wiebe est définie comme x b(γ) = 1 – exp(–a * (((===0)/Δ=)^(m+1)).
  • Pertes de transfert de chaleur[: Incorporer une corrélation de type Woschni pour le transfert convectif de chaleur aux parois du cylindre. Cela nécessite une température et une pression résolues par l'angle de vilebrequin. La corrélation donne le coefficient de transfert de chaleur en fonction de l'alésage, de la vitesse du piston, de la température et de la pression.
  • Recirculation des gaz d'échappement (EGR)[: Réglez le rapport thermique spécifique et la composition des gaz pour refléter la dilution. Le bloc MATLAB Simscape Engine Cylinder fournit un composant prêt à l'emploi pour de telles études, mais vous pouvez également mettre en œuvre un modèle de dilution simplifié en combinant linéairement les chaleurs spécifiques de l'air et du gaz inerte.

Pour une extension rapide, commencez par la fonction Wiebe. Des approximations par exemple se trouvent dans le MATLAB File Exchange. Cette simple addition capture l'effet du timing de l'étincelle sur l'efficacité et la pression maximale.

Intégration de l'analyse dans un flux de travail plus important

Dans un projet d'ingénierie, l'analyse de cycle doit alimenter les simulations de niveau véhicule ou la génération de rapports. Utilisez MATLABs pour exporter des tableaux d'état et des mesures de performance vers CSV pour un traitement ultérieur en Excel. Combinez les graphiques en PDF en utilisant la fonction .

results = table(r_values', eta', mep', ...
 'VariableNames', {'CompressionRatio','Efficiency','MEP_bar'});
writetable(results, 'otto_sweep.csv');

Vous pouvez également appeler MATLAB de Python ou Simulink pour coupler le modèle thermodynamique avec d'autres sous-systèmes. Pour une étude hybride de véhicules, la carte d'efficacité du cycle Otto peut être importée comme une table de recherche dans Simulink.

Conclusion

La conception d'un modèle de cycle Otto standard dans le MATLAB offre une plate-forme puissante mais simple pour explorer la thermodynamique des moteurs à allumage par étincelles. En définissant systématiquement les paramètres, en calculant les points d'état et en calculant les paramètres de performance tels que le travail net, l'efficacité thermique et la pression moyenne efficace, vous obtenez une compréhension immédiate de la façon dont les choix de conception affectent la sortie. L'approche modulaire présentée ici vous permet d'étendre le modèle à des chaleurs variables spécifiques, à la combustion finie et au transfert de chaleur sans réécrire la logique de base.