수소는 무게 대비 에너지 밀도가 높을 뿐 아니라 환경오염물질 배출이 없는 친환경 연료이자 최종 부산물의 재사용이 가능하여 대체에너지원으로 크게 각광을 받고 있다. 수소를 에너지원으로 활용하기 위해서는 수소의 생산, 저장과 같은 기술을 활용해야 하며 저장 및 운송에 대한 기술은 새로운 방법들이 제안되고 연구되어 왔다. 하지만 각기 다른 장단점으로 인해 활용도가 달라지고 있다. 수소 저장 기술방법에는 압축 기체 저장법, 액체수소 저장법 및 수소저장소재를 이용한 저장법 등이 있다¹⁾. 압축 기체 저장법과 액체수소 저장법은 현재 일반적으로 적용되는 방법이지만, 압축 기체 저장법의 경우 상대적으로 체적 당 저장밀도가 낮고 고압의 가스용기를 사용함으로써 위험성이 크다는 단점이 존재한다. 액체수소 저장법은 체적 당 저장밀도가 압축 기체 저장법보다 크지만 수소의 액화점의 -235℃ 이하로 유지되어야 하는 특수 단열 시스템이 요구되고 장시간의 저장에 기술적 한계가 존재한다(Table 1). 따라서 현재의 기술과 활용성 면에서 고압기체 저장법이 주로 사용되고 있으며 이를 통해 저장 운송 시스템에 적용방법이 모색되고 있다. 통상적으로 운송용으로 사용되는 가압 기체 수소저장시스템은 수소기체를 고압으로 압축하여 제한된 체적의 용기에 저장하는 방식으로 약 70 MPa 정도의 고압 수소기체를 저장하고 공급한다. 수소저장시스템은 탱크 주위의 화재와 같은 탱크의 손상이 야기되는 상황이나 정기적인 유지 보수를 수행해야 하는 상황 시 내부의 압축 수소를 외부로 안전하게 배출하는 기능을 충족해야 한다. 내부의 고압 압축 수소를 외부 대기압에 배출 시, 갑작스런 고압가스 배출로 용기 내부 온도가 급강하하는 현상이 발생할 수 있고 급강하된 온도 환경으로 인해, 고분자로 성형된 라이너 및 그 주변을 와인딩하고 있는 탄소섬유 메트릭스(matrix)는 극저온 온도에 노출된다. 이런 환경에서는 고압 기체수소로 팽창되어 있던 용기가 빠른 속도로 수축 변형이 일어날 가능성이 발생한다. 따라서 배출 조건에서도 안정적인 저장 시스템의 운용을 위해 저장 용기의 특성변화를 예측할 수 있는 수치해석 기반의 검토가 필요하다. 가스의 배출에 대한 연구는 가스터빈²⁾에서부터 자동차 배출가스³⁾, 메탄가스 배출에 대한 연구⁴⁾ 등 다방면으로 진행되었으며 배출 조절 및 제어를 위한 연구⁵⁾들이 수행되었다.

수소 저장 탱크에서 내부의 고압 압축 수소의 거동에 관한 연구는 수소 밀도의 변화에 대한 연구⁶⁾에서부터 수소 충전에 대한 연구^7,8) 들까지 다방면으로 진행되었다. 수소 배출 조건에서 수소의 거동에 따른 저장 탱크의 변화에 관한 연구는 일부가 수행되었으나 수요나 필요성에 의한 응용 연구가 더 필요한 실정이다.

본 연구에서는 Type Ⅲ 수소 저장 탱크를 대상으로 고압 수소의 배출 조건에서 수소의 거동과 저장 탱크의 변화를 수치적으로 모사하였다. 특히 수소 유체의 정확한 모사를 위해 수소의 가스 모델(gas model)에 따른 정확성과 유체 거동의 영향성에 대해 논의하고 적합한 gas model을 선정하였다.

수소 배출에 대한 모델링을 수립하고 정식화하는 과정에서는 수치 해석의 해석 기법을 검증하는 과정이 필요하다. 따라서 존재하는 실험값과 비교 검증을 통해 검증과정이 진행되었다. 비교 검증은 Type 3 형태의 고압 탱크의 실험값을 토대로 진행되었고 Type 3 탱크를 모델로 선정하였다. 해석을 위한 고압 탱크는 74 L 용량의 실린더로 내부 길이가 849 mm, 내경이 354 mm로 되어 있으며 normal working pressure 70 MPa인 고압 탱크이다. 고압의 수소 배출을 통한 내부 상태 변화를 위한 수치 해석적 접근법은 고압 기체의 배출에 따른 탱크 내부의 압력 변화와 열전달을 동시에 해결하는 문제로 정식화할 수 있다.

3.3 해석 조건

본 해석은 압축성 난류 유동에 대한 비정상상태의 유동해석을 진행하였다. 해당 기체는 수소(hydrogen)를 사용하였으며, 난류 모델은 standard k-ε model을 사용하였다. 고압 용기에서 분출하는 압축 수소를 모사하기 위해서 전산 유체역학 상용프로그램인 Ansys Fluent ver 18. 2를 사용하였으며 해석영역은 Ansys workbench의 design modeler 및 mesh를 통해 모델링하였고 모델링 및 mesh 형상을 Fig. 2에 나타내었다. 격자는 해석의 안정적인 격자인 square mesh로 구성하였다. 해석의 신뢰도를 위해 residual value 값은 1×10^-6 이하 값에서 수렴성을 판단하였으며 매 time step 당 20회 이상 반복계산을 수행하여 해석의 정확성을 얻었다. 수소 저장 탱크에서 수소 배출(discharge) 수치해석은 실험값¹⁴⁾과 비교를 통해 타당성을 검증하였다. 실험값의 배출 조건¹⁴⁾은 Table 3과 같다.

Fig. 2. 
Modeling of storage tank

Table 3. 
Hydrogen discharge condition

Initial pressure (MPa)	36.3
Final pressure (MPa)	2.5
Discharge time (s)	1,140
Ambient temperature (K)	305
Initial temperature (K)	331

주어진 압축 수소 탱크의 배출조건에 따라(Table 3) 시간의 경과에 따른 압축 수소의 압력 및 온도 변화를 예측하기 위해 수치 해석을 진행하였고 실험값¹⁴⁾의 결과와 수치 해석의 결과를 비교하여 수치 해석 방법의 타당성을 검증하였다. 36.3 MPa 초기 압력과 331 K 초기 온도 상태의 가압수소는 약 1,100초 후 배출 시간의 경과 후 내부 압력이 2.5 MPa에 도달하였고 이 조건을 토대로 gas model에 대한 수치 해석을 진행하였다. 주어진 조건에서 압축된 수소가 배출됨에 따라 내부의 압력은 일률적으로 낮아지는 것을 Fig. 3에서 확인할 수 있다. 모델별 계산의 검증을 위해 실험값의 시간별 온도값과 해석의 시간별 온도값을 비교하여 오차를 계산하였다. 약 1,100초 후 배출 시간의 경과 후 내부 압력이 2.5 MPa에 도달하는 실험값과 비교했을 때 gas model별로 실험값과 다소 오차가 발생됨을 보인다. 그중 Redich-Kwong 모델이 가장 작은 오차가 발생함을 Tables 4, 5에서 보여준다. 또한 주어진 조건에서 압축된 수소가 배출됨에 따라 내부의 온도가 감소하는 변화를 Figs. 4, 5에서 볼 수 있다. 온도의 감소변화에 따라 초기(400초 미만)에서는 상대적으로 많은 온도 감소의 양상을 보이다가 후반(500초 이후)에는 다소 일정한 온도 감소 양상을 보이는 것을 확인할 수 있다. 이와 같은 거동의 변화 때문에 gas model에 따라 초기에는 실험값과 약 0.5% 미만의 차이를 보이나 후반에는 다소 작아진 약 0.2% 미만의 차이를 보인다. 초기에는 Soave-redich-kwong 모델과 Redich-kwong 모델이 실험값과 가장 작은 오차를 보인 반면 후반에는 Redich-Kwong이 가장 작은 오차를 보인다. 실험값의 압력 및 온도 변화의 값 비교를 통해 전체적으로 Redich-Kwong 모델이 가장 작은 오차를 보였다.

Fig. 3. 
Pressure distribution of hydrogen tank with gas model

Fig. 4. 
Temperature distribution of hydrogen tank distribution with gas model (initial state)

Fig. 5. 
Temperature distribution of hydrogen tank distribution with gas model (developed state)

본 연구에서는 Type Ⅲ 수소 저장 탱크를 대상으로 고압 수소의 배출조건에서 저장 탱크의 내부 변화(압력, 온도)를 수치적으로 모사하였다. 수소 가스 모델(gas model)에 따른 수치 해석 결과값을 실험값과 비교하였고 값의 차이를 정량적으로 분석하였다. 주어진 배출조건에서 실험값(온도, 압력값)과 비교를 통해 전체적으로 Redich-Kwong의 가스 모델이 가장 작은 오차를 보임을 알 수 있었다. 본 해석 방법론은 향후 진행될 다양한 타입의 수소 저장 탱크별 수소 주입 및 배출조건별 수소 거동 분석 연구에 이용될 것으로 보인다. 나아가 수소 배출조건에 따른 수소의 거동 분석과 이에 따른 저장 탱크의 구조적 영향성 연구에 활용될 것으로 예상된다.