선택한 단어 수는 입니다.
최소 단어 이상 선택하여야 합니다.
최소 단어 이상 선택하여야 합니다.
선택한 단어 수는 30입니다.
최대 10 단어까지만 선택 가능합니다.
최대 10 단어까지만 선택 가능합니다.
문제 정의
- 본 연구에서는 수원시 인근 부지에서 PHC 파일을 이용한 현장 열성능 실험을 수행하였다. 에너지 파일은 PHC 파일 내부에 PB(Polybutylene) 소재의 파이프(외경/내경: 20mm/16mm)를 W자 및 코일 형태로 결속한 후 내부에 시멘트 그라우트를 채운 형태로 구성되어 있으며 열교환기 배치 상태는 Fig.
- 또한 3차원 유한요소해석을 수행하였으며 그로부터 예측된 순환수 온도 변화를 이용하여 열교환율을 산정하고 이를 실험값과 상호 비교하였다. 수치해석과 현장 열성능 실험을 통해 PHC 에너지 파일 내 W자형과 코일형 지중 열교환기의 열효율을 분석하여 이의 적용성에 대해 논의해 보았다.
가설 설정
- 코일형 파이프의 공간적인 위치는 원통형 좌표계에서 일정한 반경 r을 갖는 대신 회전각 θ 및 깊이 z가 변화하는 3차원 특성을 갖는다. 따라서 본 연구에서는 이를 반영하여 Fig. 6과 같이 Spiral 코일형 열교환기를 Ring-코일 형태로 가정하고, 파이프-지반 경계에서 축대칭 조건에 대한 1차원 열평형방정식 식 (8)을 적용하였다(Nakayama et al., 2002).
- 에너지 파일로 유입되는 순환수 온도는 냉방 운영조건을 가정하여 30°C로 설정하였다.
제안 방법
- (1) 현장 열성능 실험을 통해 냉방 가동 조건을 가정하여 지중으로 투입되는 순환수 온도를 30도로 제어하였으며 8시간 운전, 16시간 휴지 운영 조건으로 4일간 열성능 실험을 수행하였다. 4일간 계측된 열교환율을 근거로 추세선 식을 산출하여 3개월 가동시 열교환율 값을 추정하였다.
- (1) 현장 열성능 실험을 통해 냉방 가동 조건을 가정하여 지중으로 투입되는 순환수 온도를 30도로 제어하였으며 8시간 운전, 16시간 휴지 운영 조건으로 4일간 열성능 실험을 수행하였다. 4일간 계측된 열교환율을 근거로 추세선 식을 산출하여 3개월 가동시 열교환율 값을 추정하였다. W자형의 경우 4일간 실험값을 근거로 추세선을 이용하여 산정된 3개월 동안의 파일 길이당 열교환율은 73.
- 코일형 열교환기를 이용한 4일간의 수치해석의 경우 ABAQUS 및 CFD 수치해석을 각각 수행하였으며 CFD 해석이 조금 더 실험값과 근사한 결과를 보였으나 큰 우수성을 보이지는 않았다. CFD 해석의 경우 8시간 가동 16시간 휴지의 부분운전 조건을 3개월 장기간 시뮬레이션 하기에는 시간이 과다하게 소요되어 비교 차원에서 4일동안만 시뮬레이션을 수행하였다. 코일형 열교환기 이용시 열교환율을 산정한 결과 실험과 수치해석 결과는 10% 가량의 오차가 발생하였다.
- Fig. 2와 같이 W자형과 코일형 열교환기가 설치된 에너지 파일에서 현장 열성능 실험을 수행하였다. 실제 상업용 건물의 냉난방 운영을 모사하기 위해 하루 8시간 가동-16시간 정지를 반복하여 4일간 열성능 실험을 수행하였다.
- 7에서의 수치해석에 사용된 열물성치는 Table 3과 같다. Layer 1은 지하수위 위의 건조시료를 나타내며 현장의 시료를 채취하여 Hukseflux사의 비정상 탐침(non-steady state probe) 장비인 TP-08을 이용하여 열전도도를 측정하였다. 또한 Layer 2는 지하수위 아래로써 열전도도는 Park et al.
- 일반적으로 지중열교환기 해석에 적용되는 전산유체동역학 해석코드(computational fluid dynamics code, CFD code)는 순환 유체에 대한 유체방정식과 파이프에 대한 매우 상세한 격자망을 요구하며 특히 순환 유체의 난류흐름(turbulent flow) 조건 및 코일형 파이프가 고려될 경우 계산량이 증가한다. 따라서 본 연구에서는 에너지 파일 내의 W자형과 코일형에 대해 순환유체와 주변 매질 간의 열전달 거동을 효율적으로 계산하기위한 수치해석기법을 적용하였다. 적용된 해석기법은 ABAQUS Standard에 사용자 서브루틴(user subroutine) URDFLD 및 FILM이 결합된 방법과 CFD 수치해석 프로그램인 Comsol Multiphysics로 구현하였으며 현장 열성능 실험과 비교하여 그 적절성을 살펴보았다.
- 따라서 본 연구에서는 화강풍화토로 이루어진 현장에 PHC 에너지 파일을 설치하고 W자형과 코일형 지중 열교환기를 파일 내 설치한 후 현장 열성능 실험을 수행하여 에너지 파일의 열교환율을 측정하였다. 열성능 실험은 부분운영조건(8시간 운용, 16시간 대기)으로 4일간 진행하였으며 4일간 계측된 값을 근거로 추세선 분석(regression analysis)을 통하여 3개월 동안의 열성능값을 예측하였다.
- 본 논문에서는 화강풍화토 지층으로 구성된 현장에 PHC 에너지 파일을 설치하고 W자형과 코일형 지중 열교환기를 설치한 후 현장 열성능 실험을 수행하여 열교환율을 산정하였다. 또한 3차원 유한요소해석을 수행하여 현장 실험을 그대로 모사하였다. 본 연구를 통해 다음과 같은 결론을 도출하였다.
- 본 논문에서는 화강풍화토 지층으로 구성된 현장에 PHC 에너지 파일을 설치하고 W자형과 코일형 지중 열교환기를 설치한 후 현장 열성능 실험을 수행하여 열교환율을 산정하였다. 또한 3차원 유한요소해석을 수행하여 현장 실험을 그대로 모사하였다.
- 은 순환수의 질량 흐름 속도(mass flow rate,kg/s)이고 c는 순환수의 비열 그리고 L은 파일의 길이이다. 본 연구에서 제작된 열성능 실험기는 지중의 보어홀에 매설된 지중 열교환기 내부를 순환수로 공급하고 유체의 온도변화를 감지하는 시스템으로 구성되어 있다.지반의 유효 열전도도를 측정할 수 있는 열응답 실험장치와 같이 순환수를 공급하는 수조 내에 히터를 통하여 일정한 전력을 공급할 수 있다.
- 그리고 에너지 파일의 열교환율을 산정하기 위해서는 열전달 매개체의 열 저항을 산출하여 열교환율을 산정하여야 한다. 본 연구에서는 실험을 통하여 열교환기로 유입되는 온도를 일정하게 제어한 후 유출수 온도를 측정한 후 식 (1)을 이용하여 열교환율을 산정하였다.
- 2와 같이 W자형과 코일형 열교환기가 설치된 에너지 파일에서 현장 열성능 실험을 수행하였다. 실제 상업용 건물의 냉난방 운영을 모사하기 위해 하루 8시간 가동-16시간 정지를 반복하여 4일간 열성능 실험을 수행하였다. 온도 제어 전 10분간 무부하 운전을 실시하였으며 지반의 초기 온도는 17°C 였다.
- 따라서 본 연구에서는 화강풍화토로 이루어진 현장에 PHC 에너지 파일을 설치하고 W자형과 코일형 지중 열교환기를 파일 내 설치한 후 현장 열성능 실험을 수행하여 에너지 파일의 열교환율을 측정하였다. 열성능 실험은 부분운영조건(8시간 운용, 16시간 대기)으로 4일간 진행하였으며 4일간 계측된 값을 근거로 추세선 분석(regression analysis)을 통하여 3개월 동안의 열성능값을 예측하였다. 또한 3차원 유한요소해석을 수행하였으며 그로부터 예측된 순환수 온도 변화를 이용하여 열교환율을 산정하고 이를 실험값과 상호 비교하였다.
- 절점 기반 대류성 경계조건을 실제 3차원 파이프 형상에 대해 적용하기 위하여 코일형 열교환기의 경우 3.3절의 수치해석 외에도 CFD 기반의 3차원 유한요소모델 해석 프로그램인 COMSOL Multiphysics(Verson 4.3)을 이용하여 열교환기 형상 및 배치에 따른 에너지파일의 열적 거동을 해석하였다. 대류와 전도에 의한 수치 모델의 지배 방정식은 식 (10)과 같다.
- 지반은 지하수위 위치에 따라 건조-포화 영역으로 나뉘었으며 파이프 주변의 큰 온도 구배를 고려하여 주변 영역의 요소를 조밀하게 분할하였다.
- 2와 같으며 두 파일 간 간격은 10m 이다. 코일형 열교환기의 양 끝에서의 피치 간격은 5cm 정도로 유지되었으며 W자형 열교환기는 시공 여건을 고려하여 PHC 파일 벽면에 설치하였으며 열교환기 입출구 부위에 전기식 센서인 RTD(Resistance Temperature Detector) 센서를 설치하였다. 실험지역에 대해 지반조사를 수행한 결과 풍화토와 풍화암으로 구성되어 있었다.
대상 데이터
- 그 중 풍화토는 19/30∼50/12의 N값을 가지고 추정할 때 중간내지 조밀한 상대밀도를 나타냈다. 본 현장의 지반은 실트 및 세립질과 조립질 모래로 분류되어 채취되었고 지하수위는 파일 두부에서 GL(-)4.5m 심도에 위치하였다. 또한 GL(-)11.
- 코일형 열교환기의 양 끝에서의 피치 간격은 5cm 정도로 유지되었으며 W자형 열교환기는 시공 여건을 고려하여 PHC 파일 벽면에 설치하였으며 열교환기 입출구 부위에 전기식 센서인 RTD(Resistance Temperature Detector) 센서를 설치하였다. 실험지역에 대해 지반조사를 수행한 결과 풍화토와 풍화암으로 구성되어 있었다. 그 중 풍화토는 19/30∼50/12의 N값을 가지고 추정할 때 중간내지 조밀한 상대밀도를 나타냈다.
- 본 연구에서는 수원시 인근 부지에서 PHC 파일을 이용한 현장 열성능 실험을 수행하였다. 에너지 파일은 PHC 파일 내부에 PB(Polybutylene) 소재의 파이프(외경/내경: 20mm/16mm)를 W자 및 코일 형태로 결속한 후 내부에 시멘트 그라우트를 채운 형태로 구성되어 있으며 열교환기 배치 상태는 Fig. 2와 같으며 두 파일 간 간격은 10m 이다. 코일형 열교환기의 양 끝에서의 피치 간격은 5cm 정도로 유지되었으며 W자형 열교환기는 시공 여건을 고려하여 PHC 파일 벽면에 설치하였으며 열교환기 입출구 부위에 전기식 센서인 RTD(Resistance Temperature Detector) 센서를 설치하였다.
- 7에는 2(가로)×2(세로)m 부분만 나타내었다. 해석 모델은 그라우트와 PHC 파일, 그리고 지반으로 구성되어 있으며 6절점 및 8절점 열전달요소(DC3D6 및 DC3D8)가 적용되었다. COMSOL을 이용한 코일형의 경우 Free tetrahedral 격자망이 사용되었으며 최대 요소 크기는 0.
데이터처리
- 각 일에 대한 열교환율은 8시간 가동시 얻어진 열교환율을 단순 산술 평균하여 계산된 값이다. 4일간의 계측된 열교환율을 근거로 Fig. 10과 같이 회귀분석(regression analysis)을 수행하여 추세선 식을 산출하였으며 산출된 식을 근거로 3개월 가동시 열교환율 값을 예측한 후 이를 수치 해석 결과와 비교하였다(Fig. 10). 4일간 계측된 실험값을 근거로 산출된 추세선 식의 R2값은 거의 1에 가까웠기에 본 회귀분석을 통한 3개월 동안의 열교환율 값을 예측하는데 큰 오차 없이 사용될 수 있는 것으로 판단하였다.
- 열성능 실험은 부분운영조건(8시간 운용, 16시간 대기)으로 4일간 진행하였으며 4일간 계측된 값을 근거로 추세선 분석(regression analysis)을 통하여 3개월 동안의 열성능값을 예측하였다. 또한 3차원 유한요소해석을 수행하였으며 그로부터 예측된 순환수 온도 변화를 이용하여 열교환율을 산정하고 이를 실험값과 상호 비교하였다. 수치해석과 현장 열성능 실험을 통해 PHC 에너지 파일 내 W자형과 코일형 지중 열교환기의 열효율을 분석하여 이의 적용성에 대해 논의해 보았다.
- 에너지 파일로 유입되는 순환수 온도는 냉방 운영조건을 가정하여 30°C로 설정하였다. 위의 운전 조건과 현장 조건을 그대로 시뮬레이션 하여 수치해석을 수행하였으며 이를 실험결과와 비교 분석하였다. Fig.
- 따라서 본 연구에서는 에너지 파일 내의 W자형과 코일형에 대해 순환유체와 주변 매질 간의 열전달 거동을 효율적으로 계산하기위한 수치해석기법을 적용하였다. 적용된 해석기법은 ABAQUS Standard에 사용자 서브루틴(user subroutine) URDFLD 및 FILM이 결합된 방법과 CFD 수치해석 프로그램인 Comsol Multiphysics로 구현하였으며 현장 열성능 실험과 비교하여 그 적절성을 살펴보았다.
이론/모형
- 또한 Layer 2는 지하수위 아래로써 열전도도는 Park et al.(2012)이 제안한 현장 함수비와 단위중량에 따른 화강풍화토의 열전도 도식을 이용하여 도출하였다. 시멘트 그라우트, PHC 파일, PB 파이프 그리고 순환수의 물성값은 관련 문헌을 참고하였다(Jeong et al.
- 본 연구에서는 Fig. 4에 도시된 순환수-파이프-그라우트(또는 지반) 간의 에너지 평형 관계를 파이프의 열특성이 등가적으로 고려된 절점 기반(node-based), 유체그라우트(지반)간 대류성 열전달 경계조건이 식 (4)로 고려하는 해석기법을 적용하였다(Choi et al., 2011).
- (2012)이 제안한 현장 함수비와 단위중량에 따른 화강풍화토의 열전도 도식을 이용하여 도출하였다. 시멘트 그라우트, PHC 파일, PB 파이프 그리고 순환수의 물성값은 관련 문헌을 참고하였다(Jeong et al., 2010). 여기서 시멘트 그라우트의 열전도도는 시멘트 배합조건에 따라 달라질 것으로 판단되나 본 연구에서는 시멘트 양생 후의 열전도도를 측정하지 못하여 Jeong et al.
성능/효과
- (2) 본 연구에서 적용한 3차원 유한요소 수치해석은 현장 실험 결과를 비교적 정확하게 모사하는 것을 알수 있었다. 코일형 교환기 실험과 수치해석 간에는 10% 가량의 오차가 발생하였는데 이는 실제 시공시 코일 피치가 일정하게 유지되지 못한 시공상 오차인 것으로 판단된다.
- (3) 일반적으로 수직 밀폐형의 열교환율이 40∼70W/m인 것을 감안할 때(European Geothermal Energy Council, 2008), 에너지 파일은 충분한 열교환율을 가지는 것으로 판단되며 본 현장처럼 파일의 근입 깊이가 짧게 고정된 경우 U 또는 W자형 보다는 열교환 면적이 넓은 코일형 열교환기를 적용하는 것이 바람직할 것으로 판단된다.
- 10). 4일간 계측된 실험값을 근거로 산출된 추세선 식의 R2값은 거의 1에 가까웠기에 본 회귀분석을 통한 3개월 동안의 열교환율 값을 예측하는데 큰 오차 없이 사용될 수 있는 것으로 판단하였다. Fig.
- 7W/m)였다. 또한 3개월 동안의 코일형 열교환기 이용시 실험을 통한 추정값은 84.8W/m(파이프 길이당 열교환율 : 5.5W/m), 수치해석을 통한 예측값은 94W/m(파이프 길이당 열교환율 : 6.1W/m)값을 보였다.
- 실험 및 수치해석 결과 3개월 동안 코일형 열교환기 이용시 W자형을 이용했을 때보다 파일 깊이당 열교환율은 약 10∼15% 상승되는 것을 알 수 있다.
- 12는 코일형 열교환기의 3개월 간 부분 운전시에 따른 열교환율 값을 나타낸다. 실험 및 수치해석 결과는 상당히 일치하였으며 이는 수치해석 모델이 열교환율 값을 비교적 정확하게 모사하는 것을 알 수 있다. 코일형 열교환기를 이용한 4일간의 수치해석의 경우 ABAQUS 및 CFD 수치해석을 각각 수행하였으며 CFD 해석이 조금 더 실험값과 근사한 결과를 보였으나 큰 우수성을 보이지는 않았다.
- 하지만 파이프 길이당 열교환율은 반대로 W자형이 코일형에 비해 350∼400% 높음을 알 수 있다. 전체 시공비 중에 열교환기 파이프 비용이 차지하는 비중은 극히 미미하기에 경제성을 고려해도 코일형 열교환기의 열교환 성능이 W자형에 비해 우수할 것으로 판단된다. 이는 코일형 열교환기 적용시 그라우트 및 지반과의 넓은 열교환 면적으로 순환수가 머무르는 시간이 길기 때문에 지반으로의 방출열량이 크기 때문인 것으로 판단된다.
- 또한 히터 내의 온도를 일정하게 제어할 수 있는 온도 제어기가 설치되어 있기에 열교환기 유입수 온도를 일정하게 유지하고 유출수 온도를 계측하면 지반과의 열교환율을 측정할 수 있는 기능도 탑재되어 있다. 즉 본 장비를 이용하면 지반의 열전도도를 측정할 수 있는 열응답 실험과 열효율을 측정할 수 있는 열성능 실험이 가능하다고 할 수 있다. Table 2는 열성능 실험기의 제원을 나타낸다.
- CFD 해석의 경우 8시간 가동 16시간 휴지의 부분운전 조건을 3개월 장기간 시뮬레이션 하기에는 시간이 과다하게 소요되어 비교 차원에서 4일동안만 시뮬레이션을 수행하였다. 코일형 열교환기 이용시 열교환율을 산정한 결과 실험과 수치해석 결과는 10% 가량의 오차가 발생하였다. 이는 수치해석시 코일형 열교환기의 모델링에 대한 오차가 어느정도 발생한 것으로 판단되며 또한 코일형 열교환기의 경우 일반 W자형에 비해 수두 손실이 많이 발생되기에 실험시 유량이 일정하게 제어되지 못한 점 등에 의한 차이인 것으로 판단된다.
- 코일형 열교환기 적용시 W자형을 이용했을 때보다 파일 깊이당 열 교환율은 10∼15%가량 열교환율이 상승되는 것을 알 수 있다.
- 실험 및 수치해석 결과는 상당히 일치하였으며 이는 수치해석 모델이 열교환율 값을 비교적 정확하게 모사하는 것을 알 수 있다. 코일형 열교환기를 이용한 4일간의 수치해석의 경우 ABAQUS 및 CFD 수치해석을 각각 수행하였으며 CFD 해석이 조금 더 실험값과 근사한 결과를 보였으나 큰 우수성을 보이지는 않았다. CFD 해석의 경우 8시간 가동 16시간 휴지의 부분운전 조건을 3개월 장기간 시뮬레이션 하기에는 시간이 과다하게 소요되어 비교 차원에서 4일동안만 시뮬레이션을 수행하였다.
- 실험 및 수치해석 결과는 상당히 일치하였으며 이는 수치해석 모델이 열교환율 값을 비교적 정확하게 모사하는 것을 알 수 있다. 코일형 열교환기를 이용한 4일간의 수치해석의 경우 ABAQUS 및 CFD 수치해석을 각각 수행하였으며 CFD 해석이 조금 더 실험값과 근사한 결과를 보였으나 큰 우수성을 보이지는 않았다. CFD 해석의 경우 8시간 가동 16시간 휴지의 부분운전 조건을 3개월 장기간 시뮬레이션 하기에는 시간이 과다하게 소요되어 비교 차원에서 4일동안만 시뮬레이션을 수행하였다.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.