Un modello combinato di resistenza-capacità termica ed elementi finiti per simulazioni molto veloci e accurate a breve e medio termine di scambiatori di calore in pozzi geotermici a tubo singolo

Informazioni Fondamentali

• ID Articolo: 2510.14421
• Titolo: Un modello combinato di resistenza-capacità termica ed elementi finiti per simulazioni molto veloci e accurate a breve e medio termine di scambiatori di calore in pozzi geotermici a tubo singolo
• Autori: Enzo Zanchini, Francesco Zanchini, Claudia Naldi (Università di Bologna)
• Classificazione: physics.app-ph
• Data di Pubblicazione: 2024
• Link Articolo: https://arxiv.org/abs/2510.14421

Riassunto

La progettazione accurata dei sistemi di pompe di calore geotermiche richiede la comprensione della temperatura del fluido in uscita degli scambiatori di calore in pozzi geotermici (BHE) a breve e lungo termine. Questo articolo si concentra sull'analisi a breve e medio termine. In questo intervallo temporale, è possibile utilizzare simulazioni 3D agli elementi finiti o modelli di resistenza-capacità termica (TRCM). Il primo offre elevata precisione ma richiede tempi di calcolo lunghi, il secondo è veloce ma ha precisione limitata a causa delle ipotesi semplificate. Questo articolo propone un nuovo metodo che combina la velocità del TRCM con la precisione della simulazione agli elementi finiti per simulazioni a breve e medio termine di scambiatori di calore in pozzi geotermici a tubo singolo. Il metodo utilizza il TRCM per stimare la risposta termica del BHE, quindi corregge mediante interpolazione utilizzando un dataset generato da 54 simulazioni agli elementi finiti di diverse configurazioni. Il modello è implementato come programma C++ ed è disponibile nel database open-source dell'Università di Bologna, fornendo in due secondi l'evoluzione temporale delle temperature di ingresso, uscita e media del fluido, della temperatura media della superficie BHE, della resistenza termica 3D e della resistenza termica efficace del pozzo.

Contesto di Ricerca e Motivazione

Descrizione del Problema

I sistemi di pompe di calore geotermiche, come tecnologia efficiente di riscaldamento e raffreddamento degli edifici, richiedono una previsione accurata delle caratteristiche di risposta termica degli scambiatori di calore in pozzi geotermici. Le sfide principali sono:

1. Contraddizione tra precisione computazionale ed efficienza: le simulazioni 3D agli elementi finiti hanno elevata precisione ma richiedono ore; il TRCM è veloce ma ha precisione limitata
2. Esigenza di simulazioni a breve e medio termine: i metodi esistenti presentano problemi di precisione o efficienza nell'intervallo temporale da 10⁻²·⁵ a 10³ ore
3. Esigenze di applicazioni ingegneristiche pratiche: è necessario uno strumento di simulazione sia veloce che accurato per la progettazione e l'ottimizzazione del sistema

Limitazioni dei Metodi Esistenti

1. Metodi di soluzione analitica: sostituiscono il fluido con un solido riscaldante o una superficie, non possono considerare gli effetti dello scambio termico tra i flussi ascendenti e discendenti
2. Simulazione numerica 3D: elevata precisione ma ogni simulazione richiede ore, non adatta alla progettazione rapida in ingegneria
3. Modelli TRCM esistenti: adottano l'ipotesi di temperatura superficiale uniforme orizzontale (HUST), portando a sottostime della risposta termica a medio termine

Motivazione della Ricerca

Sviluppare un metodo di simulazione ibrida che combini la velocità del TRCM con la precisione degli elementi finiti, soddisfacendo le esigenze pratiche dell'ingegneria.

Contributi Principali

1. Identificazione e quantificazione degli effetti dell'ipotesi HUST: attraverso il confronto tra simulazioni agli elementi finiti con e senza strato ad alta conduttività termica, si dimostra che l'ipotesi HUST sottostima T_fm, T_in, T_out, R_b3D e R_beff
2. Sviluppo di un modello TRCM migliorato: adotta una rete termica stratificata che include nodi di fluido, tubo, malta e strato, migliorando la precisione della simulazione
3. Creazione di un database di coefficienti di correzione: attraverso 54 simulazioni agli elementi finiti di diverse configurazioni, costruisce un dataset di coefficienti di correzione dipendenti dal tempo
4. Implementazione di un metodo di simulazione ibrida efficiente: combinando il calcolo di base TRCM e la correzione mediante interpolazione, raggiunge la precisione della simulazione agli elementi finiti in 2 secondi
5. Fornitura di uno strumento di simulazione open-source: il programma C++ può essere utilizzato per simulazioni BHE a breve e medio termine veloci e accurate

Dettagli del Metodo

Definizione del Compito

Input: parametri geometrici BHE (lunghezza L, raggio r_b, distanza tra tubi s, ecc.), parametri termofisici, portata, carico termico Output: evoluzione temporale della temperatura di ingresso e uscita (T_in, T_out), temperatura media del fluido (T_fm), temperatura della superficie BHE (T_b), resistenza termica 3D (R_b3D), resistenza termica efficace (R_beff) Intervallo temporale: da 10⁻²·⁵ a 10³ ore

Architettura del Modello

1. Modello di Riferimento agli Elementi Finiti

• Utilizza l'interfaccia di trasferimento di calore solido di COMSOL Multiphysics per simulare la conduzione termica nel dominio solido
• Adotta il modulo di flusso in tubo per simulare l'equilibrio energetico del fluido
• Verifica dell'indipendenza della mesh per garantire la precisione
• Ridimensionamento della coordinata verticale (fattore 20) per migliorare la qualità della mesh

2. Modello TRCM Migliorato

Struttura stratificata: BHE e strato sono divisi in m strati orizzontali, con altezza di ogni strato l = L/m

Rete termica: ogni strato contiene i seguenti nodi

• Nodo fluido: f1 (discendente), f2 (ascendente)
• Nodo tubo: p1, p2
• Nodo malta: gt1, gt2
• Nodo pozzo: b
• Nodi strato: g1, g2, ..., gn

Calcolo della resistenza termica:

R₁^Δ = 2R_b                    (fluido a parete pozzo)
R₁₂^Δ = 4R_b R_a/(4R_b - R_a)   (resistenza termica tra fluidi)
R_conv = 1/(2πr_pi h)           (resistenza termica convettiva)

Equazioni di equilibrio energetico: Per ogni nodo viene stabilita un'equazione di equilibrio energetico transitorio, formando un sistema di equazioni lineari da risolvere.

3. Metodo dei Coefficienti di Correzione

Definizione del coefficiente di correzione:

C_coeff = R_b3D,COMSOL / R_b3D,TRCM

Processo di correzione:

1. Il TRCM calcola la risposta termica di base
2. Ottenere il coefficiente di correzione mediante interpolazione in base ai parametri
3. Correggere R_b3D e ricalcolare il campo di temperatura

Punti di Innovazione Tecnica

1. Analisi quantitativa dell'effetto HUST: attraverso l'introduzione di uno strato ad alta conduttività termica che forza il confine isotermico, si analizza quantitativamente l'effetto dell'ipotesi HUST sulla precisione della simulazione
2. Progettazione della rete termica stratificata: rispetto al TRCM tradizionale, aggiunge nodi di tubo e malta, descrivendo più accuratamente il trasferimento di calore interno al BHE
3. Strategia di correzione a posteriori: sfruttando la caratteristica del TRCM di prevedere accuratamente T_b, corregge indirettamente altri parametri correggendo R_b3D
4. Tecnica di interpolazione multi-parametro: basata su quattro parametri chiave - raggio del pozzo, distanza tra tubi, conduttività termica della malta, diametro esterno del tubo - mediante interpolazione parabolica

Configurazione Sperimentale

Configurazione della Simulazione

Intervallo di parametri geometrici:

• Raggio del pozzo r_b: 64-76 mm
• Distanza tra tubi s: 27-47 mm
• Diametro esterno del tubo r_pe: 16, 20 mm
• Conduttività termica della malta k_gt: 1,0-2,2 W/(m·K)

Parametri fissi:

• Lunghezza BHE: 100 m
• Profondità di interramento: 1,8 m
• Conduttività termica dello strato: 1,8 W/(m·K)
• Carico termico: 50 W/m

Costruzione del Dataset

• 54 combinazioni di parametri BHE di diverse configurazioni
• Simulazione comparativa COMSOL e TRCM per ogni gruppo
• Punti temporali: 56 punti equidistanziati logaritmicamente (da 10⁻²·⁵ a 10³ ore)
• Generazione di un database di coefficienti di correzione per l'interpolazione

Casi di Verifica

Selezione di 3 configurazioni BHE tipiche per la verifica:

• BHE 1: configurazione standard, utilizzata per il calcolo dei coefficienti di correzione
• BHE 2: configurazione con parametri che richiedono interpolazione
• BHE 3: configurazione con piccolo raggio che richiede estrapolazione

Risultati Sperimentali

Risultati Principali

Analisi dell'Effetto dell'Ipotesi HUST

Attraverso il confronto tra risultati di simulazione con e senza strato ad alta conduttività termica:

• T_fm, T_in, T_out sono sottostimate, T_b è quasi non influenzata
• R_b3D e R_beff sono significativamente sottostimate
• Il grado di sottostima aumenta nel tempo, particolarmente evidente nel medio termine (>4 ore)

Verifica dell'Effetto di Correzione

Risultati BHE 1:

• Errore quadratico medio T_fm: 0,037°C
• Errore quadratico medio T_in, T_out: <0,038°C
• Errore R_beff al momento finale: 0,066%

Risultati BHE 2 (richiede interpolazione):

• Errore quadratico medio T_fm: 0,051°C
• Errore R_b3D al momento finale: 0,13%
• Errore R_beff al momento finale: 0,23%

Risultati BHE 3 (richiede estrapolazione):

• Errore quadratico medio T_fm: 0,052°C
• Errore R_b3D al momento finale: 0,20%
• Errore R_beff al momento finale: 0,40%

Efficienza Computazionale

• TRCM corretto: completato in 2 secondi
• Simulazione COMSOL: ore
• Miglioramento della precisione: raggiunge la precisione della simulazione agli elementi finiti

Analisi Comparativa

Confronto con la rete neurale artificiale di Pasquier-Marcotte:

• Entrambi i metodi presentano problemi di sottostima nel medio termine
• Il metodo proposto risolve efficacemente questo problema attraverso coefficienti di correzione

Lavori Correlati

Metodi di Soluzione Analitica

• Lamarche & Beauchamp: soluzione analitica a breve termine, ma non può considerare lo scambio termico tra fluidi
• Beier: modello analitico che considera l'equilibrio energetico del fluido, richiede trasformata inversa di Laplace

Metodi di Simulazione Numerica

• Li & Zheng: dominio solido 3D + equilibrio energetico fluido 1D
• Marcotte & Pasquier: simulazione completamente accoppiata di dominio solido e fluido 3D

Metodo TRCM

• Bauer et al.: rete termica triangolare standard, adotta ipotesi HUST
• Pasquier & Marcotte: TRCM migliorato e sviluppo di modello quasi-3D
• Ruiz-Calvo et al.: modello pozzo-strato, evita ipotesi isotermica ma ha precisione limitata

Metodi di Apprendimento Automatico

• Pasquier et al.: rete neurale artificiale addestrata su 15.000 simulazioni TRCM

Conclusioni e Discussione

Conclusioni Principali

1. L'ipotesi HUST è la causa fondamentale della limitazione di precisione del TRCM, causando una sottostima sistematica della risposta termica a medio termine
2. La strategia di correzione a posteriori è efficace: sfruttando il vantaggio del TRCM di prevedere accuratamente T_b, realizza il miglioramento complessivo della precisione correggendo R_b3D
3. Il metodo ibrido combina velocità e precisione: raggiunge la precisione della simulazione agli elementi finiti in 2 secondi, soddisfacendo le esigenze delle applicazioni ingegneristiche
4. L'effetto di interpolazione e estrapolazione è buono: mantiene elevata precisione anche per parametri al di fuori dell'intervallo di interpolazione

Limitazioni

1. Applicabile solo a BHE a tubo singolo: i tipi a doppio tubo o coassiale richiedono la ricostruzione del database di coefficienti di correzione
2. Limitazione dell'intervallo di parametri: i parametri al di fuori dell'intervallo di interpolazione potrebbero influenzare la precisione
3. Interpretazione fisica dei coefficienti di correzione: manca un'analisi profonda dell'associazione tra coefficienti di correzione e meccanismi fisici
4. Precisione a lungo termine non verificata: la precisione della simulazione a lungo termine oltre 1000 ore rimane da verificare

Direzioni Future

1. Estensione ad altri tipi di BHE (doppio tubo, coassiale)
2. Considerazione di meccanismi di trasferimento di calore complessi come il flusso di acque sotterranee
3. Accoppiamento con il metodo della funzione g a lungo termine, realizzando simulazioni su scala temporale completa
4. Sviluppo di strategie di correzione adattive, riducendo la dipendenza da database pre-costruiti

Valutazione Approfondita

Vantaggi

1. Identificazione accurata del problema: analisi approfondita del meccanismo fisico dell'ipotesi HUST e del grado di influenza
2. Progettazione del metodo ingegnosa: sfrutta pienamente i vantaggi del TRCM, compensando le carenze attraverso correzione a posteriori
3. Progettazione sperimentale completa: 54 configurazioni coprono lo spazio dei parametri principale, con verifica sufficiente
4. Significativo valore ingegneristico: il programma C++ open-source ha valore di applicazione pratica
5. Verifica della precisione rigorosa: più casi di verifica provano l'efficacia del metodo

Insufficienze

1. Analisi teorica insufficiente: manca un'interpretazione profonda del significato fisico dei coefficienti di correzione
2. Analisi di sensibilità ai parametri limitata: non analizza sufficientemente l'influenza di ogni parametro sull'effetto di correzione
3. Limitazione dell'intervallo di applicabilità: applicabile solo a BHE a tubo singolo, universalità limitata
4. Mancanza di verifica a lungo termine: la precisione della simulazione a lunghissimo termine rimane da verificare

Impatto

1. Contributo accademico: fornisce un nuovo approccio di modellazione ibrida al campo della simulazione BHE
2. Valore ingegneristico: migliora significativamente l'efficienza e la precisione della simulazione BHE a breve e medio termine
3. Contributo open-source: il programma C++ fornito promuove la diffusione e l'applicazione del metodo
4. Ispirazione metodologica: la strategia di correzione a posteriori può essere estesa ad altri problemi di simulazione ingegneristica

Scenari Applicabili

1. Progettazione di sistemi di pompe di calore geotermiche: valutazione rapida e accurata delle prestazioni termiche a breve e medio termine del BHE
2. Analisi di test di risposta termica: assistenza nell'interpretazione dei dati di test sul campo e nell'identificazione dei parametri
3. Ricerca di ottimizzazione del sistema: fornisce uno strumento di simulazione efficiente per l'ottimizzazione dei parametri BHE
4. Insegnamento e ricerca: come strumento standard per l'analisi del trasferimento di calore BHE

Bibliografia

L'articolo cita 45 articoli correlati, coprendo i principali metodi e i progressi recenti nella simulazione BHE, fornendo una base teorica solida per questa ricerca. Le referenze chiave includono:

• Claesson & Javed (2011): metodo di analisi della funzione g
• Bauer et al. (2011): metodo standard TRCM
• Pasquier & Marcotte (2014, 2020): metodo TRCM migliorato e metodo di apprendimento automatico
• Zanchini (2024): analisi comparativa dei metodi di simulazione BHE

Valutazione Complessiva: Questo è un articolo di alta qualità di fisica applicata che, affrontando le esigenze ingegneristiche pratiche della simulazione BHE per pompe di calore geotermiche, propone un metodo di modellazione ibrida innovativo. L'articolo ha un'analisi teorica profonda, una verifica sperimentale completa e possiede significativo valore accademico e prospettive di applicazione ingegneristica. La fornitura del programma open-source ne aumenta ulteriormente l'influenza e il valore pratico.