Introduction

최근 전 세계적으로 화석연료의 고갈로 인해 유가 상승과 수급 불안이 동반되고 있는 가운데, 특히 우리나라의 경우 화석 에너지 자원이 절대 부족하여 해외 의존도가 매우 높아 이를 경감시키고 에너지 수급을 안정시키기 위해 국내 부존 에너지를 활용하는 것이 매우 중요하다. 신재생 에너지는 화석 에너지 고갈문제와 국제 환경규제에 대비해 적극 개발해야할 에너지 자원이며, 신재생 에너지 자원 중 수력 에너지는 국내의 자연적인 지역 조건과 조화를 이루고 부존 잠재량이 많아 효율적으로 이용이 가능하다. 수자원을 이용한 수력 발전에 적용되는 수차 중 하나인 프란시스 수차(Francis turbine)는 폭넓은 낙차(H)와 비속도(Ns) 범위에 적용이 가능하고, 상대적으로 높은 효율과 함께 쉬운 구조적 설계로 다른 수차 타입보다 적용에 유리한 장점을 가지고 있다.^(1~2)

저비속도 영역의 원심식 타입 유체기계의 구성요소 중, 프란시스 수차의 스파이럴 케이싱(spiral casing)은 입구에서 유입되는 유체를 모아서 스테이 베인(stay vane)으로 유도하는 역할을 하며, 형상에 따라 유입되는 유체의 흐름과 내부 유동 특성은 달라진다. 케이싱의 형상에 따라 유체가 유도되는 과정에서 케이싱 내 손실을 유발하기도 하며, 부적절하게 유도된 유체는 스테이 베인으로 유입되는 유동에도 영향을 주어 추가적인 손실을 야기한다. 따라서 스파이럴 케이싱 내 손실을 감소시키고, 유체 흐름의 적절한 유도를 위하여 케이싱의 끝부분 설부에 위치하여 입구쪽으로 유입되는 유체 흐름을 차단하여 순환유동을 막아주는 누설벽 영향을 파악하는 것은 필수적이다.

이러한 연구 동향의 일환으로, Kim 등⁽³⁾은 삼차원 수치해석을 통해 주거환기용 시로코홴 내 스크롤 케이싱의 내부 유동 및 소음원을 파악하여 내부 유동 특성을 규명하였고, Heo 등⁽⁴⁾은 시로코홴의 스크롤 케이싱 설부 형상과 관련된 변수에 대하여 최적화 설계릍 통해 공력 성능을 향상시켰다. 더불어 Zhen 등⁽⁵⁾은 프란시스 수차의 러너 후연 형상, 케이싱 설부 형상 및 스테이 베인 수 등 내부 유로 형상에 따른 내부 유동 특성을 파악하였다.

본 연구에서는 프란시스 수차의 스파이럴 케이싱 내 설부에서의 누설벽에 따른 수력학적 특성을 파악하고자 삼차원 정상상태 Reynolds-averaged Navier-Stokes (RANS) 방정식에 의한 수치해석을 수행하였다. 실물 프란시스 수차에 대한 수치해석의 타당성을 검증하기 위해 성능시험을 토대로 해석결과와 비교하였고, 수치해석을 통한 프란시스 수차의 스파이럴 케이싱 내 설부의 누설벽에 따른 유동구조 및 수력학적 성능을 파악하여 내부 유동 특성을 규명하고자 하였다. 더불어 실물 프란시스 수차를 토대로 상사된 모델 수차에 대해서도 누설벽의 영향을 비교 분석하기 위해, 모델 수차에 대해 상사된 정격 운전조건으로 수치해석을 수행하여 성능 및 내부유동을 분석하고자 하였다.

본 연구에서는 비속도 330급 프란시스 수차를 대상으로 삼차원 수치해석을 수행하였으며, Fig. 1에 스파이럴 케이싱 및 러너 등 주요 유동 영역에 대하여 삼차원 형상을 도시하였다. 실물 프란시스 수차 설계점에서의 낙차와 유량은 IEC 60193 규격의 낙차 계수와 유량 계수인 식 (1)~(3)를 사용하여 나타내었고, 모델 수차의 경우에는 실물 수차를 토대로 상사를 수행하였으므로 속도, 낙차 및 유량 계수는 같은 값을 나타낸다.⁽⁶⁾ 주요제원은 Table 1에 나타냈으며, 정격 운전조건을 토대로 내부 유동 특성을 분석하였다.

Fig. 1 
3D modeling view of a Francis hydro turbine

여기서, g, H, n, D 및 Q는 각각 중력가속도, 낙차, 회전속도, 수차의 대표직경 및 유량을 나타낸다.

Fig. 2는 프란시스 수차의 스파이럴 케이싱 내 누설벽의 위치와 형상을 도식화 하고 있다. Fig. 2에 나타난 바와 같이 설부에서의 누설벽은 스파이럴 케이싱의 끝단과 입구의 케이싱 부분 사이에 위치하고, 설부의 누설벽은 케이싱 내에서 유체의 순환을 막아주는 역할을 한다.

Fig. 2 
The leakage wall of the tongue in a spiral casing

본 연구에서는 ANSYS 사의 상용 소프트웨어인 ANSYS CFX-17.1을 사용하여 프란시스 수차의 내부 유동장을 정상 상태로 가정하여 해석을 수행하였다.⁽⁷⁾ 프란시스 수차의 격자생성을 위해 러너, 스테이 베인, 가이드 베인 및 그 외 영역에서 각각 Turbo-Grid와 CFX-ICEM을 사용하였다. 수치해석의 경계조건 정의, 유동해석 및 결과정리는 각각 CFX-Pre, CFX-Solver 및 CFX-Post를 사용하여 수행하였다.

정상상태의 수치해석 시 사용된 지배방정식들은 유한체적법으로 이산화 되었으며, 프란시스 수차 해석의 작동유체는 25℃ 물로 하였다. 입구의 경계조건은 균일한 상태의 전압력을 주었으며, 출구조건으로는 정압력 조건을 주었다.

수치해석을 위해 해석대상인 프란시스 수차에 대해 동일한 낙차 조건과 여러 가이드 베인 개도를 적용하여, 정격 유량조건 외 저유량 조건에 대한 해석을 수행하였다. 대부분의 유체기계는 해석 시 계산시간의 단축과 수렴성 향상을 위해 인접하는 두 블레이드 사이의 유동장이 회전방향에 대해 주기적으로 형성되어 있다고 가정하여 하나의 유로에 대해서 주기조건을 적용하여 해석을 수행하지만, 본 해석대상인 프란시스 수차는 스테이 베인과 확장되어지는 스파이럴 케이싱 사이에서 발생하는 내부 유동의 상호작용에 의해 주기적이라고 가정될 수 없기 때문에 러너, 스테이 베인 및 가이드 베인 내 모든 블레이드를 포함한 전 영역에 대하여 해석을 수행하였다.

정상상태의 수치해석 시 회전하는 러너와 정지상태의 가이드 베인 및 드래프트 튜브(draft tube) 사이의 경계면은 격자면의 계산값을 그대로 정지좌표계에 입력값으로 사용하는 방식인 Frozen-rotor 방식을 적용하였다. 해석 시 사용된 난류모델로는 SST (Shear Stress Transport)모델을 사용하였다. 이 모델은 유동박리 등의 예측에 적절하며, 자유 흐름 영역뿐만 아니라 벽면 인근 영역 해석에 정확성을 보이는 모델로서 다른 2-방정식 난류모델들에 비해 유체기계해석에 적절한 모델이다.⁽⁸⁾

회전하는 러너 블레이드와 스테이 베인 및 가이드 베인 표면 인근에는 O형 격자계를 사용하였고, 그 외의 스파이럴 케이싱 및 드래프트 튜브 영역에는 벽면 근처에는 H형 격자계를 사용하였다. 수치해석 시 프란시스 수차 모델 러너의 벽근처 영역에서 저 레이놀즈 전단이송모델을 적용하기 위해 첫 번째 격자점의 y+≤2를 유지하도록 하였다. 프란시스 수차 전 영역에 대하여 육면체 격자계로 구성하였으며, 격자계 구성은 Fig. 3에 나타난 바와 같다.

Fig. 3 
Computational grids of Francis hydro turbine

본 연구에서는 스파이럴 케이싱 설부에서의 누설벽이 프란시스 수차의 수력학적 성능에 미치는 영향을 파악하고자 정상상태 삼차원 RANS 해석을 수행하였다. 성능시험을 토대로 수치해석의 타당성을 검증하기 위하여 정상상태의 RANS 해석결과와 성능시험 결과를 비교하였다. 실물 프란시스 수차의 정격 운전조건에서 정상상태의 삼차원 해석을 수행하여 누설벽이 없을 때, 유입부인 스파이럴 케이싱 내 전압력 분포 차이를 통해 수차에 대한 전체 손실량이 증가하여 결국 성능이 감소하는 것을 볼 수 있었다. 특히, 누설벽 장착에 따라 스파이럴 케이싱 내 전반적인 압력이 상승함에 따라 손실 성분의 저감과 함께 수차의 출력이 증가하는 것을 확인할 수 있었다. 실물 수차를 토대로 상사된 모델 수차의 누설벽에 따른 영향을 확인하였고, 그에 따른 결과로서 모델 수차에서도 누설벽이 없을 때 전체 전압력 차이에 의한 손실량이 증가하고, 이에 따라 성능이 감소함을 확인 할 수 있었다. 특히, 누설벽의 영향은 실물 및 모델 수차에 대해 유사한 수력학적 영향을 미쳤으며, 이로 인해 비슷한 경향의 전압력 분포를 나타내는 것을 볼 수 있었다. 본 연구를 통해 스파이럴 케이싱 설부의 누설벽에 따른 성능 및 내부 유동 특성을 관찰 할 수 있었으며, 누설벽 장착에 따라서 전압력 분포 차이에 의해 손실이 감소하고 성능이 향상됨을 확인 하였다. 본 연구의 결과를 토대로 본다면 향후 프란시스 수차의 고효율 설계를 위해서는 스파이럴 케이싱 설부의 누설벽 유무는 설계인자로서 고려될 필요가 있다.