[論文レビュー] FastCTF: A Robust Solver for Conduction Transfer Function Coefficients and Thermal Response Factors

Khodr Jaber|arXiv (Cornell University)|Mar 30, 2021

Building Energy and Comfort Optimization参考文献 29被引用数 1

ひとこと要約

本稿では、sドメインにおける級数展開のパデ近似を用いて、伝導伝熱関数（CTF）係数および熱的応答係数を計算するオープンソースのC++ソルバー、FastCTFを提示する。この手法により、時間刻みや壁の重みが異なる状況下でも、一貫性があり、安定的かつ効率的な一時的熱伝導の時領域シミュレーションが可能となり、特に小さな時間刻みにおいて、FDRなどの既存手法よりも安定性と精度に優れる。

ABSTRACT

Conduction transfer functions (CTF) are commonly used in the building services to quickly estimate hourly conduction heat loads through multilayered walls without resorting to expensive, time-consuming solutions of the heat equation. It is essential for any software developed for this purpose to be able to simulate walls of varying weight with a high degree of accuracy. A robust algorithm for computing CTF coefficients and thermal response factors based on power series expansions of solutions of the governing equations in the complex s-domain is presented and validated. These series expansions are used to to construct Pad\'e approximants of the system's transfer functions, which greatly simplifies the inversion of the solution from the complex domain to the time domain, and allows for an easy recovery of a time series representation via the Z-transform. The algorithm is also implemented in an open-source C++ code. Its performance is validated with respect to exact theoretical frequency characteristics and its results are compared with data generated by previously established methods for computing CTF coefficients / response factors.

研究の動機と目的

• 多層壁におけるCTF係数および熱的応答係数を計算する、強固で正確かつ効率的なアルゴリズムの開発。
• 既存手法の限界（特に小さな時間刻みにおける不安定性や高い計算コスト）を、sドメインにおける級数展開とパデ近似を活用することで克服する。
• FDR や DRF といった既存手法を凌駕する、オープンソースで計算コストが低い実装を提供し、多様な壁構成および時間刻みで優れた性能を発揮する。
• 理論的周波数特性およびASHRAEハンドブックやFDRから発表されたデータと照合して、本手法の妥当性を検証する。
• 時間離散化が異なる状況下でも、軽量壁および重厚壁の両方において、高い精度と安定性を維持することを保証する。

提案手法

• 本手法は、複素sドメインにおける伝達関数解の級数展開を用いて、多層壁における一時的熱伝導をモデル化する。
• これらの級数からパデ近似を構築し、有理関数伝達関数を近似することで、逆ラプラス変換のプロセスを簡略化する。
• z伝達関数はZ変換を用いて導出され、離散畳み込みを用いた効率的な時領域シミュレーションを可能にする。
• アルゴリズムはオープンソースのC++コードとして実装されており、線形代数演算には2つのFORTRANルーチンに依存している。
• パデ近似から得られる再帰関係を用いて、CTF係数 {b_k}、{c_k}、{d_k} および応答係数 {X_k}、{Y_k}、{Z_k} を計算する。
• 大規模な連立一次方程式の解法や反復回帰を避けることで、周波数ドメイン回帰や直接根の探索よりも高速かつ安定した計算が可能となる。

実験結果

リサーチクエスチョン

• RQ1 sドメイン級数展開のパデ近似は、正確な理論的解と比較して、CTF係数および熱的応答係数をどの程度正確に計算できるか？
• RQ2 特に従来のソルバーが困難とする小さな時間刻みにおいて、本手法はどのような性能を示すか？
• RQ3 軽量壁および重厚壁の両方において、FDR や DRF といった既存手法と比較して、本手法の精度と安定性はどの程度高いか？
• RQ4 熱伝導率や密度などのパrameterの不確実性が、計算されたCTF係数に与える影響はどの程度で、本手法はその誤差をどのように軽減するか？
• RQ5 異なる熱的質量と時間刻みサイズを有する壁のシミュレーションにおいて、本手法は高精度かつ強健性を維持できるか？

主な発見

• 時間刻みΔt = 3600 sではL2誤差が1.174%に留まり、Δt = 60 sでは3.457e-4%にまで低下し、時間刻みにわたる高い精度を示した。
• 60 s未満の時間刻み、例えばΔt = 10 sでは誤差が1.819%に上昇し、数値的不安定性が顕在化する閾値を示した。
• FDR や DRF よりも小さな時間刻みにおいて優れた安定性を示し、周波数特性は理論値と非常に近い。
• FastCTF係数から算出されたU値は、壁群の真のU値と一致しており、エネルギーバランスの正確性を確認した。
• FDRとほぼ同一のz伝達関数周波数応答を生成し、b_k係数にわずかな差異を除いて一貫性と強健性を示した。
• オープンソースのC++実装は計算的に効率的であり、線形代数処理にわずか2つのFORTRANルーチンのみを用いるため、移植性と統合性が向上している。