Комбинированная модель теплового сопротивления-емкости и конечных элементов для быстрого и точного моделирования скважинных теплообменников с одиночной U-образной трубой в краткосрочном и среднесрочном периодах

Основная информация

• ID статьи: 2510.14421
• Название: A combined thermal-resistance-capacity and finite-element model for very fast and accurate short- and medium-term simulations of single U-tube borehole heat exchangers
• Авторы: Enzo Zanchini, Francesco Zanchini, Claudia Naldi (Болонский университет)
• Классификация: physics.app-ph
• Дата публикации: 2024
• Ссылка на статью: https://arxiv.org/abs/2510.14421

Аннотация

Точное проектирование систем геотермальных тепловых насосов требует понимания температуры выходящей жидкости скважинного теплообменника (БХЭ) в краткосрочном и долгосрочном периодах. Данная работа сосредоточена на краткосрочном и среднесрочном анализе. В этом временном диапазоне можно применять трехмерное моделирование методом конечных элементов (МКЭ) или модели теплового сопротивления-емкости (МТСЕ). Первый подход обеспечивает высокую точность, но требует значительных вычислительных затрат, второй — быстрый, но имеет ограниченную точность из-за упрощающих предположений. В настоящей работе предложен новый метод, сочетающий скорость МТСЕ и точность моделирования конечными элементами для краткосрочного и среднесрочного моделирования одиночного U-образного скважинного теплообменника. Метод использует МТСЕ для оценки теплового отклика БХЭ, а затем применяет интерполяционную коррекцию на основе набора данных, полученного из 54 конфигураций конечно-элементного моделирования. Модель реализована в виде программы на C++ и доступна в открытой базе данных Болонского университета. Она обеспечивает временную эволюцию входной и выходной температур жидкости, средней температуры жидкости, средней температуры поверхности БХЭ, трехмерного и эффективного теплового сопротивления скважины за две секунды.

Исследовательский контекст и мотивация

Описание проблемы

Системы геотермальных тепловых насосов как эффективная технология отопления и охлаждения зданий требуют точного прогнозирования характеристик теплового отклика скважинного теплообменника. Ключевые вызовы включают:

1. Противоречие между точностью и эффективностью вычислений: трехмерное моделирование МКЭ обеспечивает высокую точность, но требует часов вычислений; МТСЕ быстра, но имеет ограниченную точность
2. Потребность в краткосрочном и среднесрочном моделировании: существующие методы имеют проблемы с точностью или эффективностью в диапазоне времени от 10⁻²·⁵ до 10³ часов
3. Требования практического инженерного применения: необходим инструмент моделирования, который одновременно быстр и точен для проектирования и оптимизации систем

Ограничения существующих методов

1. Аналитические решения: заменяют жидкость нагревающимся твердым телом или поверхностью, не могут учитывать теплообмен между восходящим и нисходящим потоками жидкости
2. Трехмерное численное моделирование: высокая точность, но одно моделирование требует часов, непригодно для быстрого инженерного проектирования
3. Существующие модели МТСЕ: используют предположение о горизонтально равномерной температуре поверхности (ГРВТ), что приводит к недооценке теплового отклика в среднесрочном периоде

Исследовательская мотивация

Разработка гибридного метода моделирования, сочетающего скорость МТСЕ и точность конечных элементов, для удовлетворения практических инженерных потребностей.

Основные вклады

1. Выявление и количественная оценка влияния предположения ГРВТ: путем сравнения моделирования конечными элементами с высокопроводящим слоем и без него доказано, что предположение ГРВТ недооценивает T_fm, T_in, T_out, R_b3D и R_beff
2. Разработка улучшенной модели МТСЕ: применение многоуровневой тепловой сети с узлами жидкости, трубопровода, раствора и геологического слоя повышает точность моделирования
3. Создание базы данных коэффициентов коррекции: построена база данных коэффициентов коррекции, зависящих от времени, на основе 54 конфигураций конечно-элементного моделирования
4. Реализация эффективного гибридного метода моделирования: сочетание базовых расчетов МТСЕ и интерполяционной коррекции достигает точности конечно-элементного моделирования за 2 секунды
5. Предоставление инструмента моделирования с открытым исходным кодом: программа на C++ может использоваться для быстрого и точного краткосрочного и среднесрочного моделирования БХЭ

Подробное описание методов

Определение задачи

Входные данные: геометрические параметры БХЭ (длина L, радиус r_b, расстояние между трубами s и т.д.), теплофизические параметры, расход, тепловая нагрузка Выходные данные: временная эволюция входной и выходной температур (T_in, T_out), средней температуры жидкости (T_fm), температуры поверхности БХЭ (T_b), трехмерного теплового сопротивления (R_b3D), эффективного теплового сопротивления (R_beff) Временной диапазон: от 10⁻²·⁵ до 10³ часов

Архитектура модели

1. Базовая модель конечных элементов

• Использование интерфейса теплопередачи в твердых телах COMSOL Multiphysics для моделирования теплопроводности в твердой области
• Применение модуля потока в трубах для моделирования энергетического баланса жидкости
• Проверка независимости сетки для обеспечения точности
• Масштабирование вертикальной координаты (коэффициент 20) для улучшения качества сетки

2. Улучшенная модель МТСЕ

Многоуровневая структура: БХЭ и геологический слой разделены на m горизонтальных слоев, каждый высотой l = L/m

Тепловая сеть: каждый слой содержит следующие узлы

• Узлы жидкости: f1 (нисходящий), f2 (восходящий)
• Узлы трубопровода: p1, p2
• Узлы раствора: gt1, gt2
• Узел скважины: b
• Узлы геологического слоя: g1, g2, ..., gn

Расчет теплового сопротивления:

R₁^Δ = 2R_b                    (жидкость к стенке скважины)
R₁₂^Δ = 4R_b R_a/(4R_b - R_a)   (тепловое сопротивление между жидкостями)
R_conv = 1/(2πr_pi h)           (сопротивление конвекции)

Уравнения энергетического баланса: Для каждого узла составляется нестационарное уравнение энергетического баланса, формирующее систему линейных уравнений для решения.

3. Метод коэффициентов коррекции

Определение коэффициента коррекции:

C_coeff = R_b3D,COMSOL / R_b3D,TRCM

Процесс коррекции:

1. МТСЕ вычисляет базовый тепловой отклик
2. Коэффициент коррекции получается путем интерполяции в зависимости от параметров
3. Коррекция R_b3D и пересчет температурного поля

Технические инновации

1. Количественный анализ влияния ГРВТ: путем введения высокопроводящего слоя, обеспечивающего граничное условие равномерной температуры, проведен количественный анализ влияния предположения ГРВТ на точность моделирования
2. Проектирование многоуровневой тепловой сети: по сравнению с традиционной МТСЕ добавлены узлы трубопровода и раствора для более точного описания теплопередачи внутри БХЭ
3. Стратегия апостериорной коррекции: использование преимущества МТСЕ в точном прогнозировании T_b, косвенная коррекция других параметров путем коррекции R_b3D
4. Техника многопараметрической интерполяции: парабольная интерполяция на основе четырех ключевых параметров: радиуса скважины, расстояния между трубами, теплопроводности раствора и наружного диаметра трубы

Экспериментальная установка

Конфигурация моделирования

Диапазон геометрических параметров:

• Радиус скважины r_b: 64-76 мм
• Расстояние между трубами s: 27-47 мм
• Наружный диаметр трубы r_pe: 16, 20 мм
• Теплопроводность раствора k_gt: 1,0-2,2 Вт/(м·К)

Фиксированные параметры:

• Длина БХЭ: 100 м
• Глубина залегания: 1,8 м
• Теплопроводность геологического слоя: 1,8 Вт/(м·К)
• Тепловая нагрузка: 50 Вт/м

Построение набора данных

• 54 комбинации различных конфигураций параметров БХЭ
• Для каждой комбинации проведено сравнительное моделирование COMSOL и МТСЕ
• Временные точки: 56 логарифмически равномерно распределенных точек (от 10⁻²·⁵ до 10³ часов)
• Создание базы данных коэффициентов коррекции для интерполяции

Верификационные случаи

Выбраны 3 типичные конфигурации БХЭ для верификации:

• БХЭ 1: стандартная конфигурация, используется для расчета коэффициентов коррекции
• БХЭ 2: конфигурация, требующая интерполяции параметров
• БХЭ 3: конфигурация с малым радиусом, требующая экстраполяции

Результаты экспериментов

Основные результаты

Анализ влияния предположения ГРВТ

Путем сравнения результатов моделирования с высокопроводящим слоем и без него:

• T_fm, T_in, T_out недооцениваются, T_b практически не влияется
• R_b3D и R_beff значительно недооцениваются
• Степень недооценки увеличивается со временем, особенно заметна в среднесрочном периоде (>4 часов)

Верификация эффекта коррекции

Результаты БХЭ 1:

• Среднеквадратичная ошибка T_fm: 0,037°C
• Среднеквадратичная ошибка T_in, T_out: <0,038°C
• Ошибка R_beff в конечный момент времени: 0,066%

Результаты БХЭ 2 (требуется интерполяция):

• Среднеквадратичная ошибка T_fm: 0,051°C
• Ошибка R_b3D в конечный момент времени: 0,13%
• Ошибка R_beff в конечный момент времени: 0,23%

Результаты БХЭ 3 (требуется экстраполяция):

• Среднеквадратичная ошибка T_fm: 0,052°C
• Ошибка R_b3D в конечный момент времени: 0,20%
• Ошибка R_beff в конечный момент времени: 0,40%

Вычислительная эффективность

• Скорректированная МТСЕ: завершение за 2 секунды
• Моделирование COMSOL: несколько часов
• Повышение точности: достижение точности конечно-элементного моделирования

Сравнительный анализ

Сравнение с искусственной нейронной сетью Pasquier-Marcotte:

• Оба метода имеют проблему недооценки в среднесрочном периоде
• Предложенный метод эффективно решает эту проблему через коэффициенты коррекции

Связанные работы

Аналитические методы решения

• Lamarche & Beauchamp: краткосрочное аналитическое решение, но не учитывает теплообмен между жидкостями
• Beier: аналитическая модель, учитывающая энергетический баланс жидкости, требует обратного преобразования Лапласа

Методы численного моделирования

• Li & Zheng: трехмерная твердая область + одномерный энергетический баланс жидкости
• Marcotte & Pasquier: полностью связанное моделирование трехмерной твердой области и жидкости

Методы МТСЕ

• Bauer и др.: стандартная треугольная тепловая сеть сопротивлений, использует предположение ГРВТ
• Pasquier & Marcotte: улучшенная МТСЕ и развитие в квази-трехмерную модель
• Ruiz-Calvo и др.: модель скважина-геологический слой, избегает предположения равномерной температуры, но имеет ограниченную точность

Методы машинного обучения

• Pasquier и др.: искусственная нейронная сеть, обученная на 15000 моделированиях МТСЕ

Выводы и обсуждение

Основные выводы

1. Предположение ГРВТ является фундаментальной причиной ограничения точности МТСЕ, приводит к систематической недооценке теплового отклика в среднесрочном периоде
2. Стратегия апостериорной коррекции эффективна: использование преимущества МТСЕ в точном прогнозировании T_b, достижение повышения общей точности путем коррекции R_b3D
3. Гибридный метод сочетает скорость и точность: достижение точности конечно-элементного моделирования за 2 секунды, удовлетворение требованиям инженерного применения
4. Интерполяция и экстраполяция показывают хорошие результаты: сохранение высокой точности даже при параметрах вне диапазона интерполяции

Ограничения

1. Применимо только к одиночным U-образным БХЭ: двойные U-образные или коаксиальные типы требуют переустановления базы данных коррекции
2. Ограничение диапазона параметров: параметры вне диапазона интерполяции могут влиять на точность
3. Недостаток физической интерпретации коэффициентов коррекции: отсутствует глубокий анализ связи коэффициентов коррекции с физическими механизмами
4. Долгосрочная точность не верифицирована: точность долгосрочного моделирования свыше 1000 часов требует дальнейшей проверки

Направления будущих исследований

1. Расширение на другие типы БХЭ (двойные U-образные, коаксиальные)
2. Учет сложных механизмов теплопередачи, таких как фильтрация подземных вод
3. Связь с долгосрочным методом g-функций для моделирования на всех временных масштабах
4. Разработка адаптивной стратегии коррекции для снижения зависимости от предварительно построенной базы данных

Глубокая оценка

Преимущества

1. Точное выявление проблемы: глубокий анализ физического механизма предположения ГРВТ и степени его влияния
2. Умное проектирование метода: полное использование преимуществ МТСЕ, компенсация недостатков через апостериорную коррекцию
3. Полное проектирование экспериментов: 54 конфигурации охватывают основное пространство параметров, верификация достаточна
4. Значительная инженерная ценность: предоставленная программа на C++ имеет практическую ценность применения
5. Строгая верификация точности: несколько верификационных случаев доказывают эффективность метода

Недостатки

1. Недостаточный теоретический анализ: отсутствует глубокое объяснение физического смысла коэффициентов коррекции
2. Ограниченный анализ чувствительности параметров: недостаточный анализ влияния каждого параметра на эффект коррекции
3. Ограничение области применения: применимо только к одиночным U-образным БХЭ, ограниченная универсальность
4. Отсутствие долгосрочной верификации: точность моделирования на очень длительных периодах требует проверки

Влияние

1. Академический вклад: предоставляет новую идею гибридного моделирования для области моделирования БХЭ
2. Инженерная ценность: значительно повышает эффективность и точность краткосрочного и среднесрочного моделирования БХЭ
3. Вклад в открытый исходный код: предоставленная программа на C++ способствует распространению и применению метода
4. Методологическое вдохновение: стратегия апостериорной коррекции может быть распространена на другие задачи инженерного моделирования

Сценарии применения

1. Проектирование систем геотермальных тепловых насосов: быстрая и точная оценка краткосрочной и среднесрочной тепловой производительности БХЭ
2. Анализ тепловых ответных тестов: вспомогательная интерпретация данных полевых испытаний и идентификация параметров
3. Исследование оптимизации систем: предоставление эффективного инструмента моделирования для оптимизации параметров БХЭ
4. Преподавание и научные исследования: использование в качестве стандартного инструмента для анализа теплопередачи БХЭ

Список литературы

Статья цитирует 45 соответствующих источников, охватывающих основные методы и последние достижения в моделировании БХЭ, обеспечивая прочную теоретическую основу для данного исследования. Ключевые источники включают:

• Claesson & Javed (2011): метод анализа g-функций
• Bauer и др. (2011): стандартный метод МТСЕ
• Pasquier & Marcotte (2014, 2020): улучшенные методы МТСЕ и машинного обучения
• Zanchini (2024): сравнительный анализ методов моделирования БХЭ

Общая оценка: Это высококачественная статья по прикладной физике, которая предлагает инновационный гибридный метод моделирования для практических инженерных потребностей моделирования БХЭ геотермальных тепловых насосов. Статья содержит глубокий теоретический анализ, полную экспериментальную верификацию и имеет значительную академическую ценность и перспективы инженерного применения. Предоставление программы с открытым исходным кодом дополнительно повышает ее влияние и практическую ценность.