単一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

要約

地中熱ヒートポンプシステムの正確な設計には、ボーリング熱交換器（BHE）の短期および長期における出口流体温度の理解が必要である。本論文は短期および中期分析に焦点を当てている。この時間範囲では、3D有限要素シミュレーションまたは熱抵抗容量モデル（TRCM）を採用できる。前者は精度が高いが計算に時間がかかり、後者は高速だが簡略化の仮定により精度が限定される。本論文は、単一U字管ボーリング熱交換器の短期および中期シミュレーションのために、TRCMの速度と有限要素シミュレーションの精度を組み合わせた新しい方法を提案する。本方法はTRCMを使用してBHEの熱応答を推定し、その後54組の異なる構成の有限要素シミュレーションから生成されたデータセットによる補間修正を行う。モデルはC++プログラムとして実装され、ボローニャ大学のオープンソースデータベースで提供されており、2秒以内に入出口および平均流体温度、平均BHE表面温度、3D および有効ボーリング熱抵抗の時間発展を提供できる。

研究背景と動機

問題記述

地中熱ヒートポンプシステムは建物の暖冷房のための効率的な技術として、その設計にはボーリング熱交換器の熱応答特性の正確な予測が必要である。主な課題は以下の通りである：

1. 計算精度と効率の矛盾：3D有限要素シミュレーションは精度が高いが数時間要する；TRCMは高速だが精度が限定される
2. 短期および中期シミュレーション需要：既存方法は10⁻²·⁵時間から10³時間の時間範囲において精度または効率に問題がある
3. 実務的工学応用需要：システム設計と最適化のための高速かつ正確なシミュレーションツールが必要である

既存方法の限界

1. 解析解法：流体を発熱固体または表面に置き換え、上下行流体間の熱交換効果を考慮できない
2. 3D数値シミュレーション：精度は高いが単一シミュレーションに数時間要し、工学的な迅速な設計には不適切である
3. 既存TRCMモデル：水平均一表面温度（HUST）仮定を採用し、中期熱応答の過小評価につながる

研究動機

TRCMの速度と有限要素精度の両方を備えたハイブリッドシミュレーション方法を開発し、実務的工学需要を満たすこと。

核心的貢献

1. HUST仮定の影響の発見と定量化：高導熱層の有無による有限要素シミュレーション比較を通じて、HUST仮定がT_fm、T_in、T_out、R_b3D、およびR_beffを過小評価することを証明
2. 改善されたTRCMモデルの開発：流体、管、グラウト、および地層ノードを含む階層化熱路ネットワークを採用し、シミュレーション精度を向上
3. 修正係数データベースの構築：54組の異なる構成の有限要素シミュレーションを通じて、時間依存の修正係数データセットを構築
4. 効率的なハイブリッドシミュレーション方法の実装：TRCM基本計算と補間修正を組み合わせ、2秒以内に有限要素シミュレーション精度を達成
5. オープンソースシミュレーションツールの提供：C++プログラムはBHEの高速かつ正確な短中期シミュレーションに使用可能

方法の詳細

タスク定義

入力：BHE幾何学パラメータ（長さL、半径r_b、管間距離s等）、熱物性パラメータ、流量、熱負荷 出力：時間発展する入出口温度（T_in、T_out）、平均流体温度（T_fm）、BHE表面温度（T_b）、3D熱抵抗（R_b3D）、有効熱抵抗（R_beff） 時間範囲：10⁻²·⁵時間から10³時間

モデルアーキテクチャ

1. 有限要素基準モデル

• COMSOL Multiphysicsの伝熱固体インターフェースを使用して固体領域の熱伝導をシミュレート
• 管流モジュールを使用して流体エネルギーバランスをシミュレート
• メッシュ独立性検証により精度を確保
• 垂直座標スケーリング（係数20）によりメッシュ品質を向上

2. 改善されたTRCMモデル

階層化構造：BHEと地層を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. TRCMで基本熱応答を計算
2. パラメータに基づいて補間により修正係数を取得
3. R_b3Dを修正し温度場を再計算

技術的革新点

1. HUST影響の定量分析：高導熱層を導入して等温境界を強制し、HUST仮定がシミュレーション精度に与える影響を定量分析
2. 階層化熱路ネットワーク設計：従来のTRCMと比較して、管およびグラウトノードを追加し、BHE内部の伝熱をより正確に記述
3. 事後修正戦略：T_bの正確な予測というTRCMの特性を利用し、R_b3Dの修正を通じて他のパラメータを間接的に修正
4. 多パラメータ補間技術：ボーリング半径、管間距離、グラウト熱伝導率、管外径の4つの主要パラメータに基づいて放物線補間を実施

実験設定

シミュレーション構成

幾何学パラメータ範囲：

• ボーリング半径r_b：64～76 mm
• 管間距離s：27～47 mm
• 管外径r_pe：16、20 mm
• グラウト熱伝導率k_gt：1.0～2.2 W/(m·K)

固定パラメータ：

• BHE長さ：100 m
• 埋設深さ：1.8 m
• 地層熱伝導率：1.8 W/(m·K)
• 熱負荷：50 W/m

データセット構築

• 54組の異なるBHEパラメータ組み合わせ
• 各組についてCOMSOLおよびTRCM比較シミュレーションを実施
• 時間点：56個の対数等間隔点（10⁻²·⁵から10³時間）
• 補間用の修正係数データベースを生成

検証ケース

3つの典型的なBHE構成を選択して検証：

• BHE 1：標準構成、修正係数計算用
• BHE 2：パラメータが補間を必要とする構成
• BHE 3：外挿が必要な小半径構成

実験結果

主要結果

HUST仮定影響分析

高導熱層の有無によるシミュレーション結果の比較：

• T_fm、T_in、T_outが過小評価され、T_bはほぼ影響なし
• R_b3DおよびR_beffが著しく過小評価
• 過小評価の程度は時間とともに増加し、中期（>4時間）で特に顕著

修正効果検証

BHE 1結果：

• T_fm二乗平均平方根誤差：0.037°C
• T_in、T_out二乗平均平方根誤差：<0.038°C
• R_beff最終時刻誤差：0.066%

BHE 2結果（補間が必要）：

• T_fm二乗平均平方根誤差：0.051°C
• R_b3D最終時刻誤差：0.13%
• R_beff最終時刻誤差：0.23%

BHE 3結果（外挿が必要）：

• T_fm二乗平均平方根誤差：0.052°C
• R_b3D最終時刻誤差：0.20%
• R_beff最終時刻誤差：0.40%

計算効率

• 修正TRCM：2秒以内に完了
• COMSOLシミュレーション：数時間
• 精度向上：有限要素シミュレーション精度を達成

比較分析

Pasquier-Marcotteニューラルネットワークとの比較：

• 両方法とも中期低估問題が存在
• 本論文の方法は修正係数を通じてこの問題を効果的に解決

関連研究

解析解法

• Lamarche & Beauchamp：短期解析解、ただし流体間熱交換を考慮できない
• Beier：流体エネルギーバランスを考慮した解析モデル、ラプラス逆変換が必要

数値シミュレーション法

• Li & Zheng：3D固体領域+1D流体エネルギーバランス
• Marcotte & Pasquier：3D固体および流体領域の完全結合シミュレーション

TRCM法

• Bauer等：標準三角熱抵抗ネットワーク、HUST仮定を採用
• Pasquier & Marcotte：改善されたTRCMおよび準3Dモデルの開発
• Ruiz-Calvo等：ボーリング-地層モデル、等温仮定を回避するが精度は限定的

機械学習法

• Pasquier等：15000組のTRCMシミュレーションで訓練されたニューラルネットワーク

結論と考察

主要結論

1. HUST仮定はTRCM精度制限の根本原因であり、中期熱応答の体系的な過小評価につながる
2. 事後修正戦略は有効：T_bの正確な予測というTRCMの利点を活用し、R_b3D修正を通じて全体精度を向上
3. ハイブリッド方法は速度と精度の両立：2秒以内に有限要素シミュレーション精度を達成し、工学応用需要を満たす
4. 補間外挿効果は良好：パラメータ範囲外でも高精度を維持

限界

1. 単一U字管BHEのみに適用：双U字管または同軸型は修正データベースの再構築が必要
2. パラメータ範囲の制限：補間範囲外のパラメータは精度に影響する可能性
3. 修正係数の物理的解釈の欠如：修正係数と物理機構の深層的関連分析が不足
4. 長期精度の未検証：1000時間を超える長期シミュレーション精度は検証が必要

今後の方向性

1. 他のタイプのBHE（双U字管、同軸型）への拡張
2. 地下水浸透流など複雑な伝熱機構の考慮
3. 長期g関数法との結合による全時間スケールシミュレーションの実現
4. 事前構築データベースへの依存を減らすための適応的修正戦略の開発

深層評価

利点

1. 問題識別の正確性：HUST仮定の物理機構と影響程度を深く分析
2. 方法設計の巧妙さ：TRCMの利点を十分に活用し、事後修正により不足を補完
3. 実験設計の完備性：54組の構成が主要パラメータ空間をカバーし、検証が充分
4. 工学的価値の顕著性：オープンソースC++プログラムは実務的応用価値を有する
5. 精度検証の厳密性：複数の検証ケースが方法の有効性を証明

不足

1. 理論分析の不足：修正係数の物理的意味の深い説明が欠如
2. パラメータ感度分析の限定：各パラメータが修正効果に与える影響の充分な分析が不足
3. 適用範囲の制限：単一U字管BHEのみを対象とし、汎用性が限定的
4. 長期検証の欠失：超長期シミュレーション精度の検証が不足

影響力

1. 学術的貢献：BHEシミュレーション分野に新しいハイブリッド建模思想を提供
2. 工学的価値：短中期BHEシミュレーションの効率と精度を著しく向上
3. オープンソース貢献：提供されたC++プログラムは方法の推進と応用を促進
4. 方法論的啓発：事後修正戦略は他の工学シミュレーション問題に推広可能

適用シーン

1. 地中熱ヒートポンプシステム設計：BHEの短中期熱性能を迅速かつ正確に評価
2. 熱応答試験分析：現場試験データの解釈とパラメータ識別を補助
3. システム最適化研究：BHEパラメータ最適化のための効率的シミュレーションツールを提供
4. 教育研究：BHE伝熱分析の標準ツールとして機能

参考文献

論文は45篇の関連文献を引用しており、BHEシミュレーションの主要方法と最新進展をカバーし、本研究に堅実な理論基礎を提供している。主要文献には以下が含まれる：

• Claesson & Javed (2011)：g関数分析法
• Bauer et al. (2011)：TRCM標準法
• Pasquier & Marcotte (2014, 2020)：改善されたTRCMおよび機械学習法
• Zanchini (2024)：BHEシミュレーション方法比較分析

総合評価：これは応用物理分野の高品質論文であり、地中熱ヒートポンプBHEシミュレーションの実務的工学需要に対して、革新的なハイブリッド建模方法を提案している。論文の理論分析は深く、実験検証は充分であり、顕著な学術価値と工学応用前景を有する。オープンソースプログラムの提供はその影響力と実用価値をさらに高めている。