김진만
(Jin Man Kim)
1
김선창
(Seon-Chang Kim)
2†
-
한국생산기술연구원 산업에너지연구부문 선임연구원
(Senior Researcher, Industrial Energy R&D Department, Korea Intitute of Industrial
Technology, Cheonan, 3056, Korea)
-
한국생산기술연구원 산업에너지연구부문 수석연구원
(Principal Researcher, Industrial Energy R&D Department, Korea Intitute of Industrial
Technology, Cheonan, 31056, Korea)
Copyright © 2016, Society of Air-Conditioning and Refrigeration Engineers of Korea
키워드
대류 열전달계수, R1336mzz(Z), 증기 생산, 초임계 유체, 초고온 히트펌프
Key words
Convective heat transfer coefficient, R1336mzz(Z), Steam generation, Supercritical fluid, Very high temperature heat pump
기호설명
$d$:
직경 [m]
$h$:
열전달 계수 [W/(m2·K)]
$i$:
엔탈피 [J/kg]
$k$:
열전도도 [W/(m·K)]
$l$:
단일 시험부 길이 [m]
$\dot{m}$:
질량 유량 [kg/s]
Nu:
Nusselt 수
P:
압력 [kPa]
Pr:
Prandtl 수
$\dot{Q}$:
열전달률 [W]
Re:
Reynolds 수
T:
온도 [℃]
하첨자
$cw$:
냉각수
$i$:
입구
$no$:
평준화
$o$:
출구
$r$:
냉매
$t$:
전열관
$wi$:
전열관 내벽
$wo$:
전열관 내벽
1. 서 론
현재 각 국의 산업공정에서는 다양한 목적으로 열원을 사용하고 있으며 산업공정에 활용되는 에너지는 19% 수준이지만 산업공정의 열원으로 인한 탄소 배출량은
전체의 36.8%에 달한다.(1) 따라서, 탄소중립 정책과 함께 전 세계 산업에서는 공정에 사용하고 있는 보일러와 같은 열원기기를 히트펌프로 대체하고자 하는 움직임이 활발히 일어나고
있다.
기존 선행연구(2-6)에 의하면 압축식 히트펌프는 열원 온도와 활용 온도에 따라 히트펌프, 고온 히트펌프 및 초고온 히트펌프로 나뉜다. 이 분류의 기준은 연구자마다, 산업군마다
다르게 정의되는데, 본 연구에서 140℃ 이상의 히트펌프는 모두 초고온 히트펌프로 보고자 한다.
현재 100℃ 수준의 히트펌프 기술은 이미 성숙 단계에 이르러 산업에서 활용이 가능하지만, 초고온 히트펌프 기술은 아직 꾸준한 기술 개발이 이루어지고
있는 실정이다.(7) 하지만, 산업에서 요구하는 열원 중 200℃ 이상 열원에 대한 수요는 62%(8)로 매우 큰 비율을 차지하며 초고온 히트펌프 기술 개발의 시급성을 보여주고 있다. 이와 같은 산업의 요구에 대응하기 위해 세계 선진 기업 중 Vicking
Heating Engines AS와 Kobel Steel은 150℃ 수준의 히트펌프 제품을, Ochsner, Mayekawa, Combitherm
및 Dürr thermea 등은 130℃ 수준의 히트펌프 제품(9)을 내놓고 있지만, 아직 초고온 히트펌프 기술은 성숙도가 많이 낮은 것으로 분석된다.
이렇게 초고온 히트펌프 기술 개발이 더딘 이유는 냉매 규제와 더불어 HFO(Hydrofluoroolefin) 계열의 냉매를 이용한 요소 부품 기술
개발의 어려움 등과 맞물려 있는 것으로 보인다. 특히, 초고온 히트펌프의 경우 냉매의 열분해 문제로 인해 온도에 따라 사용할 수 있는 냉매가 제한되며,
200℃ 이상에서는 HFO 계열의 냉매 중 R1336mzz(Z) 정도만 안정적으로 사용이 가능한 것으로 보고되었다.(10-12) 뿐만 아니라 초고온 히트펌프는 높은 활용 온도로 인해 방열부에서 초임계 상태로 열교환을 하게 되는데, HFO 냉매의 초임계 유체에 대한 열전달 연구(13-17)가 일부 수행된 바 있으나 아직 다른 냉매와 비교할 때 HFO 냉매의 초임계 열전달 연구는 상대적으로 부족한 실정이다.
따라서, 본 연구에서는 160℃ 증기 생산을 위한 초고온 히트펌프 개발의 일환으로 R1336mzz(Z)을 이용한 초임계 유체 가스 냉각기의 기초 연구를
수행하였으며, 실험을 통해 단관 내 초임계 유체의 열전달 특성을 분석하고 무차원 수를 이용한 상관식을 제시하였다.
2. 연구방법
본 연구에서는 160℃ 증기 생산을 위한 히트펌프 개발을 목표로 가스 냉각기의 기초 연구를 수행하였고, 하기와 같이 냉매의 특성을 고려해 시험부와
실험장치를 구성하였다.
2.1 실험 조건 및 실험 장치 구성
실험조건은 기존에 보고된 선행연구(18)의 사이클 작동 조건을 고려해 Table 1과 같이 설정하였다. 냉매는 R1336mzz(Z)를 사용하였으며 실험을 수행한 조건을 Table 1에 명시하였다.
Table 1 Experimental conditions for a supercritical gas cooler in a single tube
|
Inlet temperature
|
Outlet temperature
|
Flow rate
|
Pressure
|
R1336mzz(Z) (tube side)
|
200℃
|
134 ~ 145℃
|
2.5 ~ 3.5 kg/min
|
3,000 ~ 3,850 kPa
|
Water (shell side)
|
42℃
|
82 ~ 93℃
|
2 LPM
|
-
|
냉매로 사용된 R1336mzz(Z)의 경우 임계 압력은 2,903 kPa, 임계 온도는 171℃로 임계점 이상의 조건에서는 초임계 유체의 특성을 갖는다.
초임계 유체는 임계 온도 부근에서 유체의 물성이 급격하게 변하기 때문에 일반 유체와는 다른 특성을 나타내고, 이로 인해 열전달이나 유동 특성을 이해하는데
어려움이 있다. 따라서, 시험부(test section)는 Fig. 1과 같이 9개의 짧은 단일 시험부를 연결해 온도 변화에 따른 열전달 특성을 상세히 분석하고자 하였다. 단일 시험부는 이중 쉘로 구성되어 있으며 냉매는
1/2 인치 전열관 내부를, 냉각수는 1 인치 쉘 쪽으로 흐르며 열교환이 이루어지게 된다. 단일 시험부는 100 mm의 길이로 제작되었고, 각 단일
시험부의 입·출구에 RTD(Resistance Temperature Detector)를 삽입해 온도를 측정함으로써 냉매 온도 변화에 따른 열전달 특성을
단일 시험부별로 분석하고자 하였다. 각 시험부의 중앙에는 외경 1 mm의 TC(Thermocouple)가 쉘을 관통해 전열관의 외표면에 용접되어 있어
열전달 특성을 분석하는데 활용하였다.
실험 수행을 위한 전체 실험장치(Test loop)는 Fig. 2와 같이 구성하였다. 시험부 전단의 히터(Electric heater \#1)를 제어함으로써 냉매의 입구 온도를 제어하였으며, 시험부 후단의 히터(Electric
heater \#2)를 이용해 전체 시스템의 압력을 제어하였다. 다만, 시험부 후단의 히터로 인해 냉매가 지나치게 높은 온도까지 가열될 경우 마그네틱
기어펌프(Magnetic gear pump)에 손상을 일으킬 수 있기 때문에, 리시버 탱크(Receiver tank)를 지난 후 판형 열교환기(P.H.E.)를
이용해 170℃ 이하로 냉각 시켜주었다.
Fig. 1 A schematic image of the test section.
Fig. 2 A schematic image of the test loop.
2.2 데이터 분석 방법론
Fig. 1과 같이 구성된 각 시험부의 열전달 특성을 분석하기 위해 우선 식(1)과 같이 냉각수의 질량 유량($\dot{m_{cw}}$)과 입․출구 엔탈피($i_{cw}$)를 이용해 n번째 시험부의 열전달률($\dot{Q}_{cw,\:
n}$)을 계산하였다.
각 시험부에서 전열관 외측에 용접된 TC를 이용해 측정된 전열관 외벽 온도($T_{wo,\: n}$)와 전열관의 열저항을 이용하면 식(2)와 같이 전열관 내벽의 온도($T_{wi,\: n}$)를 계산할 수 있다.
그리고 최종적으로 식(3)과 같이 각 시험부의 평균 냉매 온도($T_{r,\: n}$)를 이용해 냉매의 열전달계수($h_{i,\: n}$)를 계산하였다.
분석과정에서 필요한 R1336mzz(Z)와 냉각수의 물성은 NIST(National Institute of Standards and Technology)의
REFPROP V10을 이용하였으며, 상관식들의 계산은 공학용 계산 프로그래밍 소프트웨어인 EES(Engineering Equation Solver)
V11.731을 이용하였다.
추가적으로, 실험 데이터 및 계산 결과에 대한 신뢰성을 평가하고자 불확도(Uncertainty) 분석을 수행하였다. 온도, 압력, 냉매의 질량 유량(Mass
flow rate of the refrigerant) 및 냉각수의 체적 유량(Volume flow rate of the cooling water)에
대해서는 B형 표준 불확도(Standard uncertainty)를 산출하였으며, 열전달률(Heat transfer rate)과 열전달계수(Heat
transfer coefficient)의 합성 표준 불확도(Combined standard uncertainty)는 불확도 전파(Uncertainty
propagation) 법칙(19)을 이용해 산출하였다. Table 2는 분석된 최대 불확도를 나타낸 표이다.
Table 2 Maximum uncertainties of major parameters
|
Parameter
|
The maximum uncertainty
|
Standard uncertainty
|
Temperature [℃]
|
0.06
|
Pressure [kPa]
|
1.7
|
Mass flow rate of the refrigerant [kg/min]
|
0.002
|
Volume flow rate of the cooling water [L/min]
|
0.003
|
Combined standard uncertainty
|
Heat transfer rate [W]
|
11
|
Heat transfer coefficient [W/(m2K)]
|
68
|
3. 실험 결과 및 고찰
3.1 실험장치의 열손실 분석
구축된 실험장치의 기초 검증을 위해 시험부의 열 평형(heat balance) 분석을 수행하였고, 그 결과를 Fig. 3에 나타내었다. 전반적으로 평균 약 7% 수준의 열 손실이 발생했으며, 이 결과에 따르면 냉매가 가진 에너지의 약 7% 수준은 냉각수가 아닌 외부로
빠져나가기 때문에 실질적으로 각 시험부에서 유체간의 열교환을 통해 이루어진 열전달률은 식(1)에서 구한 냉각수의 열전달률과 같다고 볼 수 있다.
Fig. 3 Heat balance evaluation of the test section.
3.2 초임계 조건에서의 열전달계수 변화
Table 1은 냉매 유량 2.5 kg/min ~ 3.5 kg/min인 범위에서 수행한 실험 결과 그래프로, Fig. 4(a), Fig. 4(b) 및 Fig. 4(c)는 각각 압력이 3,000 kPa 일 때, 3,400 kPa 일 때 및 3,850 kPa 일 때 냉매 온도에 따른 대류 열전달계수를 나타낸 그래프이다.
각 압력 조건에서 유량에 따른 열전달계수를 비교하면, 국부적으로는 조금씩 차이가 있으나 전체적인 경향은 냉매의 유량이 커질수록 열전달계수가 높아지는
방향으로 분석되었다. 이는 유량이 증가할수록 Reynolds 수가 커지게 되는데, 일반적으로 Reynolds 수가 커질수록 대류 열전달계수가 증가하기
때문에 이와 같은 경향이 나타난 것으로 이해된다.
또한, 3,000 kPa, 3,400 kPa 및 3,850 kPa에서 유사 임계 온도는 각각 173℃, 181℃ 및 188℃인데, Fig. 4를 보면 유사 임계 온도 이하에서는 냉매의 온도가 증가하면서 열전달계수가 서서히 높아진다. 그리고 냉매의 온도가 유사 임계 온도에 이르면 열전달계수의
급격한 상승을 보이다가 유사 임계점을 지난 이후로는 다시 열전달계수가 감소하는 경향을 보이고 있다. 이와 같은 열전달계수의 변화 경향은 유사 임계
온도에서 냉매의 급격한 물성 변화와 관련이 있는 것으로 판단된다.
Fig. 4 Change in heat transfer coefficient according to the temperature of the refrigerant at each flow rate condition.
Fig. 5는 압력 3,000 kPa인 조건에서 냉매의 평준화 온도(Tnorm) 증가에 따른 Pr의 변화를 나타낸 그래프이다. Tnorm는 냉매의 온도를 유사
임계 온도로 나눈 값으로 유사 임계 온도에서 1이 되는 무차원 변수이다. 그래프에서 나타난 것과 같이 Tnorm가 1일 때 Pr가 증가했다가 이후에는
급격히 감소하는 경향을 확인할 수 있다. 즉, 유사 임계 온도에 도달하게 되면 Pr의 증가로 대류 열전달 효과가 촉진되어 열전달계수가 증가했다가 Tnorm가
1보다 커지게 되면 Pr의 급격한 감소로 인하여 대류 열전달 효과가 상대적으로 저하되며 열전달계수가 다시 감소하게 되는 것으로 이해된다.
동일 유량 조건에서 압력에 따른 열전달계수를 비교하면 Fig. 6과 같으며, Fig. 6(a), Fig. 6(b) 및 Fig. 6(c)는 각각 냉매 유량이 2.5 kg/min, 3.0 kg/min 및 3.5 kg/min 일 때 냉매 온도에 따른 열전달계수의 변화를 나타낸 그래프이다.
각 그래프를 보면 압력이 낮을수록 열전달계수가 더 높게 나타나고 있으며, 압력에 따른 이러한 경향은 냉매 유량과 상관없이 나타나는 것을 알 수 있다.
그리고 압력이 낮을수록 열전달계수가 높게 나타나는 이유는 압력에 따른 물성 변화로 인해 낮은 압력에서 Re가 더 높아지기 때문으로 분석된다.
Fig. 5 Change in Prandtl number according to the normalized temperature of the refrigerant(@ P = 3,000 kPa).
Fig. 6 Change in heat transfer coefficient according to the temperature of the refrigerant at each pressure condition.
3.3 무차원 상관식 도출
앞서 분석된 데이터를 바탕으로 Re, Pr 그리고 Nu를 정의하였고, 이들 간의 관계를 정의함으로써 무차원 열전달 상관식을 도출하고자 하였다. Fig. 7은 무차원 수들 간의 관계를 나타낸 그래프로 도출된 상관식을 식(4)에 명시하였다. 제안한 상관식의 $\overline{R}^{2}$는 약 0.85 이며, 전체 135개의 데이터 중 93%인 126개의 데이터가 오차범위
25% 이내에 들어오는 것을 확인하였다.
그리고 이렇게 도출된 상관식과 기존 문헌에 제시된 상관식의 상대적인 비교를 위해 본 연구의 실험 데이터와 함께 상관식들을 비교하였다. 압력이 3,850
kPa인 경우 열전달계수와 기존 문헌에 제시된 다른 상관식들과 열전달계수를 비교하면 Fig. 8과 같이 나타난다. 원형 전열관 내의 열전달계수를 예측하는데 가장 널리 활용되는 Dittus and Boelter(20) 식은 가장 낮은 열전달계수 예측값을 나타냈고, Gnielinski(21) 식과 Petukhov(22) 식을 이용했을 경우 Dittus and Boelter(20) 식보다는 큰 값을 나타내지만 여전히 실험값을 과소예측 하는 것으로 나타났다. 반면, 본 실험을 통해 개발한 상관식의 경우 상대적으로 더 근사한 예측값을
보여주는 것을 확인했다.
Fig. 7 Correlation between dimensionless numbers.
Fig. 8 Comparison of heat transfer coefficients between experimental data and correlations (@ P = 3,850 kPa).
4. 결 론
본 연구에서는 160℃ 증기 생산을 위한 히트펌프 개발의 일환으로 초임계 유체 가스 냉각기의 기초 데이터 확보를 위해 R1336mzz(Z)를 이용한
열전달 연구를 수행하였고, 다음과 같은 결론을 내릴 수 있었다.
(1) R1336mzz(Z)의 온도가 증가함에 따라 열전달계수가 높아지다가 유사 임계 온도 부근에서는 급격한 상승을 보였고, 유사 임계 온도를 넘어가면서
열전달계수가 점차 감소하는 경향을 보이고 있다.
(2) 전반적으로는 냉매의 질량 유량이 증가할수록, 압력이 낮을수록 열전달계수가 높게 나타나는 경향을 보인다.
(3) Nusselt 수를 Reynolds 수와 Prandtl 수의 함수로 나타내는 상관식을 도출하였고, 상관식 도출에 사용된 데이터 중 약 93%의
데이터가 25% 오차범위 이내로 들어오는 것을 확인하였다.
후 기
본 연구는 2024년도 산업통상자원부 기계․장비산업기술개발사업 “350 kW급 160℃ 증기(steam) 생산용 무급유 원심식 히트펌프 개발”의
연구비 지원에 의한 연구임(No. 20018456).
References
Hamid, K., Sajjad, U., Ahrens, M. U., Ren, S., Ganesan, P., Tolstorebrov, I., Arshad,
A., Said, Z., Hafner, A., and Wang, C.-C., 2023, Potential Evaluation of Integrated
High Temperature Heat Pumps: A Review of Recent Advances, Applied Thermal Engineering,
p. 120720.
![DOI](../../Resources/images/icon_crossref.png)
Jakobs, R. and Laue, H. J., 2015, Application of Industrial Heat Pumps IEA HPP IETS
Annex 35-13 IEA HPC: Workshop Regarding Heat Pumps and IEA Projects.
![URL](../../Resources/images/icon_url.png)
Laue, H. J., Application of Industrial Heat Pumps: Final Report-Part 2.
![URL](../../Resources/images/icon_url.png)
Bobelin, D. and Bourig, A., 2012, Experimental Results of a Newly Developed Very High
Temperature Industrial Heat Pump (140 C) Equipped with Scroll Compressors and Working
with a New Blend Refrigerant.
![URL](../../Resources/images/icon_url.png)
Peureux, J. L., Sicard, F., and Bobelin, D., 2014, French Industrial Heat Pump Developments
Applied to Heat Recovery.
![URL](../../Resources/images/icon_url.png)
Peureux, J. L., Sapora, E., and Bobelin, D., 2012, Very High-Temperature Heat Pumps
Applied to Energy Efficiency in Industry, Frankfurt: Achema, Application of Industrial
Heat Pumps, pp. 1-23.
![URL](../../Resources/images/icon_url.png)
Arpagaus, C., Bless, F., Uhlmann, M., Schiffmann, J., and Bertsch, S. S., 2018, High
Temperature Heat Pumps: Market Overview, State of the Art, Research Status, Refrigerants,
and Application Potentials, Energy, Vol. 152, pp. 985-1010.
![DOI](../../Resources/images/icon_crossref.png)
Rehfeldt, M., Fleiter, T., and Toro, F., 2018, A Bottom-up Estimation of the Heating
and Cooling Demand in European Industry, Energy Efficiency, Vol. 11, pp. 1057-1082.
![DOI](../../Resources/images/icon_crossref.png)
Arpagaus, C., Bless, F., Schiffmann, J., and Bertsch, S. S., 2017, Review on High
Temperature Heat Pumps, Proceedings of the European Heat Pump Summit.
![URL](../../Resources/images/icon_url.png)
Kontomaris, K., 2014, HFO-1336mzz-Z: High Temperature Chemical Stability and Use as
a Working Fluid in Organic Rankine Cycles.
![URL](../../Resources/images/icon_url.png)
Juhasz, J. R. and Simoni, L. D., 2015, A Review of Potential Working Fluids for Low
Temperature Organic Rankine Cycles in Waste Heat Recovery.
![URL](../../Resources/images/icon_url.png)
Mateu-Royo, C., Navarro-Esbrí, J., Mota-Babiloni, A., Molés, F., and Amat-Albuixech,
M., 2019, Experimental Exergy and Energy Analysis of a Novel High-Temperature Heat
Pump with Scroll Compressor for Waste Heat Recovery, Applied Energy, Vol. 253, p.
113504.
![DOI](../../Resources/images/icon_crossref.png)
He, J., Dang, C., and Hihara, E., 2018, Supercritical Heat Transfer Characteristics
of R1233zd (E) in Vertically Upward Flow, International Journal of Heat and Mass Transfer,
Vol. 127, pp. 497-505.
![DOI](../../Resources/images/icon_crossref.png)
Lee, C., Ko, J.-W., Ji, H. Y., and Kim, S.-C., 2020, Experimental Assessment of Convective
Heat Transfer and Pressure Drop Correlation of R1234ze (E) for a Supercritical Heat
Exchanger in the Organic Rankine Cycle, Journal of Mechanical Science and Technology,
Vol. 34, pp. 4809-4818.
![DOI](../../Resources/images/icon_crossref.png)
Ibrahim, A., Hu, P., Jiang, Y., Riaz, A., Saleem, F., Dong, Y., Jia, L., and Zhao,
P., 2021, Numerical Investigations on Supercritical Heat Transfer Characteristics
of Environmental Friendly Refrigerants, International Journal of Clean Coal and Energy,
Vol. 10, No. 2, pp. 21-40.
![DOI](../../Resources/images/icon_crossref.png)
Elbarghthi, A. F. A., Hdaib, M. Y., and Dvořák, V., 2022, Heat Transfer Analysis between
R744 and HFOs inside Plate Heat Exchangers, Entropy, Vol. 24, No. 8, p. 1150.
![DOI](../../Resources/images/icon_crossref.png)
Jiang, Y.-R., Hu, P., Chen, Q., Jia, C.-Q., Zhao, P.-P., and Jia, L., 2022, Experimental
Investigation On Heat Transfer Enhancement And New Correlation Of Supercritical R1234ze
(E) In Horizontal Helically Coiled Tube.
![URL](../../Resources/images/icon_url.png)
Kim, J. M. and Kim, S.-C., 2023, Optimization of High-Temperature Heat Pump Cycle
for Producing Steam of 160℃, Korean Journal of Air-Conditioning and Refrigeration
Engineering, Vol. 35, No. 8, pp. 371-380.
![DOI](../../Resources/images/icon_crossref.png)
Taylor, B. N. and Kuyatt, C. E., 1994, Guidelines for Evaluating and Expressing the
Uncertainty of NIST Measurement Results.
![URL](../../Resources/images/icon_url.png)
Dittus, F. W. and Boelter, L. M. K., 1985, Heat Transfer in Automobile Radiators of
the Tubular Type, International Communications in Heat and Mass Transfer, Vol. 2,
No. 1, pp. 3-22.
![URL](../../Resources/images/icon_url.png)
Gnielinski, V., 1976, New Equations for Heat and Mass Transfer in Turbulent Pipe and
Channel Flow, International Chemical Engineering, Vol. 16, No. 2, pp. 359-368.
![URL](../../Resources/images/icon_url.png)
Petukhov, B. S., 1970, Heat Transfer and Friction in Turbulent Pipe Flow with Variable
Physical Properties, Advances in Heat Transfer, Vol. 6, No. 503, p. i565.
![DOI](../../Resources/images/icon_crossref.png)