Ein kombiniertes Wärmewiderstand-Kapazitäts- und Finite-Elemente-Modell für sehr schnelle und genaue Kurz- und Mittelfrist-Simulationen von einzelnen U-Rohr-Bohrlochwärmetauschern

Grundinformationen

• Papier-ID: 2510.14421
• Titel: Ein kombiniertes Wärmewiderstand-Kapazitäts- und Finite-Elemente-Modell für sehr schnelle und genaue Kurz- und Mittelfrist-Simulationen von einzelnen U-Rohr-Bohrlochwärmetauschern
• Autoren: Enzo Zanchini, Francesco Zanchini, Claudia Naldi (Universität Bologna)
• Klassifizierung: physics.app-ph
• Veröffentlichungsjahr: 2024
• Papier-Link: https://arxiv.org/abs/2510.14421

Zusammenfassung

Die genaue Auslegung von Erdwärmepumpensystemen erfordert die Kenntnis der Auslassfluidtemperatur von Bohrlochwärmetauschern (BHE) im Kurz- und Langzeitbereich. Dieser Artikel konzentriert sich auf die Kurz- und Mittelfristanalyse. In diesem Zeitbereich können 3D-Finite-Elemente-Simulationen oder Wärmewiderstand-Kapazitäts-Modelle (TRCM) verwendet werden. Erstere bieten hohe Genauigkeit, erfordern aber lange Rechenzeiten, letztere sind schnell, aber aufgrund vereinfachter Annahmen in ihrer Genauigkeit begrenzt. Dieser Artikel präsentiert eine neue Methode, die die Geschwindigkeit von TRCM mit der Genauigkeit von Finite-Elemente-Simulationen für Kurz- und Mittelfrist-Simulationen einzelner U-Rohr-Bohrlochwärmetauscher kombiniert. Die Methode nutzt TRCM zur Schätzung der thermischen Reaktion des BHE und korrigiert diese dann durch Interpolation mit einem Datensatz aus 54 Finite-Elemente-Simulationen verschiedener Konfigurationen. Das Modell wurde als C++-Programm implementiert und in der Open-Source-Datenbank der Universität Bologna bereitgestellt. Es liefert die zeitliche Entwicklung von Ein- und Auslasstemperatur sowie durchschnittlicher Fluidtemperatur, durchschnittlicher BHE-Oberflächentemperatur, 3D- und effektivem Bohrlochwärmewiderstand innerhalb von zwei Sekunden.

Forschungshintergrund und Motivation

Problembeschreibung

Erdwärmepumpensysteme sind eine effiziente Technologie für Gebäudeheizung und -kühlung. Ihre Auslegung erfordert eine genaue Vorhersage der thermischen Reaktionseigenschaften von Bohrlochwärmetauschern. Die Hauptherausforderungen sind:

1. Widerspruch zwischen Rechengenauigkeit und Effizienz: 3D-Finite-Elemente-Simulationen sind genau, benötigen aber Stunden; TRCM ist schnell, aber in der Genauigkeit begrenzt
2. Anforderungen an Kurz- und Mittelfrist-Simulationen: Bestehende Methoden weisen im Zeitbereich von 10⁻²·⁵ bis 10³ Stunden Genauigkeits- oder Effizienzprobleme auf
3. Anforderungen praktischer Ingenieuranwendungen: Es wird ein Simulationswerkzeug benötigt, das sowohl schnell als auch genau ist, für Systemauslegung und -optimierung

Einschränkungen bestehender Methoden

1. Analytische Lösungsmethoden: Ersetzen das Fluid durch einen wärmeerzeugenden Feststoff oder eine Oberfläche, können aber den Wärmeeintrag zwischen auf- und absteigendem Fluid nicht berücksichtigen
2. 3D-Numerische Simulation: Hohe Genauigkeit, aber einzelne Simulationen benötigen Stunden und sind für schnelle Ingenieurentwürfe ungeeignet
3. Bestehende TRCM-Modelle: Verwenden die Annahme einer horizontal einheitlichen Oberflächentemperatur (HUST), was zu einer Unterschätzung der Mittelfrist-Wärmereaktionen führt

Forschungsmotivation

Entwicklung einer hybriden Simulationsmethode, die die Geschwindigkeit von TRCM mit der Genauigkeit von Finite-Elementen kombiniert und praktische Ingenieuranforderungen erfüllt.

Kernbeiträge

1. Identifikation und Quantifizierung der Auswirkungen der HUST-Annahme: Durch Vergleich von Finite-Elemente-Simulationen mit und ohne hochleitfähige Schicht wird nachgewiesen, dass die HUST-Annahme T_fm, T_in, T_out, R_b3D und R_beff unterschätzt
2. Entwicklung eines verbesserten TRCM-Modells: Verwendung eines geschichteten Wärmewiderstandsnetzwerks mit Knoten für Fluid, Rohre, Verfüllung und Gestein, erhöht die Simulationsgenauigkeit
3. Erstellung einer Korrekturfaktordatenbank: Durch 54 Finite-Elemente-Simulationen verschiedener Konfigurationen wird ein zeitabhängiger Korrekturfaktordatensatz erstellt
4. Implementierung einer effizienten hybriden Simulationsmethode: Kombination von TRCM-Basisberechnungen und Interpolationskorrektur erreicht Finite-Elemente-Simulationsgenauigkeit in 2 Sekunden
5. Bereitstellung eines Open-Source-Simulationswerkzeugs: C++-Programm für schnelle und genaue BHE-Kurz- und Mittelfrist-Simulationen

Methodische Details

Aufgabendefinition

Eingaben: BHE-Geometrieparameter (Länge L, Radius r_b, Rohrabstand s usw.), thermophysikalische Parameter, Durchfluss, Wärmeleistung Ausgaben: Zeitliche Entwicklung von Ein- und Auslasstemperatur (T_in, T_out), durchschnittliche Fluidtemperatur (T_fm), BHE-Oberflächentemperatur (T_b), 3D-Wärmewiderstand (R_b3D), effektiver Wärmewiderstand (R_beff) Zeitbereich: 10⁻²·⁵ bis 10³ Stunden

Modellarchitektur

1. Finite-Elemente-Referenzmodell

• Verwendung der Wärmetransport-Festkörper-Schnittstelle von COMSOL Multiphysics zur Simulation der Wärmeleitung im Festkörperbereich
• Anwendung des Rohrdurchströmungsmoduls zur Simulation des Fluidenergiehaushalts
• Netzunabhängigkeitsverifizierung gewährleistet Genauigkeit
• Vertikale Koordinatenskalierung (Faktor 20) verbessert die Netzqualität

2. Verbessertes TRCM-Modell

Geschichtete Struktur: BHE und Gestein werden in m horizontale Schichten unterteilt, jede mit Höhe l = L/m

Wärmewiderstandsnetzwerk: Jede Schicht enthält folgende Knoten

• Fluidknoten: f1 (absteigend), f2 (aufsteigend)
• Rohrknoten: p1, p2
• Verfüllungsknoten: gt1, gt2
• Bohrlochnoten: b
• Gesteinsknoten: g1, g2, ..., gn

Wärmewiderstandsberechnung:

R₁^Δ = 2R_b                    (Fluid zu Bohrlochwand)
R₁₂^Δ = 4R_b R_a/(4R_b - R_a)   (Wärmewiderstand zwischen Fluiden)
R_conv = 1/(2πr_pi h)           (Konvektionswärmewiderstand)

Energiebilanzgleichungen: Für jeden Knoten werden instationäre Energiebilanzgleichungen aufgestellt, die ein lineares Gleichungssystem bilden.

3. Korrekturfaktormethode

Korrekturfaktordefinition:

C_coeff = R_b3D,COMSOL / R_b3D,TRCM

Korrekturprozess:

1. TRCM berechnet die grundlegende thermische Reaktion
2. Interpolation der Korrekturfaktoren basierend auf Parametern
3. Korrektur von R_b3D und Neuberechnung des Temperaturfeldes

Technische Innovationen

1. Quantitative Analyse der HUST-Auswirkungen: Durch Einführung einer hochleitfähigen Schicht zur Erzwingung einer isothermen Grenze wird die Auswirkung der HUST-Annahme auf die Simulationsgenauigkeit quantifiziert
2. Entwurf eines geschichteten Wärmewiderstandsnetzwerks: Im Vergleich zu traditionellem TRCM werden Rohr- und Verfüllungsknoten hinzugefügt, um die Wärmeleitung innerhalb des BHE genauer zu beschreiben
3. A-posteriori-Korrekturbstrategie: Nutzt die Genauigkeit von TRCM bei der Vorhersage von T_b, korrigiert indirekt andere Parameter durch Korrektur von R_b3D
4. Mehrparameter-Interpolationstechnik: Parabolische Interpolation basierend auf vier Schlüsselparametern: Bohrlochhalbmesser, Rohrabstand, Verfüllungswärmeleitfähigkeit und Rohraußendurchmesser

Experimentelle Einrichtung

Simulationskonfiguration

Geometrieparameterbereiche:

• Bohrlochhalbmesser r_b: 64-76 mm
• Rohrabstand s: 27-47 mm
• Rohraußendurchmesser r_pe: 16, 20 mm
• Verfüllungswärmeleitfähigkeit k_gt: 1,0-2,2 W/(m·K)

Feste Parameter:

• BHE-Länge: 100 m
• Einbautiefe: 1,8 m
• Gesteins-Wärmeleitfähigkeit: 1,8 W/(m·K)
• Wärmeleistung: 50 W/m

Datensatzerstellung

• 54 verschiedene BHE-Parameterkombinationen
• Vergleichssimulationen mit COMSOL und TRCM für jede Gruppe
• Zeitpunkte: 56 logarithmisch gleichmäßig verteilte Punkte (10⁻²·⁵ bis 10³ Stunden)
• Erstellung einer Korrekturfaktordatenbank für Interpolation

Validierungsfälle

Drei typische BHE-Konfigurationen zur Validierung:

• BHE 1: Standardkonfiguration, verwendet für Korrekturfaktorberechnung
• BHE 2: Konfiguration mit zu interpolierenden Parametern
• BHE 3: Konfiguration mit kleinem Radius, erfordert Extrapolation

Experimentelle Ergebnisse

Hauptergebnisse

Analyse der HUST-Annahmeauswirkungen

Durch Vergleich von Simulationen mit und ohne hochleitfähige Schicht:

• T_fm, T_in, T_out werden unterschätzt, T_b ist kaum beeinflusst
• R_b3D und R_beff werden erheblich unterschätzt
• Das Ausmaß der Unterschätzung nimmt mit der Zeit zu, besonders im Mittelfristbereich (>4 Stunden)

Validierung der Korrektureffektivität

BHE 1 Ergebnisse:

• T_fm Quadratwurzelmittelfehler: 0,037°C
• T_in, T_out Quadratwurzelmittelfehler: <0,038°C
• R_beff Fehler zum Endzeitpunkt: 0,066%

BHE 2 Ergebnisse (erfordert Interpolation):

• T_fm Quadratwurzelmittelfehler: 0,051°C
• R_b3D Fehler zum Endzeitpunkt: 0,13%
• R_beff Fehler zum Endzeitpunkt: 0,23%

BHE 3 Ergebnisse (erfordert Extrapolation):

• T_fm Quadratwurzelmittelfehler: 0,052°C
• R_b3D Fehler zum Endzeitpunkt: 0,20%
• R_beff Fehler zum Endzeitpunkt: 0,40%

Recheneffizienz

• Korrigiertes TRCM: Abschluss in 2 Sekunden
• COMSOL-Simulation: Mehrere Stunden
• Genauigkeitssteigerung: Erreicht Finite-Elemente-Simulationsgenauigkeit

Vergleichende Analyse

Vergleich mit dem künstlichen Neuronalen Netz von Pasquier-Marcotte:

• Beide Methoden zeigen Mittelfrist-Unterschätzungsprobleme
• Die vorliegende Methode löst dieses Problem effektiv durch Korrekturfaktoren

Verwandte Arbeiten

Analytische Lösungsmethoden

• Lamarche & Beauchamp: Analytische Kurzzeitlesung, kann aber Wärmeeintrag zwischen Fluiden nicht berücksichtigen
• Beier: Analytisches Modell mit Fluidenergiebilanz, erfordert inverse Laplace-Transformation

Numerische Simulationsmethoden

• Li & Zheng: 3D-Festkörperbereich + 1D-Fluidenergiebilanz
• Marcotte & Pasquier: Vollständig gekoppelte 3D-Fest- und Fluidbereichssimulation

TRCM-Methoden

• Bauer et al.: Standarddreieck-Wärmewiderstandsnetzwerk mit HUST-Annahme
• Pasquier & Marcotte: Verbessertes TRCM und Entwicklung zu Quasi-3D-Modell
• Ruiz-Calvo et al.: Bohrloch-Gesteins-Modell, vermeidet Isotherm-Annahme, aber begrenzte Genauigkeit

Maschinenlernmethoden

• Pasquier et al.: Künstliches Neuronales Netz trainiert auf 15.000 TRCM-Simulationen

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

1. Die HUST-Annahme ist die grundlegende Ursache der TRCM-Genauigkeitsbegrenzung, führt zu systematischer Unterschätzung der Mittelfrist-Wärmereaktionen
2. A-posteriori-Korrekturbstrategie ist effektiv: Nutzt die Vorteile von TRCM bei der genauen Vorhersage von T_b, erreicht Gesamtgenauigkeitssteigerung durch Korrektur von R_b3D
3. Hybridmethode kombiniert Geschwindigkeit und Genauigkeit: Erreicht Finite-Elemente-Simulationsgenauigkeit in 2 Sekunden, erfüllt Anforderungen technischer Anwendungen
4. Interpolations- und Extrapolationseffekte sind gut: Behält hohe Genauigkeit auch bei Parametern außerhalb des Bereichs

Einschränkungen

1. Gilt nur für einzelne U-Rohr-BHE: Doppel-U-Rohr oder koaxiale Typen erfordern neu erstellte Korrekturdatenbanken
2. Parameterbereiche sind begrenzt: Parameter außerhalb des Interpolationsbereichs können die Genauigkeit beeinflussen
3. Physikalische Interpretation der Korrekturfaktoren: Mangelnde tiefere Analyse der Beziehung zwischen Korrekturfaktoren und physikalischen Mechanismen
4. Langzeitgenauigkeit nicht validiert: Genauigkeit von Simulationen über 1000 Stunden hinaus bleibt zu überprüfen

Zukünftige Richtungen

1. Erweiterung auf andere BHE-Typen (Doppel-U-Rohr, koaxial)
2. Berücksichtigung komplexer Wärmeübertragungsmechanismen wie Grundwasserströmung
3. Kopplung mit Langzeit-g-Funktionsmethoden für Vollzeitskalen-Simulation
4. Entwicklung adaptiver Korrekturstrategien zur Verringerung der Abhängigkeit von vorgebauten Datenbanken

Tiefgreifende Bewertung

Stärken

1. Genaue Problemidentifikation: Tiefgehende Analyse des physikalischen Mechanismus und der Auswirkungen der HUST-Annahme
2. Geschickter Methodenentwurf: Nutzt TRCM-Vorteile vollständig, kompensiert Mängel durch A-posteriori-Korrektur
3. Umfassende Experimentalplanung: 54 Konfigurationen decken den Hauptparameterraum ab, ausreichende Validierung
4. Bedeutsamer Ingenieurwert: Bereitgestelltes Open-Source-C++-Programm hat praktischen Anwendungswert
5. Strenge Genauigkeitsvalidierung: Mehrere Validierungsfälle beweisen Methodeneffektivität

Mängel

1. Unzureichende theoretische Analyse: Mangelnde tiefgehende Erklärung der physikalischen Bedeutung von Korrekturfaktoren
2. Begrenzte Parametersensitivitätsanalyse: Unzureichende Analyse der Auswirkungen verschiedener Parameter auf Korrektureffekte
3. Begrenzte Anwendungsbereiche: Nur für einzelne U-Rohr-BHE, begrenzte Universalität
4. Fehlende Langzeitvalidierung: Genauigkeit von Ultralangzeit-Simulationen bleibt zu überprüfen

Einfluss

1. Akademischer Beitrag: Bietet neue hybride Modellierungsideen für BHE-Simulationsbereich
2. Ingenieurwert: Verbessert Effizienz und Genauigkeit von BHE-Kurz- und Mittelfrist-Simulationen erheblich
3. Open-Source-Beitrag: Bereitgestelltes C++-Programm fördert Methodenverbreitung und -anwendung
4. Methodologische Inspiration: A-posteriori-Korrekturbstrategie kann auf andere Ingenieur-Simulationsprobleme übertragen werden

Anwendungsszenarien

1. Erdwärmepumpensystem-Auslegung: Schnelle und genaue Bewertung der Kurz- und Mittelfrist-Wärmeleistung von BHE
2. Thermische Reaktionstestanalyse: Unterstützt Interpretation von Feldtestdaten und Parameteridentifikation
3. Systemoptimierungsforschung: Bietet effizientes Simulationswerkzeug für BHE-Parameteroptimierung
4. Lehre und Forschung: Dient als Standardwerkzeug für BHE-Wärmeübertragungsanalyse

Literaturverzeichnis

Das Papier zitiert 45 relevante Literaturquellen, die die Hauptmethoden und neuesten Fortschritte in der BHE-Simulation abdecken und eine solide theoretische Grundlage für diese Forschung bieten. Wichtige Literaturquellen umfassen:

• Claesson & Javed (2011): g-Funktionsanalysemethode
• Bauer et al. (2011): TRCM-Standardmethode
• Pasquier & Marcotte (2014, 2020): Verbessertes TRCM und Maschinenlernmethoden
• Zanchini (2024): Vergleichende Analyse von BHE-Simulationsmethoden

Gesamtbewertung: Dies ist ein hochqualitatives Papier der angewandten Physik, das eine innovative hybride Modellierungsmethode für die praktischen Ingenieuranforderungen der BHE-Simulation von Erdwärmepumpensystemen vorschlägt. Das Papier bietet tiefgehende theoretische Analysen und umfassende experimentelle Validierung mit bedeutsamen akademischen Werten und Perspektiven für technische Anwendungen. Die Bereitstellung des Open-Source-Programms verstärkt seinen Einfluss und praktischen Wert weiter.