Mobile QR Code QR CODE
Export citation EndNote

  1. 국립군산대학교 토목공학과 박사과정 (Graduate Student, Department of Civil Engineering, Kunsan National University, Kunsan 54150, Rep. of Korea)
  2. 국립군산대학교 토목공학과 석사과정 (Graduate Student, Department of Civil Engineering, Kunsan National University, Kunsan 54150, Rep. of Korea)
  3. 국립군산대학교 토목공학과 교수 (Professor, Department of Civil Engineering, Kunsan National University, Kunsan 54150, Rep. of Korea)



graphite, concrete pebbles, thermal energy storage, thermal cycling
그라파이트, 콘크리트 페블, 열에너지 저장, 열사이클

1. 서 론

현재 전 세계 여러 국가에서는 에너지 소비의 절반 이상이 건물의 공간 난방에 사용되고 있다. 이러한 높은 난방 에너지 수요로 인한 환경적 문제를 해결하기 위한 방안으로써 히트 펌프(heat pump) 활용 기술이 주목받고 있다. 히트 펌프는 이산화탄소 배출을 줄이는 동시에 저탄소 에너지의 활용 효율을 높일 수 있는 기술로써 기존의 화석연료 연소 기반 난방 시스템을 대체할 수 있는 효과적인 대안으로 평가되고 있다(Adamson et al. 2022; Jaiswal et al. 2022).

히트펌프는 산업 현장에서 발생하는 폐열, 지열, 태양열 등 다양한 열원을 활용하여 기존 시스템보다 더 높은 온도로 열매체를 가열하여 저장함으로써 에너지 밀도를 높이고, 산업 공정에서 필요한 열을 공급하는 장치이다(Kannan and Vakeesan 2016). 이러한 방식은 화석연료의 직접적인 사용을 줄이고, 그 결과로 탄소 배출량을 효과적으로 감소시킬 수 있다(Zhou et al. 2026).

그러나 히트펌프에 활용되는 열원들은 일정한 시간 동안만 생산되는 간헐적 특성을 지니고 있어, 열을 효율적으로 저장하고 필요할 때 공급하기 위한 기술적 대안이 필요하다. 이에 따라 열에너지 저장(thermal energy storage, TES) 시스템의 개발과 적용이 활발히 이루어지고 있다. TES 시스템은 특정 시간 동안 축적된 열에너지를 저장해 두었다가 필요 시 방출함으로써, 에너지 공급의 불균형 문제를 보완하고 전체 시스템의 안정성을 향상시키는 역할을 한다(Mawire et al. 2010; Park and Yang 2023; Wullhorst et al. 2024; Georgousis et al. 2025).

TES 시스템 설계는 주로 현열 저장, 잠열 저장, 열화학 저장 등 세 가지 방식으로 이루어지고 있다(Kim et al. 2013; Alnaimat and Rashid 2019). 이 중 현열 저장 방식인 팩-베드(packed-bed) 열저장 시스템은 일정한 부피의 저장 용기 내부에 고체 형태의 열저장 매체를 채워 넣고, 열저장 매체 사이를 고온의 열전달 유체가 통과하면서 열을 저장하거나 저온의 열전달 유체를 통해 열을 방출하는 방식으로 작동한다(Laing and Zunft 2015; Erregueragui et al. 2016; Gautam and Saini 2020; Tiwari and Saini 2020). 이러한 시스템은 잠열 및 화학적 방식에 비해 구조가 단순하여 설계 용이성이 높으며, 대규모 열에너지 저장에 적합하다는 특징을 가진다(Strasser and Selvam 2014; Abera et al. 2024). 또한, 열저장 매체로 활용되는 콘크리트 페블은 기존 세라믹 매체 대비 재료비가 약 1/10 수준으로 저렴하며, 국부적인 파손 발생 시 해당 페블만 교체 가능하여 유지관리 비용 측면에서도 경제적 이점이 크다(Singh et al. 2013; Berrhazi et al. 2020).

여러 선행 연구에서 저장층 내부의 열전달 성능을 향상시키기 위한 다양한 접근이 시도되어 왔다(Rahjoo et al. 2023; Tagle-Salazar et al. 2026). Kruger et al.(2018)은 고온 현열 저장 시스템에서 콘크리트 기반 고체 저장재의 열전도도와 열확산 특성이 팩-베드 내부의 온도 분포 및 축・방열 효율에 미치는 영향을 분석하여, 열저장 소재의 열전도도가 주요 인자임을 확인하였다. 열저장 소재의 열전도도가 낮을 경우 열저장 시스템 내부에서 큰 온도 구배가 형성되어 열저장 효율이 저하되었다(Zavattoni et al. 2014; Diaz-Heras et al. 2020). 또한, Shah and Lee(2011)는 팩-베드 입자 배열과 유동 경로의 최적화를 통해 저장층 전체의 열전달 성능이 개선될 수 있음을 제시하였고, Abedin and Mahmud(2014)은 저장재 입자 크기 및 형상의 변화가 열전달 거동에 중요한 영향을 미친다는 사실을 제시하였다(Mawire and Taole 2011; Panchal et al. 2020; Martin-Alcantara and Fernandez-Feria 2025).

한편, 콘크리트의 열전달 성능을 개선하기 위한 다양한 방안 중 탄소 기반의 고전도성 재료인 그라파이트가 주목받고 있다. Chen et al.(2024)은 그라파이트를 콘크리트에 혼입함으로써 콘크리트 내부에 연속적인 열전달 경로가 형성되고, 이에 따라 일반 OPC 콘크리트 대비 열전도도가 유의하게 증가함을 실험으로 확인하였다. Rao et al.(2019)은 그라파이트의 높은 고유 열전도도를 바탕으로, 그라파이트를 포함한 콘크리트 복합체의 열전도도가 일반 콘크리트에 비해 향상됨을 이론적 모델과 실험 결과를 통해 제시하였다. 반면에, 팩-베드 열저장 시스템에서 그라파이트를 혼입한 콘크리트를 활용한 열전달 특성과 열저장 성능에 관한 연구는 미비한 상황이다.

따라서, 본 연구에서는 팩-베드 시스템을 모사한 열저장 시스템에서 그라파이트 혼입 콘크리트 페블의 열전달 및 열저장 특성에 관한 실험연구를 수행하였다. 열전달 유체의 질량속도와 열사이클 시나리오의 온도상승 속도를 주요 실험변수로 설정하여 열전달 실험을 수행하였다. 열전달 실험을 통해 팩-베드 형태의 열저장 용기 내부의 콘크리트 온도 분포를 측정하여 온도 경사와 열전달률을 파악하고, 이를 바탕으로 열저장 성능을 분석하였다.

2. 열저장 매체 제작 및 열전달 실험

2.1 콘크리트 배합설계 및 열저장 매체 제작

시멘트, 그라파이트, 폴리프로필렌(polypropylene, PP) 섬유, 잔골재를 사용하여 Table 1과 같이 물-바인더 비(W/B)가 0.45인 콘크리트 배합을 설계하고 그라파이트 콘크리트 페블을 제작하였다. 선행 연구(Jeong et al. 2025)를 통해 열전도도 향상과 강도 발현을 위해 그라파이트를 시멘트 중량 대비 10 %로 치환하여 제작하였다. 바인더로 보통포틀랜드시멘트(OPC)를 사용하였으며, 시멘트는 수화를 통해 높은 압축강도를 발현하여 콘크리트의 기초 구조를 형성한다.

콘크리트의 열적 성능 향상을 위해 사용된 그라파이트(graphite)는 99 % 이상의 탄소(C) 성분으로 이루어져 있다. 해당 그라파이트의 평균 입자 크기는 약 50∼100 μm 수준으로, 매우 고운 분말 형태를 띠고 있어 시멘트 페이스트 내에 균질하게 분산될 수 있는 크기를 선정하였다. 그라파이트의 밀도는 약 2.28 g/cm3이며, 입자 자체의 열전도도는 약 100∼200 W/m·K 수준으로 높아 콘크리트 내부에서 열전달 경로를 형성하는 역할을 한다. 이에 따라 그라파이트는 높은 열적 안정성, 열전도도 및 비열을 가지고 있어 열저장 성능의 핵심 요소로 작용한다(Yuan et al. 2012; Warkhade et al. 2016). P.P 섬유는 고온 노출 시 용융되어 그 자리에 공극을 형성하여 폭렬을 예방하는 역할을 하며, 지름 21 μm, 비중 0.91, 길이 12.7 mm이며, 인장강도는 400 MPa인 P.P 섬유를 2.0 kg/m3 사용하였다(Han et al. 2005; Yang and Park 2019). 잔골재는 단위중량을 높이고 강도와 안정성을 높이며 성형 과정의 효율성 향상 및 워커빌리티 개선을 고려하였다. 또한, 콘크리트 제조 시 그라파이트의 미세입자 특성으로 인해 발생할 수 있는 유동성 저하를 방지하기 위하여 AE 고성능 감수제를 바인더 중량 대비 1.0 % 첨가하여 워커빌리티를 확보하였다.

본 연구에서 콘크리트 페블 제작에 사용된 잔골재의 입도 분포 특성은 Table 2Fig. 1에 나타냈다. 체가름 시험 결과, 잔골재의 조립률(Fineness Modulus, F.M.)은 3.08로 계산되어 잔골재 조립률 표준 범위인 2.3∼3.1 사이의 값을 만족하는 것으로 나타났다. 또한, Fig. 1에 나타낸 입도분포 곡선이 KS F 2527 (KATS 2022)에서 규정하는 표준 입도분포의 상・하한선 내에 안정적으로 위치하고 있어, 실험에 사용된 잔골재가 적절한 입도 분포를 갖는 것을 확인하였다.

Table 1 Mixing proportions (kg/m3)

W/B Water OPC Graphite Fine agg. Water reducing agent (%) P.P fiber
0.45 360 720 80 1,000 1.0 2.0

Table 2 Particle size distribution of the fine aggregate

Nominal diameter (mm) Cumulative remaining ratio (%) Passing ratio (%)
10 0.00 100
5 1.45 98.55
2.5 17.51 82.49
1.2 41.84 58.16
0.6 65.61 34.39
0.3 85.52 14.48
0.15 96.30 3.70
Pan 100.00 0.00
Fineness Modulus (F.M.) 3.08

Fig. 1 Particle size distribution of the fine aggregate

../../Resources/KCI/JKCI.2026.38.4.463/fig1.png

직경 30 mm의 구형 몰드를 활용하여 Fig. 2와 같이 콘크리트 페블을 제작하였다. 몰드 상단의 주입구에서 콘크리트를 주입한 후 몰드 하부까지 원활하게 채워질 수 있도록 하였다. Fig. 2(b)에 나타낸 바와 같이 콘크리트 페블 내부의 온도 변화를 파악하기 위해 K-type 온도 센서를 삽입한 페블을 제작하였다. 콘크리트 페블 타설 전 거푸집 내부에 센서를 미리 배치하여 콘크리트 내에 온도센서가 삽입되는 방식으로 제작되었다.

이후, 열전달 시나리오에 따른 콘크리트 페블 온도를 측정하였다. 제작 완료된 콘크리트 페블은 28일 수중 양생 후에 실험을 진행하여 환경적 요인에 따른 변수를 제거하였다.

Fig. 2 Graphite concrete pebbles

../../Resources/KCI/JKCI.2026.38.4.463/fig2.png

2.2 콘크리트 페블 열전달 실험 기법

Fig. 3(a)에 열저장 실험 장치를 나타내었다. 열저장 실험 장치는 열전달 유체 공급장치(블로우어), 유체 가열장치, 열전달 유체 배출관, 유체 온도 및 열사이클 조절장치로 구성하였다.

열전달 장치 하단에 위치한 유체 공급장치는 일정한 유량으로 열전달 유체(공기)를 공급하며, 공급된 유체는 유체 가열장치로 이동하여 온도 조절장치에서 설정한 목표 온도로 가열된다. 가열된 열전달 유체는 Fig. 3(b)의 페블을 담고 있는 강재 용기로 유입되어 용기를 지나가면서 콘크리트를 충열하거나 방열한 후, 유체 배출관을 통해 외기로 배출된다.

Fig. 3 Heat transfer experimental setup

../../Resources/KCI/JKCI.2026.38.4.463/fig3.png

콘크리트 페블 온도 분포를 측정하기 위하여, Fig. 4(a), (b)와 같이 강재 용기의 길이 방향(800 mm)을 따라 5개의 측정 단면을 설정하였다. 각 단면에는 Fig. 4(c)와 같이 하단과 상단에 온도 센서가 삽입된 콘크리트 페블을 배치하여 페블 자체의 온도 변화를 파악하였다.

각 측점 지점은 설치 단면과 위치에 따라 N(단면 번호)-P1(하단) 및 N-P2(상단)로 표기하였다.

Fig. 4 Installation of thermo-couples in the thermal storage vessel

../../Resources/KCI/JKCI.2026.38.4.463/fig4.png

2.3 열사이클 시나리오 설정

서론에서 언급한 바와 같이, 최근 산업용 히트펌프 및 열에너지 저장(TES) 시스템은 기존의 저온 저장 방식을 넘어 고온의 열원을 확보함으로써 시스템의 전체적인 에너지 밀도와 열역학적 효율을 극대화하는 방향으로 발전하고 있다.

일반적인 저온 저장 시스템과 달리, 400 °C 수준의 고온 저장은 히트펌프의 작동 온도 범위를 확장시켜 산업 공정열로 재활용할 수 있는 가능성을 보여준다.

이에 따라, 본 연구에서는 고온 히트펌프 시스템과의 연계 및 열저장 용량 확대를 고려하여 400 °C를 목표 온도로 설정하였다. 이를 통해, 그라파이트를 혼입한 콘크리트 페블이 400 °C 고온 환경에서도 안정적인 열적 성능과 구조적 안정성을 유지할 수 있는지 파악하고자 하였다.

열전달 유체의 유량(즉, 유체 질량 속도)과 열사이클(thermal cycling)의 온도상승 단계에서 온도상승 속도를 주요 실험변수로 고려하여 총 5가지 경우의 열전달 실험계획을 수립하였다. 열전달 유체는 뜨거운 공기(hot air)를 적용하였으며, 열전달 유체의 최대 가열온도는 400 °C이다.

실험 결과 간의 신뢰성을 확보하기 위하여 본실험 전 총 2회의 예비 열사이클을 실시하였다. 400 °C까지의 예비 가열 과정을 통해 콘크리트 페블의 열적 안정성을 확보하였으며, 내부 수분 증발 등 실험 결과에 부정적인 영향을 줄 수 있는 환경적 변수들을 제거한 후 본 실험을 진행하였다.

Fig. 5에 열전달 유체의 열사이클 시나리오를 나타내었다. 열전달 유체의 열사이클 시나리오는 열유체의 온도상승 단계, 온도가 일정하게 유지되는 단계 및 온도하강 단계 등 3단계로 구성된다.

Fig. 5 Thermal cycling scenario

../../Resources/KCI/JKCI.2026.38.4.463/fig5.png

실험 경우 1∼3의 열사이클 시나리오는 동일하지만 열전달 유체의 질량속도는 서로 상이하다. 열전달 유체의 열사이클은 2시간의 상승단계, 2시간의 일정온도 유지단계 및 5시간의 온도 구간하강 단계로 구성된다. 이때, 온도상승단계 및 온도유지단계는 열전달 장비의 온도제어장치를 통하여 자동으로 수행하며, 온도하강은 자연냉각으로 수행하였다. 또한, 실험 경우 1∼3의 열전달 유체의 유량을 서로 다른 3가지로 계획하였다. 각 실험경우의 열전달 유체의 유량은 각각 135, 270 및 405 LPM이며, 이는 0.0026, 0.0053 및 0.0079 kg/s의 열전달 유체의 질량 속도에 상응한다.

한편, 실험 경우 3∼5의 경우, 열전달 유체의 질량속도는 405 LPM, 즉 유체의 질량속도는 0.0079 kg/s로써 동일하지만, 열사이클 시나리오의 온도상승 단계의 온도상승 속도는 서로 상이하다. 각 실험 경우, 열전달 유체의 최대 목표온도 400 °C에 도달할 때까지 시간은 각각 120, 90 및 60분이며, 이에 따라 온도상승단계의 열전달 유체의 승온 속도는 각각 3.33, 4.44 및 6.67 °C/min로 구분된다.

3. 열전달 실험 결과 및 분석

3.1 열전달 유체 유량에 따른 온도 분포

열전달 시나리오의 목표온도 도달시간이 120분일 때 서로 다른 유량을 갖는 세 가지 실험 경우(Case 1∼3)의 각 단면 내 Fig. 6의 P2 측점 위치에서의 콘크리트 페블 온도변화를 비교하여 나타내었다.

열저장 용기의 열유체 입구 단면(단면 1)에서 열유체 유량에 따른 콘크리트 페블의 온도 변화 영향은 미미하다. 반면에, 열유체 출구 단면(단면 5)에서 열유체 유량에 따른 콘크리트 페블의 온도 변화 영향은 뚜렷하다.

구체적으로, 열전달 유체의 유량이 작은 경우(실험 경우 1)에 온도상승 종료시점(t=120 min)에서 입구 단면과 출구 단면 사이 즉, 단면 1과 단면 5에서 콘크리트 온도는 각각 344 °C 및 235 °C이고 두 단면의 콘크리트 페블 온도 차이는 109 °C로 비교적 크게 나타난다. 반면에, 열전달 유체의 유량이 큰 경우(실험 경우 3)에 동일한 시점에서 단면 1과 단면 5에서 콘크리트 온도는 각각 365 °C 및 318 °C이고 콘크리트 페블의 온도 차이는 47 °C로 비교적 작게 나타난다. 이는 열전달 유체의 유량이 증가할 때, 열전달 유체와 콘크리트 페블 사이의 열교환이 더욱 증가하여 열저장 용기 내 콘크리트 페블의 온도가 상승하기 때문이다.

결과적으로, 열전달 유체의 유량이 증가할 때 열전달 유체의 열에너지가 용기 입구에서 출구 방향 콘크리트로 더욱 빠르게 전달되어 단면 간 온도 차이가 감소하는 것으로 나타났다. 또한, 열전달 유체의 유량 증가는 열저장 용기 내부의 열확산을 촉진하고 온도 분포의 균일성을 향상시킨다.

Fig. 6 Temperature of concrete pebbles at various flow rates of the heat transfer fluid (Cases 1∼3)

../../Resources/KCI/JKCI.2026.38.4.463/fig6.png

3.2 온도상승 속도에 따른 온도 분포

서로 다른 열사이클 온도상승 단계 시간이 120, 90 및 60분일 때 즉, 이에 상응하는 서로 다른 3개 온도상승 속도 조건이 3.33, 4.44 및 6.67 °C/min일 때, 콘크리트 온도변화를 비교하여 Fig. 7에 나타내었다. 세 가지 실험 경우(Case 3∼5)의 P2 측점 위치에서의 콘크리트 페블의 온도 측정값을 제시하여 비교하였다.

열사이클 온도상승단계 종료 직후(120, 90 및 60분) 단면 1의 콘크리트 페블 온도는 각각 365, 360 및 346 °C로 차이가 나타난다. 반면에, 이후 온도유지단계 120분 경과 후 즉, 240, 210 및 180분일 때 콘크리트 페블의 온도는 385 °C로 서로 동일하다.

마찬가지로, 단면 5의 경우, 3개 실험 경우의 열사이클 온도상승단계 종료 시각인 120, 90 및 60분일 때, 콘크리트 페블의 온도는 각각 318, 303 및 269 °C로 서로 차이가 나지만, 온도유지 단계 120분 경과 후 즉, 240, 210 및 180분일 때 콘크리트 페블의 온도는 모두 382 °C로 서로 동일하다.

따라서, 열사이클 온도상승단계 시간이 서로 다를 때 온도상승 단계 종료 후 페블의 온도는 서로 차이를 나타내며, 추가로 120분의 일정온도유지단계가 지난 후에 세 가지 실험 경우의 콘크리트 페블은 거의 동일한 온도에 도달한다. 이러한 본 연구의 실험 결과는 온도상승 단계 시간의 차이가 열저장 시스템의 초기 충열 단계 열저장 성능에 영향을 미치며, 초기 단계 페블의 열저장량 차이를 유발할 수 있음을 의미한다.

Fig. 7 Temperature of concrete pebbles during various heating processes (Cases 3∼5)

../../Resources/KCI/JKCI.2026.38.4.463/fig7.png

3.3 온도경사 특성

3.3.1 열전달 유체 유량의 영향

온도상승 속도가 3.33 °C/min(실험 경우 1∼3)일 때, 열전달 유체의 유량에 따른 온도경사 변화를 Fig. 8에 나타내었다.

열전달 유체의 유량이 가장 낮은 경우(실험 경우 1), 온도상승이 종료되는 시점(t=120 min)에서 단면 1과 단면 5에서 콘크리트 페블의 온도는 각각 345 °C 및 235 °C로 110 °C의 온도 차가 발생한다. 열전달 유체의 유량이 가장 높은 경우(실험 경우 3), 두 단면에서 온도는 각각 365 °C 및 318 °C로 47 °C의 온도 차가 발생한다.

따라서, 실험 결과는 열전달 유체의 유량이 증가할 때 동일한 시간 동안 고온의 열전달 유체가 강재 용기 출구까지 빠르게 이동하여 페블 온도가 빠르게 상승하며, 이에 따라, 단면 간 온도경사가 작아지는 결과를 나타낸다. 이러한 현상은 유량이 높을수록 동일 시간 동안 열전달 유체가 출구까지 더욱 효과적으로 열을 전달하고 페블에 전달되는 총 열에너지의 양이 증가하기 때문이다. 결과적으로, 열전달 유체의 유량 증가는 강재 용기 내부의 온도 분포의 균일성을 향상시킨다.

한편, 온도상승 단계에서 온도경사는 뚜렷이 나타나지만, 열사이클 시나리오의 일정한 온도가 유지되면서 강재 용기 내부의 콘크리트 페블로 열전달이 점차 증진되면서 단면 간 온도 차가 감소하는 경향을 보인다. 실험 경우 3의 경우, 일정온도 유지단계 종료시간에 페블 온도는 최대 목표온도 400 °C에 근접하게 도달하며, 단면 1과 단면 5의 페블의 온도차이는 각각 2 °C이다.

따라서, 일정온도 유지단계 종료시간에 용기 내부는 두 단면 간의 온도가 일정하게 유지되는 열평형 상태에 도달하는 것을 의미하고, 열전달 유체의 유량 증가는 열평형 상태에 도달하는 시간을 단축시키고 유체 흐름방향 온도 분포를 균일하게 형성하는 데 영향을 주었다고 판단된다.

Fig. 8 Temperature gradient at various HTF flow rates (Cases 1∼3)

../../Resources/KCI/JKCI.2026.38.4.463/fig8.png

3.3.2 열전달 유체의 온도상승 속도의 영향

동일한 유량 조건(유속=0.0079 kg/s)에서 온도상승 속도를 3.33 °C/min에서 6.67 °C/min까지 변화시켰을 때의 온도경사를 비교하여 Fig. 9에 나타내었다.

온도상승 속도가 6.67 °C/min일 때(실험경우 5), Fig. 9(a)와 같이 열사이클의 온도상승 종료단계(t=60분)에서 단면 1과 단면 5의 콘크리트 페블 온도는 각각 346 °C 및 269 °C로써 77 °C의 온도차가 발생한다. 반면에, 온도상승 속도가 3.33 °C/min일 때(실험경우 3), 열사이클의 온도상승 종료단계(t=120분), 두 단면에서 콘크리트 페블의 온도는 각각 365 °C 및 318 °C로 온도차는 47 °C이다.

열전달 유체 온도상승 속도가 최대인 실험 경우 5의 경우, 열사이클 시나리오의 온도상승단계 종료 시 열교환 시간이 충분하지 않아서 온도분포가 제일 작게 나타난다. 또한, 열전달 유체 유입구와 가까운 콘크리트 페블의 온도상승에 비해 출구와 가까운 콘크리트 페블의 온도상승은 작아서 온도경사는 크게 발생한다. 따라서, 본 실험 결과는 열저장 용기 입구의 콘크리트 페블에서 출구 페블까지 충열이 이루어지기 위해 일정한 시간이 필요함을 의미한다.

한편, 온도상승 속도에 따라 온도상승 종료 시점에서 큰 온도경사의 차이가 발생하였을지라도, Fig. 9(b)Fig. 9(c)에 나타낸 바와 같이, 시간(t)이 90분 및 120분일 때, 실험 경우 5의 온도가 가장 크게 분포한다. 또한, 전반적으로 세 가지 실험경우의 온도경사는 점진적으로 감소한다.

Fig. 9(d)에 나타낸 열사이클의 일정온도 유지단계가 종료되는 시각에서 콘크리트는 전체적으로 최대목표온도에 도달하여 온도 분포는 일정하게 되고 결과적으로 온도경사는 사라진다. 이는 열사이클 온도상승단계에서 열확산에 필요한 시간이 충분하지 않아 온도경사가 크게 발생하였지만, 온도유지단계를 통해 열전달 유체와 콘크리트 페블 간의 열교환이 충분히 발생하고, 이에 따라 강재 용기 내부의 콘크리트 페블이 최대온도 400 °C에 근접하는 열적 평형 상태에 도달하기 때문이다.

따라서, 열저장 장치의 충열 단계에서 열전달 유체의 온도상승 속도의 증가는 콘크리트 열저장 매체의 온도경사 증가와 최대 충열 목표온도 도달시간을 감소시켜 열저장 효율성을 증가시킬 것으로 판단된다.

한편, 본 연구에서는 그라파이트를 혼입한 콘크리트 페블의 성능을 평가하기 위해 기존 열저장 매체인 자갈(Rock-bed)과 비교분석하였다. Zavattoni et al.(2014)에 따르면, 자갈을 이용한 팩-베드 시스템은 입자의 불규칙한 형상과 초기 열충격에 의한 입자 파쇄로 인해 공극률의 불균일성이 발생하며, 이는 채널링 현상을 유발하여 시스템의 성능을 저하시킨다. 반면, 본 연구에서 제안하는 그라파이트 혼입 콘크리트 페블은 일정한 구형 형상을 유지함으로써 유체 흐름에 대한 균일성을 확보할 수 있다. 또한, 콘크리트 특성상 고온 노출 후 일정 수준의 강도 감소가 발생함에도 불구하고, 400 °C 고온 노출 후 콘크리트 페블의 압축 하중은 약 3 kN을 유지하는 것으로 나타났다. 이는 상부 페블의 자중에 의해 하단부 페블이 파손되거나 분진이 발생할 가능성이 낮음을 의미한다.

Fig. 9 Temperature gradient at various heating rates (Cases 3∼5, fluid flow velocity = 0.0079 kg/s)

../../Resources/KCI/JKCI.2026.38.4.463/fig9.png

4. 열저장 시스템의 열전달률 및 열저장량

4.1 콘크리트 페블의 열에너지 산정

콘크리트 페블의 열에너지를 산정하기 위해, Fig. 10과 같이 강재 용기의 체적을 총 5개의 구간으로 구분하였다. 5개 체적 중 V1, V3, V5 구간 길이는 40 mm이며, V2와 V4의 구간 길이는 340 mm이다. 열저장 용기에 들어가는 페블 전체의 중량을 5개 구간의 부피 비에 따라 각 구간에 차지하는 콘크리트 페블의 질량을 산정하였다.

Fig. 10 Identification of volumes for thermal energy estimations

../../Resources/KCI/JKCI.2026.38.4.463/fig10.png

열저장 용기 내 구간별 콘크리트 페블의 열전달률은 식 (1)을 이용하여 산정하였다.

(1)
$\dot{Q}_{c,V_{n}}=c_{p,c}m_{c}(T_{c,V_{n}}(j)-T_{c,V_{n}}(i))/(j-i)$

여기서, $\dot{Q}_{c,V_{n}}$: 콘크리트 각 구간에서의 열전달률

$c_{p,c}$: 콘크리트의 비열($\text{J/kg}\cdot\text{K}$)

$m_{c}$: 콘크리트의 질량($\text{kg}$)

$T_{c,V_{n}}(i)$: $i$단계 일 때 $V_{n}$영역에서 콘크리트 온도($\text{K}$)

$T_{c,V_{n}}(j)$: $j$단계 일 때 $V_{n}$ 영역에서 콘크리트 온도($\text{K}$)

$V_{n}$: $n$번째 영역의 부피($n=1, 2, 3, 4, 5$)

열저장 용기 내 각 구간 별 콘크리트 페블의 열에너지는 다음 식 (2)를 이용하여 계산하였다.

(2)
$Q_{cn}=\int_{0}^{t}\dot{Q}_{c,V_{n}}dt\quad(n=1, 2, 3, 4, 5)$

따라서, 콘크리트 페블의 전체 열저장량($Q_{c}$)은 다음 식 (3)과 같이 각 구간의 열저장량을 합산하여 계산하였다.

(3)
$Q_{c}=Q_{c1}+Q_{c2}+Q_{c3}+Q_{c4}+Q_{c5}$

4.2 열전달 유체 유량에 따른 열전달률 및 열저장 특성

유속이 각각 0.0026, 0.0053 및 0.0079 kg/s일 때(실험 경우 1∼3) 열전달 유체의 질량 유속에 따른 콘크리트 페블의 열전달률과 열저장량을 Fig. 11에 나타내었다. Fig. 11(a)에서 양수의 열전달률은 콘크리트 페블의 열에너지 증가 상태를 의미하고, 음수의 열전달률은 콘크리트 페블의 열에너지 감소 상태를 의미한다.

온도상승 단계(t=0∼120분)에서 유속이 증가할 때, 콘크리트 페블의 열전달률은 양수로 급격히 증가한다. 이는 유속이 증가함에 따라 단위 시간당 강재 용기 내부로 유입되는 열전달 유체의 질량이 증가하므로, 이에 상응하여 콘크리트로 전달되는 열전달률이 증가하기 때문이다. 이에 따라, t=120분에서 열유체 유속이 제일 클 때 콘크리트 페블의 열저장량이 제일 크다.

온도유지단계(t=120∼240분)에서 콘크리트 페블은 열평형에 가까워지므로 열전달률은 점진적으로 감소하며, 온도유지 종료시점에서 거의 0에 가까워진다. 열전달률이 0에 수렴할 때, 세 가지 실험경우의 콘크리트 페블 열저장량은 거의 동일한 3,000 kJ에 도달한다.

한편, 열사이클의 온도하강단계(t=240∼420분)에서 콘크리트 페블로부터 열유체로 열전달이 발생하므로 열전달률은 음수로 급격히 증가한다. 이에 따라, 이 단계에서 유량이 제일 클 때, 콘크리트 페블의 열에너지가 빠르게 방열되어 궁극적으로 열저장량이 제일 작다.

콘크리트 페블의 열전달 성능은 열전달 유체와 콘크리트 페블 간의 열교환 특성에 의해 영향을 받는다. 열전달 유체 유속의 증가는 열전달 유체와 콘크리트 간의 열교환 증가를 유발하고, 결과적으로 콘크리트 페블의 충열과 방열 속도를 증가시킨다.

따라서, 본 실험결과는 열저장 시스템의 충열과 방열 속도 관점에서 열전달 유체의 유속은 열저장 운용 효율성에 영향을 미치는 것을 나타낸다.

Fig. 11 Heat transfer rate and thermal energy at different flow rates in the TES system (Heating rate = 3.33 ℃/min, Cases 1∼3)

../../Resources/KCI/JKCI.2026.38.4.463/fig11.png

4.3 온도상승 속도에 따른 열전달률 및 열저장 특성

Fig. 12는 열전달 유체의 유속을 0.0079 kg/s로 일정하게 유지한 상태에서 온도상승속도에 따른 열전달률 및 열저장량 변화를 비교한 결과이다.

열전달 유체의 온도상승 속도가 증가할수록 온도상승 단계에서 짧은 시간 내에 단위시간당 많은 열에너지가 페블에 공급되므로, 콘크리트 페블의 열전달률이 더욱 빠르게 증가한다. 이에 따라, 최대 열저장량에 도달하는 시간이 짧아진다.

일정온도 유지단계에서 콘크리트 페블의 열전달률은 양수를 유지하지만, 열전달 값은 서서히 감소한다. 일정온도 유지단계 종료시간에서 각 실험의 경우 열전달률은 0에 수렴하고, 결과적으로 최대 열저장량에 도달한다. 비교 대상 모든 실험 경우에서 콘크리트 페블의 열저장량은 거의 비슷한 값에 도달한다. 이는 실험에 사용된 콘크리트의 질량과 비열이 동일한 조건에서 각 열사이클 시나리오가 최종적으로 동일한 목표 온도인 400 °C에 도달하도록 계획되었기 때문이다. 마지막으로, 열사이클의 온도하강 단계에서 콘크리트 페블은 음수의 열전달률을 나타내고, 이에 따라 온도하강 속도가 빠를수록 콘크리트로부터 방열은 더욱 빠르게 진행된다.

본 실험 결과는 콘크리트 페블을 적용한 열저장 시스템 운용 시, 열전달 유체의 온도상승 속도에 따라 시스템 운전 속도를 유연하게 결정할 수 있음을 시사한다. 즉, 높은 온도상승 및 하강 속도를 열저장 시스템에 적용하여 단주기 사이클의 열저장 시스템 운용을 도모할 수 있고, 낮은 온도상승 및 하강 속도를 적용하여 장주기 사이클의 열저장 시스템 운용을 도모할 수 있다.

Fig. 12 Heat transfer rate and thermal energy of concrete pebbles at different heating rates in the TES system (fluid velocity=0.0079 kg/s, Cases 3∼5)

../../Resources/KCI/JKCI.2026.38.4.463/fig12.png

5. 결 론

본 연구에서는 열전달 유체의 질량유속과 열사이클 온도상승 속도를 실험변수로 설정하여 그라파이트를 혼입한 콘크리트 페블의 열전달 특성과 열저장 성능을 분석하였다. 주요 연구결과는 다음과 같다.

1) 열전달 유체의 유속이 낮을 때 열저장 용기 입구와 출구 단면 간의 온도 차가 크게 발생하고 온도 경사가 크게 나타난다. 반면에, 열유체의 유속이 증가할수록 단면 간 온도 차가 감소하여 온도 분포의 균일성이 향상되고 열저장 시스템이 고르게 목표 온도에 도달하는 경향을 나타낸다.

2) 열저장 시스템 충열 단계의 열유체 온도상승속도가 클수록 콘크리트 페블의 온도는 더욱 빠르게 상승하며, 목표 온도에 도달하는 시간은 단축된다. 한편, 일정온도유지 단계가 지난 후에 콘크리트 페블은 거의 동일한 충열 온도에 도달한다.

3) 열저장 시스템의 충열 단계에서 열전달 유체의 온도상승 속도의 증가는 콘크리트 열저장 매체의 온도경사 증가와 최대 충열 목표온도 도달시간을 감소시켜 열저장 효율성을 증가시킬 것으로 판단된다.

4) 콘크리트 페블의 열전달 성능은 열전달 유체와 콘크리트 페블 간의 열교환 특성에 의해 영향을 받는다. 열전달 유체 유속의 증가는 열전달 유체와 콘크리트 간의 열교환 증가를 유발하고, 결과적으로 콘크리트 페블의 충열과 방열 속도를 증가시킨다.

5) 본 실험 결과는 콘크리트 페블을 적용한 열저장 시스템 운용 시, 열저장 시스템에 공급되는 열유체의 온도상승 및 하강 속도를 조절하여 효과적인 충・방열 주기를 결정하고 열저장 시스템 설계를 도모할 수 있음을 나타낸다.

감사의 글

본 연구는 한국기계연구원 기관 기본사업(Project ID : NK262E)의 지원을 받아 수행된 결과입니다.

References

1 
Abedin, A. H., and Mahmud, S. (2014) Effect of Particle Size and Shape on Heat Transfer Characteristics in Packed-Bed Thermal Energy Storage Systems. Energy Conversion and Management 79, 189-198. Google Search
2 
Abera, M. M., Ancha, V. R., Amare, B., Sundar, L. S., Mouli, K. V. V. C., and Sangaraju, S. (2024) Simulation and Optimization of Energy Efficiency and Total Enthalpy Analysis of Sand-Based Packed Bed Solar Thermal Energy Storage. Frontiers in Heat and Mass Transfer 22(4), 1043-1070. DOI
3 
Adamson, K., Walmsley, T. G., Carson, J. K., Chen, Q., Schlosser, F., Kong, L., and Cleland, D. J. (2022) High-Temperature and Transcritical Heat Pump Cycles and Advancements: A Review. Renewable and Sustainable Energy Reviews 167, 112798. DOI
4 
Alnaimat, F., and Rashid, Y. (2019) Thermal Energy Storage in Solar Power Plants: A Review of the Materials, Associated Limitations, and Proposed Solutions. Energies 12(21), 4164. DOI
5 
Berrhazi, S., Ouammi, A., and Benchrifa, R. (2020) Thermo-Physical Effect of Solid Filler on the Performance of a Packed-Bed Thermal Storage. Thermal Science and Engineering Progress 20, 100716. DOI
6 
Chen, C., Liu, X., Liu, Z., Shao, L., Chang, H., Liu, Q., Miao, C., and Feng, P. (2024) Developing Heat-Conductive Concrete with Graphite-Modified Recycled Aggregates. Composites Part B: Engineering 284, 111721. DOI
7 
Diaz-Heras, M., Belmonte, J. F., and Almendros-Ibanez, J. A. (2020) Effective Thermal Conductivities in Packed Beds: Review of Correlations and Its Influence on System Performance. Applied Thermal Engineering 171, 115048. DOI
8 
Erregueragui, Z., Boutammachte, N., Bouatem, A., Merroun, O., and Zemmouri, E. L. M. (2016) Packed-Bed Thermal Energy Storage Analysis: Quartzite and Palm-Oil Performance. Energy Procedia 99, 370-379. DOI
9 
Gautam, A., and Saini, R. P. (2020) A Review on Sensible Heat Based Packed Bed Solar Thermal Energy Storage System for Low Temperature Applications. Solar Energy 207, 937-956. DOI
10 
Georgousis, N., Diriken, J., Speetjens, M., and Rindt, C. (2025) Comprehensive Review on Packed-Bed Sensible Heat Storage Systems. Journal of Energy Storage 121, 116516. DOI
11 
Han, C. G., Hwang, Y. S., Yang, S. H., and Gowripalan, N. (2005) Performance of Spalling Resistance of High Performance Concrete with Polypropylene Fiber Contents and Lateral Confinement. Cement and Concrete Research 35(9), 1747-1753. DOI
12 
Jaiswal, K. K., Chowdhury, C. R., Yadav, D., Verma, R., Dutta, S., Jaiswal, K. S., Sangmesh, B., and Karuppasamy, K. S. K. (2022) Renewable and Sustainable Clean Energy Development and Impact on Social, Economic, and Environmental Health. Energy Nexus 7, 100118. DOI
13 
Jeong, S. T., Park, J. H., Tran, T. K., and Yang, I. H. (2025) Thermal Performance of Concrete Containing Graphite at High Temperatures for the Application in a TES. Energies 18(17), 4685. DOI
14 
Kannan, N., and Vakeesan, D. (2016) Solar Energy for Future World: A Review. Renewable and Sustainable Energy Reviews 62, 1092-1105. DOI
15 
KATS (2022) Aggregates for Concrete (KS F 2527) Seoul, Korea: Korea Agency for Technology and Standards (KATS), Korea Standard Association (KSA). (In Korean) Google Search
16 
Kim, S. H., Jeon, H. K., Hwang, I. D., Seo, C. H., and Kim, S. H. (2013) The Properties of Concrete with Lightweight Aggregate Impregnated by Phase Change Material. Journal of the Korea Concrete Institute 25(3), 331-338. (In Korean) DOI
17 
Kruger, M., Haunstetter, J., Knodler, P., and Zunft, S. (2018) Slag as Inventory Material for a Thermal Energy Storage (TES): Material Investigation and Thermo-Mechanical Consideration. Energy Procedia 155, 454-463. DOI
18 
Laing, D., and Zunft, S. (2015) Advances in Thermal Energy Storage Systems 65-86. Woodhead Publishing. Google Search
19 
Mawire, A., and Taole, S. H. (2011) A Comparison of Experimental Thermal Stratification Parameters for an Oil/Pebble-Bed Thermal Energy Storage (TES) System During Charging. Applied Energy 88(12), 4766-4778. DOI
20 
Mawire, A., McPherson, M., and van den Heetkamp, R. R. J. (2010) Discharging Simulations of a Thermal Energy Storage (TES) System for an Indirect Solar Cooker. Solar Energy Materials and Solar Cells 94(6), 1100-1106. DOI
21 
Rahjoo, M., Rojas, E., Goracci, G., Gaitero, J. J., Martauz, P., and Dolado, J. S. (2023) A Numerical Study of Geopolymer Concrete Thermal Energy Storage: Benchmarking TES Module Design and Optimizing Thermal Performance. Journal of Energy Storage 74, 109389. DOI
22 
Shah, R. K., and Lee, J. H. (2011) Thermal Performance Enhancement of Packed-Bed Thermal Energy Storage Systems Through Optimization of Particle Arrangement and Flow Configuration. Applied Thermal Engineering 31, 345-354. Google Search
23 
Park, J. H., and Yang, I. H. (2023) The Material Characteristics of Ultra-High Performance Concrete for Thermal Energy Storage Media at High Temperature. Journal of the Korea Concrete Institute 35(4), 421-431. (In Korean) DOI
24 
Singh, H., Saini, R. P., and Saini, J. S. (2013) Performance of a Packed Bed Solar Energy Storage System Having Large Sized Elements with Low Void Fraction. Solar Energy 87, 22-34. DOI
25 
Strasser, M. N., and Selvam, R. P. (2014) A Cost and Performance Comparison of Packed Bed and Structured Thermocline Thermal Energy Storage Systems. Solar Energy 108, 390-402. DOI
26 
Tagle-Salazar, P. D., Cabeza, L. F., and Prieto, C. (2026) Advancing Sensible Heat Storage: A Novel Transient Heat Transfer Model for Concrete-Based TES Modules for CSP Applications. Renewable Energy 256, 124558. DOI
27 
Tiwari, S., and Saini, R. P. (2020) Experimental Investigation of Packed-Bed Solar Thermal Storage. Proceedings of the Institution of Civil Engineers - Energy 173(4), 157-166. DOI
28 
Panchal, M., Saraswat, A., and Chaudhuri, P. (2020) Experimental Measurements of Gas Pressure Drop of Packed Pebble Beds. Fusion Engineering and Design 160, 111836. DOI
29 
Wullhorst, F., Reuter-Schniete, J., Maier, L., and Muller, D. (2024) Heat Pump and Thermal Energy Storage: Influences of Photovoltaic, the Control Strategy, and Price Assumptions on the Optimal Design. Renewable Energy 236, 121409. DOI
30 
Martin-Alcantara, A., and Fernandez-Feria, R. (2025) Optimum Design and Operation of an Internally-Heated Packed-Bed Tank for High-Temperature Solid-Air Energy Storage: Practical Application from Theoretically-Derived Relations. Journal of Energy Storage 132, 117844. DOI
31 
Yuan, H., Lu, C., Xu, Z., Ni, Y., and Lan, X. (2012) Mechanical and Thermal Properties of Cement Composite Graphite for Solar Thermal Storage Materials. Solar Energy 86(11), 3227-3233. DOI
32 
Zavattoni, S. A., Barbato, M. C., Pedretti, A., Zanganeh, G., and Steinfeld, A. (2014) High Temperature Rock-Bed TES System Suitable for Industrial-Scale CSP Plant: CFD Analysis Under Charge/Discharge Cyclic Conditions. Energy Procedia 46, 124-133. DOI
33 
Zhou, J., Wu, W., Bellamy, L., and Bishop, D. (2026) Thermal Energy Storage-Coupled Heat Pump Systems: Review of Configurations and Modelling Approaches. Renewable and Sustainable Energy Reviews 226, 116226. DOI
34 
Rao, R., Wang, H., Wang, H., Tuan, C. Y., and Ye, M. (2019) Models for Estimating the Thermal Properties of Electric Heating Concrete Containing Steel Fiber and Graphite. Composites Part B: Engineering 164, 116-120. DOI
35 
Warkhade, G. S., Babu, A. V., Mane, S., and Ganesh Babu, K. (2016) Experimental Investigation of Sensible Thermal Energy Storage in Small Sized, Different Shaped Concrete Material Packed Bed. World Journal of Engineering 13(5), 386-393. DOI
36 
Yang, I. H., and Park, J. H. (2019) Mechanical and Thermal Properties of UHPC Exposed to High-Temperature Thermal Cycling. Advances in Materials Science and Engineering 2019(1), 9723693. DOI
37 
Yang, I. H., Kim, K. C., and Choi, Y. C. (2016) Effect of Cementitious Composite on the Thermal and Mechanical Properties of Fiber-Reinforced Mortars for Thermal Energy Storage. Journal of the Korea Concrete Institute 28(4), 395-405. (In Korean) DOI
38 
Mawire, A., McPherson, M., Heetkamp, R. R. J. v. d., and Mlatho, S. J. P. (2009) Simulated Performance of Storage Materials for Pebble Bed Thermal Energy Storage (TES) Systems. Applied Energy 86(7), 1246-1252. DOI