Un modelo combinado de resistencia-capacidad térmica y elementos finitos para simulaciones muy rápidas y precisas de intercambiadores de calor de sondeos U simples a corto y medio plazo

Información Básica

• ID del Artículo: 2510.14421
• Título: Un modelo combinado de resistencia-capacidad térmica y elementos finitos para simulaciones muy rápidas y precisas de intercambiadores de calor de sondeos U simples a corto y medio plazo
• Autores: Enzo Zanchini, Francesco Zanchini, Claudia Naldi (Universidad de Bolonia)
• Clasificación: physics.app-ph
• Fecha de Publicación: 2024
• Enlace del Artículo: https://arxiv.org/abs/2510.14421

Resumen

El diseño preciso de sistemas de bombas de calor geotérmicas requiere comprender la temperatura del fluido de salida de los intercambiadores de calor de sondeos (BHE) a corto y largo plazo. Este artículo se enfoca en análisis a corto y medio plazo. Dentro de este rango temporal, se pueden emplear simulaciones de elementos finitos 3D o modelos de resistencia-capacidad térmica (TRCM). El primero ofrece alta precisión pero requiere tiempo computacional prolongado, mientras que el segundo es rápido pero tiene precisión limitada debido a suposiciones simplificadas. Este artículo propone un nuevo método que combina la velocidad del TRCM con la precisión de la simulación de elementos finitos para simulaciones a corto y medio plazo de intercambiadores de calor de sondeos U simples. El método utiliza TRCM para estimar la respuesta térmica del BHE, luego se corrige mediante interpolación con un conjunto de datos generado a partir de 54 simulaciones de elementos finitos con diferentes configuraciones. El modelo se ha implementado como un programa en C++ y está disponible en la base de datos de código abierto de la Universidad de Bolonia, proporcionando en dos segundos la evolución temporal de las temperaturas de entrada, salida y promedio del fluido, la temperatura promedio de la superficie del BHE, la resistencia térmica 3D y la resistencia térmica efectiva del sondeo.

Antecedentes de Investigación y Motivación

Descripción del Problema

Los sistemas de bombas de calor geotérmicas, como tecnología eficiente de calefacción y refrigeración de edificios, requieren predicciones precisas de las características de respuesta térmica de los intercambiadores de calor de sondeos. Los desafíos clave incluyen:

1. Contradicción entre precisión computacional y eficiencia: Las simulaciones 3D de elementos finitos ofrecen alta precisión pero requieren horas; el TRCM es rápido pero tiene precisión limitada
2. Necesidad de simulaciones a corto y medio plazo: Los métodos existentes presentan problemas de precisión o eficiencia en el rango temporal de 10⁻²·⁵ a 10³ horas
3. Requisitos de aplicaciones de ingeniería práctica: Se necesitan herramientas de simulación que sean tanto rápidas como precisas para el diseño y optimización de sistemas

Limitaciones de Métodos Existentes

1. Métodos de soluciones analíticas: Reemplazan el fluido con un sólido generador de calor o una superficie, sin poder considerar los efectos de intercambio de calor entre fluidos ascendentes y descendentes
2. Simulación numérica 3D: Alta precisión pero cada simulación requiere horas, inadecuada para diseño rápido de ingeniería
3. Modelos TRCM existentes: Adoptan la suposición de temperatura de superficie horizontalmente uniforme (HUST), lo que conduce a subestimación de la respuesta térmica a medio plazo

Motivación de la Investigación

Desarrollar un método de simulación híbrida que combine la velocidad del TRCM con la precisión de elementos finitos, satisfaciendo las necesidades prácticas de ingeniería.

Contribuciones Principales

1. Identificación y cuantificación del impacto de la suposición HUST: Mediante comparación de simulaciones de elementos finitos con y sin capas de alta conductividad térmica, se demuestra que la suposición HUST subestima T_fm, T_in, T_out, R_b3D y R_beff
2. Desarrollo de un modelo TRCM mejorado: Adopta una red de resistencia térmica estratificada que incluye nodos de fluido, tuberías, lechada y estrato, mejorando la precisión de simulación
3. Establecimiento de una base de datos de coeficientes de corrección: Mediante 54 simulaciones de elementos finitos con diferentes configuraciones, se construye un conjunto de datos de coeficientes de corrección dependientes del tiempo
4. Implementación de un método de simulación híbrida eficiente: Combinando cálculos base TRCM con corrección por interpolación, alcanza la precisión de simulación de elementos finitos en 2 segundos
5. Provisión de herramienta de simulación de código abierto: El programa en C++ puede utilizarse para simulaciones rápidas y precisas de BHE a corto y medio plazo

Explicación Detallada del Método

Definición de la Tarea

Entrada: Parámetros geométricos del BHE (longitud L, radio r_b, espaciamiento entre tuberías s, etc.), parámetros de propiedades térmicas, caudal, carga térmica Salida: Evolución temporal de temperaturas de entrada y salida (T_in, T_out), temperatura promedio del fluido (T_fm), temperatura de la superficie del BHE (T_b), resistencia térmica 3D (R_b3D), resistencia térmica efectiva (R_beff) Rango temporal: 10⁻²·⁵ a 10³ horas

Arquitectura del Modelo

1. Modelo de Referencia de Elementos Finitos

• Utiliza la interfaz de transferencia de calor en sólidos de COMSOL Multiphysics para simular conducción de calor en el dominio sólido
• Emplea el módulo de flujo en tuberías para simular el balance de energía del fluido
• Verificación de independencia de malla para asegurar precisión
• Escalado de coordenadas verticales (factor 20) para mejorar la calidad de la malla

2. Modelo TRCM Mejorado

Estructura estratificada: El BHE y el estrato se dividen en m capas horizontales, cada una con altura l = L/m

Red de resistencia térmica: Cada capa contiene los siguientes nodos

• Nodos de fluido: f1 (descendente), f2 (ascendente)
• Nodos de tubería: p1, p2
• Nodos de lechada: gt1, gt2
• Nodo de sondeo: b
• Nodos de estrato: g1, g2, ..., gn

Cálculo de resistencias térmicas:

R₁^Δ = 2R_b                    (fluido a pared de sondeo)
R₁₂^Δ = 4R_b R_a/(4R_b - R_a)   (resistencia térmica entre fluidos)
R_conv = 1/(2πr_pi h)           (resistencia térmica convectiva)

Ecuaciones de balance de energía: Se establece una ecuación de balance de energía transitoria para cada nodo, formando un sistema de ecuaciones lineales para resolver.

3. Método de Coeficientes de Corrección

Definición de coeficientes de corrección:

C_coeff = R_b3D,COMSOL / R_b3D,TRCM

Proceso de corrección:

1. TRCM calcula la respuesta térmica base
2. Se obtienen coeficientes de corrección mediante interpolación según parámetros
3. Se corrige R_b3D y se recalcula el campo de temperaturas

Puntos de Innovación Técnica

1. Análisis cuantitativo del impacto HUST: Mediante la introducción de capas de alta conductividad térmica que fuerzan límites isotérmicos, se analiza cuantitativamente el impacto de la suposición HUST en la precisión de simulación
2. Diseño de red de resistencia térmica estratificada: En comparación con TRCM tradicional, añade nodos de tubería y lechada, describiendo más precisamente la transferencia de calor dentro del BHE
3. Estrategia de corrección a posteriori: Aprovecha la característica de que TRCM predice con precisión T_b, corrigiendo indirectamente otros parámetros mediante la corrección de R_b3D
4. Técnica de interpolación multiparamétrica: Basada en cuatro parámetros clave: radio del sondeo, espaciamiento entre tuberías, conductividad térmica de la lechada y diámetro exterior de la tubería, utilizando interpolación parabólica

Configuración Experimental

Configuración de Simulación

Rango de parámetros geométricos:

• Radio del sondeo r_b: 64-76 mm
• Espaciamiento entre tuberías s: 27-47 mm
• Diámetro exterior de tubería r_pe: 16, 20 mm
• Conductividad térmica de lechada k_gt: 1.0-2.2 W/(m·K)

Parámetros fijos:

• Longitud del BHE: 100 m
• Profundidad de enterramiento: 1.8 m
• Conductividad térmica del estrato: 1.8 W/(m·K)
• Carga térmica: 50 W/m

Construcción del Conjunto de Datos

• 54 combinaciones de parámetros BHE con diferentes configuraciones
• Simulaciones comparativas COMSOL y TRCM para cada grupo
• Puntos temporales: 56 puntos logarítmicamente equiespaciados (10⁻²·⁵ a 10³ horas)
• Generación de base de datos de coeficientes de corrección para interpolación

Casos de Verificación

Se seleccionan 3 configuraciones típicas de BHE para verificación:

• BHE 1: Configuración estándar, utilizada para cálculo de coeficientes de corrección
• BHE 2: Configuración con parámetros que requieren interpolación
• BHE 3: Configuración de pequeño radio que requiere extrapolación

Resultados Experimentales

Resultados Principales

Análisis del Impacto de la Suposición HUST

Mediante comparación de resultados de simulación con y sin capas de alta conductividad térmica:

• T_fm, T_in, T_out se subestiman, T_b prácticamente sin efecto
• R_b3D y R_beff se subestiman significativamente
• El grado de subestimación aumenta con el tiempo, siendo especialmente notable a medio plazo (>4 horas)

Verificación del Efecto de Corrección

Resultados de BHE 1:

• Error cuadrático medio de T_fm: 0.037°C
• Error cuadrático medio de T_in, T_out: <0.038°C
• Error de R_beff en el instante final: 0.066%

Resultados de BHE 2 (requiere interpolación):

• Error cuadrático medio de T_fm: 0.051°C
• Error de R_b3D en el instante final: 0.13%
• Error de R_beff en el instante final: 0.23%

Resultados de BHE 3 (requiere extrapolación):

• Error cuadrático medio de T_fm: 0.052°C
• Error de R_b3D en el instante final: 0.20%
• Error de R_beff en el instante final: 0.40%

Eficiencia Computacional

• TRCM corregido: Completado en 2 segundos
• Simulación COMSOL: Horas
• Mejora de precisión: Alcanza la precisión de simulación de elementos finitos

Análisis Comparativo

Comparación con la red neuronal artificial de Pasquier-Marcotte:

• Ambos métodos presentan problemas de subestimación a medio plazo
• El método propuesto resuelve efectivamente este problema mediante coeficientes de corrección

Trabajos Relacionados

Métodos de Soluciones Analíticas

• Lamarche & Beauchamp: Soluciones analíticas a corto plazo, pero sin poder considerar intercambio de calor entre fluidos
• Beier: Modelo analítico considerando balance de energía del fluido, requiere transformada inversa de Laplace

Métodos de Simulación Numérica

• Li & Zheng: Dominio sólido 3D + balance de energía de fluido 1D
• Marcotte & Pasquier: Simulación completamente acoplada de dominio sólido y fluido 3D

Métodos TRCM

• Bauer et al.: Red de resistencia térmica triangular estándar, adopta suposición HUST
• Pasquier & Marcotte: TRCM mejorado y desarrollo de modelo cuasi-3D
• Ruiz-Calvo et al.: Modelo sondeo-estrato, evita suposición isotérmica pero precisión limitada

Métodos de Aprendizaje Automático

• Pasquier et al.: Red neuronal artificial entrenada con 15000 simulaciones TRCM

Conclusiones y Discusión

Conclusiones Principales

1. La suposición HUST es la causa fundamental de las limitaciones de precisión del TRCM, causando subestimación sistemática de la respuesta térmica a medio plazo
2. La estrategia de corrección a posteriori es efectiva: Aprovechando la ventaja de que TRCM predice con precisión T_b, se logra mejora general de precisión mediante corrección de R_b3D
3. El método híbrido combina velocidad y precisión: Alcanza la precisión de simulación de elementos finitos en 2 segundos, satisfaciendo requisitos de aplicaciones de ingeniería
4. El efecto de interpolación y extrapolación es bueno: Incluso con parámetros fuera del rango de interpolación se mantiene alta precisión

Limitaciones

1. Aplicable solo a BHE tipo U simple: Los tipos U doble o coaxial requieren establecer nuevas bases de datos de coeficientes de corrección
2. Limitación del rango de parámetros: Los parámetros fuera del rango de interpolación pueden afectar la precisión
3. Falta de interpretación física de coeficientes de corrección: Carencia de análisis profundo de la asociación entre coeficientes de corrección y mecanismos físicos
4. Precisión a largo plazo no verificada: La precisión de simulaciones a largo plazo superiores a 1000 horas requiere verificación

Direcciones Futuras

1. Extensión a otros tipos de BHE (U doble, coaxial)
2. Consideración de mecanismos complejos de transferencia de calor como flujo de agua subterránea
3. Acoplamiento con métodos de función g a largo plazo, realizando simulaciones en toda la escala temporal
4. Desarrollo de estrategias de corrección adaptativa, reduciendo la dependencia de bases de datos preconstruidas

Evaluación Profunda

Ventajas

1. Identificación precisa del problema: Análisis profundo del mecanismo físico de la suposición HUST y su grado de impacto
2. Diseño de método ingenioso: Aprovecha completamente las ventajas del TRCM, compensando deficiencias mediante corrección a posteriori
3. Diseño experimental completo: 54 configuraciones cubren el espacio de parámetros principal, verificación suficiente
4. Valor de ingeniería significativo: El programa C++ de código abierto posee valor de aplicación práctica
5. Verificación de precisión rigurosa: Múltiples casos de verificación demuestran la efectividad del método

Insuficiencias

1. Análisis teórico insuficiente: Falta análisis profundo del significado físico de los coeficientes de corrección
2. Análisis limitado de sensibilidad de parámetros: No se analiza suficientemente el impacto de cada parámetro en el efecto de corrección
3. Rango de aplicabilidad limitado: Solo se aplica a BHE tipo U simple, universalidad limitada
4. Falta de verificación a largo plazo: La precisión de simulaciones a muy largo plazo requiere verificación

Impacto

1. Contribución académica: Proporciona nuevas ideas de modelado híbrido para el campo de simulación de BHE
2. Valor de ingeniería: Mejora significativamente la eficiencia y precisión de simulación de BHE a corto y medio plazo
3. Contribución de código abierto: El programa C++ proporcionado promueve la difusión y aplicación del método
4. Inspiración metodológica: La estrategia de corrección a posteriori puede generalizarse a otros problemas de simulación de ingeniería

Escenarios Aplicables

1. Diseño de sistemas de bombas de calor geotérmicas: Evaluación rápida y precisa del desempeño térmico a corto y medio plazo del BHE
2. Análisis de pruebas de respuesta térmica: Asistencia en la interpretación de datos de pruebas de campo y identificación de parámetros
3. Investigación de optimización de sistemas: Proporciona herramienta de simulación eficiente para optimización de parámetros de BHE
4. Enseñanza e investigación: Sirve como herramienta estándar para análisis de transferencia de calor de BHE

Referencias

El artículo cita 45 referencias relacionadas, cubriendo los principales métodos y avances recientes en simulación de BHE, proporcionando una base teórica sólida para esta investigación. Las referencias clave incluyen:

• Claesson & Javed (2011): Método de análisis de función g
• Bauer et al. (2011): Método estándar TRCM
• Pasquier & Marcotte (2014, 2020): Métodos TRCM mejorados y de aprendizaje automático
• Zanchini (2024): Análisis comparativo de métodos de simulación de BHE

Evaluación General: Este es un artículo de alta calidad en física aplicada que, respondiendo a las necesidades prácticas de ingeniería en simulación de BHE de bombas de calor geotérmicas, propone un método de modelado híbrido innovador. El artículo presenta análisis teórico profundo, verificación experimental completa, y posee valor académico significativo y perspectivas de aplicación práctica de ingeniería. La provisión del programa de código abierto mejora aún más su influencia y valor práctico.