Introduction

과거 양수발전시스템은 전력수요가 낮은 경부하 시간대의 유휴전력을 이용하여 하부 댐의 물을 상부 댐으로 이송시켜 에너지를 저장한 후, 전력수요가 높은 첨두부하 시간에 방류하여 발전하는 방식으로 운전되었다. 특히, 빠른 기동 및 신속한 부하 증감발 운전 특성으로 비상전력을 신속히 공급하는 역할로 운전되고 있었으나, 최근에는 간헐성 신재생 에너지원인 태양광 및 풍력 발전 등의 계통 유입으로 안정적인 전력 공급을 위해 전력 계통 백업과 주파수 조정이 더 큰 역할로 대두되고 있다. 이러한 양수발전의 운전 환경 변화로 펌프터빈은 복잡한 유동 및 압력이 유발되는 탈설계점의 운전이 지속되고 있으며, 이는 잦은 고장으로 인한 시스템의 비상정지 발생 등 설비의 수명감소로 이어지고 있다.^(1~3)

양수발전용 펌프터빈의 터빈 모드는 탈설계점의 운전 조건에서 복잡한 유동 구조 및 비정상 압력 특성을 나타내며, 특히 저유량 운전 조건에서는 흡출관(Draft tube) 내 세차 와류(Precessional vortex)가 발생하기도 한다. 이러한 세차 특성을 가진 와류는 흡출관 내 러너의 출구에서 발생하며, 러너 콘의 중심선 축으로 그 둘레를 회전하는 와류로, 러너 회전 주파수의 약 0.2~0.4배의 주파수를 나타낸다.⁽⁴⁾ 또한, 압력 유발의 원천으로 시스템과 공진을 일으키는 경우, 토크와 출력의 맥동을 유발시켜 성능을 저하시킬 뿐 만 아니라 진동 및 소음을 발생시켜 펌프터빈 시스템의 운전을 불안정하게 한다. 따라서 최근 요구되어지고 있는 펌프터빈의 안정적인 운전 영역 확장 및 신뢰성 향상을 위해, 저유량 운전 조건에서 발생하는 유동 및 압력 특성에 대한 이해와 검토는 필수적이다.

이와 관련된 연구로서, Wang 등⁽⁵⁾은 펌프터빈의 저부하 조건에서 흡출관 내 와류에 대한 특성을 일직선 및 엘보(elbow)타입 형상에 대해 Very Large Eddy Simulation (VLES) 및 Shear Stress Transport (SST) 난류 모델을 사용한 Reynolds-averaged Navier Stokes (RANS) 해석에 대한 수치해석적 분석을 통해 비교 및 검토하였으며, 또한 와류 억제를 위한 흡출관 내 물 제트 분사에 대한 영향성도 함께 규명하였다. Liao 등⁽⁶⁾은 펌프터빈의 정상 및 비정상 상태 수치해석 시 작동 유체인 물의 압축성 고려 여부의 비교를 통해, 터빈 모드의 저유량 운전 조건에서 러너 및 흡출관의 내부 유동 및 압력 맥동 특성의 차이를 분석하였다. Kirschner 등⁽⁷⁾은 펌프터빈 모델의 저유량부터 고유량 운전 조건에서 흡출관 내부의 와류 구조의 유동 가시화 및 압력 특성에 대해 모델 시험을 통해 확인하였으며, 수치해석적 분석으로 각 운전조건에서 나타나는 흡출관 내 와류에 따른 압력 맥동 특성 및 반경방향 힘에 대한 연관성을 규명하였다. 이와 같이 최근 펌프터빈의 터빈 모드 내 저유량 운전 조건에서 흡출관 내 발생하는 와류의 특성에 대한 다양한 연구가 수행되고 있지만, 펌프터빈의 저유량 범위 내에서 유량에 따른 와류의 발생 범위와 그에 따른 내부유동 특성, 와류의 크기 및 압력 맥동 등을 체계적으로 분석하는 연구는 수행되고 있지 않다.

본 연구에서는 양수발전용 펌프터빈의 터빈 모드 내 저유량 운전 범위에서 선회수(Swirl number, S) 정의를 적용하여, 선회수에 따른 세차 와류의 특성을 파악하기 위해 2상 유동(물과 증기) 3차원 정상 및 비정상 상태에 대한 RANS 해석을 수행하였다. 수치해석은 펌프터빈의 터빈 모드 최고효율점(Best efficiency point, BEP)의 낙차 조건과 가이드 베인 각도에 따른 여러 유량 조건을 적용하여 흡출관 유로 내 유동 구조 및 비정상 압력 특성을 규명하고자 하였다.

본 연구에서는 터빈 모드에서 비속도 120급 펌프터빈에 대해 삼차원 수치해석을 수행하였으며, Fig. 1에 펌프터빈의 주요 유동 영역 및 유동 방향에 대해 삼차원 형상을 도식화 하였다. Table 1에 펌프터빈의 주요 제원을 나타내었고, 비속도는 식 (1)을 사용하여 나타냈으며, 에너지, 회전속도 및 유량 계수는 IEC60193 규격에 제시된 무차원 계수인 식(2) ~ (4)를 사용하여 나타내었다.⁽⁸⁾

Fig. 1 
Three-dimensional (3D) model view of a pump-turbine

여기서, n, P, H, E, D 및 Q는 각각 회전속도, 출력, 낙차, 에너지, 대표직경 및 유량을 나타낸다. 펌프터빈 러너의 날개는 7매이며, 가이드 베인(Guide vane) 및 스테이 베인(Stay vane)은 각각 20매로 구성되어있다.

본 연구에서는 펌프터빈의 터빈 모드 내 비압축성 내부 유동장에 대한 정상 및 비정상 수치해석을 수행하기 위해, ANSYS사의 상용 소프트웨어인 ANSYS CFX-18.2를 사용하였다.⁽⁹⁾ 펌프터빈의 러너 블레이드에 대한 형상 정의는 Blade-Gen을 사용하였으며, 격자 생성은 Turbo-Grid 및 ICEM-CFD, 수치해석을 위한 경계 조건 설정은 CFX-Pre, 그리고 CFX-Solver와 CFX-Post를 사용하여 각각 유동 해석 및 결과 분석을 수행하였다.

펌프터빈의 내부유동은 3차원 비압축성 유체로 가정하고 정상 및 비정상 상태에 대한 RANS 방정식을 이용하여 유동해석이 수행되었고, 유한체적법(FVM)으로 이산화된 지배방정식을 사용하여 수치해석을 계산하였다.

일반적으로 축을 중심으로 회전하며 축대칭 구조를 가진 유체기계는 해석 시간 단축과 수렴성의 향상을 위해 날개 1매의 유로에 주기조건(Periodic condition)을 설정하여 해석을 수행하지만, 펌프터빈의 스파이럴 케이싱(Spiral casing)은 터빈모드 유동방향으로 입구에서 점차 축소되어지는 유동 단면 특성과 케이싱 설부(Casing tongue)의 형상이 스테이 베인과 스테이베인 사이에서 내부유동의 상호작용이 주기적이라고 가정될 수 없는 이유로 러너 블레이드 등 전 영역을 수치해석에 적용하였다.⁽¹⁰⁾

경계조건 설정을 위한 작동 유체는 2상 유동 해석을 위해 25℃ 물과 증기(Vapor)를 사용하였고, 캐비테이션 모델은 증기의 발생과 응축을 제어하는 반응 속도식을 기반으로, 물의 기포(Bubble)의 발생 과정을 묘사하는 Reyleigh Plesset 모델을 사용하였다.

펌프터빈의 입구와 출구는 상부 및 하부 저수지의 수위를 고려한 전압력과 정압력 조건을 각각 부여하였다. 난류모델은 유동박리를 잘 예측하며 또한 자유 흐름 영역과 벽면 근처 영역의 해석에 정확성을 보이는 SST 모델을 사용하였다.⁽¹¹,¹²⁾ Fig. 2는 펌프터빈의 격자계 구성을 나타내며, 스파이럴 케이싱과 스테이 베인은 사면체 격자로 이루어져 있고, 그 외의 영역은 육면체 격자로 이루어져 있다. 상대적으로 복잡한 유동구조를 가지는 러너 블레이드 날개 표면 인근에는 벽함수 처리 조건을 적용하기 위해 O형 격자계를 사용하여 첫 번째 격자점에서 y+≤100을 유지하도록 하였다.⁽³⁾

Fig. 2 
Computational grids of a pump-turbine

펌프터빈의 정상 상태 수치해석 시 회전자와 고정자 사이의 경계면 조건은 경계면을 통해 원주방향으로 평균된 유동 값을 입력하는 Stage 방식을 사용하였고, 비정상 상태 해석에는 회전자와 고정자의 유로사이의 실제 순간적인 유동에 대한 상호작용을 예측하는 Transient-Rotor-Stator 방식을 적용하여 해석을 진행하였다. 비정상 상태 수치해석의 수렴성 향상 및 해석 시간의 단축을 위해, 정상 상태의 결과를 초기 값으로 부여하였다. 러너 1회전 시간은 t_rev=0.16667s로, 비정상 수치해석을 위한 전체 시간(Total time)과 타임 스텝(Time step)은 각각 5t_rev 및 0.000833t_rev로 설정하여 해석을 수행하였으며, 수렴성의 향상을 위한 반복횟수(Loops coefficient)는 5회로 설정하였다.

Fig. 3은 펌프터빈의 비정상 상태 해석 시 나타나는 내부유동 및 세차 와류 특성에 따른 흡출관 벽면의 비정상 압력특성을 확인하기 위한 압력 측정 위치를 나타내고 있다. 흡출관 cone 벽면에 일정한 간격으로 4개의 방향(L1~L4)에 대해 측정점을 적용하였으며, 본 연구에서는 p1~p4에 대한 비정상 압력 특성을 비교 및 분석 하였다

Fig. 3 
Unsteady pressure measuring points on the draft tube cone

본 연구에서는 양수발전용 펌프터빈의 터빈 모드 내 선회강도에 따른 흡출관 내 세차 와류에 대한 내부 유동 및 비정상 압력 특성을 파악하기 위해 2상 유동 3차원 정상 및 비정상 상태에 대한 해석을 수행하였다. 선회 강도를 나타내는 선회수의 정의를 적용하여 유량에 따라 선회 강도 및 흡출관 내 세차 와류의 유동 구조 및 발달 범위를 확인하였다. 또한, 속도삼각형 분포를 통해 저유량 영역에서 러너 출구의 유동 특성을 통회 선회 강도 발생의 원인을 확인하였으며, 흡출관 내 유선분포 비교를 통해 매우 강한 선회 강도를 가진 저유량 영역에서는 세차 와류가 발생하는 위치에 매우 복잡한 유동 구조를 나타내어 선회 강도가 매우 높음에도 불구하고 와류가 발생하지 않음을 확인하였다. 이러한 흡출관 내 유동 특성에 대한 비정상 압력 특성을 FFT 분석을 통하여 확인하였고, 저유량 영역에서 상대적으로 낮은 주파수 영역의 높은 압력 특성을 볼 수 있었다. 특히, 세차 와류가 발생했을 때 상대적으로 더욱 낮은 주파수 영역인 0.2~0.4f_n에서 우세한 압력 특성을 나타내며, 흡출관에서 유동 방향으로 진행할 때 와류가 발생하지 않았을 때 보다 상대적으로 높은 압력 특성을 보이는 것을 확인하였다. 이와 같은 결과로 비추어 볼 때, 펌프터빈의 터빈 모드 내 저유량 영역 중, 특정 범위인 선회수 2.26 이내에서 세차 와류가 발생하고, 발생된 세차 와류는 흡출관 내에 높은 비정상 압력 특성으로 흡출관 내 파속 감소 및 시스템 공진의 발생으로 인한 압력 맥동의 증폭으로 시스템 운전에 영향을 미칠 수 있으므로, 펌프터빈의 탈설계점 내 보다 안정적인 운전을 위해 흡출관 내 발생하는 세차 와류 특성을 고려할 필요가 있다.