PCHE 내 열교환 유체의 경우 극저온 유체로 액체질소를 사용하였으며 열원 유체로는 글리콜 워터를 사용하였다. 글리콜 워터의 경우 물과 에틸렌글리콜을 질량 분율 50%로 섞어 사용하였다. 해당 유체들의 물성치는 Table 1과 같다.

대류로 인한 열손실을 최소화하기 위하여 진공 체임버 내에서 실험을 수행하였으며 전도로 인한 열손실을 최소화하기 위하여 PEEK block을 제작하여 점접촉 상태에서 PCHE가 지지되도록 하였다. 시험 설비의 개략도는 Fig. 1과 같다.

Fig. 1. 
Schematic diagram of experimental setup

PCHE의 극저온 채널의 입/출구와 열원 채널의 입/출구에서의 압력 및 온도의 경우 진공 체임버 내에서 측정하였다. 액체질소의 경우 개방 루프 형태로 구성하였다. 액체질소는 액체질소 봄베로부터 공급되어 극저온 레귤레이터(B-type; Cash Valve, McKinney, TX, USA)를 이용하여 압력 및 유량이 조절된 뒤 진공 체임버 내 PCHE에서 열교환한 후에 외부로 토출되었다. 토출 전 질량 유량계(CMF010 Micromotion; Emerson Electric, St. Louis, MO, USA)를 사용하여 액체질소의 유량을 측정하였으며 질량 유량계의 측정 정확도를 유지하기 위하여 액체질소가 PCHE에서 열교환한 뒤 히터를 거쳐 상온의 기체질소가 되어 질량 유량계를 지나도록 하였다. 열원 채널의 경우 폐루프 형태로 구성하였다. 순환 수조(Dyneo DD-1001F; Julabo, Seelbach, Germany)를 이용하여 일정한 온도로 유지되도록 하였으며 마그네틱 기어 펌프(DGS 2.3; ThrillCorp, Orlando, FL, USA)를 이용하여 일정한 유량으로 순환될 수 있도록 하였다. 사용한 요소 기기들과 계측기들의 사양은 Tables 2-4와 같다. 추가적으로 진공 체임버 내 이중관 열교환기를 구성하여 액체질소가 가능한 낮은 quality를 유지하여 PCHE 내부로 들어갈 수 있도록 구성하였다.

PCHE의 경우 기보유한 PCHE⁹⁾를 활용하였으며 해당 PCHE의 상세 규격은 Table 5와 같다. 채널 직경의 경우 극저온 채널 및 열원 채널 모두 1 mm의 반원으로 설계되었으며 각각 10장의 전열판 수로 제작되었다. 전열판의 경우 극저온 채널의 경우 I type 형태로 제작하였으며 열원 채널의 경우 N type 형태로 제작하였다. 확산 접합을 이용하여 PCHE의 코어 부분을 접합하였으며 추후 용접을 통하여 헤더를 코어 부분에 결합하였다. 제작된 PCHE는 Fig. 2와 같다.

Fig. 2. 
(a) Fabricated PCHE and (b) configuration of PCHE

본 시험에서는 PCHE 내 극저온 채널 및 열원 채널에서의 열전달량을 측정하였으며 각각 채널에서의 차압 역시 측정하였다. 모든 시험의 경우 글리콜 워터가 40℃인 조건에서 수행되었으며 열원 채널의 유량뿐만 아니라 극저온 채널의 유량 역시 변경하면서 시험을 수행하였다. 열원 채널의 유량을 줄여가며 PCHE 내 동결 조건이 발생할 경우 열전달량이 감소하는 것을 확인한 후에는 PCHE를 해동시킨 후 시험을 진행하였다. 동결을 회피하기에 유리한 형태를 취하기 위하여 counter flow가 아닌 co-flow로¹⁰⁾ 시험을 수행하였으며 극저온 유체는 I type에 흐르도록 하였고 열원 유체는 N type에 흐르도록 하였다.

PCHE 입/출구 온도의 경우 T-type 열전대를 이용하여 측정하였고 압력의 경우 동일한 위치에서 압력계(3051TG/2051CA; Rosemount Inc., Shakopee, MN, USA)와 차압계(2051CD2; Rosemount Inc.)를 이용하여 측정하였다. 온도 측정에 의한 불확도를 최소화하기 위하여 정밀한 온도 교정을 통하여 열전대의 온도 측정 불확도가 0.1 K 내로 수렴하도록 하였다. 측정한 온도 및 압력을 이용하여 극저온 채널과 열원 채널의 입/출구 엔탈피를 산출하였으며 이를 바탕으로 극저온 채널과 열원 채널에서의 열전달량을 다음과 같이 산출하였다.

여기서 Q_c, Q_h는 각각 극저온 채널과 열원 채널에서의 열전달량이고 m˙c, m˙h는 각각 극저온 채널과 열원 채널 내 유체의 유량, h_c,in, h_c,out는 극저온 채널 내 유체의 입/출구 엔탈피, h_h,in, h_h,out는 열원 채널 내 유체의 입/출구 엔탈피를 나타낸다. 열전달량의 경우 PCHE가 정상 상태를 20분 이상 유지하였을 때 평가하였다. 본 연구에서 저온 유체 및 열원 유체의 모든 열전달 물성의 경우 각각 REFPROP와 COOLPROP을 이용하여 계산하였으며 이를 반영하여 열전달량을 산출하였다.

PCHE 내에서 극저온 유체와 열원 유체가 열교환 할 경우 시간에 대한 온도 경향은 Fig. 3과 같다. 각각의 그래프의 경우 극저온 채널 입/출구 온도와 열원 채널 입/출구 온도에 해당한다. 시간이 지남에 따라 PCHE 입구로 공급되는 배관의 온도가 낮아지며 배관 내 액체질소로의 열 침입이 감소하게 된다. 최종적으로 질소는 액체 조건으로 PCHE 내 극저온 채널로 공급되며 열원 채널 내 글리콜 워터와 열교환하여 기체 상태로 PCHE에서 토출된다. 글리콜 워터의 경우 항온 순환 수조를 이용하였기 때문에 40℃의 일정한 온도 조건으로 PCHE 내 열원 채널로 공급되었으며 극저온 채널 내 질소와 열교환하여 비교적 낮은 온도 상태로 항온 수조로 재순환되었다.

Fig. 3. 
Temporal data on inlet and outlet temperatures of LN2 and GW

Fig. 4는 시간에 따른 극저온 채널과 열원 채널에서의 열전달량을 나타낸다. Fig. 4(a)에서 볼 수 있듯이 시간이 지남에 따라 극저온 채널의 열전달량과 열원 채널의 열전달량이 증가하여 정상 상태에 도달하였으며 그 뒤로는 차이가 유지된다. 극저온 채널의 열전달량과 열원 채널의 열전달량이 정상 상태에 도달했음에도 차이를 가지는 이유는 극저온 채널의 입구 건도가 일부 존재하기 때문으로 파악된다. 극저온 채널의 입구 조건을 완벽한 단상의 액체질소로 가정하였으나 열 침입으로 인하여 이상유동 상태의 입구조건이 되었을 경우 항상 가정한 열전달량보다 실제 열전달량이 작은 값을 가지게 된다. 하지만 Fig. 4(b)에서 볼 수 있듯이 극저온 채널과 열원 채널의 열전달량 차이는 미비하였다. 본 연구에서 수행한 다양한 실험 조건에서 극저온 채널과 열원 채널의 열전달량 차이는 Fig. 5와 같이 약 15% 수준으로 유지되었다. 따라서 해당 연구에서는 극저온 채널의 입구 조건을 단상의 액체질소 조건으로 가정하여 분석하였다. 추가적으로 기존 연구들과 비교하여 타당한 열전달량 실험이 수행된 것을 확인하였다(Table 6).

Fig. 4. 
Temporal data of (a) heat transfer rate and (b) heat balance error

Fig. 5. 
Heat balance test

Fig. 6은 동결이 발생하지 않는 조건에서 액체질소의 유량과 글리콜 워터의 유량을 변화시켰을 때 발생하는 PCHE의 열전달량에 관한 등고선이다. 하얀색 점들의 경우 실험을 통해 측정한 결과이며 나머지 데이터의 경우 외삽법을 이용하여 추정하였다. 이를 통해 액체질소의 유량이 증가하거나 글리콜 워터의 유량이 증가함에 따라 PCHE의 열전달량은 증가하는 것을 확인할 수 있다. 그러나 글리콜 워터의 유량 증가에 비해 액체질소의 유량 증가가 PCHE의 열전달량 증가에 효과적인 것으로 파악되었다. 이는 글리콜 워터에 비해 액체질소의 열전달계수가 유량에 따라 크게 변하기 때문으로 예상된다. 글리콜 워터의 높은 점도로 인하여 실험을 수행한 모든 영역에서 열원 채널은 층류 조건에서 운전된다. 하지만 액체질소가 흐르는 저온 채널의 경우 상변화가 일어나서 기체 상태로 존재할 뿐 아니라 액체질소의 점도가 글리콜 워터에 비해 상대적으로 낮아 비교적 높은 Reynolds number (Re) 값을 가진다. 따라서 Fig. 7에서 볼 수 있듯이 저온 채널의 Re 값이 열원 채널에 비해 굉장히 높으며 대부분 난류 조건에서 운전됨을 확인할 수 있다. 그렇기 때문에 글리콜 워터의 유량 변화에 비해 액체질소의 열전달계수가 유량에 따라 크게 변한 것으로 예상되며 이로 인해 액체질소의 유량이 증가할 때 PCHE 열전달량이 크게 증가한 것으로 판단된다. 즉 열전달계수 및 열전달량 변화 폭이 적은 글리콜 워터의 유량을 변화시키며 Fig. 8과 같이 동결 현상을 관측하였다.

Fig. 6. 
Heat transfer rate of the fabricated PCHE for various mass flow rates of LN2 and that of GW

Fig. 7. 
Reynolds number of each channels for various temperature conditions

Fig. 8. 
Temporal data on temperatures and mass flow rates. (a) Non-freezing and (b) freezing condition

Fig. 8은 극저온 채널로 유입되는 액체질소의 질량 유량이 동일한 조건에서 열원 채널로 유입되는 글리콜 워터의 질량 유량을 점차 감소시켰을 때 발생하는 동결 현상을 나타낸다. Fig. 3과 유사하게 극저온채널과 열원 채널의 입/출구 온도를 나타낼 뿐 아니라 액체질소와 글리콜 워터의 유량을 나타낸다. 동결이 일어나지 않는 조건의 경우 Fig. 8(a)와 같이 온도, 유량 모두 일정하게 정상 상태로 진입하는 것을 확인할 수 있다. 하지만 글리콜 워터의 유량이 감소하여 동결이 발생하는 조건에서는 저온 채널과 열원 채널이 열교환을 하지 못하게 되고 Fig. 8(b)와 같이 저온 채널의 출구 온도가 점차 감소함을 확인할 수 있다. 또한 동결이 발생하는 영역이 증가하게 되면서 글리콜 워터의 유량이 점진적으로 하락하게 된다. 글리콜 워터의 유량을 열원 채널의 압력 강하량과 비교하여 보게 되면 Fig. 9와 같은 경향을 확인할 수 있다. 동결이 발생하지 않을 경우 Fig. 9(a)와 같이 열원 채널의 압력 강햐량이 일정하게 유지되는 반면 동결이 발생할 경우 Fig. 9(b)와 같이 글리콜 워터의 유량 감소와 함께 열원 채널의 압력 강하량이 증가함을 확인하였다.

Fig. 9. 
Temporal data on pressures and mass flow rates. (a) Non-freezing and (b) freezing condition

따라서 글리콜 워터의 동결 현상을 대표할 수 있는 가장 간단한 인자로써 열원 채널의 차압을 다양한 조건에서 비교하였다. Fig. 10은 다양한 액체 질소 유량과 글리콜 워터 유량 조건에서 열원 채널의 차압 등고선이다. 비어 있는 점들의 경우 실험적으로 얻은 열원 채널의 차압 데이터이며 이를 바탕으로 외삽법을 이용하여 나머지 영역에서의 열원 채널의 차압을 추정하였다. Fig. 10에서 볼 수 있듯이 액체 질소의 유량이 증가함에 따라 동결이 발생될 가능성이 증가하게 되어 열원 채널의 차압이 증가하게 된다. 글리콜 워터의 유량이 증가할 경우 열원 채널 내 유속이 증가하여 열원 채널의 차압은 증가하게 된다. 하지만 특정 글리콜 워터의 유량이 특정한 값 이하로 감소하게 될 경우 마찬가지로 동결이 발생할 가능성이 증가하게 되어 열원 채널의 차압은 증가하게 된다. 따라서 PCHE 내 동결을 회피하기 위해서는 액체질소의 유량에 따라 적절한 글리콜 워터의 유량을 선정하여 운전해야 할 것으로 판단된다.

Fig. 10. 
Pressure drop of GW channels for various mass flow rates of LN2 and that of GW

Fig. 11은 채널 내 동결 현상을 회피하기 위하여 시험 가능한 액체질소와 글리콜 워터의 유량 범위를 나타낸 차압 등고선이다. 각각 점선으로 되어 있는 등고선의 경우 실험적으로 동결 현상이 발생하기 직전의 등고선과 안정적인 열교환기 운영을 위하여 제시하는 영역을 의미한다. 해당 등고선보다 낮은 열원 채널 차압을 유지하는 영역에서 열교환기를 작동할 경우 동결 위험을 회피하며 안정적으로 설비를 운영할 수 있을 것으로 예상된다.

Fig. 11. 
Allowable mass flow rates of LN2 and that of GW for avoiding freezing phenomena

본 연구에서 수행한 PCHE와 유사한 열교환기를 비슷한 영역에서 운영하게 되면 Fig. 11의 차압 등고선을 참고하여 동결을 회피할 수 있도록 글리콜 워터의 유량을 선정할 수 있을 것이다. 하지만 본 연구에서 수행한 실험과 비교하였을 때 사용된 유체의 종류가 상이하거나 PCHE 내부의 채널 직경이 크게 차이가 난다면 내부 채널에서의 열전달 현상이 변화하게 되어 동결을 회피할 수 있는 적정 글리콜 워터 유량을 재설계해야 할 것이다. Fig. 12와 같이 질소와 수소의 물성치를 다양한 온도에서 비교할 경우 큰 차이가 난다는 것을 확인할 수 있다. 다른 물성치에 비해 정압 비열 값과 열전도도 값이 두 유체 간 큰 차이가 나게 된다. 두 물성치의 경우 열전달 현상에 직간접적으로 크게 영향을 주기 때문에 다양한 유체를 사용할 때에 동결 현상을 회피하기 위해서는 추가적인 연구가 필요할 것으로 판단된다. 특히 동결현상을 나타낼 수 있는 무차원수를 이용하여 무차원화하여 해석할 수 있다면 다양한 유체, 입구 조건 등을 반영하여 동결 현상을 회피할 수 있는 디자인 가이드라인을 제시할 수 있을 것이다.

Fig. 12. 
Thermodynamic properties of nitrogen and that of hydrogen for various temperature