Comprendere il modello del ciclo di Ottone Air-Standard

Il ciclo Otto è l'idealizzazione termodinamica del motore a quattro tempi di accensione a scintilla, che sostituisce i complessi processi di assunzione, compressione, combustione e scarico con quattro passaggi reversibili che isolano i meccanismi di conversione dell'energia essenziale.

I quattro processi che formano il ciclo sono:

  • Compressione isttropica (1→2): Il pistone comprime la carica senza trasferimento di calore, aumentando la temperatura e la pressione.
  • Aggiunta di calore a volume costante (2→3)[: Un rilascio istantaneo di energia (rappresentante combustione) al centro morto superiore solleva la temperatura e la pressione ai loro valori di picco.
  • L'espansione isentropica (3→4)[: I gas ad alta pressione spingono il pistone verso il basso, fornendo il lavoro netto all'albero motore.
  • Rifiuti termici involume (→1)[]: Il contenuto del cilindro perde il calore allo stato iniziale, chiudendo il ciclo.

Poiché il modello utilizza un rapporto termico specifico costante e ignora lo scambio di gas, fornisce un punto di riferimento rapido che può essere sistematicamente raffinato. L'implementazione di questo ciclo in MATLAB ti dà uno strumento flessibile e interattivo per esplorare come rapporto di compressione, aggiunta di calore e proprietà fluide influiscono sulle prestazioni.

Impostazione di un ambiente MATLAB strutturato

Iniziare creando un nuovo script o file di funzione. Uno spazio di lavoro pulito rende debug e future modifiche semplici. Utilizzare un intestazione descrittiva che documenta lo scopo, i parametri di input e le uscite. Pre-allocate array per le variabili di stato quando si spazzano oltre i parametri per evitare la resissione dinamica overhead.

% 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

Utilizzando gli array di MATLAB o ] per memorizzare i quattro stati migliora la leggibilità quando si tracciano o esportano i risultati. Ad esempio, creare un array di struttura con i campi , , ] indicizzato da 1 a 4.

Definizione dei parametri del motore critico

L'accuratezza della simulazione dipende dai dati di ingresso. Al minimo, specificare il rapporto di compressione, la pressione di assunzione e la temperatura, il rapporto termico specifico e l'energia aggiunta per massa unitaria.

  • Rapporto di compressione r[: 8–12 per i motori a accensione a scintilla moderni; i rapporti più elevati migliorano l'efficienza ma aumentano il rischio di bussare.
  • Temperatura di immissione T1[[: 300 K (ambient) o 350 K se si considera il riscaldamento dell'aria di assunzione dai gas collettori e residui.
  • Pressione di immissione P1[[: 101325 Pa (1 atm) per aspirare naturalmente; fino a 200 kPa o più per i motori turbocompresso.
  • Rapporto termico specifico γ[: 1.4 per aria secca; 1.3–1.35 per miscele di aria combustibile e prodotti di combustione a causa della presenza di molecole triatomiche come CO2 e H2O.
  • Aggiunta termica specifica Q in[[]: circa 2,5–3 MJ/kg per benzina basata sul valore di riscaldamento inferiore e sul rapporto tra aria e gasolio stoichiometrico. Per una tipica benzina, il valore di riscaldamento inferiore è di circa 44 MJ/kg, ma miscelato con aria (circa 14.7:1 rapporto aria-fuel), la densità di energia miscela diventa di circa 2.9 MJ/kg.

Definire queste variabili MATLAB nella parte superiore dello script in modo che possano essere facilmente modificate per gli studi di 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)

Si può anche includere la massa molare d'aria (28,97 g/mol) se si passa poi a calcoli basati su mole, mantenendo il gas costante esplicite evita errori di unità.

Computing dei quattro Stati d'angolo

Con i parametri impostati, le relazioni standard dell'aria danno i punti di stato direttamente. L'analisi è per massa unitaria, quindi lavoriamo con volume specifico. Impostare il volume di clearance all'unità (v2 = 1) semplifica la scala relativa; quindi v1 = r perché il rapporto di compressione v1/v2 = r. Questa normalizzazione è conveniente perché tutti i volumi sono espressi rispetto al volume di clearance.

Stato 1 (Impresa, Centro Morto inferiore)

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

Stato 2 (Compressione isttropica al centro superiore morto)

Per un processo isentropico con costante γ:

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

Queste relazioni derivano dalla legge del gas ideale e dalla relazione politropica [P v^γ = costante[[]. L'aumento della temperatura attraverso la compressione è significativo: per r=9 e γ=1.35, T2 ≈ 300 × 9^(0.35) ≈ 690 K.

Stato 3 (Aggiunta di calore di corrente)

Prima calcolare il calore specifico a volume costante:

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

La temperatura a stato 3 può superare 2500 K per la combustione tipica della benzina. La pressione di picco P3 raggiunge spesso 60–80 bar, che imposta i requisiti strutturali per il blocco motore e la testa.

Stato 4 (espansione isttropica al volume originale)

Il rapporto di espansione da TDC a BDC è r, ma in termini di volumi specifici v3 = v2 = 1 e v4 = v1 = r, quindi il rapporto di volume è 1/r.

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

Questi quattro punti definiscono completamente il ciclo in termini di temperatura, pressione e volume specifico. Conservarli in array per la successiva trama e calcolo. Una validazione utile è controllare che P4 è leggermente superiore a P1; la differenza guida il soffiaggio di scarico.

Calcolo del lavoro netto, efficienza termica e pressione efficace media

Il lavoro specifico netto è la differenza tra il calore aggiunto e il calore rifiutato. Per il modello standard dell'aria, il calore rifiutato nel processo 4→1 è:

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

L'efficienza termica è:

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 pressione effettiva media (MEP)[[]] è una metrica di prestazione che confronta il lavoro netto per ciclo al volume di spostamento. Per il ciclo standard dell'aria, il volume di spostamento specifico è v1 – v2 = r – 1. Pertanto:

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

L'euro è indipendente dalle dimensioni del motore ed è un indicatore conveniente della qualità del motore. I valori tipici per i motori a benzina naturalmente aspirati sono 8-12 bar (0.8–1.2 MPa).

Compresi i riscaldamenti specifici variabili per una migliore precisione

Per questo motivo, sostituire il γ costante con una correlazione dipendente dalla temperatura. Un metodo robusto è quello di utilizzare i coefficienti polinomiali della NASA per c]p[FIST:2] di specie principali (O2, N

L’applicazione della variabile γ richiede una soluzione iterativa perché l’esponente isentropico stesso dipende dalla temperatura sconosciuta.

% 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 risorsa polinomiale [NASA[[[[]] offre coefficienti scaricabili per le specie di aria e combustione. Modificare lo script per caricare questi coefficienti e chiamarli durante i calcoli di stato. Osservate che le temperature specifiche variabili riducono l'efficienza prevista da 2–5 punti percentuali a rapporti di compressione tipici, che si allineano strettamente con i dati reali del motore.

Per una rapida approssimazione senza integrazione polinomiale, utilizzare un rapporto termico specifico medio che varia linearmente con la temperatura: γ(T) = 1.4 – 0.00005*(T – 300). Questa semplice correlazione cattura la tendenza e può essere implementata in un loop di tempo fino alla convergenza.

Costruire un'analisi funzionale modulare

Piuttosto che uno script monolitico, incapsula la logica del ciclo in una funzione MATLAB, che lo rende riutilizzabile per le spazzature dei parametri e riduce la duplicazione del codice.

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

Per le heats specifiche variabili, aggiungere una bandiera fluida:

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

Con la funzione in atto, uno script del driver può passare attraverso i rapporti di compressione:

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--');

Questo approccio modulare mostra immediatamente i rendimenti diminutivi del rapporto di compressione crescente oltre i 12 anni. La trama evidenzia un trade-off chiave: l'efficienza continua ad aumentare ma ad un ritmo decrescente, mentre l'eurodeputato spesso raggiunge i 10 giri per gli input di calore tipici.

Generando diagrammi diagnostici P‐V e T‐S

La visualizzazione aiuta a verificare il modello e comunicare i risultati. Le funzioni di tramamento di MATLAB producono diagrammi di qualità della pubblicazione.

Diagramma di pressione-Volume (P‐V)

Le curve di compressione e di espansione seguono P v^γ = costante[. Generare una serie di punti intermedi per curve lisce:

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;

L'area chiusa rappresenta il lavoro netto. È possibile integrare numericamente il loop utilizzando [] per verificare la W net calcolata. Ad esempio, integrare la curva superiore (compressione + espansione) e sottrarre la curva inferiore (linea di rifiuto del calore).

Diagramma di temperatura-introduzione (T‐S)

Per un gas ideale con costanti terme specifiche, il cambiamento entropia rispetto allo stato 1 è:

% 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);

Per la variabile γ, utilizzare l'integrazione numerica di dQ/T per calcolare l'entropia. Il diagramma T‐S rivela la natura isentropica dei processi 1‐2 e 3‐4 (linee verticali) e il riscaldamento/rigettazione costante-volume (percorsi curvati). L'area sotto la curva di riscaldamento rappresenta il calore aggiunto, mentre l'area sotto la curva di raffreddamento rappresenta il lavoro netto.

Per migliorare la trama T‐S, sovrapporre le linee di pressione costante (isobar) dalla relazione ideale del gas, che contribuisce a contestualizzare il ciclo in un quadro termodinamico più ampio.

Studi sulla sensibilità

Una volta che il modello di base funziona, varia sistematicamente i parametri di input per capire la loro influenza.

  • Compressione rapporto spazza[[[]: efficienza del lotto e europ versus r. Efficienza aumenta rapidamente a bassa r e poi altipiani; L'eurodeputato spesso raggiunge una velocità di rotazione moderata a causa del trade-off tra aumento della pressione e aggiunta di calore.
  • Variazione di aggiunta di calore : Cambiare i modelli Q in di diverse energie di combustibile o di aria-fuel ratio. Aumentare Q in aumenta la pressione di picco e il lavoro ma può superare i limiti materiali; includere un vincolo che P3 < 150 bar (15 MPa) per riflettere la forza tipica del motore.
  • Effetto del rapporto termico specifico[[]: Il basso γ riduce l'efficienza.

Per una spazzata a due parametri, utilizzare i loop nidificati e i risultati del negozio in una matrice. Quindi creare un diagramma di contorno o una mappa di colore per visualizzare l'effetto combinato:

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 (%)');

Tali mappe sensibili sono inestimabili per l'analisi del commercio-off nella progettazione del motore, che aiutano a identificare le regioni operative in cui l'efficienza è elevata mentre la pressione di picco rimane accettabile.

Convalida dei risultati

Prima di utilizzare il modello per le decisioni di ingegneria, controllare contro i benchmark noti. La formula di efficienza standard dell'aria η = 1 − 1/r^(γ−1) fornisce una validazione diretta per il caso costante-γ. Inoltre, confrontare la pressione di picco calcolata P3 contro i dati delle mappe del motore pubblicate o degli esempi di libro di testo (ad esempio, Stone, Introduzione ai motori di combustione interna

Un'altra tecnica di validazione è quella di controllare il [ equilibrio energetico[[]: il lavoro netto più il calore rifiutato dovrebbe pari al calore di ingresso all'interno della precisione a punto variabile. Verificare inoltre che i processi di compressione e di espansione soddisfino la relazione isentropica []T v^(γ‐1) = costante]]]]] per rilevare errori di codifica.

Estendere il Modello per Effetti Realistici

Il quadro di base può essere ampliato per includere fenomeni del mondo reale:

  • Durata della combustione dei grassi[[]: Sostituire l’aggiunta di calore istantanea con una funzione Wiebe che distribuisce il rilascio di energia su un intervallo di angolo della manovella.
  • Perdite di trasferimento di calore[[]: Incorpora una correlazione di tipo Woschni per il trasferimento di calore convettivo alle pareti del cilindro. Ciò richiede temperatura e pressione risolte a manovella. La correlazione dà il coefficiente di trasferimento di calore come funzione di foro, velocità del pistone, temperatura e pressione.
  • Ricircolo gas di scarico (EGR)[]: Regolare il rapporto termico specifico e la composizione del gas per riflettere la diluizione. Il blocco MATLAB Simscape Engine Cylinder[ fornisce un componente pronto per tali studi, ma è anche possibile implementare un modello di di diluizione semplificato combinando linearmente le heat specifiche di gas.

Per una rapida estensione, inizia con la funzione Wiebe. Esempio approssimazioni si possono trovare nella MALAB File Exchange[]]. Questa semplice aggiunta cattura l'effetto della tempistica della scintilla sull'efficienza e sulla pressione di picco.

Integrare l'analisi in un flusso di lavoro più ampio

In un progetto di ingegneria, l’analisi del ciclo deve essere utilizzata per l’esportazione di tabelle di stato e metriche di performance in CSV per un’ulteriore elaborazione in Excel. Combinare i grafici in un PDF usando la funzione .

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

Per uno studio ibrido, la mappa dell'efficienza del ciclo di Otto può essere importata come una tabella di ricerca a Simulink.

Conclusioni

Con la definizione sistematica dei parametri, dei punti di stato di calcolo e delle metriche di prestazione calcolanti come il lavoro netto, l'efficienza termica e la pressione effettiva media, si ottiene immediato intuito su come le scelte di progettazione influiscono sull'output. L'approccio modulare presentato qui consente di estendere il modello a intervalli specifici variabili, la combustione finita e il trasferimento di calore senza riscrivere i dati di base.