Flow and Heat Transfer Characteristics of Forced Convection Past an Isoflux Circular Cylinder in Galinstan for Reynolds Numbers up to 600

Sintesi Tecnica: Caratteristiche di Flusso e Trasferimento di Calore per Convezione Forzata su un Cilindro Circolare Isoflux in Galinstan per Numeri di Reynolds fino a 600

Definizione del Problema
Questo studio presenta un'indagine numerica sulla convezione termica forzata stazionaria da un cilindro circolare sottoposto a una condizione al contorno di flusso termico costante (isoflux). Il fluido in esame è il Galinstan, una lega metallica liquida caratterizzata da un basso numero di Prandtl ($Pr = 0,025$) e una conducibilità termica eccezionalmente elevata. La ricerca si concentra sul comportamento idrodinamico e termico del flusso su un intervallo di numeri di Reynolds ($Re$) compreso tra $1 \le Re \le 600$. Sebbene la letteratura precedente abbia ampiamente trattato cilindri isotermi e fluidi con numeri di Prandtl moderati o elevati, questo lavoro affronta le sfide specifiche della simulazione di metalli liquidi a basso numero di Prandtl in condizioni isoflux, un regime rilevante per i sistemi avanzati di gestione termica, come il raffreddamento elettronico e i sistemi energetici basati su metalli liquidi.

Metodologia
Le equazioni governanti — che comprendono le equazioni di portata, vorticità ed energia in coordinate cilindriche — vengono risolte utilizzando uno schema di differenze finite compatte del quarto ordine (FOCS–CC). Per garantire stabilità e convergenza, questa discretizzazione spaziale è accoppiata a una tecnica di iterazione in pseudo-tempo (PTI).

Caratteristiche metodologiche chiave:

• Trasformazione delle Coordinate: Il dominio fisico è mappato in un dominio computazionale utilizzando la trasformazione $r = e^{\pi\lambda}$ e $\theta = \pi\eta$, consentendo una risoluzione della griglia efficiente vicino alla superficie del cilindro.
• Accuratezza di Ordine Elevato: Lo schema utilizza formulazioni compatte del quarto ordine per le derivate spaziali. Una caratteristica degna di nota evidenziata dagli autori è la facilità di implementazione; il metodo estende gli schemi esistenti del secondo ordine incorporando coefficienti di correzione specifici ($\alpha, b_1, \Gamma_1$, ecc.) senza alterare la struttura numerica sottostante.
• Algoritmo di Soluzione: Il sistema di equazioni risultante viene integrato utilizzando un approccio fattorizzato per ridurre le matrici a larga banda in sistemi tridiagonali, che vengono risolti efficientemente tramite l'algoritmo della matrice tridiagonale (TDMA).
• Validazione: Il codice numerico è validato attraverso studi di indipendenza dalla griglia e confrontando il numero di Nusselt medio ($Nu$) calcolato e il coefficiente di resistenza totale ($C_D$) con i risultati numerici consolidati di Dennis e Chang (1970), Fornberg (1980), Ahmad (1996) e Morgan (1975).

Risultati Chiave e Osservazioni
Lo studio analizza sistematicamente il campo di flusso e il campo termico attraverso l'intervallo specificato di numeri di Reynolds:

1. Struttura del Flusso: A bassi numeri di Reynolds ($Re=1$), il flusso rimane completamente attaccato alla superficie del cilindro senza separazione, dominato dalle forze viscose. All'aumentare di $Re$ fino a 200 e 600, si verifica la separazione del flusso, che forma una regione di scia dietro il cilindro. La dimensione di questa scia aumenta con l'aumentare di $Re$, indicando una crescente dominanza delle forze d'inerzia.
2. Campo Termico: Per $Re=1$, le isoterme formano anelli quasi concentrici, segnalando che il trasferimento di calore è governato principalmente dalla diffusione termica. All'aumentare di $Re$, il trasporto convettivo si intensifica, causando un impacchettamento delle isoterme nel punto di ristagno a monte e un loro allungamento a valle, formando una piuma termica. Nonostante questo allungamento convettivo, il campo termico rimane relativamente diffuso rispetto ai fluidi convenzionali a causa dell'elevata diffusività termica del Galinstan.
3. Metriche di Trasferimento di Calore: Il numero di Nusselt locale presenta una variazione angolare minima a bassi $Re$, ma mostra una variazione significativa ad alti $Re$ a causa dell'influenza della scia separata. Il numero di Nusselt medio aumenta monotonicamente con $Re$ in tutto l'intervallo investigato, sebbene il tasso di incremento diminuisca ai numeri di Reynolds più elevati.
4. Correlazione Empirica: Sulla base dei dati numerici, viene proposta una nuova correlazione empirica per il numero di Nusselt medio per il Galinstan nell'intervallo $1 \le Re \le 600$:
   $$Nu = \alpha + \beta\sqrt{Re} + \gamma Re$$
   dove $\alpha = 0,26133$, $\beta = 0,13302$ e $\gamma = -5,42758 \times 10^{-4}$. Questa correlazione dimostra un eccellente accordo con i dati numerici, fornendo un coefficiente di determinazione ($R^2$) di 0,99939.

Significato e Rivendicazioni
Gli autori affermano che lo schema FOCS–CC proposto, integrato con la tecnica PTI, fornisce un quadro numerico robusto e accurato per simulare la convezione forzata attorno a cilindri circolari in metalli liquidi a basso numero di Prandtl. Lo studio stabilisce che:

• Le soluzioni stazionarie per le equazioni di Navier–Stokes incomprimibili possono essere computate con successo fino a $Re = 600$ per questa configurazione, estendendosi oltre l'insorgenza dell'instabilità del flusso instazionario.
• L'interazione tra il trasporto convettivo potenziato e la diffusione termica dominante è una caratteristica definente del Galinstan, che distingue il suo comportamento termico da quello dei fluidi con numeri di Prandtl moderati o elevati.
• L'elevata accuratezza dello schema compatto del quarto ordine consente la previsione affidabile sia delle caratteristiche idrodinamiche (resistenza) che termiche (numero di Nusselt), come confermato dalla validazione rispetto alla letteratura esistente.

L'articolo conclude che la correlazione sviluppata e il metodo numerico offrono strumenti preziosi per l'analisi del trasferimento di calore nelle applicazioni con metalli liquidi, evidenziando in particolare i rendimenti decrescenti nel potenziamento del trasferimento di calore convettivo ad alti numeri di Reynolds a causa della persistente influenza dell'elevata diffusività termica.