23.27 프로펠러 비정상 공력(Unsteady Aerodynamics)
1. 비정상 공력의 범주
프로펠러의 비정상 공력(unsteady aerodynamics)은 시간 변화하는 유동 조건에서 발생하는 블레이드 단면 및 전체 프로펠러의 공력 현상을 총칭하는 학문 분야이다. 주요 원인은 다음과 같다. 첫째, 자유 흐름의 시간 변화(돌풍, 경사 유입). 둘째, 블레이드의 회전에 따른 방위각 의존적 유입 분포. 셋째, 블레이드 피치 또는 조종면 변위의 과도 입력. 넷째, 블레이드 변형에 의한 공탄성 응답. 다섯째, 블레이드-와류 상호작용에 의한 순간 유입 교란. 이 비정상 공력의 체계적 이해는 Leishman이 Principles of Helicopter Aerodynamics(2nd ed., Cambridge University Press, 2006)에서 종합적으로 제시한 바 있다.
2. 감소 주파수와 상사성
비정상 공력의 무차원 매개변수로 감소 주파수(reduced frequency) k가 사용된다.
k = \dfrac{\omega c}{2 U}
여기서 \omega는 변화 주파수, c는 시위, U는 국부 상대 속도이다. k \ll 1에서는 준정상(quasi-steady) 근사가 유효하며, 공력 응답이 순시 유동 조건으로 결정된다. k \sim 0.1 영역에서는 비정상 효과가 유의해지고, k \geq 0.3에서는 비정상 효과가 지배적이 된다. 프로펠러 블레이드의 전형적 감소 주파수는 회전 주파수와 블레이드 시위의 관계에 의해 k = 0.05 \sim 0.5 범위에 분포한다.
3. 비정상 양력의 구성 성분
비정상 조건에서 블레이드 단면의 양력은 다음 세 성분으로 분해된다.
L_{\text{unsteady}}(t) = L_{\text{circulatory}}(t) + L_{\text{non-circulatory}}(t)
첫째, 순환 성분(circulatory lift)은 블레이드 단면 주위 순환의 변화와 후류에 방출된 와류에 의해 발생한다. 둘째, 비순환 성분(non-circulatory lift) 또는 추가 질량 성분(added mass lift)은 주변 유체 질량이 블레이드와 함께 가속되는 관성력으로, 순시 가속도에 비례한다. 이 분해는 Theodorsen이 General Theory of Aerodynamic Instability and the Mechanism of Flutter(NACA Report No. 496, 1935)에서 정식화하였다.
4. Theodorsen 함수
조화(harmonic) 받음각 변화에 대한 2차원 비정상 양력 응답은 Theodorsen 함수 C(k)로 표현된다. C(k) = F(k) + i G(k)는 복소 함수이며, k \to 0에서 C \to 1, k \to \infty에서 C \to 0.5로 수렴한다. Theodorsen 함수는 주파수 영역에서 양력의 크기 감소와 위상 지연을 정량화한다. 블레이드 요소 해석에서 각 반경의 비정상 양력을 주파수 영역으로 분해하여 적용할 수 있다.
5. Wagner 함수와 시간 영역 응답
시간 영역 비정상 응답은 Wagner 함수 \Phi(s)로 표현된다. 단위 계단 받음각 입력에 대한 양력 응답은 다음과 같이 표현된다.
\dfrac{L(t)}{L_{\text{ss}}} = \Phi(s)
여기서 s = 2 U t / c는 무차원 시간이며, L_{\text{ss}}는 정상 상태 양력이다. R. T. Jones가 제시한 근사식 \Phi(s) \approx 1 - 0.165 e^{-0.0455 s} - 0.335 e^{-0.3 s}는 수치적 구현에 편리하다. Duhamel 적분을 통해 임의의 시간 입력에 대한 응답을 산출할 수 있다.
6. Küssner 함수와 돌풍 응답
자유 흐름 수직 성분의 급격한 변화, 즉 돌풍(sharp-edged gust)에 대한 비정상 양력 응답은 Küssner 함수 \Psi(s)로 표현된다. \Psi(s) \approx 1 - 0.5 e^{-0.13 s} - 0.5 e^{-1.0 s}와 같은 근사식이 사용된다. 프로펠러 블레이드가 수직 돌풍이나 난류를 관통할 때의 양력 응답을 해석하는 데 유용하며, 비행 안정성 해석에도 적용된다.
7. 동적 실속의 비정상 거동
받음각이 빠르게 정적 실속각을 초과할 때 발생하는 동적 실속은 비정상 공력의 비선형 현상이다. 순간적 양력이 정적 C_{L,\max}를 상회하고, 이후 동적 실속 와류의 방출에 의해 급감한다. Leishman-Beddoes 반경험 모형은 이러한 비정상 응답을 상태 공간 형식으로 표현하며, 프로펠러 및 회전익 시뮬레이션에 널리 사용된다.
8. 블레이드-와류 상호작용의 비정상성
블레이드-와류 상호작용(Blade-Vortex Interaction, BVI)은 고립된 와류가 블레이드에 유도하는 순간적 받음각 변화에 의해 비정상 양력 변동을 발생시킨다. 와류 통과 시 블레이드 단면 압력 분포가 급격히 변화하며, 이는 소음의 주요 원천이다. BVI 해석에는 자유 후류 모형과 고해상도 비정상 공력 해석이 필수이다.
9. 비정상 CFD 해석
비정상 공력의 고충실도 해석에는 비정상 Navier-Stokes 방정식 기반의 CFD가 사용된다. 해석 시 고려 사항은 다음과 같다. 첫째, 시간 스텝이 BPF와 비정상 변화의 최대 주파수를 해상해야 한다. 둘째, 격자 이동 기법(sliding mesh, overset mesh)이 회전 블레이드의 운동을 처리한다. 셋째, Detached Eddy Simulation 또는 Large Eddy Simulation이 박리와 BVI를 재현한다. Yoon 외의 Computational Aerodynamic Analyses of Multirotors in Forward Flight(5th Asian/Australian Rotorcraft Forum, 2016)와 Ventura Diaz와 Yoon의 High-Fidelity Computational Aerodynamics of Multi-Rotor Unmanned Aerial Vehicles(56th AIAA Aerospace Sciences Meeting, AIAA Paper 2018-1266, 2018)가 멀티로터 비정상 CFD 해석의 대표 사례이다.
10. 비정상 공력의 실측
비정상 공력의 실험적 검증은 풍동 내 조화 진동 시험(harmonic pitching test), 돌풍 발생기(gust generator)를 이용한 시험, 고속 PIV를 활용한 비정상 유동장 측정으로 수행된다. 이러한 시험은 주로 블레이드 단면에 해당하는 2차원 익형 모형 또는 축소 프로펠러 모형을 대상으로 진행되며, Theodorsen 함수 및 Wagner 함수의 검증, 동적 실속의 특성 파라미터 식별에 활용된다.
11. 로봇공학적 적용
비정상 공력 모형은 자율 비행 로봇의 정밀 제어와 시뮬레이션에 중요하다. 주요 활용은 다음과 같다. 첫째, 급기동 비행의 응답 예측. 둘째, 돌풍 및 난류 환경에서의 제어 강건성 확보. 셋째, 블레이드 구조 피로 해석을 위한 하중 스펙트럼 생성. 넷째, 소음 예측 모형의 입력 자료. 다섯째, 학습 기반 공력 잔차 보상 모형의 물리적 구조 제공. 이러한 활용은 비정상 공력을 단순한 보정 요소를 넘어 자율 비행 시스템의 핵심 설계 요소로 격상시킨다.
12. 출처
- Leishman, J. G. Principles of Helicopter Aerodynamics, 2nd ed. Cambridge University Press, 2006.
- Theodorsen, T. General Theory of Aerodynamic Instability and the Mechanism of Flutter. NACA Report No. 496, 1935.
- Jones, R. T. The Unsteady Lift of a Wing of Finite Aspect Ratio. NACA Report No. 681, 1940.
- Bisplinghoff, R. L., Ashley, H., and Halfman, R. L. Aeroelasticity. Addison-Wesley, 1955.
- Leishman, J. G., and Beddoes, T. S. “A Semi-Empirical Model for Dynamic Stall.” Journal of the American Helicopter Society, vol. 34, no. 3, 1989.
- Ventura Diaz, P., and Yoon, S. “High-Fidelity Computational Aerodynamics of Multi-Rotor Unmanned Aerial Vehicles.” 56th AIAA Aerospace Sciences Meeting, AIAA Paper 2018-1266, 2018.
13. 버전
v1.0 (2026-04-17)