[논문 리뷰] FastCTF: A Robust Solver for Conduction Transfer Function Coefficients and Thermal Response Factors

Khodr Jaber|arXiv (Cornell University)|2021. 03. 30.

Building Energy and Comfort Optimization참고 문헌 29인용 수 1

한 줄 요약

이 논문은 s-도메인에서의 멱급수 전개에 대한 Padé 근사법을 사용하여 열전도 전이함수(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-도메인의 급수 전개와 Padé 근사법을 활용하여 해결하는 것.
• 기존의 FDR 및 DRF와 같은 방법보다 우수한 성능을 보이며 다양한 벽 구성과 시간 간격에서 작동하는 오픈소스이자 계산적으로 경량화된 구현을 제공하는 것.
• 이론적 주파수 특성과 ASHRAE 핸드북 및 FDR의 공개된 자료를 기반으로 방법을 검증하는 것.
• 다양한 시간 이산화 수준에서 경량 및 중량 벽 모두에 대해 높은 정확도와 안정성을 확보하는 것.

제안 방법

• 이 방법은 복소 s-도메인에서 전이함수 해의 멱급수 전개를 사용하여 다층 벽에서의 일시적 열전도를 모델링한다.
• 이 급수에서 유리 전이함수를 근사하기 위해 Padé 근사법을 구성하여 역 라플라스 변환 과정을 단순화한다.
• Z-변환을 통해 z-전이함수를 유도함으로써 이산 컨볼루션을 이용한 효율적인 시간 도메인 시뮬레이션을 가능하게 한다.
• 알고리즘은 오픈소스 C++ 코드로 구현되었으며, 선형 대수 연산을 위해 두 개의 FORTRAN 루틴을 활용한다.
• Padé 근사법에서 유도된 재귀 관계를 통해 CTF 계수 {b_k}, {c_k}, {d_k}와 응답 인자 {X_k}, {Y_k}, {Z_k}를 계산한다.
• 대규모 선형 방정식계의 해법이나 반복적 회귀를 피함으로써 주파수 도메인 회귀나 직접 근 찾기 방법보다 더 빠르고 안정적인 계산이 가능하다.

실험 결과

연구 질문

• RQ1 s-도메인 급수 전개의 Padé 근사법은 정확한 이론적 해와 비교해 볼 때 CTF 계수와 열응답 인자를 얼마나 정확하게 계산할 수 있는가?
• RQ2 특히 전통적 솔버가 도전적인 작은 시간 간격에서도 이 방법의 성능은 어떠한가?
• RQ3 경량 및 중량 벽 모두에서 기존의 FDR 및 DRF와 비교해 FastCTF 방법의 정확도와 안정성은 어떠한가?
• RQ4 열전도도, 밀도 등의 매개변수 불확실성이 계산된 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++ 구현은 계산적으로 효율적이며, 선형 대수 연산을 위해 오직 두 개의 FORTRAN 루틴만을 사용하여 이식성과 통합성을 향상시켰다.