Analisi comparativa di schemi risolutivi comprimibili per l’aerodinamica di getti sottosfusi con convalida Schlieren

Perché i getti di gas veloci sono importanti

Dai becchi dei razzi alle prese d’aria industriali e perfino alle eruzioni vulcaniche, i getti di gas ad alta velocità determinano come spinta, rumore e calore si propagano nell’aria circostante. Prevedere con precisione il comportamento di questi getti è complesso, perché onde di gas compresso ed espanso formano modelli intricati che gli strumenti computazionali standard fanno fatica a seguire. Questo studio pone una domanda pratica: un tipo di solver meno costoso, solitamente riservato a flussi più lenti, può affrontare questi problemi impegnativi senza perdere accuratezza?

Due modi diversi di risolvere lo stesso flusso

Gli ingegneri si affidano spesso alla fluidodinamica computazionale (CFD) per esplorare i progetti prima della costruzione. Per i flussi gassosi ad alta velocità, un solver “basato sulla densità” è stata la scelta tradizionale, perché è tarato per gestire rapidi cambiamenti di pressione e densità attorno alle onde d’urto. Un solver “basato sulla pressione”, invece, è ampiamente usato per flussi a bassa velocità, quasi incomprimibili, ed è noto per essere più economico. Negli ultimi anni alcuni ricercatori hanno riportato che i solver basati sulla pressione possono funzionare bene anche a velocità maggiori, ma le loro prestazioni per getti fortemente sottosfusi – dove il gas esce da un ugello a pressione molto superiore all’aria circostante – sono rimaste incerte.

Costruire un getto in laboratorio e fotografare onde invisibili

Per testare equamente entrambi gli approcci, gli autori hanno realizzato un esperimento semplice ma accuratamente controllato. Aria compressa è stata rilasciata da un serbatoio cilindrico attraverso un piccolo ugello convergente in aria ferma ambiente. Regolando la pressione di alimentazione, hanno impostato quattro rapporti di pressione tra il serbatoio e l’atmosfera, da poco sopra la condizione di strozzamento, in cui il gas esce dall’ugello alla velocità del suono, fino a getti chiaramente sottosfusi con forti pattern di onde d’urto. Per rendere visibili queste strutture altrimenti invisibili, hanno utilizzato un apparato ottico Schlieren: un punto luce brillante, uno specchio concavo, una lamina da rasoio e una fotocamera digitale disposti in modo che piccole variazioni di densità dell’aria appaiano come bande chiare e scure. Queste immagini hanno fornito una mappa visiva delle celle d’urto che si formano a valle dell’ugello.

Mettere alla prova i modelli al computer

Parallelamente, lo stesso ugello e le stesse condizioni di flusso sono stati ricreati nel codice CFD ANSYS Fluent. I ricercatori hanno usato mesh identiche, condizioni al contorno e modelli di turbolenza per entrambi i solver, cambiando solo lo schema di soluzione di base. Hanno monitorato pressione, velocità e turbolenza lungo e intorno al getto, e verificato quanto rapidamente ogni solver convergesse a una soluzione stabile. Entrambi gli strumenti hanno riprodotto caratteristiche chiave: la formazione di un nucleo potenziale nella condizione di strozzamento, l’apparizione di celle d’urto a forma di diamante una volta che il getto è diventato sottosfuso, e il modo in cui queste celle si allungano e si intensificano con l’aumento del rapporto di pressione. Il primo shock principale, noto come disco di Mach, è apparso in posizione quasi identica in entrambe le simulazioni e nei dati sperimentali pubblicati, con solo poche percentuali di differenza nella velocità di picco.

Dove i solver concordano e dove differiscono

Lungo la linea centrale del getto, i due solver hanno prodotto oscillazioni quasi identiche di pressione e numero di Mach, rispecchiando il pattern ripetuto di compressione ed espansione. Le differenze sono emerse più a valle, dove il getto si era rallentato a velocità subsoniche e dominava la turbolenza. Lì, il solver basato sulla densità tendeva a prevedere una maggiore energia cinetica turbolenta rispetto al solver basato sulla pressione, in linea con i limiti noti degli schemi basati sulla densità a basse velocità. Nonostante ciò, entrambi i solver hanno fornito quasi lo stesso coefficiente di scarico, una misura di quanto efficacemente l’ugello trasferisce massa, con differenze inferiori allo 0,2 percento. In termini di sforzo computazionale, il solver basato sulla pressione ha raggiunto lo stesso obiettivo di convergenza stringente in circa il 22 percento di tempo CPU in meno.

Cosa significa per i progetti reali

Per gli ingegneri che devono simulare getti ad alta velocità provenienti da ugelli relativamente semplici, questo studio offre notizie rassicuranti. Un solver basato sulla pressione, pur essendo stato originariamente progettato per flussi più lenti e incomprimibili, può catturare i principali pattern di onde d’urto e le metriche di prestazione dei getti supersonici sottosfusi con un’accuratezza comparabile a un solver più specializzato basato sulla densità, impiegando meno tempo di calcolo. Gli autori avvertono che gli schemi basati sulla densità mantengono vantaggi in alcuni regimi, specialmente dove sono presenti onde d’urto molto forti e chimica complessa. Ma per molti problemi pratici di ugelli, l’approccio basato sulla pressione, a costo inferiore, può essere una scelta affidabile, aiutando i progettisti a esplorare più opzioni senza una bolletta di simulazione proibitiva.

Citazione: Alsaedi, S.S., Al-Sadawi, L.A., Al-Haddad, L.A. et al. Comparative analysis of compressible solver schemes for underexpanded jet aerodynamics with Schlieren validation. Sci Rep 16, 15724 (2026). https://doi.org/10.1038/s41598-026-44651-w

Parole chiave: getto sottosfuso, ugello supersonico, imaging Schlieren, solver CFD, onde d’urto