Numerical Study on the Effects of Wall Wettability on Heat Transfer and Flow Characteristics during Flow Boiling in Microchannels

Sintesi Tecnica: Studio Numerico sugli Effetti della Bagnabilità della Parete sul Trasferimento di Calore e sulle Caratteristiche del Flusso durante l'Ebollizione in Flusso in Microcanali

Definizione del Problema
La miniaturizzazione dei chip ad alte prestazioni ha portato a un aumento significativo della densità di flusso termico locale, che spesso supera i 1000 W/cm², valori che i metodi di raffreddamento tradizionali faticano a dissipare. L'ebollizione in flusso in microcanali rappresenta una soluzione promettente grazie agli elevati coefficienti di scambio termico e alle capacità di assorbimento del calore latente. Tuttavia, le prestazioni di questi sistemi sono fortemente influenzate dalla dinamica delle bolle e dall'evoluzione dei regimi di flusso. Un fattore critico, sebbene complesso, in questo processo è la bagnabilità della parete. Sebbene la letteratura esistente abbia esplorato superfici idrofile, idrofobe e a bagnabilità mista, i meccanismi accoppiati che governano il modo in cui le condizioni estreme di bagnabilità influenzano l'interazione tra il comportamento delle bolle, l'efficienza del trasferimento di calore, l'instabilità del flusso e l'irreversibilità termodinamica (generazione di entropia) rimangono insufficientemente compresi. Nello specifico, è necessario quantificare i compromessi tra il miglioramento del trasferimento di calore e la stabilità del flusso sotto diversi regimi di bagnabilità.

Metodologia
Questo studio impiega un modello numerico tridimensionale basato sul metodo Volume of Fluid (VOF) per simulare l'ebollizione in flusso in un microcanale scanalato. La simulazione utilizza il modello Continuum Surface Force (CSF) per catturare gli effetti della tensione superficiale e un modello di evaporazione-condensazione (basato su Lee et al.) per gestire il cambiamento di fase.

• Modello Fisico: La geometria consiste in un microcanale con scanalature triangolari laterali, progettato per studiare la nucleazione e la crescita delle bolle. Il fluido di lavoro è acqua deionizzata e il dominio solido è il silicio.
• Condizioni di Bagnabilità: Vengono confrontate due condizioni superficiali estreme: una parete idrofila (angolo di contatto $\theta = 30^\circ$) e una parete idrofoba ($\theta = 150^\circ$).
• Approccio Numerico: Le equazioni governanti (conservazione di massa, quantità di moto ed energia) sono risolte utilizzando il metodo dei volumi finiti con l'algoritmo PISO per l'accoppiamento pressione-velocità. Uno studio di indipendenza dalla griglia ha confermato che una mesh di circa 275.000 elementi fornisce un'accuratezza sufficiente. Il modello è stato validato rispetto a dati sperimentali, mostrando una discrepanza massima inferiore all'8% nell'innalzamento della temperatura della base del chip.
• Metriche di Analisi: Lo studio valuta la dinamica delle bolle, i regimi di flusso, le fluttuazioni della temperatura in uscita, il numero di Nusselt (Nu), la caduta di pressione e la generazione di entropia (sia componenti termiche che frizionali).

Contributi Chiave e Risultati

1. Dinamica delle Bolle ed Evoluzione del Regime di Flusso:

   • Superfici Idrofile ($\theta = 30^\circ$): Forti interazioni liquido-solido causano l'adesione delle bolle alla parete e la loro coalescenza, formando strutture a film liquido. Ciò promuove la transizione verso un flusso anulare ma ostacola il distacco delle bolle.
   • Superfici Idrofobe ($\theta = 150^\circ$): Le bolle si staccano rapidamente e vengono trasportate nel flusso principale. Ciò accelera l'accumulo di vapore e porta a una transizione precoce dal flusso a bolle (bubbly flow) al flusso a slug (slug flow) e al flusso anulare.

2. Prestazioni Termiche:

   • Efficienza del Trasferimento di Calore: La superficie idrofoba dimostra prestazioni di trasferimento termico superiori. Il rapido distacco delle bolle facilita la rimozione del calore latente, mantenendo temperature della parete più basse e numeri di Nusselt più elevati. Al contrario, la superficie idrofila soffre di un accumulo di calore dovuto alle bolle aderenti, che risulta in temperature della parete più elevate e in un declino continuo del numero di Nusselt.
   • Stabilità della Temperatura: Sebbene le superfici idrofobe migliorino il trasferimento di calore medio, esse presentano un'instabilità della temperatura in uscita inferiore rispetto alle superfici idrofile. La superficie idrofila mostra fluttuazioni di temperatura più ampie a causa dell'ostacolo alla dissipazione del calore causato dai persistenti film liquidi.

3. Caratteristiche del Flusso e Instabilità:

   • Caduta di Pressione: Il canale idrofobo sperimenta cadute di pressione significativamente più elevate e oscillazioni di pressione più pronunciate. Ciò è attribuito alla maggiore frazione volumetrica di vapore e allo scarico periodico delle bolle che altera la struttura locale del flusso.
   • Campi di Velocità: Il canale idrofobo presenta disturbi del flusso locale e gradienti di velocità più intensi a causa della dinamica attiva delle bolle, mentre il canale idrofilo mostra variazioni di velocità più fluide e inferiori.

4. Analisi della Generazione di Entropia:

   • Termica vs Frizionale: Lo studio distingue tra la generazione di entropia termica (guidata dai gradienti di temperatura) e la generazione di entropia frizionale (guidata dalla dissipazione viscosa).
   • Meccanismo Dominante: La generazione di entropia termica è il contributore principale alla generazione totale di entropia. La superficie idrofila genera un'entropia termica significativamente più alta a causa dei grandi gradienti di temperatura causati dal deterioramento del trasferimento di calore.
   • Compromesso (Trade-off): Sebbene la superficie idrofoba generi un'entropia frizionale maggiore a causa degli intensi gradienti di velocità, la sua generazione totale di entropia rimane inferiore a quella della superficie idrofila grazie alla sostanziale riduzione dell'irreversibilità termica.

Significatività e Rivendicazioni
Gli autori affermano che questo studio fornisce una guida teorica preziosa per la progettazione della bagnabilità e l'ottimizzazione dei sistemi termoidraulici in microcanali dissipatori di calore sotto elevati flussi termici. I risultati evidenziano un compromesso fondamentale: mentre le superfici idrofobe migliorano significativamente il trasferimento di calore per ebollizione e riducono l'irreversibilità termica, esse introducono una maggiore resistenza al flusso e fluttuazioni di pressione. Al contrario, le superfici idrofile, pur offrendo cadute di pressione più stabili, soffrono di un grave deterioramento del trasferimento di calore e di maggiori perdite termodinamiche.

Lo studio conclude che l'ottimizzazione della bagnabilità richiede un bilanciamento tra il potenziamento del trasferimento di calore per cambio di fase e la penalità derivante dall'aumento dell'instabilità del flusso. I risultati suggeriscono che, per applicazioni che danno priorità alle prestazioni termiche e alla minimizzazione dell'entropia, le superfici idrofobe sono vantaggiose, a condizione che il sistema possa gestire la relativa caduta di pressione e le oscillazioni del flusso. Lo studio non propone specifiche applicazioni future o validazioni sperimentali oltre l'ambito della presente investigazione numerica, ma sottolinea la necessità di comprendere tali meccanismi accoppiati per un design efficace della gestione termica.