Numerical homogenization of a heat transfer problem with phase transitions in heterogeneous frozen soils

Nel vasto e gelido nord, il terreno non è semplicemente terra; è una fondazione strutturale che sostiene le più grandi città, oleodotti e strade del mondo. Per secoli, questo permafrost è rimasto solido, il suo stato congelato che fungeva da letto roccioso rigido per le costruzioni umane. Tuttavia, con il riscaldamento del clima, questo terreno ghiacciato sta iniziando a scongelarsi. Quando il ghiaccio all'interno del suolo si scioglie, il terreno perde la sua forza, portando a cedimenti irregolari e alla fessurazione degli edifici. Prevedere esattamente quanto sarà profondo questo scongelamento è una questione di sicurezza e ingegneria, ma il suolo stesso presenta un enorme complicazione. Il vero suolo non è un blocco di ghiaccio uniforme; è un miscuglio caotico di strati, lenti e sacche di materiali diversi, con proprietà che cambiano su distanze piccole anche di pochi centimetri. Per simulare lo scioglimento di un terreno così complesso con perfetta accuratezza, un computer dovrebbe tracciare ogni singolo strato, un compito che richiede una potenza di calcolo superiore a quella attualmente disponibile per previsioni su larga scala e a lungo termine.

Per risolvere questo problema, i ricercatori ricorrono a una tecnica chiamata omogeneizzazione numerica. Immaginate di cercare di prevedere il flusso dell'acqua attraverso una spugna. Invece di modellare ogni minuscolo foro e fibra, trattate la spugna come un unico materiale con una capacità "media" di lasciar passare l'acqua. Questa è l'essenza dell'omogeneizzazione: sostituire una realtà disordinata e dettagliata con un modello più semplice e mediato, abbastanza veloce da essere eseguito da un computer. Tuttavia, quando il terreno si sta sciogliendo, questa media diventa una trappola. Il confine tra il terreno ghiacciato e quello scongelato — il "fronte" — non si ferma nettamente al bordo del blocco di calcolo del computer. Esso taglia proprio attraverso il mezzo. Questo crea un dilemma: se un blocco di terreno è metà congelato e metà scongelato, come dovrebbe decidere il computer le sue proprietà? Dovrebbe comportarsi più come il ghiaccio, più come l'acqua, o qualcosa nel mezzo? La scelta fatta qui determina se il modello prevede che un edificio affondi o rimanga stabile.

In uno studio recente, Timur Zakharov della North-Eastern Federal University ha affrontato questo specifico dilemma testando tre diversi modi per gestire questi blocchi di terreno "misti". La ricerca si è concentrata su un edificio riscaldato situato su un terreno ghiacciato, uno scenario in cui il calore proveniente dalla struttura scioglie lentamente il suolo sottostante. Il team ha creato un gemello digitale di questo ambiente, utilizzando un pattern realistico e casuale di conduttività del suolo per imitare il caos naturale del terreno. Hanno poi eseguito la simulazione su una griglia grossolana, dove ogni blocco di terreno era molto più grande degli strati effettivi di terra e ghiaccio, costringendo il fronte di cambiamento di fase a tagliare attraverso il centro di questi blocchi. Hanno confrontato tre metodi distinti per decidere cosa accade all'interno di quei blocchi tagliati.

Il primo metodo, noto come approccio a tabella (table approach), pre-calcola un elenco di proprietà per ogni possibile stato di un blocco di suolo e cerca la risposta in base alla temperatura media del blocco. Il secondo metodo, uno schema nodale (nodal scheme), decide lo stato del suolo agli angoli dei blocchi e media quelle decisioni. Il terzo metodo, chiamato h-smoothing, calcola la frazione esatta del blolo che è sciolta basandosi su un gradiente di temperatura lineare, misurando essenzialmente la fetta del blocco che si è trasformata in acqua. I ricercatori hanno eseguito queste simulazioni sia in due che in tre dimensioni, confrontando i risultati con una soluzione di "riferimento" altamente dettagliata che tracciava ogni strato di suolo.

I risultati hanno rivelato un compromesso sorprendente che sfida il modo in cui gli ingegneri potrebbero scegliere i loro strumenti. Quando l'obiettivo è prevedere la temperatura del suolo, il metodo h-smoothing si è dimostrato il più accurato. Ha tracciato il flusso di calore con la massima precisione, specialmente sulle griglie grossolane dove gli altri metodi faticavano. Tuttavia, quando i ricercatori hanno osservato la posizione del fronte di fusione — il confine effettivo tra terreno ghiacciato e scongelato — la storia è cambiata completamente. Il metodo h-smoothing, pur essendo buono per la temperatura, ha costantemente sovrastimato quanto sarebbe profondo lo scongelamento. Ha previsto che il terreno si sarebbe sciolto più in profondità di quanto non sia in realtà.

La ragione di questa discrepanza risiede nelle leggi della conservazione dell'energia. Il metodo h-smoothing, nel suo sforzo di essere preciso sulla forma della zona di fusione, ha involontariamente aggiunto una piccola quantità di calore extra al sistema ad ogni passaggio temporale. Nel corso di un anno, questo surplus minimo si è accumulato, causando lo scioglimento del terreno simulato più di quanto non farebbe il terreno reale. Al contrario, gli altri due metodi, in particolare l'approccio a tabella e lo schema nodale, erano "conservativi". Non hanno aggiunto né perso energia e, di conseguenza, hanno previsto la profondità dello scongelamento con molta più precisione, anche se le loro letture della temperatura erano leggermente meno precise. Lo studio ha dimostrato che se un ricercatore ha bisogno di conoscere la temperatura esatta del suolo, l'h-smoothing è la scelta superiore. Ma se l'obiettivo è sapere quanto affonderà la fondazione a causa dello scongelamento, i metodi conservativi sono molto più affidabili.

La ricerca ha anche esplorato cosa accade quando la temperatura alla quale il suolo si scioglie non è la stessa ovunque, ma varia tra i diversi strati. In questo scenario più complesso, il metodo basato sulla tabella, che utilizza una ricerca dettagliata delle proprietà, ha superato i metodi a due tensori più semplici nel prevedere la posizione del fronte di scongelamento. Ciò suggerisce che per condizioni del terreno altamente variabili, è necessaria una mappa più dettagliata e pre-calcolata degli stati del suolo. In definitiva, lo studio fornisce una guida chiara per ingegneri e scienziati: non esiste un unico metodo "migliore". La scelta dipende interamente dalla domanda che viene posta. Se la priorità è l'accuratezza termica, un metodo vince; se la priorità è il confine fisico dello scongelamento, un altro prende il comando. Chiarificando queste differenze, il lavoro aiuta a garantire che i modelli utilizzati per proteggere gli edifici nell'Artico non siano solo veloci, ma anche affidabili.