단일 U자형 보어홀 열교환기의 초단기 및 중기 시뮬레이션을 위한 열저항-용량 및 유한요소 결합 모델

기본 정보

• 논문 ID: 2510.14421
• 제목: 단일 U자형 보어홀 열교환기의 초단기 및 중기 시뮬레이션을 위한 매우 빠르고 정확한 열저항-용량 및 유한요소 결합 모델
• 저자: Enzo Zanchini, Francesco Zanchini, Claudia Naldi (볼로냐 대학교)
• 분류: physics.app-ph
• 발표 시간: 2024년
• 논문 링크: https://arxiv.org/abs/2510.14421

초록

지열 열펌프 시스템의 정확한 설계는 보어홀 열교환기(BHE)의 단기 및 장기 출구 유체 온도 이해가 필요합니다. 본 논문은 초단기 및 중기 분석에 중점을 둡니다. 이 시간 범위 내에서 3D 유한요소 시뮬레이션 또는 열저항-용량 모델(TRCM)을 채택할 수 있습니다. 전자는 높은 정확도를 제공하지만 계산 시간이 오래 걸리고, 후자는 빠르지만 단순화된 가정으로 인해 정확도가 제한됩니다. 본 논문은 TRCM의 속도와 유한요소 시뮬레이션의 정확도를 결합한 새로운 방법을 제시하며, 이는 단일 U자형 보어홀 열교환기의 초단기 및 중기 시뮬레이션에 사용됩니다. 이 방법은 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. 사후 수정 전략: TRCM이 T_b를 정확히 예측하는 특성을 활용하여 R_b3D 수정을 통해 간접적으로 다른 매개변수 수정
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 등: 15,000개 TRCM 시뮬레이션으로 훈련한 인공신경망

결론 및 토의

주요 결론

1. HUST 가정은 TRCM 정확도 제한의 근본 원인으로, 중기 열응답의 체계적 저평가 초래
2. 사후 수정 전략이 효과적: TRCM이 T_b를 정확히 예측하는 장점을 활용하여 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 열전달 분석의 표준 도구로 활용

참고문헌

논문은 BHE 시뮬레이션의 주요 방법과 최신 진전을 포괄하는 45편의 관련 문헌을 인용하여 본 연구에 견고한 이론적 기초를 제공합니다. 핵심 문헌은 다음을 포함합니다:

• Claesson & Javed (2011): g 함수 분석 방법
• Bauer 등 (2011): TRCM 표준 방법
• Pasquier & Marcotte (2014, 2020): 개선된 TRCM 및 기계학습 방법
• Zanchini (2024): BHE 시뮬레이션 방법 비교 분석

종합 평가: 이는 지열 열펌프 BHE 시뮬레이션의 실제 공학 요구에 대응하여 혁신적인 혼합 모델링 방법을 제시한 고품질의 응용물리 논문입니다. 논문은 이론 분석이 심층적이고 실험 검증이 충분하며, 현저한 학술 가치와 공학 응용 전망을 보유합니다. 오픈소스 프로그램의 제공은 그 영향력과 실용 가치를 더욱 향상시킵니다.