Numerical simulations of transition and long-term response of a wind turbine airfoil

본 논문은 저레이놀즈수에서 FFA-W3 익형의 천이 및 장기 유동 응답을 분석하기 위해 Nek5000과 EllipSys를 이용한 수치 시뮬레이션을 제시하며, 켈빈-헬름홀츠 불안정성을 포착하는 EllipSys 코드의 능력을 검증하고 저주파 진동 및 잠재적인 버블 버스팅(bubble bursting)과 연관된 법선력 계수의 느리고 주기적인 변조를 식별한다.

핵심

풍력 터빈 블레이드를 거대하게 회전하는 날개라고 상상해 보세요. 비행기 날개와 마찬가지로, 효율적으로 작동하려면 매끄러운 기류가 필요합니다. 하지만 바람이 특정 각도로 날개에 부딪히면, 공기가 표면에서 떨어져 나가며 "걸려버릴" 수 있고, 이는 혼란스럽고 소용돌이치는 난기류를 만들어냅니다. 이 논문은 마치 첨단 풍동 실험과 같습니다. 다만 물리적인 모델을 사용하는 대신, 연구진은 슈퍼컴퓨터 안에 가상 모델을 구축하여 공기가 정확히 어떻게 움직이는지 관찰했습니다.

다음은 그들의 발견을 쉬운 개념으로 나누어 설명한 이야기입니다.

1. 가상 풍동

연구진은 거대한 풍력 터빈 블레이드(10MW급 터빈)의 특정 단면을 연구하고자 했습니다. 그들은 NEK5000과 ELLIPSYS라는 두 가지 서로 다른 컴퓨터 프로그램을 사용하여 블레이드 위로 흐르는 공기를 시뮬레이션했습니다.

NEK5000을 모든 미세한 디테일을 포착하지만 매우 느리고 비용이 많이 드는 '하이엔드 초정밀 카메라'라고 한다면, ELLIPSYS는 그보다 약간 더 빠르고 효율적인 '카메라'와 같습니다. 연구팀은 먼저 이 '더 빠른' 카메라(ELLIPSYS)가 '하이엔드' 카메라와 동일한 것을 볼 수 있는지 증명해야 했습니다. 그 결과, ELLIPSYS가 매끄러운 공기 속의 아주 작고 희미한 물결은 일부 놓치기도 했지만(너무 많이 매끄럽게 처리했기 때문), 성능에 실제로 중요한 크고 혼란스러운 소용돌이를 포착하는 데는 탁월하다는 것을 발견했습니다.

2. 터널의 폭은 얼마나 넓어야 하는가?

장기 시뮬레이션을 실행하기 전에, 그들은 가상의 "풍동"이 얼마나 넓어야 하는지 결정해야 했습니다. 터널이 너무 좁으면 공기를 압착하여 가짜 결과를 만들 수 있고, 너무 넓으면 컴퓨터 자원을 낭비하게 됩니다.

그들은 "좁은" 터널(날개 너비의 10%)을 "넓은" 터널(날개 너비의 20%)과 비교 테스트했습니다.

• 비유: 강물이 흐르는 모습을 관찰한다고 상상해 보세요. 만약 강의 좁은 구역만 본다면, 큰 파도를 놓치게 될까요?
• 결과: 그들은 좁은 터널로도 충분하다는 것을 발견했습니다. 큰 파도와 소용돌이가 좁은 공간에서도 완벽하게 형성되었습니다. 이는 더 작은, 더 좁은 시뮬레이션 박스를 사용함으로써 컴퓨터 시간을 크게 절약할 수 있음을 의미합니다.

3. "거품"과 "플랩"

연구의 가장 흥สนใจ로운 부분은 날개의 윗부분(흡입면)에서 일어났습니다.

• 박리 버블 (Separation Bubble): 공기가 날개 위를 흐를 때, 잠시 동안 공기가 떨어져 나가며 '층류 박리 버블(Laminar Separation Bubble, LSB)'이라고 불리는 작은 재순환 포켓을 만듭니다. 이것은 날개 표면에 생기는 아주 작은, 일시적인 소용돌이통과 같습니다.
• 불안정성: 이 버블 내부에서 공기는 단순히 정지해 있는 것이 아니라, 진동하며 파동(연못 위의 물결처럼)으로 말려 올라갑니다. 연구진은 이 파동이 어떻게 성장하는지 관찰했습니다. 그들은 이 버블 안의 주요 '롤러'가 **켈빈-헬름홀츠 모드(Kelvin-Helmholtz mode)**라고 불리는 일종의 불안정성임을 확인했습니다.
• 발견: 그들은 이 '더 빠른' 컴퓨터 프로그램(ELLIPSYS)이 이러한 파동이 어떻게 성장하고 버블이 어떻게 행동하는지를 정확하게 예측할 수 있으며, 초정밀 프로그램의 결과와 일치한다는 것을 확인했습니다.

4. 느리고 리드미컬한 맥동 (거대한 놀라움)

도구들을 검증한 후, 그들은 시뮬레이션을 매우 오랫동안(공기가 날개를 지나가는 것의 50배에 해당하는 시간 동안) 실행했습니다. 여기서 그들은 특별한 것을 발견했습니다.

공기가 빠르게 요동치고 있는 와중에도, 날개를 밀어내는 힘에서 매우 느리고 리드미컬한 맥동을 포착했습니다.

• 비유: 드럼 연주를 상상해 보세요. 공기의 빠른 소용돌이는 빠른 드럼 롤과 같습니다. 그들이 발견한 느린 맥동은 시뮬레이션 시간 기준으로 48초마다 한 번씩 발생하는 깊고 느린 심장 박동과 같습니다.
• 효과: 이 느린 심장 박동은 날개에 가해지는 힘을 약 10.5% 정도 위아래로 흔들리게 만들었습니다.
• 실제 터빈과의 연결: 이를 실제 회전하는 풍력 터빈으로 환산했을 때, 이 느린 맥동은 블레이드가 7.7번 완전히 회전할 때마다 한 번씩 발생한다는 것을 깨달았습니다.

5. 왜 이런 현상이 발생하는가?

연구진은 이 느린 맥동이 공기가 날개에서 "실속(stalling, 걸림)"했다가 다시 "탈실(unstalling, 놓아줌)"하는 주기에 의해 발생한다고 믿습니다.

• 주기: 공기가 걸리면서 커다란 버블을 만듭니다. 그 후, 무언가가 이 버블을 터뜨리는 트리거 역할을 하여 공기가 다시 매끄럽게 달라붙습니다. 그러면 다시 압력이 쌓이고, 버블이 형성되며, 이 주기가 반복됩니다.
• 트리거: 그들은 공기가 날개 위에서 매우 강력하게 역방향으로 소용돌이치기 때문에, 공기가 스스로 진동할 수밖에 없는 상태인 "절대적 불안정성(absolute instability)" 상태가 되기 때문이라고 추측합니다.

6. 결론

이 연구는 컴퓨터 모델링의 성공 사례입니다. 연구진은 (먼저 표준 모델과 대조하여 검증을 거친다면) 더 빠르고 효율적인 컴퓨터 프로그램(ELLIPSYS)이 복잡한 풍력 터빈 물리학을 연구하는 데 신뢰할 수 있다는 것을 입증했습니다.

그들은 두꺼운 풍력 터빈 블레이드 위에서도 약 8회전마다 한 번씩 발생하는 매우 느리고 리드미컬한 공기의 "호흡"이 존재한다는 것을 발견했습니다. 이 호흡은 양력(터빈을 돌리는 힘)을 크게 오르내리게 합니다. 이 느린 리듬을 이해하는 것은 중요합니다. 당장 터빈을 파손시키지는 않더라도, 이러한 크고 느린 힘의 변화는 수년간의 운용 과정에서 재료에 피로를 누적시킬 수 있기 때문입니다.

요약하자면: 그들은 가상 풍동을 구축했고, 빠른 컴퓨터가 그 역할을 수행할 수 있음을 증명했으며, 풍력 터빈 블레이드 표면에서 공기 거품이 형성되고 터지는 현상으로 인해 발생하는 느리고 리드미컬한 "심장 박동"을 발견했습니다.

기술 요약

기술 요약: 풍력 터빈 에어포일의 천이 및 장기 응답에 관한 수치 시뮬레이션

문제 정의
높은 유동 입사각과 두꺼운 에어포일 형상으로 인해 발생하는 풍력 터터빈 블레이드의 유동 박리는 성능과 구조적 하중에 상당한 영향을 미친다. 선행 연구들은 와류 방출(wake shedding), 켈빈-헬름홀츠(Kelvin-Helmholtz, KH) 모드, 저주파 진동(LFO)을 포함하여 박리된 유동에서 나타나는 다수의 고유 진동수를 확인하였다. LFO는 강한 양력 변화를 특징으로 하며, 종종 층류 박리 버블(laminar separation bubble, LSB) 내의 절대 불안정성(absolute instability)과 연관된 주기적인 스톨(stalling) 및 탈스톨(unstalling) 현상에 기인하는 것으로 간주된다. LFO는 다양한 에어포일 연구에서 관찰되어 왔으나, 두꺼운 풍력 터빈 에어포일의 낮은 받음각(AoA)에서 이러한 현상이 발생하는지, 그리고 해당 조건에서 이를 구동하는 구체적인 메커니즘이 무엇인지에 대해서는 추가적인 조사가 필요하다. 본 연구는 현 현상(chord Reynolds number, $Re_c$)이 $1 \times 10^5$인 DTU 10-MW 참조 풍력 터빈 에어포일(FFA-W3 시리즈) 단면에 초점을 맞춘다.

방법론
본 연구는 서브그리드 스케일 모델 없이 비압축성 나비에-스토크스 방정식을 해결하는 벽면 해상 대형 에디 시뮬레이션(wall-resolved large eddy simulations, LES)을 채택한다. 두 가지 솔버가 사용된다:

1. NEK5000: 최소한의 수치 소산과 높은 정확도로 알려진 스펙트럴 엘리먼트 방법(spectral element method, SEM) 코드이다.
2. ELLIPSYS: DTU/Risø에서 개발된 유한 체적법 코드로, 전통적으로 RANS에 사용되지만 여기서는 암시적 LES(implicit LES)에 적용되었다.

본 연구는 ELLIPSYS를 NEK5000에 대해 체계적으로 검증하는 과정을 포함한다. 주요 방법론적 단계는 다음과 같다:

• 도메인 민감도: 관련 유동 구조의 해상도를 보장하기 위해 $L_z = 0.1c$와 $L_z = 0.2c$에서의 시뮬레이션을 비교하여 스팬 방향 도메인 폭($L_z$)을 평가한다.
• 솔버 검증: 두 솔버 간의 평균 유동장, 섭동 진화 및 스펙트럼 특성을 비교한다.
• 안정성 분석: 파라볼릭 안정성 방정식(Parabolized Stability Equations, PSE)과 스펙트럴 적절 직교 분해(Spectral Proper Orthogonal Decomposition, SPOD)를 활용하여 불안정성(Tollmien-Schlichting 및 KH 모드)의 성장과 결합 모드의 공간적 구조를 분석한다.
• 장기 적분: 저주파 현상을 포착하기 위해 계산 효율이 더 높은 ELLIPSYS 솔버를 사용하여 50 유동 통과 시간($T = 50c/U_\infty$) 동안의 유동 진화를 시뮬레이션한다.

주요 결과

• 도메인 폭: 주요 교란의 진화를 포착하고 시간 평균 유동 통계치를 재현하는 데 $L_z = 0.1c$의 스팬 폭이 충분하다는 것이 밝혀졌다. 이는 더 넓은 도메인($L_z = 0.2c$)의 결과와 일치한다.
• 솔버 검증: ELLIPSYS는 평균 유동장과 가장 증폭된 섭동의 진화에 대해 NEK5000과 밀접한 일치를 보였다. 그러나 ELLIPSYS는 높은 수치 소산으로 인해 부착된 경계층 내의 톨미엔-슈리히팅(Tollien-Schlichting, TS) 파동의 진폭을 과소 예측하였다. 반대로, ELLIPSYS는 층류 박리 버블(LSB) 내의 KH 모드 진화를 정확하게 예측하였으며, 이는 PSE 예측과 잘 일치하였다.
• 천이 메커니로: 흡입면(suction side)의 천이는 박리된 전단층에서 발달하는 KH 불안정성에 의해 구동된다. 가장 증폭된 KH 모드의 주파수는 $f \approx 15$이며, 그 하류 진화는 두 솔버 및 PSE에 의해 잘 포착되었다.
• 저주파 진동(LFO): 장기 시뮬레이션 결과, 법선력 계수($C_N$)가 10.5%의 진폭과 48 유동 통과 시간의 주기를 가지고 느리게 변조됨을 발견하였다. 이는 주파수 $f = 0.021$ 및 스트로우할 수(Strouhal number, $St$)$St = 0.0012$에 해당한다.
  • 이 주파수는 일차 KH 불안정성 주파수보다 약 720배 낮다.
  • DTU 10-MW 터빈의 맥락에서, 이 주기는 약 7.7회의 블레이드 회전에 해당한다.
  • 관찰된 $St$값은 다른 에어포일에서 보고된 값들보다 낮으며, 상대적으로 작은 유효 받음각($3.1^\circ - 3.3^\circ$)에서 발생한다.

의의 및 주장
본 논문은 ELLIPSYS가 천이 유동의 벽면 해상 LES를 위한 실행 가능한 도구임을 입증하며, 장기 적분을 위한 계산 효율적인 대안을 제공한다. 본 연구는 높은 역류(흡입면에서 최대 -17%까지)로 인해 발생할 수 있는 절대 불안정성에 의해 트리거되는, 유동의 주기적인 스톨 및 탈스톨 현상이 관찰된 LFO의 원인임을 확인하였다.

중요한 발견은 두꺼운 에어포일(두께 24%)에서 낮은 받음각 조건에서도 LFO가 식별되었다는 점이며, 이는 문헌에서 흔히 다뤄지는 얇은 에어포일의 스톨 현상과는 구별되는 조건이다. 저자들은 에어포일 양측의 역류가 절대 불안정성을 허용할 만큼 충분하며, 이것이 LSB가 완전한 스톨 유동으로 급격히 붕괴(bursting)되는 것을 구동할 수 있다고 언급하였다. 본 연구는 두꺼운 에어포일 형상이 저주파 진동의 발생과 특성에 어떻게 영향을 미치는지에 대한 이해에 기여하며, 이러한 현상이 기존의 표준 에어포일에서 기록된 것보다 더 낮은 입사각과 낮은 스트로우할 수에서 발생할 수 있음을 강조한다.