23.19 와류 모델링 기법(Vortex Lattice Method)

23.19 와류 모델링 기법(Vortex Lattice Method)

1. 와류 격자법의 개념

와류 격자법(Vortex Lattice Method, VLM)은 비점성·비회전·비압축성 유동 가정하에 양력면을 이산 와류 격자로 근사하여 3차원 공력 해석을 수행하는 포텐셜 유동 기법이다. 날개, 블레이드, 프로펠러의 양력 분포, 유도 속도, 유도 항력 등의 해석에 사용되며, 패널법(panel method)과 유사하지만 두께가 없는 얇은 양력면(lifting surface)에 초점을 맞춘 단순화된 모형이다. Falkner가 The Calculation of Aerodynamic Loading on Surfaces of Any Shape(Aeronautical Research Council R&M 1910, 1943)에서 이론적 기초를 정립한 이래, VLM은 공력 해석의 표준 기법 중 하나로 자리 잡았다.

2. 수학적 기반

VLM은 라플라스 방정식 \nabla^2 \phi = 0의 기본 해 중 와류 요소의 중첩으로 유동장을 구성한다. 핵심 와류 요소는 수평 말굽 와류(horseshoe vortex)로, 한 쌍의 자유 와류 꼬리를 가진 유한 길이 구속 와류로 구성된다. Biot-Savart 법칙에 따라 각 말굽 와류는 공간 내 임의 점에서 다음과 같은 유도 속도를 제공한다.

\mathbf{v}_{\text{induced}}(\mathbf{x}) = \dfrac{\Gamma}{4 \pi} \oint \dfrac{d\boldsymbol{\ell} \times (\mathbf{x} - \mathbf{x}')}{|\mathbf{x} - \mathbf{x}'|^3}

여기서 \Gamma는 와류 순환, d\boldsymbol{\ell}은 와류 필라멘트의 미소 길이, \mathbf{x}'는 필라멘트 위치이다. 양력면 전체는 많은 수의 말굽 와류를 격자 형식으로 배치하여 모형화된다.

3. 양력면의 이산화

VLM에서 양력면은 시위 방향과 스팬 방향으로 분할된 사각형 패널의 격자로 구성된다. 각 패널에는 하나의 구속 와류(bound vortex)가 1/4 시위선에 배치되며, 패널 후연으로부터 무한 후방으로 두 개의 자유 와류 꼬리가 연장된다. 해의 미지수는 각 패널의 와류 순환 \Gamma_i이며, 이는 각 패널의 3/4 시위점(control point)에서의 불투수 조건(impermeability condition)에 의해 결정된다. 이 조건은 다음과 같이 표현된다.

\hat{\mathbf{n}}_i \cdot (\mathbf{V}_\infty + \sum_{j=1}^{N} \Gamma_j \mathbf{v}_{ij}) = 0

여기서 \hat{\mathbf{n}}_i는 패널 i의 외향 법선 벡터, \mathbf{v}_{ij}는 패널 j의 단위 순환이 패널 i의 제어점에서 유도하는 속도이다. 이 조건은 선형 대수 방정식 [A]\{\Gamma\} = \{b\}로 정리된다.

4. 양력과 항력의 산출

패널별 순환 \Gamma_i가 결정되면, Kutta-Joukowski 정리로부터 각 구속 와류에 작용하는 양력 벡터가 산출된다.

d\mathbf{L}_i = \rho \, (\mathbf{V}_\infty + \mathbf{v}_{\text{induced},i}) \times \Gamma_i \, d\boldsymbol{\ell}_i

전체 양력은 모든 패널에 대한 합산으로 얻어진다. 유도 항력은 Trefftz 평면(Trefftz plane)에서의 후류 와류 시트의 운동 에너지 유속으로부터 계산된다. VLM은 익형의 두께 효과와 점성 효과를 포함하지 않으므로, 결과는 이상 유도 항력 성분에 해당하며, 점성 항력은 별도의 경계층 해석기로 보완해야 한다.

5. 프로펠러 해석에의 확장

프로펠러 해석에서 VLM은 회전하는 블레이드를 양력면으로 표현하고, 블레이드 회전을 고려한 좌표계에서 격자를 구성한다. 자유 후류는 블레이드 회전과 축방향 이동에 의해 나선형 경로를 따르며, 각 블레이드의 후연에서 방출된 와류 시트가 이 경로에 분포된다. 이 접근은 Hough와 Ordway가 The Generalized Actuator Disk Theory(Developments in Theoretical and Applied Mechanics, vol. 2, 1965)에서 제시한 일반 액추에이터 디스크 해석의 후속으로 확장되었으며, 회전익 VLM은 현재까지 초기 설계 단계의 주요 해석 도구 중 하나이다.

6. 자유 후류와 강제 후류

VLM의 후류는 두 가지 접근으로 처리된다. 첫째, 강제 후류(prescribed wake) 접근은 후류 와류 시트의 기하를 사전에 결정한 형태(예를 들어 일정 피치의 나선)로 고정하고, 블레이드 양력 분포만 해석한다. 둘째, 자유 후류(free wake) 접근은 후류 와류 시트 자체가 유도 속도에 의해 이동하며 스스로 정합되는 형태로 변형되도록 반복 계산한다. 자유 후류 접근은 더 정확한 유도 속도 예측을 제공하지만, 계산 비용과 수렴성 관리가 복잡해진다. 이러한 구분은 회전익 VLM 해석의 기본 분류이다.

7. VLM의 장단점

항목장점단점
계산 비용낮음
설계 유연성신속한 형상 변형 및 최적화
양력 분포 예측3차원 양력 분포와 유도 항력 정확도
두께 효과무시됨
점성 효과별도 해석 필요
고받음각 박리예측 불가
압축성약한 압축성 보정 가능충격파 포함 불가

이 표는 VLM의 일반적 장단점을 요약한 것이며, 구체적 성능은 구현 상세와 적용 대상에 따라 달라진다.

8. 대표 구현과 적용 사례

AVL(Athena Vortex Lattice)은 Drela와 Youngren이 MIT에서 개발한 대표적 VLM 구현이며, 고정익 무인기, 소형 항공기, 안정성 계수 해석에 폭넓게 사용된다. 또한 Tornado, XFLR5 등의 공개 소프트웨어는 VLM에 기반한 분석 기능을 제공한다. 회전익 전용 VLM으로는 CAMRAD II, RCAS(Rotorcraft Comprehensive Analysis System) 등이 상용·연구 분야에서 활용된다.

9. 로봇공학적 적용

소형 무인기의 공력 해석에서 VLM은 다음과 같이 활용된다. 첫째, 고정익 무인기의 날개 및 꼬리 날개의 양력 분포와 유도 항력의 초기 예측. 둘째, 비행 동역학 모형에 필요한 안정성 계수(stability derivatives)의 산출. 셋째, 프로펠러 블레이드의 초기 설계 단계에서 반경 방향 양력 분포와 유도 속도의 추정. 넷째, 소형 전기 비행체의 최적 비행 속도 및 설계 진행비 결정. VLM은 계산 비용이 낮고 설계 파라미터 감도 분석이 용이하므로, 설계 초기 단계의 반복 해석에 적합하다.

10. 출처

  • Falkner, V. M. The Calculation of Aerodynamic Loading on Surfaces of Any Shape. Aeronautical Research Council R&M No. 1910, 1943.
  • Hough, G. R., and Ordway, D. E. “The Generalized Actuator Disk Theory.” Developments in Theoretical and Applied Mechanics, vol. 2, 1965.
  • Katz, J., and Plotkin, A. Low-Speed Aerodynamics, 2nd ed. Cambridge University Press, 2001.
  • Drela, M., and Youngren, H. AVL User Guide. MIT Technical Documentation, 2004.
  • Leishman, J. G. Principles of Helicopter Aerodynamics, 2nd ed. Cambridge University Press, 2006.

11. 버전

v1.0 (2026-04-17)