Stability-Based Optimization of Planar Quiet Wind Tunnel Nozzles: Two-Dimensional Design and Three-Dimensional Assessment

Sintesi Tecnica: Ottimizzazione della Stabilità per Ugelli di Tunnel del Vento Planari Silenziosi

Definizione del Problema
La progettazione di ugelli per tunnel del vento supersonici silenziosi è fondamentale per investigare la transizione naturale dello strato limite (BL) in flussi ad alta velocità, poiché le installazioni convenzionali soffrono di disturbi nel flusso indisturbato elevati che compromettono l'estrapolazione dal suolo al volo. Mentre gli ugelli assialsimmetrici sono intrinsecamente bidimensionali (2D) nella simmetria del flusso, gli ugelli planari possiedono pareti laterali che introducono effetti tridimensionali (3D) inerenti. Le metodologie di progettazione standard spesso si affidano ad assunzioni 2D (flusso omogeneo nella direzione trasversale), trascurando l'espansione anisotropa causata dall'espansione verticale senza una corrispondente espansione trasversale. Questa semplificazione genera gradienti di pressione trasversali e strutture di flusso secondario, che possono indurre meccanismi di instabilità — come il crossflow — assenti nei modelli 2D. Il problema centrale affrontato è se le strategie di ottimizzazione basate sulla stabilità 2D rimangano robuste quando applicate ad ugelli planari in cui sono presenti effetti di parete 3D.

Metodologia
Lo studio impiega un framework di ottimizzazione modulare basato sulla stabilità, integrato in un ciclo di Ottimizzazione Bayesiana (BO). Il processo prevede tre fasi principali:

1. Generazione del Contorno Inviscido: Il toolkit HYPNOZE utilizza il Metodo delle Caratteristiche (MoC) per generare i contorni dell'ugello basandosi su parametri geometrici (curvatura della gola, angolo di apertura, lunghezza di rettificazione) e condizioni operative.
2. Analisi dello Strato Limite Viscoso: Il solver DEKAF calcola i profili dello strato limite laminare e lo spessore di spostamento, fornendo correzioni viscose alla geometria effettiva.
3. Valutazione della Stabilità: Il toolkit VESTA applica la Teoria della Stabilità Lineare Locale (LST) per valutare l'amplificazione dei disturbi. L'analisi si concentra sulle instabilità di primo modo, secondo modo e Görtler.

L'ottimizzazione minimizza una funzione di costo definita come la somma lineare dei fattori N ($N_{tot} = N_{1st} + N_{2nd} + N_{G\ddot{o}}$) per evitare l'insensibilità delle formulazioni radice-somma-quadrata rispetto alle deboli amplificazioni. Le variabili di progetto includono l'angolo di apertura dell'ugello ($\omega$), la lunghezza di rettificazione ($X_{BC}$) e una distribuzione parametrizzata della temperatura della parete ($T_w(x)$).

Per valutare la robustezza, i design 2D ottimizzati vengono rivalutati utilizzando simulazioni di Fluidodinamica Computazionale (CFD) completamente 3D (utilizzando CFD++) che includono esplicitamente le pareti laterali. L'analisi di stabilità su questi flussi base 3D utilizza sia LST che le Equazioni di Stabilità Paraboliche Lineari (LPSE) per catturare effetti non paralleli e non locali, in particolare lungo le pareti laterali dove si sviluppano flussi secondari.

Contributi Chiave

• Framework di Ottimizzazione 2D: Lo studio dimostra che un'ottimizzazione 2D basata sulla stabilità, che combina il controllo geometrico e la gestione termica attiva, può attenuare con successo le instabilità di primo modo e Görtler lungo la parete superiore di un ugello planare, raggiungendo livelli di amplificazione subcritici ($N_{combined} < 7.5$) ed estendendo la regione di flusso silenzioso.
• Valutazione della Robustezza 3D: Il documento fornisce una valutazione sistematica dei design 2D in un ambiente 3D, rivelando che mentre le previsioni di stabilità della parete superiore rimangono accurate, le regioni delle pareti laterali presentano meccanismi di instabilità distinti guidati da flussi secondari.
• Controllo Termico sulle Instabilità 3D: Un risultato significativo è che una distribuzione della temperatura della parete ottimizzata per i meccanismi 2D (specificamente un brusco calo della temperatura a valle della gola) attenua anche l'amplificazione di tipo crossflow lungo le pareti laterali, suggerendo un accoppiamento tra il controllo termico e i meccanismi di instabilità 3D.
• Limiti di Validità delle Teorie di Stabilità: Il confronto tra LST e LPSE delinea il dominio di validità delle teorie di stabilità lineare nei flussi planari 3D. L'integrazione di LPSE fallisce nelle regioni di forte non-parallelismo e rapida variazione streamwise (tipicamente a valle del punto di inflessione), indicando dove sono necessari approcci a maggiore fedeltà.

Risultati

• Esiti dell'Ottimizzazione: La configurazione ottimale (IT048) presenta un angolo di apertura poco profondo ($\approx 2^\circ$), una lunghezza di rettificazione breve e un profilo termico specifico: temperature elevate vicino alla gola (350 K) per sopprimere le onde di primo modo, seguite da un brusco calo a $\approx 292$ K per smorzare ulteriormente le instabilità, e una leggera risalita verso l'uscita. Questo design mantiene il flusso laminare lungo la parete superiore, a differenza del caso di riferimento che predice una transizione precoce.
• Prestazioni della Parete Superiore: Nelle simulazioni CFD 3D, le caratteristiche di stabilità della parete superiore del design ottimizzato rimangono coerenti con le previsioni 2D, con errori relativi inferiori al 5%. Il framework 2D funge da approssimazione affidabile di primo ordine per la parete superiore.
• Prestazioni della Parete Laterale: Lo strato limite della parete laterale presenta forti effetti 3D, inclusi vortici d'angolo e componenti di velocità trasversale. Questi inducono instabilità di tipo crossflow che non sono catturate dal modello 2D. Sebbene l'ottimizzazione geometrica e il controllo termico riducano i livelli di amplificazione rispetto al caso di riferimento, la soglia del fattore N ($N_{crit} \approx 7.5$) viene comunque raggiunta a circa $x/L \approx 0.5$ lungo le pareti laterali. Ciò suggerisce che la transizione indotta dalla parete laterale può verificarsi a monte della regione uniforme prevista, anche per i design ottimizzati.
• Interazione dei Meccanismi: Il profilo termico ottimizzato disaccoppia efficacemente la sovrapposizione spaziale delle regioni di amplificazione del primo modo e di Görtler. Inoltre, la strategia di controllo termico, sebbene derivata dalla fisica 2D, fornisce un beneficio secondario riducendo l'amplificazione del crossflow lungo le pareti laterali.

Significatività e Rivendicazioni
Il documento afferma che l'ottimizzazione basata sulla stabilità è uno strumento efficace e fisicamente interpretabile per controllare la transizione sulla parete superiore di ugelli planari silenziosi. Stabilisce che, sebbene l'ottimizzazione 2D sia sufficiente per la progettazione della parete superiore, essa è insufficiente per una valutazione completa dell'ugello perché non riesce a catturare le instabilità di crossflow indotte dalle pareti laterali.

Gli autori concludono che una progettazione affidabile degli ugelli per configurazioni planari richiede la considerazione esplicita degli effetti 3D. Sebbene l'attuale ottimizzazione 2D fornisca una base robusta, la presenza di rischi di transizione indotti dalle pareti laterali richiede l'integrazione di approcci a maggiore fedeltà (come la teoria PSE a piani o la teoria di stabilità BiGlobal spanwise) nei futuri framework di progettazione per risolvere completamente la dinamica delle instabilità multidirezionali. Lo studio evidenzia il controllo della temperatura della parete come una leva di progettazione versatile che influenza sia i meccanismi di instabilità 2D che 3D, offrendo una nuova via per mitigare la transizione delle pareti laterali.