High-order gas-kinetic scheme for numerical simulations of wind turbine with nacelle and tower using ALM and IBM

본 논문은 나셀과 타워를 포함한 풍력 터빈의 정밀한 시뮬레이션을 위해 고차 가스 운동론적 scheme(GKS) 에 액추에이터 라인 모델 (ALM) 과 잠입 경계법 (IBM) 을 통합하고 GPU 병렬 계산을 적용하여, 블레이드 - 타워 상호작용 및 와류 간섭 효과를 정확히 예측하는 새로운 수치 기법을 제안합니다.

핵심

1. 문제: 거대한 풍력 터빈을 어떻게 시뮬레이션할까?

풍력 터빈은 지름이 200m 이상인 거대한 구조물입니다. 날개가 회전할 때 공기를 가르며 복잡한 소용돌이를 만들어내는데, 이 소용돌이는 매우 작고 빠르게 변합니다.

• 기존의 어려움: 컴퓨터로 이 모든 것을 계산하려면 엄청난 시간이 걸립니다. 마치 거대한 폭포를 한 방울의 물방울 단위로 세어보려는 것과 비슷하죠.
• 간단한 방법의 한계: 예전에는 터빈의 '날개'와 '기둥 (타워)', '머리 부분 (네이클)'을 모두 정밀하게 모델링하지 않고, 날개만 대충 표시하거나 아예 기둥을 무시하는 경우가 많았습니다. 하지만 실제 바다에 설치된 터빈은 기둥과 머리 부분의 영향도 중요하기 때문에, 이를 무시하면 정확한 예측이 어렵습니다.

2. 해결책: 세 가지 도구를 하나로 합치다

이 연구팀은 세 가지 강력한 도구를 하나로 섞어 **'초고해상도 풍력 터빈 시뮬레이션'**을 만들었습니다.

① 액추에이터 라인 모델 (ALM) = "날개의 유령"

실제 날개의 복잡한 모양을 다 그리지 않고, 날개 위치에 **'가상의 힘'**을 가하는 점들을 배치합니다.

• 비유: 마치 날개 대신 공기 중에 **'보이지 않는 손'**이 바람을 밀어낸다고 상상하는 것입니다. 날개의 모양을 다 그리지 않아도, 그 손이 바람에 미치는 힘만 정확히 계산하면 됩니다.

② 침수 경계법 (IBM) = "유연한 벽"

기존에는 컴퓨터 격자 (그물망) 가 터빈의 모양에 딱 맞춰져야 했습니다. 하지만 터빈은 둥글고 복잡하죠. IBM 은 격자를 그대로 두고, 그물망 위에 **'투명한 벽'**을 만들어내는 기술입니다.

• 비유: 물이 흐르는 강에 돌멩이를 넣으면 물이 돌을 돌아서 흐르죠? IBM 은 돌멩이 (타워와 네이클) 를 물속의 격자 안에 **'투명한 장벽'**처럼 배치해서, 격자를 다시 만들지 않아도 물 (바람) 이 자연스럽게 돌아서 흐르게 합니다.

③ 고차 기체 운동론적 scheme (GKS) = "초고속 카메라"

이 모든 것을 계산하는 핵심 엔진입니다. 기존에 쓰이던 2 차원 (저해상도) 방식 대신, 4 차원 (초고해상도) 방식을 사용합니다.

• 비유: 기존 방식이 1 초에 30 장의 사진을 찍는 일반 카메라라면, 이 새로운 GKS 는 1 초에 수천 장의 고화질 사진을 찍는 초고속 카메라입니다. 바람이 만들어내는 아주 작은 소용돌이 (난류) 까지 놓치지 않고 선명하게 포착합니다.

3. 실험 결과: 무엇을 발견했나?

이 연구팀은 이 새로운 방법을 두 가지 사례에 적용해 보았습니다.

사례 1: NREL 5MW 터빈 (날개만 있는 경우)

• 결과: 날개의 힘과 뒤에 남는 바람의 흐름을 기존 방법보다 훨씬 정확하게 예측했습니다.
• 발견: 고해상도 (4 차) 카메라를 쓰니, 날개 끝에서 생기는 작은 소용돌이가 어떻게 커져서 큰 난류가 되는지 그 과정을 훨씬 더 일찍, 더 선명하게 볼 수 있었습니다.

사례 2: NTNU 터빈 (날개 + 기둥 + 머리 포함)

• 결과: 이제 기둥 (타워) 과 머리 (네이클) 까지 모두 포함했습니다.
• 놀라운 발견:
  • 주기적인 리듬: 날개가 기둥 앞을 지날 때마다 바람의 흐름이 살짝 흔들립니다. 마치 시계추가 흔들리듯, 터빈의 출력과 힘이 규칙적으로 변하는 것을 포착했습니다.
  • 소용돌이의 충돌: 기둥 뒤에서 생기는 소용돌이와 날개 끝에서 생기는 소용돌이가 부딪히면서, 바람이 더 일찍 '난기류' 상태로 변했습니다.
  • 비대칭성: 바람이 터빈을 지나갈 때, 기둥 때문에 바람의 흐름이 한쪽으로 치우치는 현상 (비대칭) 을 정확히 재현했습니다.

4. 결론: 왜 이 연구가 중요한가?

이 연구는 "GPU(그래픽 카드) 를 여러 개 연결해서" 이 복잡한 계산을 빠르게 처리할 수 있게 했습니다.

• 기존: 풍력 터빈의 기둥을 무시하거나, 날개만 대충 계산하면 실제 풍력 발전량을 예측하기 어렵습니다.
• 이제: 이 새로운 방법으로 날개, 기둥, 머리를 모두 포함해서, 실제 바다나 땅에 설치된 터빈의 성능을 정확하고 빠르게 예측할 수 있게 되었습니다.

한 줄 요약:

"이 논문은 거대한 풍력 터빈이 만들어내는 복잡한 바람의 흐름을, **초고해상도 카메라 (고차 GKS)**로 찍고, **가상의 손 (ALM)**과 **투명한 벽 (IBM)**을 이용해 날개와 기둥까지 모두 포함하여 정확하게 예측할 수 있는 새로운 방법을 개발했습니다."

이 기술은 앞으로 더 크고 효율적인 풍력 터빈을 설계하고, 신재생 에너지를 더 많이 확보하는 데 큰 도움이 될 것입니다.

기술 요약

논문 요약: ALM 및 IBM 을 통합한 고차 기체 운동론적 scheme (GKS) 을 이용한 나셀과 타워가 포함된 풍력 터빈 시뮬레이션

1. 연구 배경 및 문제 제기 (Problem)

• 풍력 터빈의 대형화: 최근 해상 풍력 발전의 확대로 풍력 터빈의 크기가 급격히 증가하고 있으며 (지름 200m 이상), 이로 인해 날개 흐름의 레이놀즈 수가 1,000 만 수준에 달하게 되었습니다.
• 시뮬레이션의 난제: 이러한 대형 터빈은 와류 (vortex) 가 지배적인 다중 스케일 난류 구조를 생성하며, 이를 정확하게 포착하기 위해서는 높은 공간 및 시간 해상도가 필요합니다. 기존의 2 차 (second-order) 수치 기법으로는 이러한 난류 와류를 효율적으로 해석하기 어렵습니다.
• 기하학적 복잡성: 풍력 터빈의 날개뿐만 아니라 나셀 (nacelle) 과 타워 (tower) 의 존재는 국부 유속을 감소시키고 날개 - 타워 상호작용을 통해 와류의 전이를 유발하여 유동장에 중요한 영향을 미칩니다. 그러나 많은 기존 연구에서는 계산 비용 절감을 위해 나셀과 타워를 생략하거나 단순화했습니다.
• 목표: 복잡한 기하학적 구조 (나셀, 타워) 를 포함하는 풍력 터빈의 와류 지배적 난류 후류 (wake) 를 정확하고 효율적으로 시뮬레이션할 수 있는 새로운 수치 프레임워크 개발.

2. 방법론 (Methodology)

이 연구는 고차 기체 운동론적 scheme (High-order Gas-Kinetic Scheme, GKS) 을 기반으로 하여, 작동선 모델 (Actuator Line Model, ALM) 과 침수 경계법 (Immersed Boundary Method, IBM) 을 통합한 새로운 접근법을 제시합니다.

• 고차 GKS 의 확장:
  • 기존의 2 차 GKS 를 3 차원 약한 압축성 (weakly compressible) 등온 유동 시뮬레이션으로 확장했습니다.
  • 4 차 정확도를 갖는 2 단계 시간 적분법 (two-stage fourth-order framework) 과 5 차 WENO (Weighted Essentially Non-Oscillatory) 공간 재구성 기법을 적용하여 높은 해상도를 확보했습니다.
• ALM (날개 모델링):
  • 회전하는 날개를 메시화하지 않고, 날개 방향을 따라 분포된 작동점 (actuator points) 의 집합으로 표현합니다.
  • 유동장에서 속도를 샘플링하여 양력/항력을 계산하고, 가우시안 커널을 통해 이를 주변 유동장에 체적력 (body force) 으로 분포시킵니다.
• IBM (나셀 및 타워 모델링):
  • 복잡한 형상의 나셀과 타워를 오일러 격자 (Eulerian mesh) 와 라그랑주 점 (Lagrangian points) 의 집합으로 표현합니다.
  • 직접 힘 (direct-force) 형식을 사용하여 경계 조건을 강제하며, 격자와 입자 간의 상호작용을 통해 유동장에 외력을 부과합니다.
• 통합 및 병렬화:
  • ALM 과 IBM 에서 계산된 힘은 운동량 방정식의 외력 항으로 통합됩니다.
  • 대규모 난류 후류 시뮬레이션을 위해 다중 GPU (Multiple-GPU) 환경 (Nvidia Tesla V100, RTX A5000) 에서 MPI 와 CUDA 를 활용한 병렬 계산이 구현되었습니다.

3. 주요 기여 (Key Contributions)

1. 최초 통합: 고차 GKS 프레임워크에 ALM 과 IBM 을 동시에 통합하여, 날개, 나셀, 타워가 모두 포함된 풍력 터빈의 3 차원 시뮬레이션을 수행했습니다.
2. 고차 정확도 검증: 약한 압축성 난류 채널 유동과 원통 주위 유동 시뮬레이션을 통해 고차 GKS 의 수치적 정확도를 검증했습니다.
3. 스미어링 커널 (Smearing Kernel) 영향 분석: NREL 5MW 터빈 시뮬레이션을 통해 ALM 의 가우시안 커널 폭 ($\epsilon$) 이 날개 하중 및 후류 회복에 미치는 영향을 분석하고 최적값을 도출했습니다.
4. 고차 vs 2 차 비교: 동일한 격자 해상도에서 고차 GKS 가 2 차 GKS 보다 와류 지배적 난류 (vortex-dominated turbulence) 를 훨씬 효과적으로 해결함을 입증했습니다.
5. 타워 - 날개 상호작용 규명: NTNU "Blind Test 1" 데이터를 기반으로 나셀과 타워를 포함한 시뮬레이션을 수행하여, 타워 와류와 팁 와류의 상호작용이 후류의 비대칭성과 난류 전이 시기를 어떻게 변화시키는지 정량적으로 규명했습니다.

4. 주요 결과 (Results)

• NREL 5MW 터빈 (ALM 만 적용):
  • $\epsilon = 3\Delta x$일 때 ALM 예측값이 기존 참조 해 (Reference solutions) 와 가장 잘 일치했습니다.
  • 고차 GKS 는 2 차 GKS 에 비해 동일한 격자에서 더 작은 와류 구조를 해상도하며, 특히 과도기적 후류 (transitional wake) 와 원거리 후류 (far wake) 에서 더 빠른 속도 회복과 정확한 난류 운동 에너지 (TKE) 분포를 보였습니다.
• NTNU 터빈 (ALM + IBM 적용, 나셀 및 타워 포함):
  • 주기적 하중: 타워가 있을 때 (ALM+N+T), 블레이드가 타워를 통과할 때마다 양력과 추력 계수 ($C_P, C_T$) 가 주기적으로 변동하는 것을 포착했습니다. 타워가 없을 때는 정상 상태 값을 보였습니다.
  • 와류 상호작용: 타워 뒤에서 발생하는 카르만 와류 (Kármán street) 가 팁 와류와 상호작용하여 난류 전이 (transition) 를 더 일찍 발생시킵니다.
  • 비대칭성: 타워의 영향으로 인해 후류의 속도와 TKE 분포가 비대칭적이 되었으며, 이는 실험 데이터와 높은 일치도를 보였습니다.
  • 정확도: 고차 GKS 와 IBM 을 결합한 방법은 나셀과 타워가 포함된 복잡한 유동장의 평균 속도와 난류 운동 에너지를 실험 데이터와 매우 잘 일치하게 예측했습니다.

5. 의의 및 결론 (Significance)

• 정확성과 효율성: 본 연구에서 개발된 다중 GPU 기반의 고차 GKS-ALM-IBM 통합 방법은 복잡한 기하학적 구조를 가진 풍력 터빈의 고 레이놀즈 수 난류 후류 시뮬레이션을 위한 정확하고 효율적인 도구를 제공합니다.
• 실제 적용 가능성: 기존의 단순화된 모델 (타워/나셀 생략) 의 한계를 극복하여, 실제 풍력 터빈 설계 및 성능 평가에 필요한 더 현실적인 유동 정보를 제공할 수 있습니다.
• 향후 전망: 본 연구는 대기 경계층의 난류 유입 (turbulent inflows) 과 전단 속도 프로파일을 향후 연구에 통합하여 실제 대기 조건 하의 풍력 터빈 거동을 더 정밀하게 예측하는 기반을 마련했습니다.

요약하자면, 이 논문은 고차 수치 기법 (GKS) 과 모델링 기법 (ALM, IBM) 의 결합을 통해 풍력 터빈의 핵심 구성 요소 (날개, 나셀, 타워) 가 모두 포함된 복잡한 난류 유동을 고해상도로 해석할 수 있는 새로운 패러다임을 제시한 중요한 연구입니다.