23.20 자유 후류 해석(Free Wake Analysis)

1. 자유 후류 해석의 개념

자유 후류 해석(free wake analysis)은 프로펠러 또는 회전익 블레이드에서 방출된 후류 와류 필라멘트가 주변 유도 속도장에 의해 이동하며 스스로 정합적인 기하를 형성하도록 반복 계산을 수행하는 해석 기법이다. 강제 후류(prescribed wake)가 후류 기하를 사전에 고정하는 것과 달리, 자유 후류는 유동장의 자기 유도(self-induction)를 통해 실제 후류의 비선형 변형을 재현한다. 이 접근은 근후류와 원후류의 유도 속도 정확도를 크게 개선하며, 회전익 호버링, 전진 비행, 근접 비행 시 기체 상호작용 해석에 필수적이다. Leishman의 Principles of Helicopter Aerodynamics(2nd ed., Cambridge University Press, 2006)에서 자유 후류 해석의 체계가 상세히 정리되어 있다.

2. Biot-Savart 법칙과 자기 유도

후류 와류 필라멘트는 유한 개의 선형 와류 세그먼트로 이산화되고, 각 세그먼트는 Biot-Savart 법칙에 따라 공간 내 임의 점에서 다음의 유도 속도를 제공한다.

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

후류 내 임의 지점 \mathbf{x}에서의 유도 속도는 모든 와류 세그먼트의 기여를 합산하여 얻어진다. 각 후류 표지점(marker)의 위치 \mathbf{x}(t)는 다음의 운동 방정식을 따라 시간에 대해 전개된다.

\dfrac{d \mathbf{x}}{dt} = \mathbf{V}_\infty + \mathbf{v}_{\text{induced}}(\mathbf{x})

이 방정식의 수치 적분을 통해 후류의 시간 변화와 최종 정상 기하가 얻어진다.

3. 와류 특이점의 수치적 처리

Biot-Savart 법칙의 분모에 포함된 |\mathbf{x} - \mathbf{x}'|^{-3} 항은 와류 필라멘트 중심에서 특이성(singularity)을 가진다. 이를 수치적으로 처리하기 위해 다음과 같은 특이점 제거 기법이 적용된다. 첫째, Rankine 코어 모형은 와류 중심 영역을 유한한 코어 반경 r_c로 확장하여 특이성을 제거한다. 둘째, Scully 모형, Vatistas 모형과 같은 경험적 코어 모형은 실제 관측된 와류 속도 분포에 근사한 코어 함수를 제공한다. Vatistas 외가 A Simpler Model for Concentrated Vortices(Experiments in Fluids, vol. 11, no. 1, 1991)에서 제시한 모형이 실용적으로 자주 사용된다.

4. 수치 적분 방법

자유 후류의 시간 적분에는 다양한 수치 기법이 사용된다. 대표적으로 Euler 전진법, 개선된 Euler법, 2차 Runge-Kutta법, 4차 Runge-Kutta법 등이 있다. 회전익 해석에서는 블레이드 회전과 연동하여 시간 스텝을 방위각 증가 단위로 설정하는 경우가 일반적이다. Bagai와 Leishman이 Rotor Free-Wake Modeling Using a Pseudo-Implicit Technique—Including Comparisons with Experimental Data(Journal of the American Helicopter Society, vol. 40, no. 3, 1995)에서 가상 암시적 기법(pseudo-implicit technique)을 제시하여, 자유 후류 해석의 수치 안정성과 수렴성을 크게 개선하였다.

5. 정상 및 비정상 자유 후류

자유 후류는 해석 목적에 따라 정상(steady) 또는 비정상(unsteady) 형태로 구분된다. 정상 자유 후류는 회전익의 주기적 운동 평균에 대한 후류의 평균 기하와 유도 속도를 계산한다. 비정상 자유 후류는 블레이드의 방위각 위치와 비정상 운동(급기동, 비주기적 조종면 변화)에 따른 후류의 시간 변화를 재현한다. 비정상 해석은 BVI 소음 예측, 기동 중의 하중 해석, 천이 비행 조건의 성능 예측에 사용된다.

6. 계산 비용과 병렬화

자유 후류의 기본 계산 비용은 O(N^2)으로, N은 후류 표지점의 수이다. 긴 후류와 세밀한 이산화는 계산 비용을 크게 증가시키므로, 최근 연구에서는 빠른 다극 전개(Fast Multipole Method)와 GPU 가속을 활용하여 O(N \log N) 또는 O(N)의 계산 복잡도로 감축하고 있다. 이러한 가속 기법은 실시간 시뮬레이션, 설계 최적화, 대규모 다중 로터 해석에서 유용하다.

7. 결과의 활용

이 표는 자유 후류 해석 결과의 대표적 활용을 요약한 것이다. 실제 해석은 블레이드 요소 해석과 결합하여 통합 성능 예측을 제공한다.

8. 다른 해석 기법과의 결합

자유 후류 해석은 일반적으로 블레이드 요소 이론 또는 VLM과 결합되어 사용된다. 블레이드 양력 분포는 블레이드 요소 해석으로 계산하고, 후류 와류 필라멘트의 강도는 양력 분포로부터 결정된다. 이 후류 필라멘트들이 자유롭게 이동하면서 유도 속도장을 형성하며, 유도 속도가 블레이드의 새로운 국부 받음각 분포를 결정한다. 이와 같은 반복 해석이 수렴하면 정합적인 공력 해가 얻어진다. 이 접근은 CAMRAD II, UMARC, CHARM 등의 종합 회전익 해석 코드에서 표준적으로 구현된다.

9. 실험 자료와의 검증

자유 후류 해석은 실험 자료와의 비교를 통해 정확성이 검증된다. 대표적 실험으로는 Caradonna와 Tung이 Experimental and Analytical Studies of a Model Helicopter Rotor in Hover(Vertica, vol. 5, no. 2, 1981)에서 제시한 모델 로터 호버링 시험이 있다. 이 시험은 블레이드 표면 압력, 후류 기하, 블레이드 하중을 동시에 측정하여 자유 후류 해석의 검증 기준으로 사용된다. 또한 PIV 기반 실험은 후류의 순간 속도장을 측정하여, 자유 후류 해석의 시간 변화 예측과 비교하는 자료를 제공한다.

10. 로봇공학적 활용

자유 후류 해석은 멀티로터 무인기와 도심항공교통 기체의 정밀 공력 해석에 사용된다. 특히 다음의 응용에서 필수적이다. 첫째, 인접 로터 간 후류 간섭 해석. 둘째, 프로펠러 후류가 기체 날개 또는 수평 안정판에 유도하는 공력 변화 해석. 셋째, 지면 근접 비행 시 재순환 유동의 정량화. 넷째, BVI 소음의 방위각 분포 예측. 자유 후류 해석은 이들 응용에서 CFD 대비 계산 비용이 낮으면서도 후류 구조의 주요 특징을 정확히 재현한다.

11. 출처

• Leishman, J. G. Principles of Helicopter Aerodynamics, 2nd ed. Cambridge University Press, 2006.
• Bagai, A., and Leishman, J. G. “Rotor Free-Wake Modeling Using a Pseudo-Implicit Technique—Including Comparisons with Experimental Data.” Journal of the American Helicopter Society, vol. 40, no. 3, 1995.
• Vatistas, G. H., Kozel, V., and Mih, W. C. “A Simpler Model for Concentrated Vortices.” Experiments in Fluids, vol. 11, no. 1, 1991.
• Caradonna, F. X., and Tung, C. “Experimental and Analytical Studies of a Model Helicopter Rotor in Hover.” Vertica, vol. 5, no. 2, 1981.
• Johnson, W. Helicopter Theory. Princeton University Press, 1980.

12. 버전

v1.0 (2026-04-17)