Distributed Roughness-Induced Transition on a Blunt Body at Mach 6: a Numerical Investigation

Questo studio presenta la prima simulazione numerica diretta di rugosità distribuita su un cilindro smussato a Mach 6, rivelando che la disposizione della rugosità determina il meccanismo di transizione promuovendo o i modi a striscia sinuosa o le onde T-S 2D, con queste ultime sostenute unicamente da un ciclo di feedback acustico interno che elimina la necessità di un forzamento esterno.

L'essenziale

Immaginate un veicolo spaziale che sfreccia attraverso l'atmosfera a sei volte la velocità del suono. Per sopravvivere al calore intenso, la sua superficie è ricoperta da un materiale speciale che brucia lentamente (ablazione) per proteggere la nave. Tuttavia, mentre questo materiale brucia, non lascia una superficie perfettamente liscia, ma lascia dietro di sé una texture irregolare e rugosa, simile alla carta vetrata.

Questo articolo è una simulazione al computer ad alta velocità che pone una domanda semplice ma critica: in che modo queste piccole protuberanze sulla superficie trasformano un flusso d'aria liscio e ordinato in un flusso d'aria turbolento e caotico?

Ecco la ripartizione delle loro scoperte utilizzando analogie quotidiane:

1. L'allestimento: Il cilindro "carta vetrata"

I ricercatori hanno costruito un modello digitale di un cilindro smussato (come il muso di un razzo) che vola a Mach 6. Invece di una superficie liscia, lo hanno ricoperto di piccole "protuberanze" artificiali (rugosità) per imitare la texture sabbiosa lasciata dal materiale che brucia.

Hanno testato tre modi diversi per disporre queste protuberanze:

• Allineate: Come soldati in file e colonne perfette.
• Sfalsate: Come un muro di mattoni, dove le protuberanze in una fila sono spostate rispetto alla fila dietro di essa.
• Casuali: Come ciottoli sparsi su un marciapiede senza alcun schema.

2. La vecchia teoria vs. La nuova scoperta

Per molto tempo, gli scienziati hanno pensato che la transizione alla turbolenza fosse causata da un "accumulo lento" di energia, simile a come un'altalena guadagna altezza se la spingi nel modo giusto nel tempo. Questo è chiamato "crescita transitoria".

La scoperta del documento:
La simulazione ha dimostrato che questa teoria dell' "accumulo lento" non spiega realmente ciò che sta accadendo qui. Le protuberanze sulla superficie non si sono limitate ad amplificare lentamente l'energia; hanno agito come destabilizzatori. Hanno preso il flusso d'aria e lo hanno reso immediatamente instabile, trasformandolo in un tipo specifico di onda che cresce molto velocemente.

Pensatela in questo modo: la vecchia teoria pensava che le protuberanze stessero gentilmente dando un piccolo colpetto a un domino per farlo cadere. La nuova scoperta mostra che le protuberanze stanno in realtà dando un calcio al domino, facendolo schiantare contro quello successivo immediatamente.

3. I due tipi di "onde"

A seconda di come erano disposte le protuberanze, il flusso d'aria reagiva in due modi diversi:

• Il "Serpente" (Modo Sinuoso): Quando le protuberanze erano allineate (file perfette), il flusso d'aria iniziava a ondeggiare lateralmente come un serpente. È un'oscillazione molto specifica e organizzata.
• L' "Onda Piatta" (Onde Tollmien-Schlichting o T-S): Quando le protuberanze erano sfalsate o casuali, il flusso d'aria iniziava a incresparsi su e giù in un modello d'onda 2D piatto. Questo è un tipo classico di onda che si trova solitamente in flussi d'aria molto più lenti e a bassa velocità, il che è stato sorprendente da trovare in questo ambiente ad alta velocità.

L'intuizione chiave: La disposizione delle protuberanze dettava quale "danza" avrebbe eseguito l'aria. La danza del "serpente" accadeva con le protuberanze allineate, mentre la danza della "onda piatta" accadeva con le altre.

4. Il finale a "Forcella"

Una volta che queste onde sono diventate abbastanza forti, hanno innescato la fase finale dello schianto. Le costanti "scie" d'aria create dalle protuberanze (che sono come lunghi nastri invisibili di aria lenta) improvvisamente si sono attorcigliate e sono scattate in vortici a forcina.

Immaginate un elastico che viene teso. Improvvisamente, si attorciglia e forma un anello che sembra una forcina. Questi anelli sono la nascita della turbolenza. Una volta formate queste forcine, l'aria fluida si rompe completamente nel caos, e il calore sulla superficie del veicolo spaziale aumenta drasticamente.

5. La sorpresa della "Camera dell'Eco"

Una delle scoperte più affascinanti è stata come la turbolenza sia iniziata in primo luogo per i casi sfalsati e casuali.

Di solito, gli scienziati pensano che sia necessario un "impulso" esterno (come una raffica di vento o una vibrazione) per avviare queste onde. Ma la simulazione ha mostrato qualcosa di auto-sostenuto:

1. La turbolenza inizia in un punto sul cilindro.
2. Poiché l'aria dietro l'onda d'urto si muove più lentamente del suono (subsonica), il rumore di questa turbolenza viaggia all'indietro verso monte come un'eco.
3. Questa "eco" colpisce la parte liscia della superficie davanti alla turbolenza ed eccita l'aria lì, causando l'inizio di nuova turbolenza.
4. Questo crea un ciclo di feedback: la turbolenza produce rumore, il rumore viaggia all'indietro, e il rumore crea altra turbolenza.

È come un microfono che cattura l'uscita del proprio altoparlante e crea un feedback stridente, ma in questo caso è l'aria che diventa turbolenta.

Riassunto

Questo articolo ha utilizzato un supercomputer per osservare il flusso d'aria sopra un cilindro irregolare ad alta velocità. Hanno scoperto che:

• Il modello delle protuberanze decide esattamente come l'aria diventa turbolenta.
• La vecchia idea di "lento accumulo di energia" non è il colpevole principale; invece, le protuberanze destabilizzano direttamente onde specifiche.
• Queste onde crescono fino a torcersi in forme a "forcina", causando il caos dell'aria.
• In alcuni casi, la turbolenza crea la propria "eco" che viaggia all'indietro per avviare il processo ancora una volta, senza bisogno di alcun aiuto esterno.

Questo aiuta gli ingegneri a capire che le minuscole, casuali protuberanze lasciate dagli scudi termici in combustione non sono solo piccole imperfezioni; sono gli architetti primari di come e quando la superficie di un veicolo spaziale diventa pericolosamente calda.

Sommario tecnico

Sintesi Tecnica: Transizione Indotta da Rugosità Distribuita su un Corpo Ottuso a Mach 6

Problematica
La previsione della transizione laminare-turbolenta (LTT) su geometrie ottuse ad alta velocità rimane una sfida critica per la progettazione ingegneristica, in particolare per i veicoli di rientro atmosferico dove l'ablazione superficiale crea rugosità distribuite. Sebbene esistano correlazioni empiriche (ad esempio, Hollis 2017, 2019, 2025), esse mancano di un fondamento fisico chiaro. Sforzi teorici precedenti, come la correlazione di Reshotko & Tumin (2004), si basano su meccanismi di crescita transitoria ma presentano limitazioni: assumono una scalatura lineare delle perturbazioni con l'altezza della rugosità (valida solo per piccoli $k/\delta$), richiedono costanti di adattamento ai dati sperimentali e faticano a spiegare la transizione nei casi in cui l'altezza della rugosità eccede lo spessore dello strato limite ($k/\delta > 1$). Inoltre, le precedenti indagini numeriche spesso simulavano patch di rugosità piuttosto che rugosità completamente distribuite, o non riuscivano a osservare l'amplificazione dell'instabilità modale, lasciando poco chiari i meccanismi specifici che guidano la transizione nei flussi di corpi ottusi ipersonici con rugosità completamente distribuita.

Metodologia
Questo studio presenta la prima Simulazione Numerica Diretta (DNS) di un cilindro ottuso in cross-flow a Mach 6 con elementi di rugosità distribuiti su tutta la superficie. Le simulazioni utilizzano il solver CHAMPS per integrare le equazioni di Navier-Stokes comprimibili (CNSE) per un gas ideale. La geometria è un cilindro di 0,1524 m di diametro, modellato sulla base delle condizioni sperimentali di Hollis (2017) (Run 38, Test 6975), con un numero di Mach in corrente libera di 5,98.

La rugosità è modellata come una distribuzione sinusoidale con un'ampiezza fissa ($A \approx 35\%$ dello spessore dello strato limite) e una lunghezza d'onda corrispondente alla rugosità nominale a grana di sabbia. Sono state investigate tre diverse configurazioni di sfasamento trasversale (spanwise):

1. Staggered (Sfalsata): Le file sono sfasate di $180^\circ$.
2. Aligned (Allineata): Le file sono in fase (scostamento di $0^\circ$).
3. Random (Casuale): Le file hanno scostamenti di fase distribuiti casualmente.

Per analizzare il processo di transizione, gli autori hanno impiegato:

• DNS: Per catturare l'intera evoluzione non lineare del flusso.
• Equazioni di Disturbo Lineare (LDE): Risolte su profili mediati nel tempo derivati dalla DNS per isolare i meccanismi di instabilità mediante forcing a impulso.
• Teoria della Stabilità Lineare (LST): Applicata a profili mediati trasversalmente per il confronto.
• Analisi Spettrale: Trasformate di Fourier veloci (FFT) dei campi di velocità e pressione per identificare le frequenze e i numeri d'onda dominanti.

Risultati Chiave

• Flusso Medio e Streak: Tutte le configurazioni hanno generato "streak" di velocità longitudinale tramite un meccanismo di lift-up guidato da coppie di vortici contro-rotanti (CVP). La configurazione allineata ha prodotto gli streak stazionari più forti a causa di uno sfasamento coerente che eleva il fluido a basso flusso di massa più in alto nello strato limite. Le configurazioni sfalsata e casuale hanno prodotto pattern di streak più deboli e complessi. Contrariamente alle previsioni della teoria della crescita transitoria ottimale, gli ampi degli streak hanno mostrato una crescita relativa bassa (meno di un fattore 2) e sono rimasti non lineari in tutto il dominio, sfidando l'utilità della crescita transitoria lineare come strumento predittivo per questo regime.

• Meccanismi di Instabilità: Il processo di transizione è stato guidato da instabilità convezionive, ma il tipo di modo specifico dipendeva dalla fase della rugosità:

  • Caso Allineato: Dominato da un modo di streak sinuoso fondamentale. Le onde di Tollmien-Schlichting (T-S) sono state debolmente destabilizzate.
  • Casi Sfalsato e Casuale: Dominati da onde T-S 2D (modo varicoso). Queste configurazioni hanno mostrato una significativa destabilizzazione delle onde T-S, portando a una rapida crescita esponenziale.
  • Breakdown (Rottura): In tutti i casi, una volta che le onde di instabilità dominanti hanno raggiunto un'ampiezza sufficiente, hanno innescato la formazione di vortici a forcina (hairpin vortices) negli streak stazionari, portando alla rottura e alla turbolenza.

• Posizioni di Transizione: La DNS ha previsto posizioni di transizione ($\Theta_{tr}$) significativamente più a valle (circa $11^\circ$–$17^\circ$) rispetto alle correlazioni empiriche allo stato dell'arte (che prevedevano $6^\circ$–$13^\circ$). Il caso allineato è transitato per primo tra i casi di rugosità a causa del forcing richiesto per innescare la LTT, mentre i casi sfalsato e casuale sono transitati più tardi.

• Meccanismo di Feedback: Un risultato innovativo nei casi sfalsato e casuale è stata l'identificazione di un ciclo di feedback. La turbolenza generata nella regione subsonica del cilindro ha disperso onde acustiche che si sono propagate verso monte. Queste lente onde acustiche hanno eccitato le onde T-S sul lato opposto del cilindro (attraverso la linea di ristagno) al loro punto neutro. Questo meccanismo, che non richiede un forcing esterno, sostiene il processo di transizione.

Significatività e Rivendicazioni
Gli autori affermano che questo lavoro fornisce una nuova comprensione fisica della transizione indotta da rugosità su corpi ottusi ipersonici, sfidando la prevalente dipendenza dalla teoria della crescita transitoria. I contributi chiave sono:

1. Rifiuto della Crescita Transitoria come Driver Primario: Lo studio dimostra che per rugosità completamente distribuite con $k/\delta \approx 0,35$, lo sviluppo non lineare degli streak e i bassi rapporti di amplificazione osservati rendono la teoria della crescita transitoria ottimale un predittore insufficiente del reale processo di transizione.
2. Destabilizzazione delle Instabilità Modali: Il documento stabilisce che la rugosità distribuita può destabilizzare efficacemente le instabilità modali (specificamente le onde T-S) anche su corpi ottusi dove il flusso a parete liscia è tipicamente stabile per quanto riguarda gli autovettori. La configurazione della rugosità determina se la transizione è guidata dai modi di streak o dalle onde T-S.
3. Scoperta di un Ciclo di Feedback Acustico: L'identificazione di un meccanismo di feedback auto-sostenuto, in cui l'acustica generata dalla turbolenza eccita le onde T-S sul lato opposto del corpo, offre una potenziale spiegazione per la transizione nelle regioni subsoniche del flusso ipersonico, un fenomeno precedentemente non spiegato dai meccanismi lineari.
4. Sensibilità alla Fase: I risultati evidenziano l'estrema sensibilità della LTT alla fase relativa degli elementi di rugosità, suggerendo che topologie semplici possono delimitare il comportamento di pattern di rugosità più complessi e casuali.

Gli autori concludono che, sebbene il proposto meccanismo di feedback acustico richieda un'ulteriore prova rigorosa, l'esposizione di questi meccanismi di instabilità tramite DNS ad alta fedeltà suggerisce un cambiamento fondamentale nella comprensione dei disturbi responsabili della LTT sui corpi ottusi ipersonici.